如何在ggplot中绘制多变量时间序列及均值与置信区间

问题出在你把col=variable放到了ggplot的全局美学映射里，这会让所有图层（包括均值线和置信区间）都按照variable分组，导致置信区间被拆成每个变量的区间，而非你需要的跨变量汇总区间。

正确的做法是把变量曲线的美学映射单独放到对应的geom_line里，同时用单独的汇总数据框来绘制均值线和置信区间，具体代码如下：

library(ggplot2)
library(reshape2)

# 1. 生成模拟数据
nr <- 25
nc <- 15
dat <- data.frame(matrix(data=rnorm(nr*nc), nrow=nr, ncol = nc))
year <- seq(2000, 2000+3*(nr-1), by = 3)
dat$year <- year
dat_long <- melt(data=dat, id.vars = "year")

# 2. 计算跨变量的均值、标准差及置信区间上下限，整理成单独的数据框
summary_dat <- data.frame(
  year = year,
  mean = apply(dat[, -ncol(dat)], 1, mean), # 排除year列计算均值
  sd = apply(dat[, -ncol(dat)], 1, sd),
  upper = apply(dat[, -ncol(dat)], 1, function(x) mean(x) + 2*sd(x)),
  lower = apply(dat[, -ncol(dat)], 1, function(x) mean(x) - 2*sd(x))
)

# 3. 绘图：分层绘制，先画置信区间，再均值线，最后变量曲线
ggplot() +
  # 绘制置信区间，用汇总数据框
  geom_ribbon(data = summary_dat, aes(x = year, ymin = lower, ymax = upper), 
              alpha = 0.5, fill = "gray") +
  # 绘制均值线，用汇总数据框
  geom_line(data = summary_dat, aes(x = year, y = mean), 
            color = "black", size = 1) +
  # 绘制各变量的时间曲线，把col=variable放在局部aes里
  geom_line(data = dat_long, aes(x = year, y = value, color = variable), 
            alpha = 0.6) +
  # 可选：调整图例和主题
  labs(x = "年份", y = "数值", color = "变量") +
  theme_minimal()

关键说明： #

• 拆分数据源：汇总数据（均值、区间）和变量长数据分别对应不同的geom，各自指定data参数，避免图层间的数据源干扰。
• 局部映射变量颜色：变量曲线的color = variable仅在geom_line的局部美学映射里设置，不会影响均值线和置信区间的分组逻辑。
• 置信区间用填充色：geom_ribbon用fill设置填充色，和变量曲线的color（线条色）区分开，避免视觉混淆。

这样就能同时展示所有变量的时间趋势，以及每个时间点跨变量的均值和置信区间了。

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