用R语言绘制ARIMA拟合模型与原始序列的图形

Question

我曾使用GRETL。在那里，当我对ARIMA模型进行验证的预测时，我会得到蓝色线条表示的拟合系列和红色线条表示的原始系列。后来，我转向了R，但是我找不到任何命令可以执行同样的操作。我正在使用forecast软件包中的Arima模型。

具体细节如下：

在GRETL中，我通常会执行model->time series -> arima -> forecast操作。它会自动打印出拟合和原始系列。有没有办法在R中执行相同的操作呢？

- Pankaj Parag

Answer 1

这个问题非常开放，但下面是一个非常非常基础的答案。直接从forecast包中Arima的帮助文件示例开始：

fit <- Arima(WWWusage,order=c(3,1,0))

你说你想要原始系列为红色，拟合系列为蓝色：

plot(fit$x,col="red")
lines(fitted(fit),col="blue")

生成的图表大致如下：

- joran

---

最后一条评论有些不公平——这并不是非常基础的问题。你的答案几乎不能在介绍中找到，而ARIMA模型对象的结构相当复杂。 - naught101

@naught101 或许吧。另一方面，我成功地回答了这个问题，而且我对 Arima 模型基本上一无所知，以前甚至从未安装过 forecast 包。 - joran

也许是这样，但Arima应该是Arima，fitted（fit）返回NULL，并且fit $ x不存在。 - user1244215

@user1244215，我不知道你在说什么。这段代码运行起来就像我描述的那样。（我刚刚重新检查了一下。） - joran

@user1244215，你需要加载预测库。 - scai

@joran，你能在这个图表中添加预测线吗？ - ali srn

Answer 2

您发现的差异是因为 arima() 函数与 Arima() 函数不同。前者包含在基础的 stats 包中，而后者来自于 forecast 包，并包括 fitted() 函数以预测观察值。

- paul