Re: [R] R-help Digest, Vol 32, Issue 26

From: Alan Cobo-Lewis <>
Date: Wed 26 Oct 2005 - 22:52:36 EST on Wednesday, October 26, 2005 at 6:00 AM -0500 wrote:

Try Harold's suggestion. The df still won't agree, because lmer (at least in its current version) just puts an upper bound on the df. But that should be OK, because all those t tests are approximations anyways, and you can get better confidence intervals (credible intervals, whatever) by using the mcmcsamp() function that works with lmer() alan


"Doran, Harold" <> responded:
>>There is an issue with implicit nesting in lmer. In your lme() model you nest
>>block/irrigation/density/fertilizer. In lmer you need to do something like
>>(I dind't include all of your variables, but I think the makes the point)
>>lmer(yield~irrigation*density*fertilizer+(1|fertilizer:density)+(1|density), data)
>>Which notes that fertilizer is nested in density.
>>Try this and then compare the results.
"Ronaldo Reis-Jr." <>, wrote:
>>I make the correct model with aov, lme do compare with lmer.
>>But I cant make a correct model in lmer. Look that the aov and lme results are
>>similars, but very different from lmer. In aov and lme is used the correct DF
>>for each variable, in lmer it use a same DF for all? Denom=54.

Alan B. Cobo-Lewis, Ph.D.		(207) 581-3840 tel
Department of Psychology		(207) 581-6128 fax
University of Maine
Orono, ME 04469-5742

______________________________________________ mailing list
PLEASE do read the posting guide!
Received on Wed Oct 26 23:13:29 2005

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