R语言ggplot2调用geom_forecast()报错 绘制10年预测图问题

• 方法1错误点：直接对整个数据框调用ts()会将year和AvgTMean两列都识别为独立的时间序列值，且未指定时间序列的起始年份、采样频率，导致自动生成的坐标轴和需求不匹配。
• 方法2错误点：geom_forecast()默认不会自动继承ggplot()全局设置的美学映射，也未指定预测步长，因此无法识别AvgTMean字段。

以下两种方案均可直接运行，满足x轴映射year、y轴映射AvgTMean、预测未来10年的要求：

方案1：基于ts对象+autoplot实现（最简便） #

构造ts对象时仅传入气温列，明确指定起始年份和年尺度频率，绘图后手动匹配坐标轴标签即可：

library(tidyverse)
library(forecast)

# 正确构造年尺度时间序列：仅传入气温列，起始年1980，年数据frequency=1
AvgTMaxYearTS <- ts(
  data = AvgTMaxYear$AvgTMean,
  start = min(AvgTMaxYear$year),
  frequency = 1
)

# 绘图，h=10指定预测未来10年
autoplot(AvgTMaxYearTS) +
  geom_forecast(h = 10) +
  # 匹配要求的坐标轴名称
  scale_x_continuous(name = "year") +
  scale_y_continuous(name = "AvgTMean")

方案2：直接基于数据框+ggplot原生语法实现 #

在geom_forecast()中显式指定美学映射，不需要手动构造ts对象：

library(tidyverse)
library(forecast)

ggplot(data = AvgTMaxYear, aes(x = year, y = AvgTMean)) +
  geom_line() +
  # 显式传入映射，指定预测10年
  geom_forecast(mapping = aes(x = year, y = AvgTMean), h = 10)

补充：如果需要自定义预测模型、置信区间样式，可以在geom_forecast()中额外传参，比如传入level = c(0.8, 0.95)调整置信水平，传入method = "ets"指定使用指数平滑模型等。

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