Hi: Thank you very much,and for your case if I wirte as regress=lme(yield~nitro*Variety,data=Oats,random=~Variety|Block)
what this mean? Regards 2008/6/22 Daniel Malter <[EMAIL PROTECTED]>: > Hi, > > random=~1|B/C > C is nested in B > > ##Example > data=Oats > regress=lme(yield~nitro*Variety,data=Oats,random=~1|Block/Variety) > ##i.e. variety is nested in block > summary(regress) > ##End of example > > Best, > Daniel > > > ------------------------- > cuncta stricte discussurus > ------------------------- > > -----Ursprüngliche Nachricht----- > Von: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Im > Auftrag von Manli Yan > Gesendet: Sunday, June 22, 2008 9:30 PM > An: r-help@r-project.org > Betreff: [R] problem in R for Linear mixed model~ > > Dear R users: > I just got confused some R code used in linear mixed model~ > example,two factors,A, B,C,A is fixed ,B,C are random,and B is nested in > C,if I wannt to use linear mixed model,are the following code correct for > each case? > case1:want to know random effect of B, > case1<-lme(y~A*B*C,random=~B|C) where "B|C" stand for what?,mean B is > nested in C? > > case2: how to wirte random effect of C? > case2<-lme(y~A*B*C,random=~C)? this doesnt work out,it seem it must have > somehing like #|$ > > case3.omitting the random effect for B from case1 > case3<-update(case1,random=~1|C),so I just type 1,so the random effect of > B will be removed from the model,there only left random effect of c ,the > random effect I removed ,which include both random intercept and slope > ,correct?? > > case4:omitting the random intercept > case4<-update(case1,random=B-1|C) > this code I got from some paper,it said by inputing B-1|C,then the random > intercept is removed,so,if I want to remove random slode,I input B-2|C,it > doesnt work out. > > case5 :how to know the both random effect of B,and C,I dont know how to > wirtie this in R . > And I am a little confused of these R code,especially the #|# part,what > deos this syntax really mean in LME package, > > Thank you for your time~ > > [[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<http://www.r-project.org/posting-guide.html> > and provide commented, minimal, self-contained, reproducible code. > > [[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.
