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

Question

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

Answer 1

对于每项研究，都应该有一个指标显示它是方法 A 还是方法 B，对吗？假设变量名为 StudyType（= 16 如果方法 A，= 17 如果方法 B），然后您可以通过指定“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)