計算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

回答