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

Reply via email to