[R] extracting p-values from lmer()

From: toka tokas <tokkass_at_yahoo.com>
Date: Tue 06 Dec 2005 - 02:05:41 EST


Dear R users,

I've been struggling with the following problem: I want to extract the Wald p-value from an lmer() fit, i.e., consider

library(lme4)
n <- 120
x1 <- runif(n, -4, 4)
x2 <- sample(0:1, n, TRUE)
z <- rnorm(n)
id <- 1:n
N <- sample(20:200, n, TRUE)
y <- rbinom(n, N, plogis(0.1 + 0.2 * x1 - 0.5 * x2 + 1.5 * z))

m1 <- lmer(cbind(y, N - y) ~ x1 + x2 + (1 | id), family = binomial, method = "AGQ") m1

how to extract the p-value for 'x2' from object m1?

Thanks in advance for any hint,
tokas                 


Just $16.99/mo. or less.
dsl.yahoo.com



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 Tue Dec 06 02:14:44 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:41:28 EST