Displaying 1 result from an estimated 1 matches for "yeart".
Did you mean:
year
2004 Jul 06
2
lme: extract variance estimate
...[it's the 0.1590195 in the output below]
>print(fit)
>Linear mixed-effects model fit by maximum likelihood
> Data: datag2
> Log-likelihood: -145.0028
> Fixed: lst1 ~ lst
>(Intercept) lst
> 1.1080629 0.7582595
>
>Random effects:
> Formula: ~1 | yeart
> (Intercept) Residual
>StdDev: 0.1590195 0.3313497
However I have not been able to find that number anywhere in the
model object returned by lme. Can anyone tell me how to pull it
or compute it from a model returned by lme? Ditto for models with
random intercept fitted by glmmPQL...