見せかけの回帰のソースを表示

{{Otheruses|見せかけの回帰|見せかけの相関|擬似相関}}
[[File:Spurious Regression.svg|thumb|400px|[[ランダムウォーク]]過程の見せかけの回帰。シミュレーションで発生させた互いに[[確率論的独立性|独立]]なランダムウォークを回帰したもの。独立なのにも関わらず非常に高い[[t検定]]統計量の値となっている。]]
'''見せかけの回帰'''（みせかけのかいき、{{lang-en-short|spurious regression}}）とは、[[統計学]]や[[計量経済学]]において、統計的に[[確率論的独立性|独立]]である無関係の二つの[[時系列]]変数が[[最小二乗法]]による[[回帰分析]]において統計的に有意な係数の推定値を取ってしまうという問題である。[[クライヴ・グレンジャー]]と{{仮リンク|ポール・ニューボールド|en|Paul Newbold}}によって[[1974年]]に[[モンテカルロ法]]を用いたシミュレーションで発見され<ref name="GrangerNewbold">{{Harvnb|Granger and Newbold|(1974)|ref=Granger,Newbold1974}}</ref>、{{仮リンク|ピーター・フィリップス (統計学者)|en|Peter C. B. Phillips}}によって[[1986年]]に理論的に示された<ref name="Phillips">{{Harvnb|Phillips|(1986)|ref=Phillips1986}}</ref>。[[単位根]]過程と呼ばれる時系列変数同士の回帰分析によって起こる問題であり、単位根過程は経済データなどで頻繁にみられるため、1980年代以降の計量経済学における時系列分析では常に注意が払われる問題となっている。

==概要==
見せかけの回帰は時系列変数の回帰において説明変数と被説明変数が双方とも[[単位根]]過程である時に発生しうる。説明変数と被説明変数が共に単位根過程ならば、二つの変数は[[確率論的独立性|独立]]、つまり何の因果関係も相関関係もないとしても、その回帰係数は通常の手続き<ref>[[t検定]]や[[F検定]]</ref>において統計的に有意<ref>ここでの統計的有意性とは係数が0であるという[[帰無仮説]]が棄却されるという意味で用いている。</ref>となる場合が多い。これを見せかけの回帰と呼ぶ。単位根過程は例えば[[ランダムウォーク]]過程やドリフト付きランダムウォーク過程など時系列データにおいて頻繁に見られるため、見せかけの回帰を考慮しなければ本来無関係であるはずの二つの変数に対して誤った関係を見出しかねない。

[[擬似相関]]は観測されない、もしくは推計の対象としていない第三の変数がそれぞれの変数に対して影響を与えるという意味では、因果関係は存在しないが相関関係は存在する。つまり、説明変数は被説明変数の（無条件の）予測としては機能しうる。しかし、見せかけの回帰は因果関係も[[相関関係]]も存在しない確率論的に独立である二変数に対しあたかも関係があるかのように推定されてしまう場合があるということを述べている。

例えば、経済データにおける単位根である可能性が疑われる系列の例として[[株価]]などの資産価格、[[GDP]]や消費、雇用者数などのマクロ集計量などがある。これらの変数同士の回帰は、仮に有意な関係が得られたとしても、見せかけの回帰である可能性も否定できない。

==数式での表現==
見せかけの回帰は以下のようにして数式として表現される<ref name="Phillips" />。次のような2つの[[ランダムウォーク]]過程を考える。
:<math> x_t = x_{t-1} + w_t </math>
:<math> y_t = y_{t-1} + v_t </math>
ただし、<math>w_t,v_t</math> は[[独立同分布|IID]]（独立かつ同一分布に従う）で
:<math>\operatorname{E}[w_t] = \operatorname{E}[v_t] = 0</math>
:<math>\operatorname{Var}(w_t) = \sigma_w^2,\quad \operatorname{Var}(v_t) = \sigma_v^2</math>
であるとする。また <math>w_t</math> と <math>v_t</math> は互いに独立であると仮定する。この仮定により <math>x_t</math> と <math>y_t</math> もまた互いに独立な無関係である変数となることが分かる。

