Ronaldo,
According to Douglas Bates's paper in 'R' News, It would seem that the
correct model for nested split plot random effects with lmer , in your
example ,with x2 nested within x1, would be --
lmer(y~x1 + x2 +(1|x1)+(1|x1:x2))
Try it with your model any see how it compares with your aov
Em Seg 24 Out 2005 18:08, Doran, Harold escreveu:
Ronaldo
See the article on lmer pasted below for syntax. It is the only current
source documenting the code. In lmer(), the nesting structure for the
ranmdom effects is handled in a slightly different way. If your
observations are nested as
|density),
data)
Which notes that fertilizer is nested in density.
Try this and then compare the results.
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ronaldo Reis-Jr.
Sent: Tuesday, October 25, 2005 12:21 PM
To: R-Help
Subject: Re: [R] lme and lmer
Em Ter 25 Out 2005 15:24, Doran, Harold escreveu:
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)
Hi,
I have this:
lme(y~x1+x2,random=~1|x1/x2)
How to make this random effect using lmer?
I try this:
lmer(y~x1+x2+(1|x1/x2)
But it dont work.
Any idea?
Thanks
Ronaldo
--
System halted!
--
| // | \\ [***]
| ( õ õ ) [Ronaldo Reis Júnior
://CRAN.R-project.org/doc/Rnews/},
}
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ronaldo Reis-Jr.
Sent: Monday, October 24, 2005 3:56 PM
To: R-Help
Subject: [R] lme and lmer syntax
Hi,
I have this:
lme(y~x1+x2,random=~1|x1/x2)
How to make this random