glmmPQL对象的比值比的置信区间（逻辑回归）

Question

计算glm对象的优势比的置信区间是直截了当的：

mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv")
# set up logistic model
mylogit <- glm(admit ~ gre + gpa + rank, data = mydata, family = "binomial")
# calc odds ratios and confidence intervals
exp(cbind(OR = coef(mylogit), confint(mylogit)))

Dataset and code credit

但是，相同的方法不适用于glmmPQL返回的MASS:glmmPQL()个对象（具有随机效果）。

来自glmmPQL（）的MWE帮助：

library(nlme) # will be loaded automatically if omitted
mylogit <- summary(glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID,
+                 family = binomial, data = bacteria))
exp(cbind(OR = coef(mylogit), confint(mylogit)))

使用摘要输出的系数手动计算置信区间（系数+ 2 *标准错误）对我来说也失败了，因为我在摘要对象中找不到它们（但是当将摘要打印到控制台时它们就在那里） 。手动取出它们不是长期的方法......

这对glmmPQL个对象有效吗？