ARIMAモデルを使用した予測に足を踏み入れているとき、季節性とドリフトに適合するARIMAに基づいて予測を改善する方法を理解しようとしています。
私のデータは次の時系列です(3年以上、明確な上昇傾向と目に見える季節性があり、ラグ12、24、36の自己相関ではサポートされていないようです)。
> bal2sum3years.ts
Jan Feb Mar Apr May Jun Jul Aug
2010 2540346 2139440 2218652 2176167 2287778 1861061 2000102 2560729
2011 3119573 2704986 2594432 2362869 2509506 2434504 2680088 2689888
2012 3619060 3204588 2800260 2973428 2737696 2744716 3043868 2867416
Sep Oct Nov Dec
2010 2232261 2394644 2468479 2816287
2011 2480940 2699780 2760268 3206372
2012 2951516 3119176 3032960 3738256
によって提案されたモデルauto.arima(bal2sum3years.ts)
は私に次のモデルを与えました:
Series: bal2sum3years.ts
ARIMA(0,0,0)(0,1,0)[12] with drift
Coefficients:
drift
31725.567
s.e. 2651.693
sigma^2 estimated as 2.43e+10: log likelihood=-321.02
AIC=646.04 AICc=646.61 BIC=648.39
ただし、acf(bal2sum3years.ts,max.lag=35)
0.3を超えるacf係数は表示されません。ただし、データの季節性はかなり明白です。毎年の初めに急上昇します。これは、シリーズがグラフでどのように見えるかです。
fit=Arima(bal2sum3years.ts,seasonal=list(order=c(0,1,0),period=12),include.drift=TRUE)
関数によって呼び出されるを使用した予測によりforecast(fit)
、次の12か月の平均は、データの最後の12か月と定数に等しくなります。これはplot(forecast(fit))
、
また、自己相関ではありませんが、正の平均(非ゼロ)をもつ残差もチェックしました。
私の意見では、近似は元の時系列を正確にモデル化していません(青は元の時系列、赤はfitted(fit)
:
ゲストは、モデルが間違っていますか?何か不足していますか?モデルを改善するにはどうすればよいですか?モデルは文字通り過去12か月かかり、次の12か月を達成するために定数を追加するようです。
私は時系列予測モデルと統計の比較的初心者です。