コンテンツにスキップ
Wikipedia

ARCHモデル

出典: フリー百科事典『ウィキペディア(Wikipedia)』

ARCHモデル(アーチモデル、: autoregressive conditional heteroscedasticity model, ARCH model)とは、金融経済学統計学計量経済学などにおいて分散不均一性を示す時系列データに適用されるモデル。日本語では、「分散自己回帰モデル」「分散不均一モデル」等と称される。1982年ロバート・エングルによって提案された[1] 。特に金融時系列データへの適用事例が多い。

分散不均一性

[編集 ]

株式の収益率をプロットすると、ある時期(景気が安定して拡大している時期など)には変動の程度(ボラティリティ)が平均して小さく、別の時期(不況の直前など)にはボラティリティが平均して大きくなる傾向が観察される。このようなボラティリティが時期によって異なった水準を示すことをボラティリティ・クラスタリング(: volatility clustering)、または分散不均一性(: heteroscedasticity)と呼ぶ。分散不均一性は金融時系列データをはじめ幅広く見られる現象である。

ARCH(q)モデル

[編集 ]

時刻 t {\displaystyle t} {\displaystyle t} における時系列データ y t {\displaystyle y_{t}} {\displaystyle y_{t}} の時刻 t 1 {\displaystyle t-1} {\displaystyle t-1} までの情報による条件付き期待値を μ t {\displaystyle \mu _{t}} {\displaystyle \mu _{t}} とする。 y t {\displaystyle y_{t}} {\displaystyle y_{t}} μ t {\displaystyle \mu _{t}} {\displaystyle \mu _{t}} の差を u t = y t μ t {\displaystyle u_{t}=y_{t}-\mu _{t}} {\displaystyle u_{t}=y_{t}-\mu _{t}} とする。さらに

u t = σ t ε t {\displaystyle u_{t}=\sigma _{t}\varepsilon _{t}} {\displaystyle u_{t}=\sigma _{t}\varepsilon _{t}}

と分解できるとする。ただし ε t {\displaystyle \varepsilon _{t}} {\displaystyle \varepsilon _{t}} は平均が0、分散が1の確率変数で、 σ t {\displaystyle \sigma _{t}} {\displaystyle \sigma _{t}} はボラティリティであり、時刻 t 1 {\displaystyle t-1} {\displaystyle t-1} までの情報で確定していると考える。すなわち、時刻 t 1 {\displaystyle t-1} {\displaystyle t-1} の時点で、時刻 t {\displaystyle t} {\displaystyle t} におけるこの時系列データのボラティリティは予測できる、と考えるのである。他方、 u t {\displaystyle u_{t}} {\displaystyle u_{t}} そのものは実際に時刻 t {\displaystyle t} {\displaystyle t} になり確率変数 ε t {\displaystyle \varepsilon _{t}} {\displaystyle \varepsilon _{t}} の値が確定するまでは確定しない。よって y t {\displaystyle y_{t}} {\displaystyle y_{t}} 自体は

y t = μ t + u t = μ t + σ t ε t {\displaystyle y_{t}=\mu _{t}+u_{t}=\mu _{t}+\sigma _{t}\varepsilon _{t}} {\displaystyle y_{t}=\mu _{t}+u_{t}=\mu _{t}+\sigma _{t}\varepsilon _{t}}

と表せる。ARCH(q)モデルの下で条件付ボラティリティ σ t {\displaystyle \sigma _{t}} {\displaystyle \sigma _{t}} は以下の式で決定される。

σ t 2 = α 0 + α 1 u t 1 2 + + α q u t q 2 = α 0 + i = 1 q α i u t i 2 {\displaystyle \sigma _{t}^{2}=\alpha _{0}+\alpha _{1}u_{t-1}^{2}+\cdots +\alpha _{q}u_{t-q}^{2}=\alpha _{0}+\sum _{i=1}^{q}\alpha _{i}u_{t-i}^{2}} {\displaystyle \sigma _{t}^{2}=\alpha _{0}+\alpha _{1}u_{t-1}^{2}+\cdots +\alpha _{q}u_{t-q}^{2}=\alpha _{0}+\sum _{i=1}^{q}\alpha _{i}u_{t-i}^{2}}

