我制作了一个预测时间序列的情节。它是通过以下代码实现的:
forecast1 <- HoltWinters(ts, beta = FALSE, gamma = TRUE)
forecast2 <- forecast(forecast1, h = 60)
autoplot(forecast2)
哪里&#39; ts&#39;是一个时间序列对象。 所以我想添加另一个时间序列来比较预测值和实际值,从我上次的实际观察开始。我用经典的情节实现了它,添加了一个实际时间序列的线。这是我的情节:
如何将这一新行添加到我的第一张图中?
答案 0 :(得分:5)
这是最简单的方法:
library(ggplot2)
library(forecast)
smpl1 <- window(AirPassengers, end = c(1952, 12))
smpl2 <- window(AirPassengers, start = c(1953, 1), end = c(1953,12))
hw <- HoltWinters(smpl1, beta = FALSE, gamma = TRUE)
forecast <- forecast(hw, h = 12)
autoplot(forecast) +
autolayer(smpl2, series="Data") +
autolayer(forecast$mean, series="Forecasts")
预测包中的autolayer
命令允许您将涉及时间序列和预测的图层添加到现有图中。