Multilevel Modeling using glmmPQL
On 3/11/2009 11:29 AM, Howard Alper wrote:
Hi, I'm trying to perform a power simulation for a simple multilevel model, using the function glmmPQL in R version 2.8.1. I want to extract the p-value for the fixed-effects portion of the regression, but I'm having trouble doing that. I can extract the coefficients (summary(fit)$coeff), and the covariance matrix (summary(fit)$varFix), but I can't grab the p-value (or the t-statistic.) Could someone explain how to do this? Please send responses to halper at health.nyc.gov.
library(MASS)
fm1 <- glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID,
family = binomial, data = bacteria)
summary(fm1)$tTable[,5]
Look at str(summary(fm1)) .
Thanks in advance,
Howard
DISCLAIMER:\ Sample Disclaimer added in a VBScript.\ \ ...{{dropped:6}}
______________________________________________ R-help at r-project.org 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.
Chuck Cleland, Ph.D. NDRI, Inc. (www.ndri.org) 71 West 23rd Street, 8th floor New York, NY 10010 tel: (212) 845-4495 (Tu, Th) tel: (732) 512-0171 (M, W, F) fax: (917) 438-0894