huang min <minhuangr <at> gmail.com> writes: > I would like to extract the variance components estimation in lme function > like > > a.fit<-lme(distance~age, data=aaa, random=~day/subject) >
Try VarCor(a.fit) Using lme for the type of problems you have is just fine; in many respects, lme4 currently is a welcome complement rather than a replacement of nlme/lme. Dieter ______________________________________________ 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.