[R] Fwd: Extract fixed effects SE from lmer

From: Douglas Bates <dmbates_at_gmail.com>
Date: Thu 30 Jun 2005 - 09:29:37 EST

I forgot to cc: the list on this reply.

On 6/29/05, La Sorte, Frank A. <lasortef@missouri.edu> wrote:
> Hi,
>
> Does anyone know how to extract fixed effects SE values from generalized linear mixed models estimated using the lmer function in the lme4 library? I searched attributes and structure with no luck.
>

Try vcov(myLmerModel). That returns a variance-covariance matrix for the fixed effects. The usual manipulations are used to get the standard errors. For the binomial and Poisson families you should use the optional argument useScale=FALSE as these families do not have a separately estimated scale parameter in the variance function.

In the "show" method for the summary.lmer class the relevant piece of code is

              corF <- as(as(vcov(object, useScale = useScale), "pdmatrix"),
                         "corrmatrix")

______________________________________________
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 Jun 30 09:34:11 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:33:07 EST