「ARIMAモデル」の版間の差分
Albeit-Kun (トーク | 投稿記録) |
|||
(3人の利用者による、間の4版が非表示) | |||
1行目: | 1行目: | ||
'''【ありまもでる/えいあーるあいえむえいもでる (ARIMA (autoregressive integrated moving average) model)】''' | '''【ありまもでる/えいあーるあいえむえいもでる (ARIMA (autoregressive integrated moving average) model)】''' | ||
− | : | + | === 概要 === |
+ | |||
+ | <math>y_{t} \,</math> を非定常過程とし,<math>\varepsilon_{t} \,</math> を<math>\mbox{E}(\varepsilon_{t})=0 \,</math>,<math>\mbox{V}(\varepsilon_{t})=\sigma^{2} \,</math>,<math>\mbox{E}(\varepsilon_{t}\varepsilon_{s})=0 \,</math> <math>(t \ne s) \,</math>のホワイトノイズとする.<math>L \,</math> をラグ演算子 <math>L^{i}y_{t}=y_{t-i} \,</math>,<math>L^{i}\varepsilon_{t}=\varepsilon_{t-i} \,</math>(<math>i=1,2,\cdots \,</math>),<math>\phi(L) \,</math>, <math>\theta(L) \,</math> を <math>\textstyle \phi(L) \equiv 1-\sum_{i=1}^{p}\phi_{i}L^{i} \,</math>,<math>\textstyle \theta(L) \equiv 1+\sum_{i=1}^{p} \theta_{i}L^{i} \,</math>とする.<math>d \,</math> を自然数として, <math>y_{t} \,</math> の <math>d \,</math> 階階差 <math>(1-L)^{d}y_{t} \,</math> が定常な<math>\mbox{ARMA}(p,q) \,</math> モデルで表現できるとき, モデル <math>\phi(L)(1-L)^{d}y_{t} =\theta(L)\varepsilon_{t} \,</math>を次数 <math>(p,d,q) \,</math> の自己回帰和分移動平均(ARIMA)モデルと呼び, <math>\mbox{ARIMA}(p,d,q) \,</math> モデルと略記する. | ||
+ | |||
+ | === 詳説 === | ||
+ | |||
+ | t を時点を表わす添字 (整数) とし, <math>x_{t}\, </math> を <math>\mbox{E}(x_{t})=0\, </math> の[[弱定常過程]],<math>\varepsilon_{t}\, </math> を<math>\mbox{E}(\varepsilon_{t})=0\, </math>,<math>\mbox{V}(\varepsilon_{t})=\sigma^{2}\, </math>,<math>\mbox{E}(\varepsilon_{t}\varepsilon_{s})=0 (t \ne s)\, </math>を満たすホワイトノイズ (white noise) とする.また, <math>L\, </math> を時間を後退させる作用をもつラグ演算子 (lag operator)<math>L^{i}x_{t}=x_{t-i}\, </math>, <math>L^{i}\varepsilon_{t}=\varepsilon_{t-i}\, </math>(<math>i=1,2,\cdots\, </math>)とし, <math>\phi(L)\, </math>, <math>\theta(L)\, </math> を <math>\textstyle \phi(L)=1-\sum_{i=1}^{p}\phi_{i}L^{i}\, </math>,<math>\textstyle \theta(L)=1+\sum_{i=1}^{q} \theta_{i}L^{i}\, </math> で定義される多項式ラグ演算子とする.(ただし, 多項式 <math>\phi(z)=0\, </math>, <math>\theta(z)=0\, </math> には共通根はないものとする.)<math>\phi_{i}\, </math> (<math>i=1,2,\cdots,p\, </math>), <math>\theta_{i}\, </math> (<math>i=1,2,\cdots,q\, </math>) はパラメータ (一定) である. | ||
+ | |||
+ | 弱定常過程 <math>x_{t}\, </math> の確率的変動が <math>\phi(L)x_{t} =\theta(L)\varepsilon_{t}\, </math>(すなわち,<math>x_{t}=\phi_{1}x_{t-1}+\cdots+\phi_{p}x_{t-p}+\varepsilon_{t}+\theta_{1}\varepsilon_{t-1}+\cdots+\theta_{q}\varepsilon_{t-q}\, </math>)で表わされるとき, このモデルを次数 <math>(p,q)\, </math> の自己回帰移動平均モデル(autoregressive moving average model) と呼び, <math>ARMA(p,q)\, </math> モデルと略記する.<math>ARMA(p,q)\, </math> モデルが条件「<math>\phi(z)=0\, </math> の根はすべて単位円 <math>|z|=1\, </math> 外にある」を満たせば, これを <math>MA(\infty)\, </math> モデル :<math>x_{t}=\psi(L)\varepsilon_{t}\, </math> として表現出来る.(ただし <math>\textstyle \psi(z)=\theta(z)/\phi(z)=\sum_{i=0}^{\infty}\psi_{i}z^{i}\, </math>,<math>\psi_{0}=1\, </math>, <math>\textstyle \sum_{i=0}^{\infty}|\psi_{i}|<\infty\, </math>.)また逆に条件 (反転可能性の条件 (invertibility condition) )「<math>\theta(z)=0\, </math> の根はすべて単位円 <math>|z|=1\, </math> 外にある」を満たせば, これを <math>AR(\infty)\, </math> モデル :<math>\pi(L)x_{t}=\varepsilon_{t}\, </math> として表現出来る.(ただし <math>\textstyle \pi(z)=\phi(z)/\theta(z)=\sum_{i=0}^{\infty}\pi_{i}z^{i}\, </math>, <math>\pi_{0}=1\, </math>, \<math>sum_{i=0}^{\infty}|\pi_{i}|<\infty\, </math>.) | ||
+ | |||
+ | AR なる用語は <math>x_{t}\, </math> を自身の過去の値に回帰することに由来しており, [[ARモデル]]は理解しやすい構造を持っている.一方, [[MAモデル]]は過程の理論的性質を調べる上で重要である.また <math>ARMA(p,q)\, </math> モデルは <math>AR(p)\, </math> モデルと <math>MA(q)\, </math> モデルを混合したモデルであり, これらのモデルを単独で使用した場合に比べてより少ないパラメータで定常過程の種々の性質を表現出来る点に特徴がある. | ||
+ | |||
+ | 次に<math> d\, </math> を自然数として, 次数 <math>d\, </math> の階差演算子 (difference operator)<math>(1-L)^{d}\, </math> を<math>\textstyle (1-L)^{d}=\sum_{i=0}^{d} {}_{d}C_{i}(-1)^{i}L^{i}\, </math>(<math>d=1,2,\cdots\, </math>) と定義する.非定常過程 <math>y_{t}\, </math> の <math>d\, </math> 階階差 <math>x_{t}=(1-L)^{d}y_{t}\, </math> が弱定常となるとき,<math>y_{t}\, </math> は階差次数 <math>d\, </math> の和分モデル (integrated model) に従うと言い,このモデルを <math>I(d)\, </math> モデルと略記する.中でも <math>x_{t}=(1-L)^{d}y_{t}\, </math> が特に <math>ARMA(p,q)\, </math> モデルに従い,<math>\phi(L)(1-L)^{d}y_{t} = \theta(L)\varepsilon_{t}\, </math>と表わせるとき, このモデルを次数 <math>(p,d,q)\, </math> の[[自己回帰和分移動平均モデル]] (autoregressive integrated moving average model) と呼び,<math>ARIMA(p,d,q)\, </math> モデルと略記する.特に <math>ARIMA(p,0,q)\, </math> モデルが <math>ARMA(p,q)\, </math> モデルに他ならない.また <math>ARIMA(0,1,0)\, </math> モデルは[[ランダムウォーク]]モデル (randomwalk model) となる. | ||
+ | |||
+ | Box and Jenkins [2] が従来の研究成果をふまえて, [[ARIMAモデル]]の 1) 同定 (identification), 2) 推定 (estimation), 3) 診断 (diagnostic checking), 4) [[予測]]と制御 (forecasting and control)に関する統計的分析法を体系的に提示して以来, ARIMA モデルは[[時系列解析]]に不可欠なパラメトリックモデルとして重要な役割を果たしている. | ||
+ | |||
+ | 特に 1) モデルの同定に関しては, 多くの時系列において階差次数 <math>d\, </math> が高々 1 (ないし 2) であることが経験的に知られている.また AR 次数 <math>p\, </math>, MA 次数 <math>q\, </math> については, [[自己相関関数]]と偏自己相関関数の特徴をもとにこれらを決定する方法 [2] 以外に, AIC (Akaike's information criterion) 最小化法が用いられることも多い.また, 季節変動を含む経済時系列解析に有効なモデルとして, 周期 <math>s\, </math> (例えば4 半期データ, 月次データに応じて <math>s=4, 12\, </math> 等) の季節変動を取り扱う季節的ARIMA モデル(seasonal ARIMA model) があげられる [2]. 2) パラメータの推定に関しては最尤法をはじめとした各種の非線形推定法が提案されており, 3) モデルの診断についても, 残差系列がホワイトノイズに従うか否かを残差系列の自己相関関数にもとづき検定する方法をはじめとして幾つかの統計的仮説検定法が提案されている [1], [3], [6]. 4) 予測や制御に関しては, 例えば (単純) [[指数平滑法|指数平滑化法]]が <math>ARIMA(0,1,1)\, </math> モデルの予測に最適な方法であることが明らかにされており,さらに ARIMA モデルの状態空間表現 (state space representation) や[[カルマンフィルター]] (Kalmanfilter)との関連性が明らかにされている.また ARIMA モデルは各種の需要予測や[[経済予測]]に有効なモデルであることが知られている. | ||
+ | |||
+ | ARIMA モデルは様々な方向に拡張されている. 例えば金融時系列解析の分野における分散変動を考慮した[[ARCHモデル|ARCH, GARCH]], EGARCH (exponential GARCH), IGARCH (integrated GARCH) モデルや,階差次数 <math>d\, </math> を実数に拡張した ARFIMA モデル (AR fractionally integratedMA model) はその一例である [3], [5].また, [[計量経済モデル]]との関連では,AR(I)MA モデルを多変量化した VAR(I)MA モデル (vector AR(I)MA model),外生変数 (exogenous variables) を取り入れた AR(I)MAX モデル(AR(I)MA model with exogenous variables) 等の構築や, これらをもとにしたグレンジャー因果関係 (Granger causality) の検証がなされており,さらにランダムウォークモデルとの関連で階差次数 <math>d=1\, </math> の統計的仮説検定を取り扱う単位根検定 (unit root test) や, 複数の和分モデルの一次結合が定常モデルに従う共和分モデル (cointegrated model) の構築等がなされている [4], [5], [7]. | ||
+ | |||
+ | |||
+ | |||
+ | ---- | ||
+ | '''参考文献''' | ||
+ | |||
+ | [1] T. W. Anderson, ''The Statistical Analysis of Time Series'', John Wiley, 1971. | ||
+ | |||
+ | [2] G. E. P. Box and G. M. Jenkins, ''Time Series Analysis : Forecasting and Control'', Holden-Day, 1970. (rev. ed., Holden-Day, 1976. 3rd ed. by G. E. P. Box, G. M. Jenkins and G. C. Reinsel, Prentice-Hall, 1994.) | ||
+ | |||
+ | [3] P. J. Brockwell and R. A. Davis, ''Time Series : Theory and Methods (2nd ed.),'' Springer-Verlag, 1991. | ||
+ | |||
+ | [4] W. A. Fuller, ''Introduction to Statistical Time Series (2nd ed.)'', John Wiley, 1996. | ||
+ | |||
+ | [5] J. D. Hamilton, ''Time Series Analysis'', Princeton University Press, 1994. | ||
+ | |||
+ | [6] E. J. Hannan, ''Multiple Time Series'', John Wiley, 1970. | ||
+ | |||
+ | [7] G. C. Reinsel, ''Elements of Multivariate Time Series Analysis (2nd ed.)'', Springer-Verlag, 1997. | ||
+ | |||
+ | [[category:予測|ありまもでる/えいあーるあいえむえいもでる]] |
2008年11月5日 (水) 16:26時点における最新版
【ありまもでる/えいあーるあいえむえいもでる (ARIMA (autoregressive integrated moving average) model)】
概要
構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_{t} \,} を非定常過程とし,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \varepsilon_{t} \,} を構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{E}(\varepsilon_{t})=0 \,} ,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{V}(\varepsilon_{t})=\sigma^{2} \,} ,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{E}(\varepsilon_{t}\varepsilon_{s})=0 \,} 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (t \ne s) \,} のホワイトノイズとする.構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L \,} をラグ演算子 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L^{i}y_{t}=y_{t-i} \,} ,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L^{i}\varepsilon_{t}=\varepsilon_{t-i} \,} (構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i=1,2,\cdots \,} ),, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \theta(L) \,} を 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \phi(L) \equiv 1-\sum_{i=1}^{p}\phi_{i}L^{i} \,} ,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \theta(L) \equiv 1+\sum_{i=1}^{p} \theta_{i}L^{i} \,} とする.構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d \,} を自然数として, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_{t} \,} の 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d \,} 階階差 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (1-L)^{d}y_{t} \,} が定常な構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{ARMA}(p,q) \,} モデルで表現できるとき, モデル 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \phi(L)(1-L)^{d}y_{t} =\theta(L)\varepsilon_{t} \,} を次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (p,d,q) \,} の自己回帰和分移動平均(ARIMA)モデルと呼び, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{ARIMA}(p,d,q) \,} モデルと略記する.
詳説
t を時点を表わす添字 (整数) とし, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x_{t}\, } を 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{E}(x_{t})=0\, } の弱定常過程,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \varepsilon_{t}\, } を構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{E}(\varepsilon_{t})=0\, } ,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{V}(\varepsilon_{t})=\sigma^{2}\, } ,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{E}(\varepsilon_{t}\varepsilon_{s})=0 (t \ne s)\, } を満たすホワイトノイズ (white noise) とする.また, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L\, } を時間を後退させる作用をもつラグ演算子 (lag operator)構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L^{i}x_{t}=x_{t-i}\, } , 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L^{i}\varepsilon_{t}=\varepsilon_{t-i}\, } (構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i=1,2,\cdots\, } )とし, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \phi(L)\, } , 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \theta(L)\, } を 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \phi(L)=1-\sum_{i=1}^{p}\phi_{i}L^{i}\, } ,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \theta(L)=1+\sum_{i=1}^{q} \theta_{i}L^{i}\, } で定義される多項式ラグ演算子とする.(ただし, 多項式 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \phi(z)=0\, } , 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \theta(z)=0\, } には共通根はないものとする.)構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \phi_{i}\, } (構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i=1,2,\cdots,p\, } ), 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \theta_{i}\, } (構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i=1,2,\cdots,q\, } ) はパラメータ (一定) である.
弱定常過程 の確率的変動が 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \phi(L)x_{t} =\theta(L)\varepsilon_{t}\, } (すなわち,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x_{t}=\phi_{1}x_{t-1}+\cdots+\phi_{p}x_{t-p}+\varepsilon_{t}+\theta_{1}\varepsilon_{t-1}+\cdots+\theta_{q}\varepsilon_{t-q}\, } )で表わされるとき, このモデルを次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (p,q)\, } の自己回帰移動平均モデル(autoregressive moving average model) と呼び, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ARMA(p,q)\, } モデルと略記する.構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ARMA(p,q)\, } モデルが条件「構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \phi(z)=0\, } の根はすべて単位円 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle |z|=1\, } 外にある」を満たせば, これを モデル :構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x_{t}=\psi(L)\varepsilon_{t}\, } として表現出来る.(ただし 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \psi(z)=\theta(z)/\phi(z)=\sum_{i=0}^{\infty}\psi_{i}z^{i}\, } ,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \psi_{0}=1\, } , 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \sum_{i=0}^{\infty}|\psi_{i}|<\infty\, } .)また逆に条件 (反転可能性の条件 (invertibility condition) )「構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \theta(z)=0\, } の根はすべて単位円 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle |z|=1\, } 外にある」を満たせば, これを 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle AR(\infty)\, } モデル : として表現出来る.(ただし 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \pi(z)=\phi(z)/\theta(z)=\sum_{i=0}^{\infty}\pi_{i}z^{i}\, } , 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \pi_{0}=1\, } , \構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle sum_{i=0}^{\infty}|\pi_{i}|<\infty\, } .)
AR なる用語は 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x_{t}\, } を自身の過去の値に回帰することに由来しており, ARモデルは理解しやすい構造を持っている.一方, MAモデルは過程の理論的性質を調べる上で重要である.また モデルは 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle AR(p)\, } モデルと 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle MA(q)\, } モデルを混合したモデルであり, これらのモデルを単独で使用した場合に比べてより少ないパラメータで定常過程の種々の性質を表現出来る点に特徴がある.
次に構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d\, } を自然数として, 次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d\, } の階差演算子 (difference operator)構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (1-L)^{d}\, } を構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle (1-L)^{d}=\sum_{i=0}^{d} {}_{d}C_{i}(-1)^{i}L^{i}\, } (構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d=1,2,\cdots\, } ) と定義する.非定常過程 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_{t}\, } の 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d\, } 階階差 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x_{t}=(1-L)^{d}y_{t}\, } が弱定常となるとき,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_{t}\, } は階差次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d\, } の和分モデル (integrated model) に従うと言い,このモデルを 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle I(d)\, } モデルと略記する.中でも 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x_{t}=(1-L)^{d}y_{t}\, } が特に モデルに従い,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \phi(L)(1-L)^{d}y_{t} = \theta(L)\varepsilon_{t}\, } と表わせるとき, このモデルを次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (p,d,q)\, } の自己回帰和分移動平均モデル (autoregressive integrated moving average model) と呼び,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ARIMA(p,d,q)\, } モデルと略記する.特に モデルが 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ARMA(p,q)\, } モデルに他ならない.また 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ARIMA(0,1,0)\, } モデルはランダムウォークモデル (randomwalk model) となる.
Box and Jenkins [2] が従来の研究成果をふまえて, ARIMAモデルの 1) 同定 (identification), 2) 推定 (estimation), 3) 診断 (diagnostic checking), 4) 予測と制御 (forecasting and control)に関する統計的分析法を体系的に提示して以来, ARIMA モデルは時系列解析に不可欠なパラメトリックモデルとして重要な役割を果たしている.
特に 1) モデルの同定に関しては, 多くの時系列において階差次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d\, } が高々 1 (ないし 2) であることが経験的に知られている.また AR 次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle p\, } , MA 次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle q\, } については, 自己相関関数と偏自己相関関数の特徴をもとにこれらを決定する方法 [2] 以外に, AIC (Akaike's information criterion) 最小化法が用いられることも多い.また, 季節変動を含む経済時系列解析に有効なモデルとして, 周期 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle s\, } (例えば4 半期データ, 月次データに応じて 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle s=4, 12\, } 等) の季節変動を取り扱う季節的ARIMA モデル(seasonal ARIMA model) があげられる [2]. 2) パラメータの推定に関しては最尤法をはじめとした各種の非線形推定法が提案されており, 3) モデルの診断についても, 残差系列がホワイトノイズに従うか否かを残差系列の自己相関関数にもとづき検定する方法をはじめとして幾つかの統計的仮説検定法が提案されている [1], [3], [6]. 4) 予測や制御に関しては, 例えば (単純) 指数平滑化法が 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ARIMA(0,1,1)\, } モデルの予測に最適な方法であることが明らかにされており,さらに ARIMA モデルの状態空間表現 (state space representation) やカルマンフィルター (Kalmanfilter)との関連性が明らかにされている.また ARIMA モデルは各種の需要予測や経済予測に有効なモデルであることが知られている.
ARIMA モデルは様々な方向に拡張されている. 例えば金融時系列解析の分野における分散変動を考慮したARCH, GARCH, EGARCH (exponential GARCH), IGARCH (integrated GARCH) モデルや,階差次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d\, } を実数に拡張した ARFIMA モデル (AR fractionally integratedMA model) はその一例である [3], [5].また, 計量経済モデルとの関連では,AR(I)MA モデルを多変量化した VAR(I)MA モデル (vector AR(I)MA model),外生変数 (exogenous variables) を取り入れた AR(I)MAX モデル(AR(I)MA model with exogenous variables) 等の構築や, これらをもとにしたグレンジャー因果関係 (Granger causality) の検証がなされており,さらにランダムウォークモデルとの関連で階差次数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d=1\, } の統計的仮説検定を取り扱う単位根検定 (unit root test) や, 複数の和分モデルの一次結合が定常モデルに従う共和分モデル (cointegrated model) の構築等がなされている [4], [5], [7].
参考文献
[1] T. W. Anderson, The Statistical Analysis of Time Series, John Wiley, 1971.
[2] G. E. P. Box and G. M. Jenkins, Time Series Analysis : Forecasting and Control, Holden-Day, 1970. (rev. ed., Holden-Day, 1976. 3rd ed. by G. E. P. Box, G. M. Jenkins and G. C. Reinsel, Prentice-Hall, 1994.)
[3] P. J. Brockwell and R. A. Davis, Time Series : Theory and Methods (2nd ed.), Springer-Verlag, 1991.
[4] W. A. Fuller, Introduction to Statistical Time Series (2nd ed.), John Wiley, 1996.
[5] J. D. Hamilton, Time Series Analysis, Princeton University Press, 1994.
[6] E. J. Hannan, Multiple Time Series, John Wiley, 1970.
[7] G. C. Reinsel, Elements of Multivariate Time Series Analysis (2nd ed.), Springer-Verlag, 1997.