Displaying 1 result from an estimated 1 matches for "predictionresponse".
2008 May 28
1
confidence interval for the logit - predict.glm
...i, predallx, type="link", se.fit=TRUE)
upperlogit <- predictionlogit$fit + 2*predictionlogit$se.fit
lowerlogit <- predictionlogit$fit - 2*predictionlogit$se.fit
upperresponse[,i] <- exp(upperlogit)/(1+exp(upperlogit))
lowerresponse[,i] <- exp(lowerlogit)/(1+exp(lowerlogit))
predictionresponse <- predict(stepmodfi, predallx, type="response", se.fit=FALSE)
fishoccur[,i] <- predictionresponse
type="link" should be the equivalent to the S-Plus version of type="1", which was indicated in the online posting explaining this procedure.
So I first tried to...