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

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的制剂通过研究者提出并排序的比例的。组合置信区间的上限大于所有研究比例。我知道这个结果没有任何问题，但只需要计算组合置信区间的确切关系来表明结果是正确的。

2016-09-29 M. Er

你是什么意思“的确切关系”的意思是提取它们？ – Wolfgang

我的意思是手动找到置信区间并显示它给出了相同的包装结果。 –

使用''method =“PM”''您正在使用迭代算法。你想用手做这个？ – Wolfgang

您是否使用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]

希望它为你工作

2016-09-29 14:58:07

回答