Re: [R] anova on binomial LMER objects

From: Robert Bagchi <>
Date: Wed 28 Sep 2005 - 21:13:56 EST

Hi Patrick

thanks for your advice. I have now tried glmmPQL, and it worked fine - I'm getting consistent results between plots and models fitted by glmmPQL. Plus it allows predict() and resid() which is another advantage over lmer at present.

quick question though: why does one need to use PQL for binomial models? Is there a good reference for this?

A few of my colleagues have also had similar problems, so I'm copying this message on to R-help as it might be useful there.

Many thanks

Patrick A. Jansen wrote:

> Hi dr Bacghi,
> I ran into exactly the same problem with lmer models that had an
> rbind() response variable, as you posted to the R-list.
> The sums of squares produced by anova() seem wrong. They are almost
> identical to the mean squares, and hence F-values approach 1.
> Since you need PQL for binomial models anyway, you might want to use
> GlmmPQL instead. Seems to work fine with anova().
> Best regards,
> Patrick Jansen
> *dr Patrick A. Jansen*
> University of Groningen
> Community and Conservation Ecology group
> E
> W
> <>
> c/o:
> Instituto Smithsonian de Investigaciones Tropicales
> Att. Patrick A. Jansen – Gamboa
> Apartado 0843-03092, Balboa, Ancón, Panamá, República de Panamá
> or:
> Smithsonian Tropical Research Insititute
> Att: Patrick A. Jansen – Gamboa
> Unit 0948, APO AA, 34002-0948, U.S.A.
> T +507-212-8904 (office) / +507-6516-2008 (cell)
> mailing list PLEASE do read the posting guide! Received on Wed Sep 28 21:22:42 2005

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