[R] lme output

[Dieter Menne]

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