I am using the following code to get estimates of the variances of the levels of factor 1 (factor1).
lme(sp1~factor1-1,data=z.a.1,random=pdBlocked(list(pdIdent(~factro1-1),pdIdent(~factor2-1))), weights=varIdent(~1|factor1)) In this study factor 1 is crossed with factor 2. is this correct? [[alternative HTML version deleted]] ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html