[R] Is it possible with two random effects in lme()?

[Bettina Lado Lindner]

Hello, 
I want two random effect in R. One is BLOQUE nested in REP and the other is 
TRAT. I don’t know how put TRAT independant of the other random effect.

My model until this step is ok:
data_2010_ECCI_ensayo$le1=groupedData(REND~ REP|BLOQUE, 
data=data_2010_ECCI_ensayo$le1)
lmer_ECC2010_le1AR1BI=lme(REND ~ REP + TRAT, 
data=data_2010_ECCI_ensayo$le1,random=~1,correlation=corAR1())

If I want put TRAT as other independant random effect I don’t know how to do. 
When I did this give me an error...
lmer_ECC2010_le1AR1BI=lme(REND ~ REP + TRAT, 
data=data_2010_ECCI_ensayo$le1,random=list(~1, TRAT=~1),correlation=corAR1())
lmer_ECC2010_le1AR1BI=lme(REND ~ REP + TRAT, 
data=data_2010_ECCI_ensayo$le1,random=list(~1, ~1|TRAT),correlation=corAR1())

I know that in lmer is easier but I use error correlation so I need use lme.
Bettina
        [[alternative HTML version deleted]]


______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.