[R-sig-ME] crossed effects with lmer but correlation structure with lme
Ben Bolker
bbolker at gmail.com
Fri Mar 30 22:09:33 CEST 2012
Denis Vile <denis.vile at ...> writes:
>
> Hi all,
>
> Was my problem not sufficiently well exposed or no one could help me ?
>
> Sincerely,
>
> Denis
The canonical reference for this is p. 163ff of Pinheiro
and Bates 2000 (section 4.2.2, http://tinyurl.com/crossedRE )
You are correct that lme4 doesn't handle 'R-side' correlation
structures, nor will it in the near future ...
(I've added this information to glmm.wikidot.com/faq ...)
Ben Bolker
>
> Le 29/03/2012 10:09, Denis Vile a écrit :
> > Dear R users,
> >
> > I'm trying to fit a crossed-effects mixed model that would include a
> > spatial correlation structure..
> > The data come from four controlled experiments (control, treatment1,
> > treatment2, treatment1+treatment2) on plants grown in a growth
> > chamber. Individual replicates of different genotypes were grown
> > together and response traits were measured. A covariate X is included
> > in the model with a quadratic form.
> >
> > We fitted the following model using lmer:
> >
> > fm1 <- lmer(Y ~ 1 + Trt1*Trt2*poly(X, degree=2, raw=T) +
> > (1|idGenotype) + (1|Trt2:idGenotype) + (1|Trt1:idGenotype) +
> > (1|Trt1:Trt2:idGenotype), data=...)
> >
[snip snip snip]
More information about the R-sig-mixed-models
mailing list