[R] extracting intercepts from lme model summaries
Jarrod Hadfield
jarrod.hadfield at ic.ac.uk
Mon Sep 16 15:44:02 CEST 2002
--
Using the function lme below,
summary(lme(clip~1,random=~1|Block/Nurse.Nest/Orig.Nest, data=ringy1)))
I get the summary:
Linear mixed-effects model fit by REML
Data: ringy1
Log-restricted-likelihood: -760.1605
Fixed: clip ~ 1
(Intercept)
16.8508
Random effects:
Formula: ~1 | Block
(Intercept)
StdDev: 0.1995932
Formula: ~1 | Nurse.Nest %in% Block
(Intercept)
StdDev: 0.00551305
Formula: ~1 | Orig.Nest %in% Nurse.Nest %in% Block
(Intercept) Residual
StdDev: 0.3064555 0.5537723
Number of Observations: 803
Number of Groups:
Block Nurse.Nest %in% Block
51 99
Orig.Nest %in% Nurse.Nest %in% Block
195
I can find an extract most values except the intercept value for the
final formula; in this case 0.30694555. Is there a way of doing this?
Yours,
Jarrod Hadfield
>
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list