-- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._