如何在metafor漏斗图X轴显示原始比例并按3个年份分组着色

1. 先确认分组变量的有效性 #

首先检查你的年份分组变量是否为因子类型且包含3个明确水平。如果是字符型变量，先转换为因子并验证：

# 假设年份分组变量名为year_cat
dat$year_cat <- factor(dat$year_cat, levels = c("2010-2015", "2016-2020", "2021-2025")) # 替换为你的实际分组名称
# 检查分组水平和样本分布
levels(dat$year_cat)
table(dat$year_cat)

如果table()输出显示某分组样本量为0，那绘图自然只会显示2种颜色；若3个分组都有数据，继续下一步。

2. 完整绘图代码（含颜色映射与X轴转换） #

metafor绘制比例漏斗图时，通常先对比例做logit转换用于效应量计算，再通过参数将X轴转换回原始比例，同时按分组指定颜色：

library(metafor)

# 第一步：计算logit转换的效应量及方差（xi为事件数，ni为总样本量）
dat <- escalc(measure = "PLO", xi = round(prop * ni), ni = ni, data = dat)

# 定义对应3个分组的颜色
group_cols <- c("blue", "red", "orange")

# 绘制漏斗图
funnel(dat$yi, vi = dat$vi,
       col = group_cols[dat$year_cat], # 按分组因子映射颜色
       pch = 19, # 设置点的形状为实心圆
       xlab = "原始比例",
       atransf = transf.ilogit, # 将X轴刻度从logit转换回原始比例
       at = logit(c(0.1, 0.3, 0.5, 0.7, 0.9)) # 指定要显示的原始比例刻度对应的logit值
)

# 添加图例区分分组
legend("topright", legend = levels(dat$year_cat), 
       col = group_cols, pch = 19, bty = "n")

3. 常见问题排查 #

• 仍只显示2种颜色：检查dat$year_cat的每一行取值，确认是否存在观测被错误归类（比如字符大小写不一致、分组名称拼写错误），可以用unique(dat$year_cat)查看所有唯一值。
• X轴未显示原始比例：确保atransf = transf.ilogit参数正确添加，at参数可根据你的数据比例范围调整（比如如果比例集中在0.2-0.6，就设置at = logit(c(0.2,0.4,0.6))）。

内容的提问来源于stack exchange，提问作者Mohamed Rahouma