[R] lme object manipulation

[Richard.Cotton at hsl.gov.uk]

> >  > I have an lme object, say lme_res2,  which was
> >  > generated using the varIdent.  I'm trying to extract
> >  > the double 1.532940 from the object, but I can't
> >  > find it by attributes(lme_res2) or
> >  > attributes(summary(lme_res2)).  How can I pull it
> >  > out (so that I can save it to another variable)? 
> >  > Thanks.

> >  > Variance function:
> >  >  Structure: Different standard deviations per
> >  > stratum
> >  >  Formula: ~1 | race
> >  >  Parameter estimates:
> >  >        0        1
> >  > 1.000000 1.532940

Try:
modelint <- intervals(lme_res2)
modelint$varStruct

The help page for intervals.lme claims this is called varFunc, but I 
suspect this is a typo.

Regards,
Richie.

Mathematical Sciences Unit
HSL


------------------------------------------------------------------------
ATTENTION:

This message contains privileged and confidential inform...{{dropped:20}}