[R-sig-ME] Why se.fit differ in predict.glm and predict.glmmadmb?

Xandre alex at iim.csic.es
Tue May 3 22:51:28 CEST 2016


Dear list,

I am running a GLM (family="binomial") without random effects using both 
glm and glmmadmb.

Summaries are almost identical, however when I used the predict function 
as follows:

predict(glm1,newdatos1,type="link",se.fit =T)

predict(admb1,newdatos1,type="link",se.fit =T)

I realized that se.fit differ a lot between them, admb se.fit resulted 
much much higher (fit is almost identical). This is just and example of 
what I found:

glm1$se.fit	admb1$se.fit
0.04290869	0.2676562
0.04435600	0.2733130
0.04095631	0.2728592
0.03402992	0.2718389
0.03000669	0.2713617
0.03633637	0.2722059

Maybe I'm missing something or I am making a big mistake. Any help with 
this?


Many thanks,

Alexandre Alonso


	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list