在漏斗图上绘制两组数据点(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)

​

​