在漏斗图上绘制两组数据点（使用R中的metafor）

Question

我正在对流行病学数据进行元分析，其中研究可以按两组（例如按方法）分开。有没有办法改变漏斗图上绘制的点以识别组（例如点和叉号）？
(就像使用trimfill时pch.fill所做的那样) 谢谢。

- Cullen

---

你是在谈论funnel()命令吗？它接受pch选项。请参阅技术文档的第81页。 - Penguin_Knight

也许我没有表达清楚。我想在同一个漏斗图上绘制两个不同的点组（来自同一项元分析）。每个组代表分析中的一个子组。如何更改绘制的点以不同方式表示这些组？pch似乎会更改所有绘制的点，谢谢。 - Cullen

每个研究都应该有一个指标，显示它是否采用A方法或B方法，对吧？假设变量叫做StudyType（如果是A方法则=1，如果是B方法则=2），然后您可以通过指定“pch = StudyType”来使用不同的符号进行绘图。对于TypeA，图表将使用pch = 1，而对于TypeB，图表将使用pch = 2。 - Penguin_Knight

Answer 1

每项研究都应该有一个显示它是方法A还是方法B的指标，对吧？假设变量名为StudyType（如果采用方法A，则=16，如果采用方法B，则=17），您可以通过指定“pch=StudyType”来使用不同的符号绘制它们。绘图将使用pch = 16代表method A，pch = 17代表method B。

### 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)

- Penguin_Knight