使用 rma.glmm 计算估计值：当 2 个估计值相同时出错答案

【问题标题】：Computing estimates with rma.glmm: Error when the 2 estimates are identical使用 rma.glmm 计算估计值：当 2 个估计值相同时出错
【发布时间】：2017-12-09 00:13:52
【问题描述】：

我正在使用 rma.glmm 对来自 2 个不同研究的元比例进行分析。例如，在 6 个月时出现 (xi) 和没有 (xii) 不良事件的个体比例：

library(metafor)

#6 months
study=c("1", "2")
ni = c(41, 19)
xi = c(26, 14)
xii = c(15, 5)
NP6monthAT <- data.frame(study, xi, xii, ni)
res2 <- rma.glmm(measure="PLO", xi = xi, ni = ni, data = NP6monthAT)
predict (res2, transf = transf.ilogit, digits=2)
res2 <- rma.glmm(measure = "PLO", xi = xii, ni = ni, data = NP6monthAT)
predict (res2, transf = transf.ilogit, digits=2)

但是，偶然在特定时间点，两个不同人群中的两个比例相同 (11.1%)，我收到一个错误：

#12 months
study=c("1", "2")
ni=c(81, 45)
xi=c(9, 5)
xii=c(72, 40)
NNPNNP12monthAT<-data.frame(study, xi, xii, ni)
res4<-rma.glmm(measure="PLO", xi=xi, ni=ni, data=NNPNNP12monthAT)
predict(res4, transf=transf.ilogit, digits=2)

rma.glmm 中的错误（测量 =“PLO”，xi = xi，ni = ni，数据 = NNPNNP12monthAT）：无法拟合 ML 模型。

我知道估计值将等于 11.1（因为这在两个总体中都是如此）...但我想要具有置信区间的输出，关于我可以做些什么来获取这些信息有什么建议吗？

【问题讨论】：

如果您发布一个可重现的示例，您更有可能获得帮助，人们可以在其中运行您的代码并重现您的错误。
@jmuhlenkamp 我确实包含了代码，它在上面列出。它可以运行并产生错误。
您需要在代码中定义 study 以使其可重现。
好的，研究已定义。

标签： r metafor

【解决方案1】：

该问题源于用于拟合此模型的lme4::glmer()：

> glmer(cbind(xi,ni-xi) ~ 1 + (1 | study), family=binomial, data=NNPNNP12monthAT)
Error: Response is constant

当不同研究的对数几率相同时，显然不可能有任何研究间的异质性。一个简单的解决方案是然后拟合 FE 模型（因为 tau^2 等于 0，无论如何这就是你会得到的）。所以，使用：

res4 <- rma.glmm(measure="PLO", xi=xi, ni=ni, data=NNPNNP12monthAT, method="FE")

【讨论】：