5

私は疫学データのメタアナリシスを行っており、研究を2つのグループに分けることができます(方法論など)。ファンネルプロットにプロットされたポイントを変更して、グループ(ドットやクロスなど)を識別する方法はありますか?
(trimfillが使用される場合のpch.fillと同じように)ありがとう。

4

1 に答える 1

3

研究ごとに、方法Aか方法Bかを示す指標があるはずですよね?変数の名前が StudyType (メソッド A の場合は = 16、メソッド B の場合は = 17) であるとすると、"pch=StudyType" を指定することで別の記号でプロットできます。プロットでは、TypeA に pch=16、TypeB に pch=17 を使用します。

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

### Attach study type (I just randomly made up some here):
StudyType = sample(c(16,17),13,replace=T)

### meta-analysis of the log relative risks using a random-effects model
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg,
           data=dat.bcg, measure="RR", method="REML")

### standard funnel plot (Notice the pcd command)
funnel(res, pch=StudyType)

ここに画像の説明を入力

于 2013-03-24T14:51:08.277 に答える