根据模型自动绘制寓言预测的寓言工具图

Question

有没有办法将自动打印与寓言一起使用，而不是通过模型对其进行分面？下面的代码生成了一个很好的小图形，但将预测叠加在一起。

library(tidyverse)
library(tsibble)
library(feasts)
library(fable)
library(fabletools)

tourism_melb <- tourism %>%
    filter(Region == "Melbourne") %>% 
    filter(Purpose == "Business") %>% 
    select(-Region, -State, -Purpose)

fableModels <- tourism_melb %>% 
    model(arima = ARIMA(), tslm = TSLM(Trips ~ trend()), mean = MEAN(window = 4))

forecasts <- fableModels %>% 
    forecast(h = "2 years")

forecasts %>% autoplot(tourism_melb)

它看起来是这样的：

​

我正在寻找更像这样的东西，除了刻面之外，这样我就不必纠结于轴比例和标签之类的东西：

library(gridExtra)

arimaPlot <- forecasts %>% 
    filter(.model == "arima") %>% 
    autoplot(tourism_melb)

tslmPlot <- forecasts %>% 
    filter(.model == "tslm") %>% 
    autoplot(tourism_melb)

grid.arrange(arimaPlot, tslmPlot)

​

​

这似乎是一件很常见的事情，但我知道自动绘制这些东西的目的是非常快速和肮脏的，所以我不知道它是否有那个功能。我浏览了fabletools，但没有找到任何东西。

Answer 1

也许这就是你要找的。当autoplot返回一个ggplot对象时，您可以简单地添加facet_wrap(~.model, ncol = 1)

library(tidyverse)
library(tsibble)
library(feasts)
library(fable)
library(fabletools)

tourism_melb <- tourism %>%
  filter(Region == "Melbourne") %>% 
  filter(Purpose == "Business") %>% 
  select(-Region, -State, -Purpose)

fableModels <- tourism_melb %>% 
  model(arima = ARIMA(), tslm = TSLM(Trips ~ trend()), mean = MEAN(window = 4))

forecasts <- fableModels %>% 
  forecast(h = "2 years")

forecasts %>% autoplot(tourism_melb) + facet_wrap(~.model, ncol = 1)

​

​