自动从MASS :: glmmPQL中选择p值

Question

我想自动从这个模型中选择p值：

otu1 <- glmmPQL(X577110~treat,random=~1|time,family=poisson,data=otu)

因为我将运行1000个这样的模型，所以最好能够自动选择p值。我尝试了summary(otu1)$p-value但它不起作用。

Answer 1

你想要这样的东西吗？这些是我在摘要输出中可以找到的唯一p值......

library("MASS")
ss <- summary(glmmPQL(y ~ trt + week, random = ~ 1 | ID,
              family = binomial, data = bacteria)) 

您可以使用names(ss)，str(ss)在摘要对象中搜索您可能感兴趣的组件。

print(ss$tTable[,"p-value"])
 ##   (Intercept)      trtdrug     trtdrug+         week 
 ##  2.192186e-09 5.395735e-02 2.091910e-01 3.057753e-04