[R] anova with models from glmmPQL

From: Emmanuelle TASTARD <tastard_at_cict.fr>
Date: Wed 19 Oct 2005 - 22:19:51 EST


Hi !  

I try to compare some models obtained from glmmPQL.  

model1 <-
glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4 +I(freq8_4^2), random=~1|num, binomial); model2 <-
glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4 , random=~1|num, binomial);
anova(model1, model2)  

here is the answer :
Erreur dans anova.lme(model1, model2) : Objects must inherit from classes "gls", "gnls" "lm","lmList", "lme","nlme","nlsList", or "nls"  

Yet, according to “?glmmPQL”, glmmPQL gives an object of class lme: Value:

     A object of class '"lme"': see 'lmeObject'.  

What should I do to be able to compare my models?  

Thank you.
Emmanuelle    

Emmanuelle TASTARD
UMR 5174 'Evolution et Diversité Biologique' Université Paul Sabatier Bat 4R3
31062 TOULOUSE CEDEX 9 France
tel : 05 61 55 67 59  

        [[alternative HTML version deleted]]



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Oct 19 22:24:10 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:40:46 EST