以下の書籍
を中心に時系列解析を勉強していきます。
前回
5. VARモデル
5.1 Granger因果性
変数間の因果性の有無を判断することをデータから(-
)行うべく、予測を基準とした因果性を
は提案した:
因果性は現在と過去の変数値
が将来の
に関して
の現在と過去の値以上の情報を有しているかどうかを基準としている*1。
因果性が良く用いられるのは、
モデルの枠組みでは
因果性の有無が
検定により判断できるからである。もし
因果性が存在しないならば、
モデルにおいてある変数の項(方程式)において別変数の過去値に対する係数はすべて
になる。そこで
検定を用いればよい。
モデルにおける
のモデルを最小二乗法により推定し、その残差平方和を
とする。
モデルにおける
のモデルに制約を課したモデルを最小二乗法により推定し、その残差平方和を
とする。
検定量
で計算する。ここでは
因果性検定に必要な制約の数である。
を
%点と比較し、
の方が大きければ、ある変数(群)から
への
因果性は存在し、小さければ
因果性は存在しないと結論付ける。
因果性を理解する上では、
による分布ラグモデルからの解釈も重要である。ある時系列
が別の時系列
と期待値
の誤差項
を用いて
という形で表現できるとき、は
の分布ラグモデルに従うという。
分布ラグモデルの特徴はどの時点を考えても説明変数
が誤差項
と無相関であることである。これは
将来の
に関して現在と過去の
が持つ以上の情報を有さないことを意味する。逆に
から
への
因果性が存在しない場合、
が
の分布ラグモデルで表現できることも示すことができる。
因果性は実際の因果性とは逆の方向に観測される可能性もあり、経済理論などにより因果関係の方向性が明確な場合を除いては、定義どおりに予測に有用か否かに用いるのが望ましい
因果性は定性的概念であり、その強度に関しては一切の情報を持たない。
5.4 インパルス応答関数
因果性は定量的な解析ができない点がデメリットであった。そこを補うのがインパルス応答関数である。
一般的な変量
モデル
を考え、は対角行列でないとする。このとき、非直交化インパルス応答関数を定義できる、
5.4.1 非直交化インパルス応答関数
非直交化インパルス応答関数
数学的には、(1単位の変化に対する)非直交化インパルス応答関数は
と表すことができる。インパルス応答分析では攪乱項の確定的な変化に対する
の確定的な変化を調べている。非直交化インパルス応答関数を得るには、
から逐次的に計算していくのが1つの方法である。
非直交化インパルス応答関数には問題が1点ある。それは攪乱項が同時点で互いに相関しているにもかかわらず、その相関を無視している点である。その点を改善したのが次節で説明する直交化インパルス応答関数である。
5.4.2 直交化インパルス応答関数
非直交化インパルス応答関数の問題点である、同時点における攪乱項の相関を無視している点を解消するのに、攪乱項の分散共分散行列が修正分解可能だと仮定する方法がある。
直交化インパルス応答関数
一般にインパルス応答関数というときは、直交化インパルス応答関数を指す。
攪乱項の分散共分散行列を
とするとき、これは正定値行列であるから、
を対角成分が
に等しい下三角行列、
を対角行列として
と書ける。このとき直交化攪乱項は
と定義され、これは各変数固有の変動を表す。
から、直交化攪乱項は確かに互いに無相関である。
のショックに対する
のインパルス応答関数は
に
単位のショックを与えたときの
の反応を時間の関数として見たときで、
と表すことができる。他にも分散共分散行列を分解して表すこともできる。
直交化攪乱項の計算では、修正分解することで
が成立するから、元々の攪乱項は直交化攪乱項に対角成分が
に等しい下三角行列を掛けた形で表現できることが分かる。このため、元来の攪乱項はその添字以下の添字を持つ直交化攪乱項の線形和で表される。このことは、攪乱項に与えた添字によっては、影響を受けない直交化攪乱項が発生しなくなることを意味する。このような構造を再帰的構造という。このためインパルス応答関数を計算する場合には、変数を外生性の高い順に並べる必要があることに注意しなければならない。
参考文献
*1:このように予測という観点から評価しているもので、いわゆる「原因=結果」という意味合いでの因果関係を有するか否かを判断するものではないことはしっかりと押さえておかなければならない。