r - 如何使 autoplot.varprd 与动物园对象一起使用

Question

library(vars)
library(ggfortify)
library(zoo)
data(Canada, package = 'vars')
Canada <- zoo(Canada)
#autoplot.zoo works fine with zoo objects
autoplot(Canada)
d.var <- vars::VAR(Canada, p = 3, type = 'const')
#however autoplot.varprd does not work with zoo objects
autoplot(stats::predict(d.var, n.ahead = 50), conf.int = TRUE)

这是我上面得到的错误：

autoplot(stats::predict(d.var, n.ahead = 50), conf.int = TRUE) get.dtindex.continuous(model$model$y, length = nrow(fcst[[1]]) 中的错误， : 无法将 ts 对象索引转换为日期

如何使 autoplot.varprd 与动物园对象一起使用并显示日期？

如果我从 zoo 对象开始并将其转换为 ts，那么 autoplot.varprd 将起作用，但会显示数字而不是日期。

谢谢。

Answer 1

我们可以通过将预测转换为 zooreg 对象然后调用来避免使用 ggfortify autoplot.zoo：

library(ggplot2)

p <- predict(d.var, n.ahead = 50)
freq <- frequency(Canada)
z <- zooreg(data.frame(lapply(p$fcst, "[", TRUE, 1)), end(Canada) + 1/freq, freq = freq)
autoplot.zoo(z) + scale_x_yearqtr() + facet_free()

请注意，如果未加载 ggfortify，则autoplot.zoo可以选择将其缩短为autoplot.