1

この質問に対する答えを見つけることができませんでしたが、難しくないことを願っています。

私はメタアナリシスを実施しており、調査をフォレスト プロットに含めたいと考えていますが、効果サイズのメタアナリティック推定から除外したいと考えています。つまり、比較目的でその研究を含めたいのですが、実際のメタ分析には寄与したくありません.

これは、同じプロトコルを使用した約 30 の研究がある奇妙なケースです。それらはすべて、元の研究を直接複製したものです。比較目的で元の研究をフォレスト プロットに含めたいと考えていますが、約 30 の新しい研究のメタ分析には含めたくありません。欠損データ ( meta のフォレスト プロットの欠損変数を除く) を処理できることがわかりましたが、フォレスト プロットにデータを含めたい場合に、モデル自体から除外するデータの行を指定する方法がわかりません。 . 助言がありますか？

4

1 に答える 1

1

もう少し手作業が必要ですが、最初にforest()関数を使用してフォレスト プロットを作成し (すべてのスタディを表示し、下部に要約推定値用のスペースを残します)、次にモデルを適合させます (たとえば、rma()関数を使用して) すべてを除外します。分析に含めたくないスタディを選択し、addpoly()関数を使用して要約推定値を追加します。

次に例を示します。

library(metafor)

### load BCG vaccine data
data(dat.bcg)

### calculate log relative risks and corresponding sampling variances
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, slab=paste(dat$author, dat$year, sep=","))
dat

### forest plot of all studies
forest(dat$yi, dat$vi, xlim=c(-9.5, 7.5), ylim=c(-1.5,16), atransf=exp, at=log(c(.05, .25, 1, 4, 20)))
text(-9.5, 15, "Author(s) and Year",    pos=4)
text( 7.5, 15, "Observed RR [95% CI]",  pos=2)

### random-effects model (without Aronson, 1948)
res <- rma(yi, vi, data=dat, subset=author!="Aronson")
res

### add summary estimate to the bottom
addpoly(res, atransf=exp, row=-1, mlab="RE Model (w/out Aronson, 1948)")

### horizontal separation line
abline(h=0)
于 2014-05-28T15:53:33.583 に答える