auto.arima() のオプションを試す
auto.arima() は多くの異なるモデルを推定する必要があり、できるだけ高速に動作させるためにさまざまな近道が使われます。そのため、最小の AICc 値を持たないモデルが返されることがあります。auto.arima() により良いモデルを見つけさせるには、任意引数 stepwise = FALSE を追加して、より大きなモデル集合を探索するようにします。
ここでは、あらかじめ読み込まれている a10 データに対して ARIMA モデルを見つけてみます。a10 には、1991 年から 2008 年までのオーストラリアにおける糖尿病治療薬補助金(月次、単位は豪ドル百万)が含まれています。演習を始める前に、コンソールで内容を確認してください。
この演習はコースの一部です
Rで学ぶ予測入門
演習の手順
auto.arima()のデフォルト設定を用いてa10の ARIMA モデルを推定し、fit1に保存します。- 段階的探索を使わずに
auto.arima()を実行してa10の ARIMA モデルを求め、fit2に保存します。 - コンソールで両方の
fit1とfit2に対してsummary()を実行し、どちらがより良いモデルかを判断してください。小数点以下 2 桁で、そのモデルの AICc 値はいくつですか? その数値をAICcに代入します。 - 最後に、AICc に基づいてより良いモデルを使い、2 年先の予測をプロットします。
hをそれに合わせて設定してください。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Find an ARIMA model for a10
fit1 <- ___
# Don't use a stepwise search
fit2 <- ___
# AICc of better model
AICc <- ___
# Compute 2-year forecasts from better model
___ %>% ___ %>% ___