如何为R中时间序列分解图的X轴添加日期或年份？

你的问题本质是原始时间序列对象没有绑定日期索引，导致autoplot默认用观测序号作为X轴标签。以下是具体修正步骤：

1. 先确保时间序列携带日期信息 #

如果你的原始数据是普通数据框，先把日期列转为Date类型，再创建带周期的时间序列对象（日度数据适合用msts，支持年周期365天）：

library(forecast)
library(ggplot2)
library(lubridate) # 处理日期的工具包

# 假设你的数据框为df，包含date（日期列）和load（电力负荷列）
df$date <- as.Date(df$date)
# 创建多周期时间序列，指定年周期365天，起始时间对应第一个观测的年和年内天数
load_ts <- msts(df$load, seasonal.periods = 365, start = c(year(df$date[1]), yday(df$date[1])))

2. 分解后修改autoplot的X轴刻度 #

完成mstl分解后，直接通过scale_x_continuous将序号映射为日期/年份：

# 执行时间序列分解
decomp_result <- mstl(load_ts)

# 获取完整的日期序列
full_dates <- df$date

# 绘制分解图并替换X轴
autoplot(decomp_result) +
  scale_x_continuous(
    # 按年设置刻度（每365个观测一个刻度）
    breaks = seq(1, length(full_dates), by = 365),
    # 用年份作为标签，若要显示具体日期，换成format(full_dates[x], "%Y-%m-%d")
    labels = function(x) year(full_dates[x])
  ) +
  labs(x = "年份") +
  theme_minimal()

3. 备选方案：转成tibble后用ggplot直接绘图 #

如果觉得上面的方法麻烦，可将分解结果转为tibble，直接绑定日期列绘图，更灵活：

library(broom)
library(dplyr)

# 将分解结果转为整洁格式，并添加日期列
decomp_tidy <- tidy(decomp_result) %>%
  mutate(date = rep(full_dates, nrow(decomp_result)/length(full_dates)))

# 用ggplot绘制，直接以date为X轴
ggplot(decomp_tidy, aes(x = date, y = value)) +
  geom_line(linewidth = 0.8) +
  facet_wrap(~component, scales = "free_y", ncol = 1) +
  labs(x = "日期", y = "数值") +
  theme_bw()

内容的提问来源于stack exchange，提问作者Saïd Maanan