Can't make spagetti
So, I downloaded the glmmML R package to get p-values for both fixed and random effects but I don't know how to extract said p-values.
I've been using my usual bag of tricks with str(mod), unclass(mod), coef(mod), etc. but, although I can get the fixed effects coefficients. I'd really like to get what's on the Pr(>|z|) column and that p value of 9.567e-10 for the variance component....
id <- factor(rep(1:20, rep(5, 20))) y <- rbinom(100, prob = rep(runif(20), rep(5, 20)), size = 1) x <- rnorm(100) dat <- data.frame(y = y, x = x, id = id) mod<-glmmML(y ~ x, data = dat, cluster = id) >mod Call: glmmML(formula = y ~ x, data = dat, cluster = id) coef se(coef) z Pr(>|z|) (Intercept) -1.3840 0.9609 -1.440 0.150 x 0.6809 0.3459 1.968 0.049 Scale parameter in mixing distribution: 3.436 gaussian Std. Error: 0.6324 LR p-value for H_0: sigma = 0: 9.567e-10 Residual deviance: 96.2 on 97 degrees of freedom AIC: 102.2