時系列データ分析(2) のバックアップ差分(No.1)

バックアップ一覧
現在との差分を表示
ソースを表示
バックアップを表示
時系列データ分析(2) へ行く。
- 1 (2009-01-01 (木) 19:20:10)
- 2 (2009-01-01 (木) 22:35:09)
- 3 (2009-01-01 (木) 22:57:18)

追加された行はこの色です。
削除された行はこの色です。

[[時系列データ分析]]

*モデル選択の実際 [#c66db483]

　[[Kleiber&Zeileis[2008]>http://www.amazon.co.jp/dp/0387773169]]には、一気に24通りのSARIMAを実行して結果を出すスクリプトが例として載っています(pp.161-162)。このスクリプトはGPL(ver.2)のもとで配布されているAERパッケージに含まれていますから、スクリプトを取り上げて独自の説明をするだけなら著作権侵害にもならないし、本の売れ行きを妨げることにもならないでしょう。この本は計量経済学全体に関するRによって統一された初めての教科書ですし、売れて欲しいと思います。これが売れなければ日本語の教科書なんか出ませんからね。

　AERのインストールが終わったとしましょう(2008年以降、並河が授業で配布したR-2.8.0にはインストールしてあります)。スクリプトは

 (インストールした場所)\R-2.8.0\library\AER\demo\Ch-TimeSeries.R

　にあります。

 ###################################################
 ### chunk number 25: arima-setup
 ###################################################
 　
 nd_pars <- expand.grid(ar = 0:2, diff = 1, ma = 0:2,
   sar = 0:1, sdiff = 1, sma = 0:1)
 nd_aic <- rep(0, nrow(nd_pars))
 for(i in seq(along = nd_aic)) nd_aic[i] <- AIC(arima(nd,
   unlist(nd_pars[i, 1:3]), unlist(nd_pars[i, 4:6])),
   k = log(length(nd)))
 nd_pars[which.min(nd_aic),]