Displaying 1 result from an estimated 1 matches for "predsxx".
Did you mean:
preds
2004 Mar 01
1
glm logistic model, prediction intervals on impact af age 60 compared to age 30
...one a logistic glm using Age as explanatory variable for some
allergic event.
#the model
model2d<-glm(formula=AEorSAEInfecBac~Age,family=binomial("logit"),data=emrisk)
#predictions for age 30 and 60
preds<-predict(model2d,data.frame(Age=c(30,60)),se.fit=TRUE)
# prediction interval
predsxx<-cbind(fit=preds$fit,lower=preds$fit-1.96*preds$se,upper=preds$fit+1.96*preds$se)
#transformation
model2dres<-family(model2d)$linkinv(predsxx)
In my next step I want to know the confidence interval (CI) for the change
in risk for the allergic event to occur for age 60 compared to age 30.
T...