Hi all,
I have the output of summary() of an lme object called "lme.exp1", for
example
#############################################
> summary(lme.exp1)
Linear mixed-effects model fit by REML
Data: DATA
Log-restricted-likelihood: -430.8981
Fixed: fixed.exp1
....
Random effects:
Formula: ~-1 + mu1 + log.sig1 | animID
Structure: Diagonal
mu1 log.sig1 Residual
StdDev: 2.605223 0.1013134 3.876431
Correlation Structure: General
Formula: ~1 | animID/fieldN
Parameter estimate(s):
Correlation:
1
2 0.612
Variance function:
Structure: Different standard deviations per stratum
Formula: ~1 | type
Parameter estimates:
mu1 log.sig1
1.00000000 0.03722675
Number of Observations: 376
Number of Groups: 32
##############################################
If I want to reuse these estimates, I can use
> summary(lme.exp1)$coefficients$fixed;
to extract fixed effect estimates, and
> summary(lme.exp1)$sigma
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
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
