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

【问题标题】：Plotting two sets of data points on funnel plot (metafor for R)在漏斗图上绘制两组数据点（R 的 metafor）
【发布时间】：2013-03-24 22:51:08
【问题描述】：

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

【问题讨论】：

你说的是 funnel() 命令吗？它接受 pch 选项。请参阅technical document 的第 81 页。
也许我没有像我应该的那样清楚。我想在同一个漏斗（来自同一个荟萃分析）上绘制两组不同的点。每个组代表分析中的一个子组。如何更改绘制的点以不同地代表组？ pch 似乎会改变所有绘制的点，谢谢
对于每一项研究，都应该有一个指标显示是方法A还是方法B，对吧？假设变量名为 StudyType（= 1 如果方法 A，= 2 如果方法 B），然后您可以通过指定“pch=StudyType”用不同的符号绘制它们。该图将对 TypeA 使用 pch=1，对 TypeB 使用 pch=2。

标签： r

【解决方案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)

【讨论】：

谢谢——这很有用