Auston wrote:
Hi list,Would someone please give me a hint how to extract the between and within subject StdDev from "lme" fit? Below is the output from the fit, but I couldn't find the number 0.3927714 anywhere in the fitted object. Thanks in advance for your help! Random effects: Formula: ~1 | group (Intercept) Residual StdDev: 0.3927714 0.7857085 Best, Auston ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Hi Auston, Extracting those standard deviation values isn't so hard. If you want to extract the variance on your parameter estimates, use diag(lme.output$varFix) ; If you want to extract the variance for the random effects, use as.numeric(VarCorr(lme.output)[,"Variance"]) ; At last, if you want your within-group error standard deviation, simply try lme.output$sigma ; Is there something missing? Cheers, -- *Luc Villandré* /Biostatistician McGill University Health Center - Montreal Children's Hospital Research Institute/ ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
