如何在R语言双Y轴图中为cNDVI添加标准差误差棒？

问题根源 #

• geom_errorbar中color=target未放入aes()映射内，导致ggplot无法识别该变量，触发报错。
• 数据集sd_NDVI字段存在NA值，这些行的误差棒会被自动跳过，导致部分误差棒不显示。
• 未配置双Y轴逻辑，且误差棒的分组位置未与点/线对齐，即使不报错也会显示错位。

修正步骤及完整代码 #

核心修改点 #

1. 将颜色映射移入aes()：把color=target放到geom_errorbar的aes()参数内，让ggplot按分组变量渲染颜色。
2. 处理NA值：通过na.rm=TRUE忽略缺失值，或提前过滤掉sd_NDVI为NA的行。
3. 对齐分组位置：所有与target分组相关的图层（点、线、误差棒）使用相同的position_dodge偏移量，保证同日期下的分组元素对齐。
4. 配置双Y轴：通过scale_y_continuous设置主轴（温度）和次轴（原始NDVI值），完成数值转换的反向映射。

修正后的代码：

# 可选：提前过滤sd_NDVI为NA的行，或在geom_errorbar中用na.rm=TRUE
NDVI_clean <- NDVI[!is.na(NDVI$sd_NDVI), ]

plot <- ggplot(NDVI_clean) +
  # 温度柱状图
  geom_bar(aes(x = date, y = TempAmb_Avg), 
           stat = "identity", colour = "blue", fill = "blue", width = 0.8) +
  # NDVI点图：设置分组偏移
  geom_point(aes(x = date, y = 30*cNDVI, color = target), 
             position = position_dodge(width = 0.8)) +
  # NDVI线图：与点图使用相同偏移量
  geom_line(aes(x = date, y = 30*cNDVI, color = target), 
            position = position_dodge(width = 0.8)) +
  # 误差棒：颜色映射入aes，对齐分组位置，忽略NA
  geom_errorbar(aes(x = date, ymin = 30*cNDVI - 30*sd_NDVI, 
                    ymax = 30*cNDVI + 30*sd_NDVI, color = target),
                width = 0.2, position = position_dodge(width = 0.8),
                na.rm = TRUE) +
  # 配置双Y轴
  scale_y_continuous(
    name = "TºC",
    sec.axis = sec_axis(~./30, name = "cNDVI")
  ) +
  labs(x = "TIME") +
  theme_bw()

plot

效果说明 #

• 解决object 'target' not found报错，误差棒按target分组显示对应颜色。
• 同日期下的点、线、误差棒对齐，双Y轴分别显示温度和原始NDVI值，更符合直观阅读逻辑。
• 自动忽略sd_NDVI缺失的行，保证有效数据的误差棒正常渲染。

内容的提问来源于stack exchange，提问作者Cláudio Siva