[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