[R] Retrieving p-values and z values from lmer output

Rick Bilonick rab45+ at pitt.edu
Mon Aug 21 23:10:42 CEST 2006


On Mon, 2006-08-21 at 16:52 -0400, Rick Bilonick wrote:
> I can't find a way to retrieve z values and p-values from the output
> from lmer in the lme4 package. How is this done?
> 
> Rick B.
> 
> ______________________________________________
> R-help at 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
> and provide commented, minimal, self-contained, reproducible code.
> 

I was able to do it this way:

> slotNames(summary(fit.lmer))
 [1] "isG"       "methTitle" "logLik"    "ngrps"     "sigma"     "coefs"
 [7] "vcov"      "REmat"     "AICtab"    "flist"     "Zt"        "X"
[13] "y"         "wts"       "wrkres"    "method"    "useScale"
"family"
[19] "call"      "cnames"    "nc"        "Gp"        "XtX"       "ZtZ"
[25] "ZtX"       "Zty"       "Xty"       "Omega"     "L"         "RZX"
[31] "RXX"       "rZy"       "rXy"       "devComp"   "deviance"  "fixef"
[37] "ranef"     "RZXinv"    "bVar"      "gradComp"  "status"

> slot(summary(fit.lmer),"coefs")
               Estimate Std. Error    z value     Pr(>|z|)
(Intercept) -0.02626325 0.03114850 -0.8431628 3.991374e-01
timeAfter    0.42322569 0.02870844 14.7422023 3.453370e-49

> slot(summary(fit.lmer),"coefs")[,4]
 (Intercept)    timeAfter
3.991374e-01 3.453370e-49

Rick B.



More information about the R-help mailing list