[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