Displaying 1 result from an estimated 1 matches for "lmeres2".
Did you mean:
lme_res2
2011 Jan 17
1
Using anova() with glmmPQL()
...>78,1,0)
lmeres1=lme(fixed=stmassheavy~mnth*roostsitu,random=~1|subject/mnth,na.action=na.exclude,family=binomial)
anova(object=lmeres1,test="Chisq")
but I get errors doing that. After a certain amount of web searching, I find that I'm supposed to use glmmPQL for this so I tried:
lmeres2=glmmPQL(fixed=stmassheavy~mnth*roostsitu,random=~1|subject/mnth,na.action=na.exclude,family=binomial)
anova(object=lmeres2,test="Chisq")
The glmmPQL command runs, but I get "Error in anova.glmmPQL(object = lmeres, test = "Chisq") : 'anova' is not available for PQL...