[R] lme4 extracting individual variance components

FabbroThomas thomas.fabbro at unifr.ch
Fri Jun 24 14:47:47 CEST 2005


Hi,
For further calculations I need to extract indivdual Variances of 
different random effects from a fitted model.
I found out how to extract the correlations 
(VarCorr(m1)@reSumry$group1) but I was not able to find a way to 
extract the other components individually.

To extract the Residuals I tried: (ranef(m1)@ stdErr) which 
unfortunately did not work.
Thank you very much for your help!

Thomas


My model:

 > m1<- lmer(y ~ trtt + (trtt-1|group3) + (trtt-1|group2) + 
(trtt-1|group1), d1)
 > m1
Linear mixed-effects model fit by REML
Formula: y ~ trtt + (trtt - 1 | group3) + (trtt - 1 | group2) + (trtt - 
      1 | group1)
    Data: d1
       AIC      BIC    logLik MLdeviance REMLdeviance
  1819.454 2003.915 -874.7269   1736.421     1749.454
Random effects:
  Groups   Name                Variance Std.Dev. Corr
  group1      trtt1/TR1       0.115094 0.33926
           trtt1/TR2 0.338576 0.58187   0.177
           trtt2/TR1       0.141726 0.37647  -0.002 -0.007
           trtt2/TR2 0.327869 0.57260  -0.007 -0.002  0.321
  group2     trtt1/TR1       0.026259 0.16205
           trtt1/TR2 0.021771 0.14755  0.325
           trtt2/TR1       0.025403 0.15938  0.898  0.062
           trtt2/TR2 0.024479 0.15646  0.048  0.012  0.156
  group3   trtt1/TR1       0.028206 0.16795
           trtt1/TR2 0.133718 0.36567  0.848
           trtt2/TR1       0.057741 0.24029  0.990  0.905
           trtt2/TR2 0.110075 0.33178  0.638  0.944  0.731
  Residual                     0.021744 0.14746
# of obs: 1437, groups: group1, 1088; group2, 31; group3, 3

Fixed effects:
                        Estimate  Std. Error   DF t value  Pr(>|t|)
(Intercept)            4.875320    0.102476 1433  47.575 < 2.2e-16 ***
trtt1/TR2   -4.647211    0.149408 1433 -31.104 < 2.2e-16 ***
trtt2/TR1          0.577186    0.052293 1433  11.037 < 2.2e-16 ***
trtt2/TR2   -4.237753    0.166187 1433 -25.500 < 2.2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1



I am working with:

R 2.1.0
lme4 0.96-1
Matrix 0.96-3
lattice 0.11-6
latticeExtra 0.1-3

On Mac OS X 10.3.9




More information about the R-help mailing list