[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