AICの導出について
- ページ: The R Book ご意見
- 投稿者: 松本
- 優先順位: 重要
- 状態: 提案
- カテゴリー: 疑問
- 投稿日: 2005-05-14 (土) 13:03:32
メッセージ
先日は、decomposeについての質問にお答え頂きありがとうございました。
AICの導出プログラムについてお聞きしたいのですが、
The R Book第5章にあるAICの導出プログラムで、
order=c(j,1,0)とあるので、ARIMA(j,1,0)の次数決定に使えるとは思うのですが、ARIMA(j,k,0)の次数決定にはどのようにすればよいのでしょうか?
for(j in 1:15,k in 1:15)のように直すのかと思いためしましたが、エラーが出たので、どのようにすれば良いのか教えていただけませんか?
- forをつかった二重ループにすれば計算できます。あとでプログラム例を載せておきます。 -- 矢野
- それと第5章を見直していて、ARIMAモデルに関する記述(式(5.12)に記述漏れがあります)に間違いを見つけてしまいました。早急に訂正文をthe R bookのサポートページに掲載しておきます(このページでもお知らせします)。お買い上げいただいた皆様、本当に申し訳ありませんでした。 -- 矢野
- ご要望のようにするのは以下のようにfor文を二重に使用します。-- 矢野
for (j in 1:5){
for(d in 1:5){
tmp <- arima(lipTs, order=c(j,d,0))
cat("ORDER", j, "DIFF", d, "AIC=", tmp$aic,"\n")
}
}
- ひとつ注意点があります。たとえば、1階差分をとったARIMA(q,1,0)のAICと2階差分をとったARIMA(q,2,0)のAICを比較するのは得策ではありません。というのは1階差分をとった時系列と2階差分をとった時系列では想定される真のモデルが異なるため、それぞれの場合のAICの意味が異なるからです。そして、そのことを補足せねばならないなと思ってthe R bookを見直して自分の愚かな間違いに(1年も経ってから)気がつきました。お買い上げいただいた皆様、本当に申し訳ありません。心よりお詫び申し上げます。また、共著者の方々と出版に携わってくださった皆様にも本当に申し訳ありません。正誤表へのリンク:The R book正誤表 -- 矢野
- 矢野先生、本当にありがとうございます。 -- 松本
- Rを始めたばかりで、今、為替レートの時系列分析をしようとして奮闘していますが、分からない点がやはりあって、できるだけ自分で調べようとしているのですが、どうしてもGoogle等では、Rという名前のせいか、探しづらく、矢野先生にはお世話になりっぱなしで、本当に感謝しています。 -- 松本
- ご指摘いただいた注意点ですが、何回差分をとるかは、AICではなく、違う方法にします。私が使うデータは、そのままでは定常でないですが、一階差分をとれば、定常性が確認できたので、ARIMA(p,1,q)の次数決定を教えていただいたようにしようと思います。教えていただいたプログラムの、orderをc(j,1,d)とすれば良いと思いますので。 -- 松本
- あ、catの所のDIFFをなにか違うORDER2等に変える必要もありそうですね。 -- 松本
- 先生じゃないので先生は止めてください。 -- 矢野
- それち1階差分でデータが定常になるのであれば1階差分で使った方が良いかもしれません。何度も差分をとると分析結果を解釈しにくくなりますので。 -- 矢野
- 矢野さん、でいいんでしょうか?ありがとうございます。親切に教えていただいたので、おそらくこれで上手い事いきそうです。教えていただいたプログラムでAIC導出も上手い事行きました。 -- 松本