在漏斗图上绘制两组数据点（metafor for R）

Question

我正在进行流行病学数据的荟萃分析，其中研究可以分为两组（例如通过方法学）。有没有办法改变漏斗图上绘制的点来识别组（例如点和十字）？ （就像使用trimfill时的pch.fill一样） 感谢。

Answer 1

对于每项研究，应该有一个指标显示它是否是方法B的方法A，对吧？假设变量被称为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)