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]]

