如何使autoplot.varprd与zoo对象一起使用

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对象然后调用autoplot.zoo来避免使用ggfortify：

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。