如何使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)

这是我在上面看到的错误：

自动绘图(stats：：d.var( n.ahead = 50)，conf.int = TRUE) get.dtindex.continuous中的错误(模型$y，length = nrow(fcst[1]) )：未能转换ts对象索引

我如何让autoplot.varprd与动物园对象一起工作并显示日期？

如果我从一个动物园对象开始，并将它转换为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。