[R] [R-lme] Extract estimated variances from output of lme?

[Ben Bolker]

Tingting Zhan <tingting.zhan <at> jefferson.edu> writes:

> 
> Hi all,
> 
> I have the output of summary() of an lme object called "lme.exp1", for
> example
> 
> #############################################
> > summary(lme.exp1)
>

  [snip]

> for the common variance parameter sigma.   But if I need the covariance
> estimates 0.612 and 0.0372, do I have a function to extract these numbers
> too?   I'm looking for something similar to function ranef() [of course
> ranef() does not serve my purpose].  Any advice appreciated!
> 
> Tingting


  ?VarCorr

and

str(VarCorr(lme.exp1))

to figure out how to extract the values you want (it's
a slightly odd object, a matrix of mode 'character' with
labels interspersed with numeric values ...)

  r-sig-mixed-models is probably the best list for this
sort of query ...