Re: [R] zero random effect sizes with binomial lmer

From: Andrew Robinson <A.Robinson_at_ms.unimelb.edu.au>
Date: Tue 02 Jan 2007 - 02:14:30 GMT

Good observation, Daniel. I had not picked up the different default fitting method.

Cheers

Andrew

On Mon, Jan 01, 2007 at 08:59:55PM -0500, Daniel Ezra Johnson wrote:
> > From: Andrew Robinson <A.Robinson_at_ms.unimelb.edu.au>
> > Date: Mon 01 Jan 2007 - 19:19:29 GMT
> >
> > I tried an earlier version of R, on a different platform, and got quite
> > different results. Sadly, the *earlier* results are the ones that make
> > sense.
>
> Andrew, I tried installing R 2.3.1, it seemed to work as you said, but I
> realized it was because the default method was PQL instead of Laplace.
> Switching the method to Laplace, the zero effect came back.
>
> It seems like there's no easy way around this. I would like to compare
> random effect sizes between different subsets of my data, but this may not
> be possible using the current functions...
>
> Daniel
>
> ______________________________________________
> 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
> and provide commented, minimal, self-contained, reproducible code.

-- 
Andrew Robinson  
Department of Mathematics and Statistics            Tel: +61-3-8344-9763
University of Melbourne, VIC 3010 Australia         Fax: +61-3-8344-4599
http://www.ms.unimelb.edu.au/~andrewpr
http://blogs.mbs.edu/fishing-in-the-bay/

______________________________________________
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
and provide commented, minimal, self-contained, reproducible code.
Received on Tue Jan 02 13:22:37 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Tue 02 Jan 2007 - 02:30:28 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.