使用ggplot绘制双轴时序图：如何为不同target正确显示NDVI标准差

1. 为何grass的标准差无法显示？ #

• y轴缩放不匹配：geom_errorbar的y参数使用了原始NDVI值（0-1范围），但点线图用的是100*NDVI（0-100范围），导致误差棒被绘制在左轴（温度轴）的极低位置，肉眼无法察觉。
• sd_NDVI数值过小：grass的sd_NDVI大多在0.002-0.016之间，乘以100后误差范围仅0.2-1.6，对应右轴（NDVI轴，0-1）的长度非常短，容易被忽略。
• 存在NA值：数据集中date=19286的grass行sd_NDVI为NA，该行的误差棒无法生成。
• Dodge位置不匹配：geom_errorbar的position_dodge(width=0.05)和点线图的默认position_dodge()（宽度0.9）不一致，导致误差棒和点的位置错位。

2. 如何正确绘制两个target的标准差？ #

修正代码中的参数错误，优化显示效果：

time %>% 
  ggplot()+
  # 绘制温度柱状图，过滤NA值避免空柱子
  geom_bar(aes(x=date, y=TempAmb_Avg),
           stat="identity", na.rm=TRUE)+
  # NDVI点图，统一缩放为100*NDVI，设置dodge宽度
  geom_point(aes(x=date, y=100*NDVI , color=target), 
             position = position_dodge(width=0.5),
             size=2)+
  # NDVI线图，添加group确保按target分组连接，匹配dodge宽度
  geom_line(aes(x=date, y=100*NDVI, color=target, group=target), 
            position = position_dodge(width=0.5))+
  # 误差棒：统一y轴缩放，匹配dodge宽度，开启na.rm跳过无效值
  geom_errorbar(aes(x = date, 
                    y = 100*NDVI,
                    ymin=100*NDVI - 100*sd_NDVI, 
                    ymax=100*NDVI + 100*sd_NDVI, 
                    width=0.2,
                    color = target
  ), 
  position = position_dodge(width = 0.5),
  na.rm=TRUE) +
  labs(x="TIME", y="TºC") +
  scale_y_continuous(sec.axis=sec_axis(~./100, name="NDVI")) +
  theme_bw()

核心修正点：

• 统一所有NDVI相关图层的y轴缩放为100*NDVI，确保轴对齐
• 同步所有图层的position_dodge宽度，保证点、线、误差棒位置对应
• 开启na.rm=TRUE，跳过NA值不影响其他数据绘制
• 增大误差棒的width参数，让小范围误差棒更易识别

3. 如何确保NDVI与sd_NDVI对应正确的target？ #

• 数据层面：保证原始数据中每一行的target、NDVI、sd_NDVI是一一对应的（当前数据集已满足该条件，每行对应一个target的一组观测值）。
• 绘图层面：在所有涉及NDVI的图层（geom_point、geom_line、geom_errorbar）的aes中绑定color=target或group=target，ggplot会自动按target分组，将对应行的NDVI和sd_NDVI匹配绘制。
• 避免手动拆分数据后拼接，保持原始数据的行对应关系，确保分组逻辑统一。

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