[R] extract value from mer object ?

From: Maas James Dr (MED) <J.Maas_at_uea.ac.uk>
Date: Tue, 17 May 2011 13:20:16 +0100

What is the easiest way to extract a value from a 'mer' object from glmer? The first I need is the trtpair Std.Dev. which in this case is 0.17542?

I've managed to get the fixed effects numbers from summary(fednmaout)@coeffs

but no luck with the Random effects stuff.

TIA Jim

> fednmaout

Linear mixed model fit by REML
Formula: lor ~ as.factor(t2) + as.factor(t3) + (1 | trtpair)

   AIC BIC logLik deviance REMLdev
 97.17 107.6 -43.58 83.12 87.17
Random effects:
 Groups Name Variance Std.Dev.
 trtpair (Intercept) 0.030771 0.17542
 Residual 0.230785 0.48040
Number of obs: 60, groups: trtpair, 3

Fixed effects:

               Estimate Std. Error t value
(Intercept)      0.3752     0.3563   1.053
as.factor(t2)1  -0.1677     0.2909  -0.576
as.factor(t3)1   0.1796     0.2909   0.617

Correlation of Fixed Effects:
            (Intr) a.(2)1
as.fctr(2)1 -0.816       

as.fctr(3)1 -0.816 0.500

Dr. Jim Maas
University of East Anglia

R-help_at_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. Received on Tue 17 May 2011 - 12:26:44 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 17 May 2011 - 16:20:08 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive