Re: [R] lme output

From: Dieter Menne <dieter.menne_at_menne-biomed.de>
Date: Tue 17 Jan 2006 - 21:08:41 EST

John Logsdon <j.logsdon <at> quantex-research.com> writes:

>
> I am trying to extract the solution from a simple lme calculation.
>
> > y<-c(-1.118,-.5,.5,1.118,10)
> > gp<-factor(c(rep('one',4),'two'))
> > res<-lme(y~1,rand=~1|gp)
...
> Random effects:
> Formula: ~1 | gp
> (Intercept) Residual
> StdDev: 7.026737 0.9999747
>
> The problem is to get the random effects intercept - essentially the
> between groups standard deviation - into a scalar.

VarCorr(res)

will give you the matrix with the relevant numbers



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 Received on Tue Jan 17 22:23:46 2006

This archive was generated by hypermail 2.1.8 : Wed 18 Jan 2006 - 02:10:24 EST