mydat <- data.frame(stage1 = c(4, 3, 2, 1), n1 = c(10, 40, 30, 20))
mydat
stage1 n1
1 4 10
2 3 40
3 2 30
4 1 20
4 つのスタディ (行ごとに 1 つ) の単純なデータ セットがあり、各スタディのイベント数 (ステージ 1) とサンプル サイズ (n1) があります。研究のステージ 1 では、研究 1 の 10 人中 4 人がこの病気にかかっています。比較すると、研究 #2 では 40 人中 3 人しかこの病気にかかっていません。
library(meta)
metaprop(mydat$stage1, mydat$n1)
この関数を使用してmetaprop、4 つの研究から単一比率のメタ分析を実行します。ただし、各研究が後の段階で他の効果の推定値を報告するとします。
mydat2 = data.frame(stage1 = c(4, 3, 2, 1), n1 = c(10, 20, 30, 40), stage2 = c(7, 5, 3, 4), n2 = c(10, 20, 30, 40))
mydat2
stage1 n1 stage2 n2
1 4 10 7 10
2 3 20 5 20
3 2 30 3 30
4 1 40 4 40
ここで、ステージ 2 では、10 人の可能性のある人のうち合計 7 人が研究 1 の病気にかかっています。相関関係が正しく調整されるように、どのように正確に考慮に入れるでしょうか? data.frame でステージを作成し、インジケーターを調査する必要があります。
> mydat3 = data.frame(event = c(mydat2$stage1, mydat2$stage2), n = c(mydat2$n1, mydat2$n2), stage = c(rep(1, 4), rep(2, 4)), study = c(rep(c(1, 2, 3, 4), 2)))
> mydat3
event n stage study
1 4 10 1 1
2 3 20 1 2
3 2 30 1 3
4 1 40 1 4
5 7 10 2 1
6 5 20 2 2
7 3 30 2 3
8 4 40 2 4
metapropしかし、それを関数呼び出しにどのように含めるのでしょうか? 私は他の関数 ( だけでなくmetaprop) を使用することにもオープンです。