[R] variance component analysis for nested model
Russell Senior
seniorr at aracnet.com
Tue Oct 28 06:48:22 CET 2003
>>>>> "Pascal" == Pascal A Niklaus <Pascal.Niklaus at unibas.ch> writes:
Pascal> lme should do the job (r1,r2,r3 are your random factors):
Pascal> library(nlme) y.lme <- lme(y ~ 1,random = ~ 1 | r1/r2/r3)
Pascal> summary(y.lme)
Pascal> This is equivalent to a call to varcomp in S-Plus
Thanks! This was the clue I needed.
--
Russell Senior ``I have nine fingers; you have ten.''
seniorr at aracnet.com
More information about the R-help
mailing list