[R-sig-ME] crossed effects with lmer but correlation structure with lme

Denis Vile denis.vile at supagro.inra.fr
Sun Apr 1 22:31:48 CEST 2012


Thank you very much Ben ! I'm going to look at this attentively.
Sincerely,
Denis

Le 30/03/2012 22:09, Ben Bolker a écrit :
> 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]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

-- 

*Denis VILE*
Chargé de Recherche
Laboratoire d'Ecophysiologie des Plantes sous Stress Environnementaux 
(*LEPSE*)
UMR 759 *INRA*-SUPAGRO // Institut de Biologie Intégrative des Plantes 
(IBIP, bât 7)
2 place Pierre Viala
34060 Montpellier Cedex 2
Tel +33 (0)4 99 61 31 87
Fax +33 (0)4 67 52 21 16
http://www1.montpellier.inra.fr/ibip/lepse/




More information about the R-sig-mixed-models mailing list