ここで次のような線形回帰式を[[最小二乗法]]で推定する。
:<math>y_t = \alpha + \beta x_t + \epsilon_t </math>
<math>x_t</math> と <math>y_t</math> には何の関係もないので <math>\beta = 0</math> である。しかし、この回帰式の、<math>T</math> 個の観測値においての最小二乗法による <math>\beta</math> 推定量 <math>\widehat\beta_T</math> は以下の確率変数に[[確率変数の収束#分布収束|分布収束]]する。
:<math>\widehat\beta_T\ \xrightarrow{d}\ \frac{\sigma_v}{\sigma_w}\frac{\int_0^1V(t)W(t)d t - \int_0^1V(t)d t\int_0^1W(t)d t}{\int_0^1W(t)^2d t - \left(\int_0^1W(t)d t\right)^2} </math>
ここで <math>W(t),\ V(t)</math> は互いに独立な[[ブラウン運動]]である。そして[[帰無仮説]] <math>\beta = 0</math> についての[[t検定]]統計量 <math>t_\beta</math> はある確率変数 <math>Z</math> に対して
:<math>\frac{t_\beta}{\sqrt{T}}\ \xrightarrow{d}\ Z </math>
となる。検定統計量 <math>t_\beta</math> を <math>\sqrt{T} </math> で割ったものがある確率変数 <math>Z</math> に分布収束するので、検定統計量そのものはマイナス無限大か無限大に発散する。よって観測個数 <math>T</math> が十分大きければ検定統計量 <math>t_\beta</math> は棄却値を上回ることが予想される。

ここで問題となっているのは真の値 <math>\beta = 0</math> の下で誤差項 <math>\epsilon_t</math> が非[[定常過程]]（[[単位根]]過程）となっていることである。

==共和分==
{{Main|共和分}}
見せかけの回帰は無関係な単位根過程同士に関係が存在するような推定結果が得られてしまうという問題であるが、実際に単位根過程がある（線形）方程式を常に成立させるように変動している場合がある。このような関係を'''共和分'''（{{lang-en-short|cointegration}}）と呼ぶ。単位根過程間での関係は見せかけの回帰を生じるさせるようなものか、共和分であるかのいずれかなので、単位根過程間での有意な関係の有無を調べたい時には見せかけの回帰を起こすリスクが存在する最小二乗法ではなく、共和分関係の推定と共和分検定を用いることが推奨される。

==見せかけの回帰の回避==
見せかけの回帰を回避する方法はいくつか提案されている。
===ラグ付きの回帰===
1期前の説明変数、もしくは被説明変数を説明変数として加える事で見せかけの回帰は回避できる<ref name="Hamilton">{{Harvnb|Hamilton|(1994)|pp=561–562|ref=Hamilton1994}}</ref>。つまり、単位根過程 <math>x_t,\ y_t</math> に対し
:<math>y_t = \alpha + \beta_1 x_t + \beta_2 x_{t-1} + \beta_3 y_{t-1} + \epsilon_t </math>
という回帰式で最小二乗法を適用する。例えば、[[#数式での表現|数式での表現]]で定義されているように <math>x_t,\ y_t</math> がランダムウォークならば、<math>\beta_3 = 1,\ \beta_1 = \beta_2 = 0</math> であり、結果として誤差項が <math>v_t = \epsilon_t</math> を満たし、定常過程となるので通常の最小二乗法で係数が一致推定できる。しかし、仮説検定においては <math>y_t</math> の非定常性により通常の手続きと異なる方法を取る必要がある。

===差分を取った回帰===
説明変数と被説明変数の差分に対する関係を見る方法もある<ref name="Hamilton"/>。つまり、単位根過程 <math>x_t,\ y_t</math> に対し
:<math>y_t-y_{t-1} = \alpha + \beta (x_t-x_{t-1}) + \epsilon_t </math>
という回帰式で最小二乗法を適用する。ここで <math>x_t,\ y_t</math> が1次の単位根過程ならば説明変数と被説明変数が共に定常過程となるので通常と全く同じ方法で一致推定や仮説検定なども行える。この簡便さや扱いやすさから差分を取った回帰は多くの研究者によって用いられているが欠点もある。一つが <math>x_t,\ y_t</math> が単位根過程でない、つまり定常過程の場合である。もしそうならば最小二乗法による推定は誤った結果を導きかねない。もう一つが <math>x_t,\ y_t</math> に共和分の関係が存在する場合である。この場合も差分を取った回帰では共和分関係を特定する事はできず、間違った結果を導きかねない。

==脚注==
<references/>

==参考文献==
* {{cite journal 
|last=Granger 
|first=Clive W. J. 
|last2=Newbold 
|first2=Paul 
|authorlink1 = クライヴ・グレンジャー
|year=1974 
|title=Spurious Regressions in Econometrics 
|journal=Journal of Econometrics
|volume=2 
|issue=2 
|pages=111–120 
|url=http://wolfweb.unr.edu/homepage/zal/STAT758/Granger_Newbold_1974.pdf|accessdate=2016-02-07|format=PDF
|doi=10.1016/0304-4076(74)90034-7
|ref=Granger,Newbold1974}}
* {{Citation 
|last=Hamilton
|first=James D.
|authorlink = ジェームズ・ハミルトン (計量経済学者)
|year=1994 
|title=Time Series Analysis
|isbn = 0691042896
|publisher = [[Princeton University Press]]
|ref=Hamilton1994}}
* {{cite journal 
|last=Phillips 
|first=Peter C. B.  
|year=1986 
|title=Understanding Spurious Regressions in Econometrics 
|journal=Journal of Econometrics
|volume=33
|issue=3 
|pages=311–340
|url=http://cowles.yale.edu/sites/default/files/files/pub/d07/d0757.pdf|accessdate=2016-02-07|format=PDF
|doi=10.1016/0304-4076(86)90001-1
|ref=Phillips1986}}

==関連項目==
* [[単位根]]
* [[共和分]]

{{デフォルトソート:みせかけのかいき}}
[[Category:統計学]]
[[Category:計量経済学]]
[[Category:時系列分析]]
[[Category:数学に関する記事]]