在漏斗图上绘制两组数据点(metafor for R)

Question

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

Answer 1

对于每项研究，都应该有一个指标来显示它是方法 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)