「大人の教養・知識・気付き」を伸ばすブログ

一流の大人(ビジネスマン、政治家、リーダー…)として知っておきたい、教養・社会動向を意外なところから取り上げ学ぶことで“気付く力”を伸ばすブログです。

MENU

【数理統計】時系列解析Vol.03:AR過程

 前回は時系列解析の基礎として定常性やホワイトノイズを導入した後、いよいよMA過程を導入した。MA過程はランダムな変動の集積として現時点の値を説明するものだった。


zeitgeist.hatenablog.com

引き続き、以下の書籍を基に学んだことを整理していく。

3. ARMA過程(2)

 前回学んだことを復習しておこう。

 時系列解析で焦点を当てているのは異時点間の相関であるからそれをモデルに織り込むことが論点となる。その相関のモデル化には2つの考え方がある。1つは異なる時点の値に共通の成分を含める方法である。もう1つの方法はある時点に過去時点の値を盛り込む方法である。前者がMA過程であり後者がAR過程である。実はこれらを組み合わせたARMA過程を考えることもできる。

 このうち後者の方を今回は扱おう。

AR過程

 ある時点の値にそれ以前の値を盛り込む方法として自己回帰(AR)過程がある。n次AR過程は


\begin{aligned}
y_t=c+\sum_{j=1}^{n}\phi_{j}y_{t-j}+\varepsilon_t,\ \varepsilon_{t}\sim W.N.(\sigma^2)
\end{aligned}
と定義され、y_t\sim AR(n)と書く。確率的変動はMA過程と同様に攪乱項(\varepsilon_t)がもたらす。

f:id:suguru_125:20200506203236j:plainf:id:suguru_125:20200506203241j:plain
シミュレーション結果

また上に掲げるシミュレーションから予想されるように、定常性の有無はパラメータの値に依存して決まる。
 具体例を考えてみよう。たとえば1次AR過程は


\begin{aligned}
y_t=c+\phi_{1}y_{t-1}+\varepsilon_{t},\ \varepsilon_{t}\sim W.N.(\sigma^2)
\end{aligned}
と書ける。
 さて一般のn次AR過程について、ホワイトノイズの期待値が0であることに注意すると

\begin{aligned}
E[y_t]&=E\left[c+\sum_{j=1}^{n}\phi_{j}y_{t-j}+\varepsilon_t\right]\\
           &=c+\sum_{j=1}^{n}\phi_{j}E[y_{t-j}]+E[\varepsilon_{t}] \\
           &=c+\sum_{j=1}^{n}\phi_{j}E[y_{t-j}].
\end{aligned}
である。この過程が弱定常であると仮定すると、E[y_t]=\cdots=E[y_{t-n}]=\muであるから、

\begin{aligned}
\mu=c+\mu\sum_{j=1}^{n}\phi_{j}\ \ \ \ \therefore \mu=\frac{c}{1-\sum_{j=1}^{n}\phi_{j}}.
\end{aligned}
となる*1
 また分散は

\begin{aligned}
V[y_t]&=V\left[c+\sum_{j=1}^{n}\phi_{j}y_{t-j}+\varepsilon_t\right] \\
            &=\sum_{j=1}^{n}{\phi_{j}}^2V[y_{t-j}]+V[\varepsilon_t] \\
            &=\sum_{j=1}^{n}{\phi_{j}}^2V[y_{t-j}]+\sigma^2.
\end{aligned}
である。この過程が弱定常であると仮定すると、V[y_t]=V[y_{t-1}]=\cdots=\gamma_{0}であるから

\begin{aligned}
\gamma_{0}&=\gamma_{0}\sum_{j=1}^{n}{\phi_{j}}^2+\sigma^2 \\
\left(1-\sum_{j=1}^{n}{\phi_{j}}^2\right)\gamma_{0}&=\sigma^2.
\end{aligned}
すなわち、

\begin{aligned}
V[y_t]=\gamma_0=\displaystyle\frac{\sigma^2}{1-\sum_{j=1}^{n}{\phi_{j}}^2}
\end{aligned}
である。
 さらにラグk\leq nに対して自己共分散\gamma_kを計算する。AR過程の方程式

\begin{aligned}
y_t=c+\sum_{i=1}^{n}\phi_{i}y_{t-i}+\varepsilon_t
\end{aligned}
に対して両辺とそれぞれy_t,\cdots,y_{t-n}との共分散を取ることで


{\displaystyle 
\begin{eqnarray}
 \left\{
    \begin{array}{l}
        Cov[y_t,y_t]&=Cov\left[c+\sum_{i=1}^{n}\phi_{i}y_{t-i}+\varepsilon_t,y_t\right],\\
                         &\vdots  \\
        Cov[y_t,y_{t-n}]&=Cov\left[c+\sum_{i=1}^{n}\phi_{i}y_{t-i}+\varepsilon_t,y_{t-n}\right]
    \end{array}
  \right.
\end{eqnarray}
}

であり、これをまとめると


{\displaystyle 
\begin{eqnarray}
 \left\{
    \begin{array}{l}
         \gamma_0&=\displaystyle{\sum_{i=1}^{n}\phi_{i}\gamma_i},\\
                          &\vdots,  \\
         \gamma_n&=\displaystyle{\sum_{i=1}^{n}\phi_{i}\gamma_{n-i}}
    \end{array}
  \right.
\end{eqnarray}
}

このように自己共分散\gamma_kは上記の連立方程式の解として得られる。この連立方程式の両辺を\gamma_0で割って自己相関係数にしたものをYule-Walker方程式と呼ぶ。
 なお定常性の有無については、以下の方程式(AR方程式)


\begin{aligned}
1-\phi_{1}z-\cdots-\phi_{n}z^{n}=0
\end{aligned}

の解の絶対値がすべて1よりも大きいことと定常性をもつことが同値であることが知られている。また同方程式が共役な複素数z=a\pm{bi}をもつとき、この過程は周期が\displaystyle\frac{2\pi}{\cos^{-1}{\frac{a}{\sqrt{a^2+b^2}}}}となる循環成分をもつことも知られている。

まとめ

 以上のとおり、AR過程を導入した。AR過程はMA過程に比較して定常性を必ずしも持つわけではないという欠点をもつものの、自己共分散が解析(数値)的に計算できることや循環成分を表現することができることなど、解析において利点を有する(特にMA過程よりも直感的に理解しやすい)。

*1:ここから\sum_{j=1}^{n}\phi_{j}\neq 1でなければならないと示唆される。

プライバシーポリシー お問い合わせ