计算元包中各比例的组合置信区间

Question

考虑以下代码:

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的提法是研究者提出的，是一种比例。组合置信区间有一个上限，该上限大于所有研究的比例。我知道这个结果没有错，只是需要计算组合置信区间的确切关系来显示结果是正确的。

Answer 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]

希望它对你有用