Re: [R] lmm WITHOUT random factor (lme4)

From: Mark Difford <mark_difford_at_yahoo.co.uk>
Date: Thu, 17 Mar 2011 08:13:17 -0700 (PDT)

On Mar 17, 2011; 04:29pm Thierry Onkelinx wrote:

>> You cannot compare lm() with lme() because the likelihoods are not the
>> same. Use gls() instead of lm()

Hi Thierry,

Of course, I stand subject to correction, but unless something dramatic has changed, you can. gls() can be used if you need to accommodate a correlation structure.

The method I have outlined, i.e. anova(lme$obj, lm$obj), is detailed in Pinheiro & Bates (2000) beginning page 154. Please refer to this if you doubt me.

Regards, Mark.

--
View this message in context: http://r.789695.n4.nabble.com/lmm-WITHOUT-random-factor-lme4-tp3384054p3384802.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
R-help_at_r-project.org 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 Thu 17 Mar 2011 - 17:28:29 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 17 Mar 2011 - 17:50:22 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.

list of date sections of archive