使用ggplot的geom_ribbon绘制含分位数的时间序列带状图问题

Hey James，别慌！我来帮你搞定这个时序数据的百分位带状图问题。很多时候用geom_ribbon卡壳，都是数据格式或者映射没搞对，咱们一步步来：

第一步：确认你的数据格式 #

geom_ribbon需要宽格式数据——也就是每一行对应一个时间点，同时包含第5、50、95百分位的三列数据（比如命名为p5、p50、p95）。如果你的数据是长格式（比如一行是一个时间+一个百分位的数值），得先转成宽格式。

先给你个示例数据参考：

# 构造模拟的时序百分位数据
set.seed(123)
df <- data.frame(
  time = seq.Date(as.Date("2023-01-01"), as.Date("2023-12-31"), by = "day"),
  p5 = rnorm(365, mean = 10, sd = 2),   # 第5百分位
  p50 = rnorm(365, mean = 15, sd = 2),  # 第50百分位（中位数）
  p95 = rnorm(365, mean = 20, sd = 2)   # 第95百分位
)

如果你的数据是长格式，用tidyr的pivot_wider转成宽格式就行：

library(tidyr)

# 假设长格式数据是这样的：time, percentile, value
long_df <- df %>%
  pivot_longer(cols = c(p5, p50, p95), names_to = "percentile", values_from = "value")

# 转成宽格式
wide_df <- long_df %>%
  pivot_wider(names_from = "percentile", values_from = "value")

第二步：用geom_ribbon绘制带状图 #

数据格式没问题的话，直接调用geom_ribbon，把ymin设为第5百分位，ymax设为第95百分位，再叠加geom_line画出中间的第50百分位线：

library(ggplot2)

ggplot(wide_df, aes(x = time)) +
  # 绘制5%-95%的带状区域，用半透明填充
  geom_ribbon(aes(ymin = p5, ymax = p95), fill = "#69b3a2", alpha = 0.3) +
  # 绘制第50百分位的中线，加粗突出
  geom_line(aes(y = p50), color = "#333333", linewidth = 1) +
  # 美化图表标签和主题
  labs(
    x = "时间",
    y = "数值",
    title = "时序数据的5th-95th百分位带状图",
    subtitle = "中线为第50百分位"
  ) +
  theme_minimal()

常见坑点排查 #

• 时间轴格式错误：如果你的time列是字符串类型，一定要转成日期格式（as.Date()）或者时间格式（as.POSIXct()），不然ggplot会把它当成离散变量，图表会乱掉。
• 映射参数放错位置：ymin和ymax要放在aes()里面，这样才会和x轴的时间一一对应；如果放在aes()外面，只会生成固定的带状区域。
• 数据缺失值：如果某行的p5或p95有NA，geom_ribbon会跳过那个时间点，你可以用na.omit()清理数据，或者设置na.rm = TRUE。

这样应该就能顺利画出你想要的时序带状图啦！

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