I suggest you post this to the r-sig-mixed-models list instead. You are much more likely to get useful answers there.
Cheers, Bert On Thu, Sep 26, 2013 at 11:25 AM, Bettina Lado Lindner <betti...@hotmail.com > wrote: > 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. > > -- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm [[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.
