在 R 预测包中添加预测线的图例条目

Question

我有一张通过fcast <- forecast(Arima(data$train, order=...))获得的预测对象图：

fcast %>% autoplot(series='Forecast') +
  autolayer(fitted(fit), series='Fitted') + 
  autolayer(data$train, series = 'Train') +
  autolayer(data$test, series='Test') +
  xlab("Observation [days]") +
  ylab("Energy [Watts]") +
  guides(colour=guide_legend(title="Data series"), 
         fill=guide_legend(title="Prediction interval"))

看起来像这样：

如您所见，我在图例中添加了所有内容，除了表示预测点的深蓝色线。

如何将预测点（深蓝色线）添加到图例中？

编辑： 该解决方案基于 Rob Hyndman 的回答：

clrs <- c("blueviolet", "blue", "darkgoldenrod4", "red")

autoplot(fcast) +
  autolayer(fcast$mean, series="Forecast") +
  autolayer(fitted(fit), series='Fitted') + 
  autolayer(data$train, series = 'Train') +
  autolayer(data$test, series='Test') +
  xlab("Observation [days]") +
  ylab("Energy [Watts]") +
  guides(colour=guide_legend(title="Data series"), 
         fill=guide_legend(title="Prediction interval")) +
  scale_color_manual(values=clrs)

Answer 1

只需使用autolayer 再次添加即可：

autoplot(fcast) +
  autolayer(fcast$mean, series="Forecast") +
  autolayer(fitted(fit), series='Fitted') + 
  autolayer(data$train, series = 'Train') +
  autolayer(data$test, series='Test') +
  xlab("Observation [days]") +
  ylab("Energy [Watts]") +
  guides(colour=guide_legend(title="Data series"), 
         fill=guide_legend(title="Prediction interval"))

那么你可能需要玩弄颜色。

Answer 2

我做了以下，但我仍然无法得到传说

 autoplot(arima_kw_Pat)+
      xlab("Day")+
      ylab("Under Treatment")+
      ggtitle("Patients Under Treatment in Kuwait ARIMA")+
      guides(colour=guide_legend(title="Data series"), 
              fill=guide_legend(title="Prediction interval"))+
      scale_color_manual(values=clrs)

相反，我收到以下错误 is_missing(values) 中的错误：找不到对象“clrs”