つまりARCH(q)モデルでは、q 期前までの平均からの乖離部分 u t i {\displaystyle u_{t-i}} {\displaystyle u_{t-i}} の2乗が条件付きボラティリティに影響を与えている。仮定から v t = u t 2 E t 1 [ u t 2 ] = u t 2 σ t 2 {\displaystyle v_{t}=u_{t}^{2}-E_{t-1}[u_{t}^{2}]=u_{t}^{2}-\sigma _{t}^{2}} {\displaystyle v_{t}=u_{t}^{2}-E_{t-1}[u_{t}^{2}]=u_{t}^{2}-\sigma _{t}^{2}} であるのでARCHモデルの決定式は

u t 2 = α 0 + i = 1 q α i u t i 2 + v t {\displaystyle u_{t}^{2}=\alpha _{0}+\sum _{i=1}^{q}\alpha _{i}u_{t-i}^{2}+v_{t}} {\displaystyle u_{t}^{2}=\alpha _{0}+\sum _{i=1}^{q}\alpha _{i}u_{t-i}^{2}+v_{t}}

と書き直すことが出来る。さらに v t {\displaystyle v_{t}} {\displaystyle v_{t}} E [ u t i 2 v t ] = 0 , i = 1 , {\displaystyle E[u_{t-i}^{2}v_{t}]=0,\;i=1,\dots } {\displaystyle E[u_{t-i}^{2}v_{t}]=0,\;i=1,\dots } であることも分かる。つまり u t 2 {\displaystyle u_{t}^{2}} {\displaystyle u_{t}^{2}} から見ると q 次の自己回帰モデルと見なせる。よって u t 2 {\displaystyle u_{t}^{2}} {\displaystyle u_{t}^{2}} について自己回帰(: autoregressive)であり、条件付き(: conditional)ボラティリティ σ t {\displaystyle \sigma _{t}} {\displaystyle \sigma _{t}} が分散不均一性(: heteroscedasticity)を示すことから頭文字を取りARCHモデルと名付けられている。 u t 2 {\displaystyle u_{t}^{2}} {\displaystyle u_{t}^{2}} についての定常性条件から次の z {\displaystyle z} {\displaystyle z} についての方程式

1 i = 1 q α i z i = 0 {\displaystyle 1-\sum _{i=1}^{q}\alpha _{i}z^{i}=0} {\displaystyle 1-\sum _{i=1}^{q}\alpha _{i}z^{i}=0}

の全ての解の絶対値が1より大きくなるように係数 α i , i = 1 , , q {\displaystyle \alpha _{i},\;i=1,\dots ,q} {\displaystyle \alpha _{i},\;i=1,\dots ,q} に条件が課される場合が多い。

GARCH(p,q)モデル

[編集 ]

1986年にロバート・エングルの弟子Tim Bollerslev (英語版)はARCHモデルを一般化したGARCHモデル(がーちモデル、: Generalized ARCH model, GARCH model)を提案した[2] 。GARCHモデルでは、条件付ボラティリティ σ t {\displaystyle \sigma _{t}} {\displaystyle \sigma _{t}} は以下のように決定される。

σ t 2 = α 0 + α 1 u t 1 2 + + α q u t q 2 + β 1 σ t 1 2 + + β p σ t p 2 = α 0 + i = 1 q α i u t i 2 + i = 1 p β i σ t i 2 . {\displaystyle \sigma _{t}^{2}=\alpha _{0}+\alpha _{1}u_{t-1}^{2}+\cdots +\alpha _{q}u_{t-q}^{2}+\beta _{1}\sigma _{t-1}^{2}+\cdots +\beta _{p}\sigma _{t-p}^{2}=\alpha _{0}+\sum _{i=1}^{q}\alpha _{i}u_{t-i}^{2}+\sum _{i=1}^{p}\beta _{i}\sigma _{t-i}^{2}.} {\displaystyle \sigma _{t}^{2}=\alpha _{0}+\alpha _{1}u_{t-1}^{2}+\cdots +\alpha _{q}u_{t-q}^{2}+\beta _{1}\sigma _{t-1}^{2}+\cdots +\beta _{p}\sigma _{t-p}^{2}=\alpha _{0}+\sum _{i=1}^{q}\alpha _{i}u_{t-i}^{2}+\sum _{i=1}^{p}\beta _{i}\sigma _{t-i}^{2}.}

