Re: [R] Extracting variance components from lmer

From: Douglas Bates <bates_at_stat.wisc.edu>
Date: Thu 23 Feb 2006 - 01:35:11 EST

On 2/21/06, Christoph Buser <buser@stat.math.ethz.ch> wrote:
> Hi Rick
>
> There may be a better way, but the following should work:
>
> attributes(vc.fit)$sc

That works but a more direct way would be

attr(vc.fit, "sc")

By the way, that value is the estimated standard deviation not the estimated variance.

>
> Regards,
>
> Christoph Buser
>
> --------------------------------------------------------------
> Christoph Buser <buser@stat.math.ethz.ch>
> Seminar fuer Statistik, LEO C13
> ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND
> phone: x-41-44-632-4673 fax: 632-1228
> http://stat.ethz.ch/~buser/
> --------------------------------------------------------------
>
>
> Rick DeShon writes:
> > Hi All.
> >
> > I need a bit of help extracting the residual error variance from the VarCorr
> > structure from lmer.
> >
> > #Here's a 2-way random effects model
> > lmer.1 <- lmer(rating ~ (1|person)+(1|rater), data = dat)
> >
> > #Get the structure
> > vc.fit <- VarCorr(lmer.1)
> >
> > #results in.....
> > $person
> > 1 x 1 Matrix of class "dpoMatrix"
> > (Intercept)
> > (Intercept) 0.7755392
> >
> > $rater
> > 1 x 1 Matrix of class "dpoMatrix"
> > (Intercept)
> > (Intercept) 0.2054469
> >
> > attr(,"sc")
> > [1] 0.5051518
> >
> > #I can pull out the person and rater variance components easy enough. For
> > example...
> > vc.person <- vc.fit$person@x
> >
> > I'm sure it's simple but I have not been able to grab the residual variance
> > in the last matrix. I simply wish to asign the residual to a scalar
> > variable. Any suggestions would be appreciated!
> >
> > Thanks!
> >
> > Rick DeShon
> >
> >
> > --
> > Rick DeShon
> > 306 Psychology Building
> > Department of Psychology
> > Michigan State University
> > East Lansing, MI 48824-1116
> >
> > [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Feb 23 02:59:18 2006

This archive was generated by hypermail 2.1.8 : Thu 23 Feb 2006 - 06:23:11 EST