我在这里使用预测:原理和实践教科书中关于 VAR 模型的章节中的示例数据来fable
使用包制作预测图。
我希望能够从autoplot
. 例如,在下面的代码中,如何Income
仅选择和显示绘图?我还希望能够访问Income
作为响应变量的模型的预测点值和上限和下限。
当我尝试查看forecast()
调用的结果时,我得到一个“丰富显示对象时出错”,所以我看不到如何对这个对象进行子集化。这可能是因为我在 Deepnote(如 Jupyter 笔记本)上运行此分析,我不确定。但也许在调用之前或之后以某种方式过滤、选择或子集会forecast()
得到我正在寻找的数据和绘图。
library(fpp3)
fit <- us_change %>%
model(
aicc = VAR(vars(Consumption, Income)),
bic = VAR(vars(Consumption, Income), ic = "bic")
)
fit %>%
select(aicc) %>%
forecast() %>%
autoplot(us_change %>% filter(year(Quarter) > 2010))