R ggplot2使用geom_ribbon绘图与stat_summary效果不符如何解决

你的代码存在3个核心问题，导致输出图形和原始效果不一致：

• 统计逻辑错误：你先调用summarise()计算了每个(hour, name)分组的均值，后续再计算sd(value)时，每个分组下只有1个均值数值，计算出的标准差全为NA，完全无法生成ribbon需要的区间范围。原始代码中mean_sdl是基于分组下的所有原始观测值同时计算均值和标准差，和你先汇总均值再算标准差的逻辑完全不同。
• 映射缺失：geom_ribbon()没有添加fill = name或group = name的映射，即使区间计算正确，也会把出发延误、到达延误两个类别的区间混叠在一起，无法和点、线的颜色分组对应。
• 语法错误：geom_ribbon()代码行结尾缺少连接图层的+号，导致后续的theme_bw()根本没有被添加到绘图对象中。

library(nycflights13)
library(tidyverse)

# 按分组同时计算均值、标准差，对齐stat_summary的计算逻辑
plot_df <- flights %>%  
  select(hour, dep_delay, arr_delay) %>% 
  filter(hour > 4) %>%
  pivot_longer(!hour) %>% 
  group_by(hour, name) %>%
  summarise(
    delay_mean = mean(value, na.rm = TRUE),
    delay_sd = sd(value, na.rm = TRUE),
    .groups = "drop"
  ) %>% 
  # 按照原始代码mult=0.2的参数计算区间上下限
  mutate(
    ymin = delay_mean - delay_sd * 0.2,
    ymax = delay_mean + delay_sd * 0.2
  )

# 组合geom图层绘图
plot_df %>% 
  ggplot(aes(x = hour, color = name, fill = name)) +
  geom_point(aes(y = delay_mean), size = 3) +
  geom_line(aes(y = delay_mean), size = 1.1) +
  geom_ribbon(aes(ymin = ymin, ymax = ymax), alpha = 0.3) +
  theme_bw()

• 统计步骤对齐：分组汇总时同时基于原始观测计算均值和标准差，从根源上避免单值计算标准差得到NA的问题，和原始stat_summary的计算规则完全一致。
• 全局映射优化：在ggplot初始映射中统一指定color = name、fill = name，让点、线的颜色分组和ribbon的填充分组自动对齐，不需要在每个geom里重复写映射规则。
• 语法补全：补全图层连接的+号，保证主题正常生效。

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