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