すなわち、現在の条件付ボラティリティは p 期前までの条件付ボラティリティと q 期前までの平均からの乖離部分の2乗により決定される。Bollerslev も当該論文中の実証分析の節で述べているが、ARCHモデルを金融時系列データに適用すると分散の長期記憶性を再現する為に次数 q が大きくなる傾向があったが、GARCHモデルは比較的小さい次数(p = 1, q = 1)でも十分に分散の長期記憶性が再現されるので、ARCHモデルに比べると倹約的なモデルとなる。GARCHモデルにおいては u t 2 {\displaystyle u_{t}^{2}} {\displaystyle u_{t}^{2}}自己回帰移動平均モデルとして表され、その定常条件は

1 i = 1 max { p , q } ( α i + β i ) z i = 0 {\displaystyle 1-\sum _{i=1}^{\max\{p,q\}}(\alpha _{i}+\beta _{i})z^{i}=0} {\displaystyle 1-\sum _{i=1}^{\max\{p,q\}}(\alpha _{i}+\beta _{i})z^{i}=0}

の全ての解の絶対値が1より大きくなることである。ただし α i = 0 , i > q {\displaystyle \alpha _{i}=0,\;i>q} {\displaystyle \alpha _{i}=0,\;i>q} かつ β i = 0 , i > p {\displaystyle \beta _{i}=0,\;i>p} {\displaystyle \beta _{i}=0,\;i>p} である。

GARCHモデルの拡張

[編集 ]

GARCHモデルは様々な拡張がなされている。以下で代表的なものを述べる。

EGARCHモデル

[編集 ]

Daniel B. Nelson が1991年に提案したExponential GARCH(p,q)モデル(EGRACH(p,q)モデル)は以下のようにボラティリティが決定する[3]

log σ t 2 = ω + i = 1 p β i log σ t i 2 + i = 1 q ( α i ε t i + γ i ( | ε t i | E [ | ε t i | ] ) ) {\displaystyle \log \sigma _{t}^{2}=\omega +\sum _{i=1}^{p}\beta _{i}\log \sigma _{t-i}^{2}+\sum _{i=1}^{q}{\Big (}\alpha _{i}\varepsilon _{t-i}+\gamma _{i}{\Big (}|\varepsilon _{t-i}|-E[|\varepsilon _{t-i}|]{\Big )}{\Big )}} {\displaystyle \log \sigma _{t}^{2}=\omega +\sum _{i=1}^{p}\beta _{i}\log \sigma _{t-i}^{2}+\sum _{i=1}^{q}{\Big (}\alpha _{i}\varepsilon _{t-i}+\gamma _{i}{\Big (}|\varepsilon _{t-i}|-E[|\varepsilon _{t-i}|]{\Big )}{\Big )}}

EGARCHモデルにおいては通常のGARCHモデルと異なり、 u t i {\displaystyle u_{t-i}} {\displaystyle u_{t-i}} ではなく、それを σ t i {\displaystyle \sigma _{t-i}} {\displaystyle \sigma _{t-i}} で割った ε t i {\displaystyle \varepsilon _{t-i}} {\displaystyle \varepsilon _{t-i}} がボラティリティに影響を与える。条件付き分散の対数に対してモデル化が行われているため、通常のGARCHモデルに比べると非負性や定常性のための制約が緩くなるという利点がある。

GJR GARCHモデル

[編集 ]

Lawrence R. Glosten, Ravi Jagannathan (英語版), David E. Runkle によって1993年に提案されたGJR GARCHモデルは以下のようにボラティリティが決定する[4]

