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@(protected)
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html