考虑下面的守则:
library(metafor)
dat <- escalc(measure="PR", xi=dpos+cpos, ni=dpos+dneg+cpos+cneg, data=data)
res <- rma(yi, vi, data = dat, method="PM")习和ni的配方是由研究者提出的,并且有一定的比例。综合置信区间有一个上限,大于所有研究的比例。我知道这个结果没有什么问题,只是需要精确的关系来计算综合置信区间,以证明结果是正确的。
发布于 2016-09-29 14:58:07
你试过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]希望它对你有用
https://stackoverflow.com/questions/39762993
复制相似问题