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

のために働くホープ

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

でそれらを抽出することができますが、「正確な関係」とはどういう意味ですか？ – Wolfgang

私は信頼区間を手で見つけて、それがパッケージの結果と同じであることを示しています。 –

'' method = "PM" "'で、あなたは反復アルゴリズムを使用しています。あなたは手でこれをしたいですか？ – Wolfgang

PropCIsパッケージで試しましたか？

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

このように表示されます。しかし、あなたはそれがあなた

2016-09-29 14:58:07

答えて