σ t 2 = ω + α u t 1 2 + β σ t 1 2 + γ u t 1 2 I t 1 {\displaystyle \sigma _{t}^{2}=\omega +\alpha u_{t-1}^{2}+\beta \sigma _{t-1}^{2}+\gamma u_{t-1}^{2}I_{t-1}} {\displaystyle \sigma _{t}^{2}=\omega +\alpha u_{t-1}^{2}+\beta \sigma _{t-1}^{2}+\gamma u_{t-1}^{2}I_{t-1}}

ただし、 I t 1 {\displaystyle I_{t-1}} {\displaystyle I_{t-1}} u t 1 {\displaystyle u_{t-1}} {\displaystyle u_{t-1}} が負ならば1、正ならば0を取る変数である。株価収益率などが持つ、下落局面でボラティリティがより増加するレバレッジ効果を捉えるためのモデルである。

Heston-Nandi GARCH モデル

[編集 ]

Steven L. Heston (英語版), Saikat Nandi により2000年に提案されたHeston-Nandi GARCH(p,q)モデルは以下のようにボラティリティが決定する[5]

σ t 2 = ω + i = 1 p β i σ t i 2 + i = 1 q α i ( ε t i γ i σ t i ) 2 {\displaystyle \sigma _{t}^{2}=\omega +\sum _{i=1}^{p}\beta _{i}\sigma _{t-i}^{2}+\sum _{i=1}^{q}\alpha _{i}{\Big (}\varepsilon _{t-i}-\gamma _{i}\sigma _{t-i}{\Big )}^{2}} {\displaystyle \sigma _{t}^{2}=\omega +\sum _{i=1}^{p}\beta _{i}\sigma _{t-i}^{2}+\sum _{i=1}^{q}\alpha _{i}{\Big (}\varepsilon _{t-i}-\gamma _{i}\sigma _{t-i}{\Big )}^{2}}

Heston-Nandi GARCHモデルもEGARCHモデルと同様に u t i {\displaystyle u_{t-i}} {\displaystyle u_{t-i}} ではなく ε t i {\displaystyle \varepsilon _{t-i}} {\displaystyle \varepsilon _{t-i}} がボラティリティに影響を与える。また、このモデルもGJR GARCHモデルと同様にレバレッジ効果を捉えることができる。さらにデリバティブオプションと親和性が高く、Heston-Nandi GARCHモデルに従う株式のオプションについて、その無裁定価格が導出されている。しかし、Heston-Nandi GARCHモデルはモデルが過適合を起こしやすいという欠点もある。

多変数モデルへの拡張

[編集 ]

ここまで述べてきたGARCHモデルはいずれも単一変数の時系列データに対して適用されるものであったが、多変数の時系列データに対してその相関構造を内包しつつ適用可能なGARCHモデルも存在する。例としてBEKKモデル[6] やCCC-GARCHモデル[7] 、DCC-GARCHモデル[8] などがある。

脚注

[編集 ]
  1. ^ Engle 1982
  2. ^ Bollerslev 1986
  3. ^ Nelson 1991
  4. ^ Glosten, Jagannathan and Runkle 1993
  5. ^ Heston and Nandi 2000
  6. ^ Engle and Kroner 1995
  7. ^ Bollerslev 1990
  8. ^ Engle 2002

参考文献

[編集 ]

関連項目

[編集 ]
標本調査
記述統計学
連続データ
位置
分散
モーメント
カテゴリデータ
推計統計学
仮説検定
パラメトリック
ノンパラメトリック
その他
区間推定
モデル選択基準
その他
ベイズ統計学
確率
その他
相関
相関係数
その他
モデル
回帰
線形
非線形
時系列
分類
線形
二次
非線形
その他
教師なし学習
クラスタリング
密度推定 (英語版)
その他
統計図表
生存時間分析
歴史
応用
出版物
全般
その他
カテゴリ カテゴリ

AltStyle によって変換されたページ (->オリジナル) /