计算 metafor 包中比例的组合置信区间？答案

【问题标题】：Calculate the combined confidence interval of proportions in metafor package?计算 metafor 包中比例的组合置信区间？
【发布时间】：2016-09-29 23:32:09
【问题描述】：

考虑下面的代码：

library(metafor)
dat <- escalc(measure="PR", xi=dpos+cpos, ni=dpos+dneg+cpos+cneg, data=data)
res <- rma(yi, vi, data = dat, method="PM")

xi 和 ni 的公式是研究人员提出的，是一种比例。组合置信区间的上限大于所有研究比例。我知道这个结果没有问题，但只需要计算组合置信区间的确切关系即可表明结果正确。

【问题讨论】：

“确切的关系”是什么意思？
我的意思是手动找到置信区间并表明它给出了与包相同的结果。
使用method="PM"，您正在使用迭代算法。你想手动做吗？

标签： r

【解决方案1】：

您是否尝试过使用 PropCIs 包？

library(PropCIs)
Exact=exactci(x=data$dpos+data$cpos,n=data$dpos+data$dneg+data$cpos+data$cneg)

你可以这样看。但是你可以用

提取它们

Inferior=Exact[[1]][1]
Superior=Exact[[1]][2]

希望对你有用

【讨论】：