自己回帰和分移動平均モデル

統計学や計量経済学、特に時系列分析において、自己回帰和分移動平均（じこかいきわぶんいどうへいきん、、略称: ARIMA）モデルは、自己回帰移動平均（ARMA）モデルの一般化である。これらのモデルは、データの理解を深めるため、または将来のポイントを予測するために、時系列データに適用される。
ARIMAモデルは、データが（分散/自己共分散ではなく）平均に関して非定常性を示す場合に適用され、初期の差分ステップ（モデルの「Integrated 和分」部分に対応）を 1回以上適用して平均関数（すなわち、トレンド）の非定常性を排除することができる。時系列に季節性が見られる場合は、季節成分を除去するために季節的差分を適用することができる。ウォルドの分解定理によれば、ARMAモデルは規則的な（つまり純粋に非決定論的な）広義の定常時系列を記述するのに理論的には十分であるので、ARMAモデルを使用する前に、例えば差分を使用して非定常時系列を定常化することが主な動機となる。時系列に予測可能なサブプロセス（純粋な正弦波や複素数指数プロセス）が含まれている場合、予測可能な成分はARIMAのフレームワークでは平均非ゼロで周期的な（つまり季節的な）成分として扱われるので、季節的な差分処理によって除去されることに注意が必要である。
ARIMAのAR（autoregressive、自己回帰）の部分は、関心のある展開する変数がそれ自体の遅延した値（すなわち、以前の値）に回帰されることを示している。MA（moving average、移動平均）の部分は、回帰誤差が実際には、同時期および過去の様々な時点で発生した誤差項の線型結合であることを示している。I （integrated 和分）の部分は、データの値が過去の値との差分に置き換えられていることを示している（この差分処理は複数回行われる場合もある）。これらの特徴の目的は、モデルがデータにできるだけ適合するようにすることである。
非季節ARIMAモデルは、一般にformula_1 と表記される。パラメータ"p"、"d"、"q"は非負の整数で、pは自己回帰モデルの次数（タイムラグの数）、"d"は差分の階数（データの過去の値を差し引いた回数）、"q"は移動平均モデルの次数を表す。
季節ARIMAモデルは、通常formula_2 と表記される。ここで、"m"は各季節の期間の数を意味し、大文字の"P"、"D"、"Q"はARIMAモデルの季節部分の自己回帰項、差分項、移動平均項を意味する 。
3つの項のうち2つがゼロの場合、モデルを表す頭字語のうちゼロではないパラメータを用いてモデルを記載することができる。例えば、formula_3 は formula_4、formula_5 は formula_6、formula_7 は formula_8 と記載される。
ARIMAモデルは、ボックス・ジェンキンス法に従って推定することができる。
定義.
時系列データ"X"（"t"は整数インデックス、"Xt"は実数）が与えられたとき、formula_9 モデルは下記のようになる。
または同等に
ここで、formula_12はラグ演算子、formula_13はモデルの自己回帰部分のパラメーター、formula_14は移動平均部分のパラメータ、formula_15は誤差項である。誤差項formula_15は、一般に、平均ゼロの正規分布からサンプリングされた独立同分布の変数であるとみなされる。
ここで、多項式 formula_17 が多重度"d"の単位根（因子"1 - L"）を持つ場合、次のように書き換えることができる。
formula_1 プロセスでは、"p" = "p'-d"で用いて、次のように表す。
したがって、"d"個の単位根を持つ自己回帰多項式を持つformula_21 プロセスの特殊なケースと考えることができる。 このため、"d > 0" のARIMAモデルで正確に記述されるプロセスは広義の定常ではない。
上記は次のように一般化できる。
これは、ドリフト formula_23 を伴う "formula_1" プロセスを定義する。
差分.
定常時系列の性質は、観測された時刻に依存しない。具体的には、広義の定常時系列では、平均と分散/自己共分散は時間の経過とともに一定になる。統計における差分とは、非定常時系列を平均的な意味で定常化するために（つまり、非定常トレンドを除去するために）適用される変換であり、分散/自己共分散の非定常性とは関係がない。同様に、季節性時系列に季節差分を適用して季節成分を除去する。信号処理、特にフーリエ・スペクトル解析理論の観点からは、トレンドは非定常時系列のスペクトルにおける低周波数部分であり、季節はそのスペクトルにおける周期的な周波数部分である。したがって、差分はハイパス（つまり、ローストップ）フィルタとして、季節差分はコムフィルタとして機能し、それぞれ低周波のトレンドと周期的な周波数の季節を（時間領域で直接ではなく）スペクトル領域で抑制することができる。この観点から、差分と季節差分の哲学、数学、力、欠点を説明することができる。
データの差分を取るために、連続した観測値の差を計算する。数学的には次のようになる。
差分は時系列のレベルの変化を取り除き、トレンドと季節性を排除し、結果的に時系列の平均値を安定させる。
定常時系列を得るために、2回に渡ってデータの差分を取る必要がある場合もあり、これは 2次差分と呼ばれる。
データの差分を取るもう一つの方法として、季節差分がある。これは、観測値と前の季節（例えば1年）の対応する観測値との差を計算するものである。これは次のように示される。
そして、この差分を取ったデータを用いて、ARMAモデルを推定する。
例.
いくつかのよく知られていた特殊なケースは、自然に生じたり、他の一般的な予測モデルと数学的に同等であったりする。例えば：
次数の選択.
次数"p"および"q"は、サンプル自己相関関数（ACF）、偏自己相関関数（PACF）、拡張自己相関関数（EACF）法を用いて決定することができる。
その他の代替法として、AIC、BICなどがある。非季節性ARIMAモデルの次数を決定するためには赤池情報量規準（AIC）が有用である。AICは次のように書かれる。
ここで、"L"はデータの尤度、"p"は自己回帰モデル部分の次数、"q"は移動平均モデル部分の次数であり、"k' はARIMAモデルの切片を表す。AICでは "k" = 1の場合はARIMAモデルに切片があり（"c" ≠ 0）、"k" = 0 の場合はARIMAモデルに切片がない（"c" = 0）ことになる。
ARIMAモデルの補正AICは、次のように書くことができる。
ベイズ情報量規準（BIC）は、次のように書くことができる。
目標は、良いモデルのAIC、AICc、BICの値を最小化することである。調査するモデルの範囲でこれらの基準の一つの値が低ければ低いほど、そのモデルはデータに適している。AICとBICは 2つのまったく異なる目的で使用される。AICがモデルを現実の状況に近づけようとするのに対し、BICは完全な適合性をみつけようとする。BICのアプローチは、現実の複雑なデータに完璧にフィットすることはないと批判されることが多いが、AICに比べてパラメータが多いことでモデルに大きなペナルティを与えるため、選択のための有効な手法であることに変わりはない。
AICcは、差分の次数が等しいARIMAモデルの比較にのみ使用できる。差分の次数が異なるARIMAモデルについては、RMSEをモデルの比較に使用することができる。
ARIMAモデルを使用した予測.
ARIMAモデルは、2つのモデルの「カスケード」と見なすことができる。1つ目は非定常である：
2番目は広義の定常である：
こうして、プロセスYtを自己回帰予測法の一般化を用いて予測することができる。
予測区間.
ARIMAモデルの予測区間（予測の信頼区間）は、残差が無相関で正規分布しているという仮定に基づいている。これらの仮定のいずれかが当てはまらない場合、予測間隔が正しくない可能性がある。このため、研究者は予測区間を作成する前に、過程をチェックするために残差のACFとヒストグラムをプロットする。
95%の予測区間は
formula_43
ここで、formula_44 は formula_45 の分散である
formula_46 のとき、パラメータや次数に関係なく、すべてのARIMAモデルに適用される。対して。
formula_47 の場合、 formula_48
一般に、ARIMAモデルからの予測区間は、予測期間が長くなるにつれて広がる。
バリエーションと拡張.
ARIMAモデルにはいくつかのバリエーションがある。複数の時系列を使用する場合は formula_50 をベクトルと考えることができ、VARIMAモデルが適切な場合がある。モデルに季節効果が疑われるときは、モデルのARやMAの次数を増やすよりも、SARIMA（季節ARIMA）モデルを使用する方が一般的には良いと考えられる。時系列が長距離依存性を示すと疑われる場合、フラクショナルARIMA（FARIMAまたはARFIMA）モデルとも呼ばれる自己回帰フラクショナル和分移動平均モデルでは、"d" パラメーターに非整数値を持たせることができる。
ソフトウェアの実装.
ARIMAモデルに適切なパラメータを見つけるため、Box-Jenkinsパラメータ最適化のような方法論を適用する様々なパッケージがある。