[R] VarCorr
Stephanie Cooke
cooke.stephanie at gmail.com
Wed Jan 19 20:32:58 CET 2011
I have a loop that I would like to use to extract the "stddev" for
each itteration so I can average the "stddev" for all the runs. It
would be helpful to know how to extract the "stddev" for each run from
the VarCorr. Thanks
MCruns<-1000
sighatlvec<-rep(NA,MCruns)
sighatbvec<-rep(NA,MCruns)
sighatevec<-rep(NA,MCruns)
for(mc in 1:MCruns)
{
samples<-(sqrt(sigbsq)*rnorm(6))
labs<-(sqrt(siglsq)*rnorm(2))
errors<-(sqrt(sigesq)*rnorm(24))
y<-errors+samples[sampfac]+labs[labfac]
lmer<-lmer(y~1+(1|sampfac)+(1|labfac))
print(VC<-VarCorr(lmer))
}
More information about the R-help
mailing list