Hello R-helpers,

 

I would like to calculate the standard error for the predicted value
from geeglm.  

As an example, I would like to calculate the GEE mean of treatments and
their standard error.  I first specified the model as 

 mod <- geeglm(resp ~ trt,
data=dat,id=id,family=Gaussian,corstr="ar1",weights=weight)

 

Then I predicted the GEE mean and se using the following code

pred.mean <-
predict(mod,data.frame(trt=factor(unique(dat$trt),levels=levels(dat$trt)
)),se=T)

 

I got the error message "Error in XRinv^2 %*% rep(res.var, p) :
non-conformable arguments".  

 

Thanks for your help,

 

Sue

 

 


	[[alternative HTML version deleted]]

