如何在ggplot分组柱状图上叠加双Y轴的样本量折线图？

问题出在你一开始把fill=Type设成了全局映射，叠加折线时ggplot会默认找Type对应的值，但你的样本量数据n是按Analyte汇总的，没有Type列，所以报错。解决核心是拆分图层的数据映射，单独处理柱状图和折线图，再配置双Y轴。

步骤1：整理样本量数据 #

先把样本量n和对应的Analyte做成单独的数据框，不需要Type字段：

n_df <- data.frame(Analyte = letters[1:10], n = n)

步骤2：计算轴转换系数 #

因为Value和样本量n的数值范围差异大，需要计算一个转换系数，让折线能和柱状图在同一画布显示，右侧轴再转换回原始样本量数值：

max_value <- max(main$Value)
max_n <- max(n_df$n)
scale_factor <- max_value / max_n  # 用Value的最大值除以n的最大值得到转换比例

步骤3：绘制双轴图 #

分开设置柱状图和折线图的数据源与映射，避免全局映射冲突：

ggplot() +
  # 绘制分组柱状图，单独指定数据源和映射
  geom_bar(data = main, aes(x = Analyte, y = Value, fill = Type),
           stat = 'identity', position = position_dodge()) +
  # 绘制样本量折线，用n_df数据源，y值是转换后的数值，group=1确保离散x轴能连成线
  geom_line(data = n_df, aes(x = Analyte, y = n * scale_factor, group = 1),
            color = "darkred", size = 1) +
  # 加上折线的点标记
  geom_point(data = n_df, aes(x = Analyte, y = n * scale_factor),
             color = "darkred", size = 2) +
  # 配置双Y轴：左侧是Value，右侧是转换回原始值的样本量
  scale_y_continuous(
    name = "Value",
    sec.axis = sec_axis(~ . / scale_factor, name = "样本量(n)")
  ) +
  # 调整主题和图例
  theme_minimal() +
  labs(fill = "Type")

关键说明 #

• 不要把fill=Type放在全局ggplot(aes())里，而是只在柱状图的geom_bar中指定，这样折线图层不会被强制要求Type字段。
• 折线的group=1必须加，因为x轴是离散的Analyte，ggplot需要明确分组才能把点连成线。
• 双Y轴只是视觉上的适配，本质是通过转换系数让两个不同量级的数据在同一坐标系显示，转换系数可以根据你的数据范围调整。

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