web-dev-qa-db-ja.com

Rプロットで、元のシリーズを使用したアリマフィットモデル

GRETLを使っていました。そこで、arimaモデルの検証のために予測を行うと、フィットされたシリーズが青い線で、元のシリーズが赤い線で表示されます。その後、Rに切り替えたところ、同じことを行うコマンドが見つかりませんでした。予測パッケージの有馬モデルを使用しています。

詳細、

GRETLでは、モデル->時系列->有馬->予測に使用します。取り付けられたシリーズと元のシリーズが自動的に印刷されます。 Rで同じことをするアイデアはありますか?

20
Pankaj Parag

この質問はかなり自由回答ですが、これが非常に基本的な答えです。 Arimaパッケージのforecastのヘルプファイルで提供されている例の1つから直接始めます。

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

元のシリーズを赤で、フィットしたシリーズを青で表示するとします。

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

次のようなプロットが作成されます。

arima

このような基本的な質問については、Rの基本的な機能のいくつかに慣れるために、ある程度の時間を費やすことをお勧めします here

40
joran

あなたが見つけている違いは、arima()関数がArima()関数とは異なるためです。前者は基本の{stats}パッケージに含まれ、後者は{forecast}パッケージから提供され、観測値を予測するfit()関数が含まれています。

16
paul