一般化モーメント法のソースを表示

'''一般化モーメント法'''（いっぱんかモーメントほう、{{lang-en-short|generalized method of moments, GMM}}）とは、[[計量経済学]]において統計モデルのパラメーターを推定するための一般的な方法である。

一般化モーメント法においては、モデルについてのいくつかのモーメント条件が特定されている必要がある。これらのモーメント条件はモデルのパラメーターとデータの関数である。例えば、真のパラメーターの下で期待値が0となるようなものがある。この時、一般化モーメント法はモーメント条件の標本平均のある[[ノルム]]を最小化する。

一般化モーメント法による推定量は一致性、漸近正規性を持つことが知られ、さらにモーメント条件以外の情報を使わないすべての推定量のクラスにおいて統計的に効率的であることも知られている。

一般化モーメント法は[[ラース・ハンセン]]により[[1982年]]に、[[カール・ピアソン]]が[[1894年]]に導入したモーメント法の一つの一般化として提案された。ハンセンは一般化モーメント法とそれのファイナンスへの応用により[[2013年]]の[[ノーベル経済学賞]]を受賞した。

== 概要 ==
利用可能なデータは ''T'' 個の観測値 {{nowrap|{''Y<sub>t</sub>'' }<sub> ''t'' {{=}} 1,...,''T''</sub>}} からなると仮定する。ここでそれぞれの観測値 ''Y<sub>t</sub>'' は ''n'' 次元の多次元[[確率変数]]であるとする。ここでこのデータはある統計モデルから生成されるとし、その統計モデルは未知パラメーター {{nowrap|''θ'' ∈ Θ}} によって定義されるものとする。この推定問題の目的は真のパラメーター ''θ''<sub>0</sub> もしくは少なくとも適度に近い推定量を見つけることである。

一般化モーメント法の一般的な仮定はデータ ''Y<sub>t</sub>'' が{{仮リンク|弱定常|en|Stationary process}}かつ{{仮リンク|エルゴード|en|Ergodic process}}な[[確率過程]]であることである（[[独立同分布|独立かつ同一分布]]に従う確率変数 ''Y<sub>t</sub>'' はこの条件の特殊ケースである）。

一般化モーメント法を適用する為に、モーメント条件を特定する必要がある。つまり以下のようなベクトル値関数 ''g''(''Y'',''θ'') が既知でなくてはならない。
: <math>
    m(\theta_0) \equiv \operatorname{E}[\,g(Y_t,\theta_0)\,]=0,
  </math>
ここで E は[[期待値]]、''Y<sub>t</sub>'' は一般的な観測値を表す。加えて関数 ''m''(''θ'') は {{nowrap|''θ'' ≠ ''θ''<sub>0</sub>}} ならば0と異なる値を取らなくてはならない。そうでなければパラメーター ''θ'' は識別不可能である。

一般化モーメント法の基本的なアイデアは理論的な期待値 E[⋅] を実証的なもの、つまり[[標本平均]]に置き換えることである。
: <math>
    \hat{m}(\theta) \equiv \frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)
  </math>
そして、この時、この表現のあるノルムを ''θ'' について最小化する。ノルムを最小化する ''θ'' が ''θ''<sub>0</sub> の推定量である。

[[大数の法則]]により、十分大きな ''T'' について <math style="vertical-align:-.37em">\scriptstyle\hat{m}(\theta)\,\approx\;\operatorname{E}[g(Y_t,\theta)]\,=\,m(\theta)</math> であり、 よって <math style="vertical-align:-.3em">\scriptstyle\hat{m}(\theta_0)\;\approx\;m(\theta_0)\;=\;0</math> が成り立つことが予想される。一般化モーメント法はできるだけ <math style="vertical-align:-.3em">\scriptstyle\hat{m}(\;\!\hat\theta\;\!)</math> を0に近づけるような <math style="vertical-align:0">\scriptstyle\hat\theta</math> を探す。数学的にはこの方法は <math style="vertical-align:-.3em">\scriptstyle\hat{m}(\theta)</math> のある[[ノルム]]を最小化することと同値である（''m'' のノルムを ||''m''|| と表し、''m'' とゼロの間の距離を測るものとする)。結果として得られた推定量の持つ性質はノルム関数の選択にもよるので、ゆえに一般化モーメント法の理論はノルム全体の族を考慮する。以下を定義する。
: <math>
    \| \hat{m}(\theta) \|^2_{W} = \hat{m}(\theta)'\,W\hat{m}(\theta),
  </math>
ここで ''W'' は[[定符号二次形式|正値定符号]]である加重行列で ''m′'' は転置を表す。実践上、加重行列 ''W'' は利用可能なデータセットに基づいて計算され、そのようにして計算された加重行列を <math style="vertical-align:0em">\scriptstyle\hat{W}_T</math> とする。よって一般化モーメント法による推定量は以下のように書ける。
: <math>
    \hat\theta = \operatorname{arg}\min_{\theta\in\Theta} \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg)' \hat{W}_T \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg)
  </math>
適切な条件の下で、一般化モーメント法による推定量は一致性と漸近正規性を持つ。そして加重行列 <math style="vertical-align:0em">\scriptstyle\hat{W}_T</math> を正しく選択すれば効率的な推定量となる。

== 性質 ==
=== 一致性 ===
一致性とは、推定量の持つ統計的な性質であり、十分に多くの観測値がある場合、推定量は真の値に任意に近づいていくということである。
: <math>\hat\theta \xrightarrow{p} \theta_0\ \text{as}\ T\to\infty</math>
（[[確率変数の収束#確率収束|確率収束]]を参照）。一般化モーメント法による推定量が一致性を持つ必要十分条件は以下の通りである。
# <math>\hat{W}_T \xrightarrow{p} W</math> を満たす。ただし ''W'' は[[定符号二次形式|正値定符号]]行列である。
# <math>\,\theta=\theta_0</math> である時に限り <math>\,\operatorname{E}[\,g(Y_t,\theta)\,]=0</math> を満たす。
# パラメーターが値を取りうる集合 <math>\Theta \subset \mathbb{R}^{k}</math> は[[コンパクト (数学)|コンパクト]]集合である。
# <math>\,g(Y,\theta)</math> は ''θ'' について[[連続]]である。
# <math>\operatorname{E}[\,\textstyle\sup_{\theta\in\Theta} \lVert g(Y,\theta)\rVert\,]<\infty</math>

第二の条件（'''大域的識別'''条件と呼ばれる）は、非線形モデルにおいて確かめるのが難しい。

実証計量経済学者は実際に大域的識別条件を確かめずに、それが成立していると単に仮定することがしばしばある<ref>{{harvtxt|Newey|McFadden|1994|p=2127}}</ref>

が、あまり推奨されない。識別条件が成立しない非線形モデルの例については、Dominguez and Lobato (2004)を参照のこと。

=== 漸近正規性 ===
漸近正規性は有用な性質であり、漸近正規性により推定量の[[信頼区間]]を計算することや[[仮説検定]]を行うことができる。一般化モーメント法による推定量の漸近分布について述べる前に、以下の2つの補助的な行列を定義する。
: <math>G = \operatorname{E}[\,\nabla_{\!\theta}\,g(Y_t,\theta_0)\,], \qquad
        \Omega = \operatorname{E}[\,g(Y_t,\theta_0)g(Y_t,\theta_0)'\,]</math>
以下の1から6までの条件の下で、一般化モーメント法による推定量は漸近正規性を持つ。
: <math>\sqrt{T}\big(\hat\theta - \theta_0\big)\ \xrightarrow{d}\ \mathcal{N}\big[0, (G'WG)^{-1}G'W\Omega W'G(G'W'G)^{-1}\big]</math>
（[[確率変数の収束#分布収束|分布収束]]を参照）。条件は以下の通りである。
# <math>\hat\theta</math> は一致性を持つ。
#<math>\,g(Y,\theta)</math> は <math>\theta_0</math> のある近傍 ''N'' において連続微分可能である。
# <math>\operatorname{E}[\,\lVert g(Y_t,\theta) \rVert^2\,]<\infty</math>
# <math>\operatorname{E}[\,\textstyle\sup_{\theta\in N}\lVert \nabla_\theta g(Y_t,\theta) \rVert\,]<\infty</math>
# 行列 <math>G'WG</math> は[[正則行列]]である。

=== 効率性 ===
ここまで行列 ''W'' の選択については、それが半正値定符号で無くてはならないということを除き何も述べてこなかった。実際、どのような半正値定符号行列であっても一般化モーメント法による推定量は一致性と漸近正規性を持つ。唯一の違いはその推定量の漸近分散にある。加重行列を以下のように取る。
: <math> W \propto\ \Omega^{-1} </math>
すると、一般化モーメント法による推定量はすべての漸近正規的な推定量の中で最も効率的となる。この場合の効率性は、推定量が可能な限り最小の分散行列（ある行列 ''A'' が行列 ''B'' より小さいとは ''B-A'' が半正値低符号であるということである。）を持つという意味である。

この場合、一般化モーメント法による推定量の漸近分散についての公式は以下のように単純化される。
: <math>\sqrt{T}\big(\hat\theta - \theta_0\big)\ \xrightarrow{d}\ \mathcal{N}\big[0, (G'\,\Omega^{-1}G)^{-1}\big]</math>

このような加重行列を選ぶことが最適になるという証明は、しばしば他の推定量の効率性を証明する時の証明を少しばかり模倣したものを取り入れる。大雑把に言えば、加重行列を分散についての"サンドイッチ公式"が単純な表現になるように選べば、その加重行列は最適となる。

{| style="margin-left:30pt"
| colspan="2" | '''''証明''''' 加重行列を任意の ''W'' とした時と <math>W=\Omega^{-1}</math> とした時の漸近分散の差について考える。もし、その差がある行列 ''C'' についての対称な積の形式 ''CC''' に分解できれば、それはその差が非負値定符号であることを意味し、ゆえに定義より <math>W=\Omega^{-1}</math> は最適になる。
|-
| <math>\,V(W)-V(\Omega^{-1})</math>
|    <math>\,=(G'WG)^{-1}G'W\Omega WG(G'WG)^{-1} - (G'\Omega^{-1}G)^{-1}</math>
|-
| || <math>\,=(G'WG)^{-1}\Big(G'W\Omega WG - G'WG(G'\Omega^{-1}G)^{-1}G'WG\Big)(G'WG)^{-1}</math>
|-
| || <math>\,=(G'WG)^{-1}G'W\Omega^{1/2}\Big(I - \Omega^{-1/2}G(G'\Omega^{-1}G)^{-1}G'\Omega^{-1/2}\Big)\Omega^{1/2}WG(G'WG)^{-1}</math>
|-
| || <math>\,=A(I-B)A',</math>
|-
| colspan="2" |
ここで行列 ''A'' と ''B'' を記法の単純化のために導入している。''I'' は[[単位行列]]である。行列 ''B'' は対称かつ[[冪等]]な行列であることが分かる。これは ''I-B'' もまた対称かつ冪等であることを意味する。つまり <math>I-B=(I-B)(I-B)'</math> が成り立つ。ゆえに以前の表現を以下のように分解することが可能である。
|-
| || <math>\,=A(I-B)(I-B)'A' = \Big(A(I-B)\Big)\Big(A(I-B)\Big)' \geq 0</math>
|}

== 実装 ==
今まで述べてきた方法を実装するにあたっての一つの難しい点は {{nowrap|''W'' {{=}} Ω<sup>−1</sup>}} として加重行列を取ることである。なぜならば Ω の定義より、それを計算するためには ''θ''<sub>0</sub> の値が既知でなければならず、''θ''<sub>0</sub> はまさに未知であり、そもそも推定しようとしている量である。

この問題を解決するための方法がいくつか存在する。以下であげるものうち、2段階GMMが最も一般的である。

*'''2段階GMM'''（{{lang-en-short|Two-step GMM}}）
** ステップ1　<math style="vertical-align:-.5em">\scriptstyle\hat{W}_T\;=I</math> （[[単位行列]]）とし、事前の一般化モーメント法による推定量 <math style="vertical-align:-.5em">\scriptstyle\hat\theta_{(1)}</math> を計算する。この推定量は ''θ''<sub>0</sub> についての一致推定量ではあるが、効率的ではない。
** ステップ2
:: <math>\hat{W}_T = \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\hat\theta_{(1)})g(Y_t,\hat\theta_{(1)})'\bigg)^{-1},</math>
::とする。ただし、ステップ1における推定量 <math style="vertical-align:-.5em">\scriptstyle\hat\theta_{(1)}</math> を用いた。この行列は Ω<sup>−1</sup> に確率収束し、ゆえにこの加重行列を用いて推定量 <math style="vertical-align:0">\scriptstyle\hat\theta</math> を計算すれば、その推定量は漸近的に効率的である。

*'''繰り返しGMM'''（{{lang-en-short|Iterated GMM}}）

:行列 <math>\hat{W}_T</math> を複数回計算することを除けば、本質的には2段階GMMと同じ方法である。つまりステップ2で得た推定量を加重行列として再び用いて推定量を計算し、これを繰り返す。このような推定量は、<math style="vertical-align:-.5em">\scriptstyle\hat\theta_{(i)}</math> と記すが、以下のシステム方程式を解いた場合と同値になる<ref>{{harvtxt|Imbens|Spady|Johnson|1998|p=336}}</ref>。
: <math>
    \bigg(\frac{1}{T}\sum_{t=1}^T \frac{\partial g}{\partial\theta'}(Y_t,\hat\theta_{(i)})\bigg)' \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\hat\theta_{(i)})g(Y_t,\hat\theta_{(i)})'\bigg)^{\!-1} \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\hat\theta_{(i)})\bigg) = 0
  </math>
:このような繰り返しを行っても漸近的な改善は達成できないが、あるモンテカルロ実験では有限標本における推定量の振る舞いが若干よくなる{{要出典|date=2009年6月}}。

*'''連続更新GMM'''（{{lang-en-short|Continuously updating GMM}} CUGMM もしくは CUE）

:<math style="vertical-align:0">\scriptstyle\hat\theta</math> を加重行列 ''W'' と同時に推定する。つまり、
: <math>
    \hat\theta = \operatorname{arg}\min_{\theta\in\Theta} \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg)' \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)g(Y_t,\theta)'\bigg)^{\!-1} \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg)
  </math>
:として推定する。モンテカルロ実験において、この方法は伝統的な2段階GMMよりよいパフォーマンスを見せている。連続更新GMMは（裾が厚くなるが、）中位点のバイアスが小さくなり、そして多くの場合における過剰識別制約のためのJ検定がよりもっともらしい結果となる<ref>{{harvtxt|Hansen|Heaton|Yaron|1996}}</ref>。

最小化の手続きの実装におけるもう一つの重要な問題は、（高次元であることもありうる）パラメーター空間 ''Θ'' を探索し、目的関数を最小化する ''θ'' の値を見つけるということになっているということである。このような手続きについて一般的に推奨される方法は存在せず、それは個々の場合による問題となる（[[数理最適化]]）。

== J検定 ==
モーメント条件の数がパラメーターベクトルの次元より大きい時、そのモデルは''過剰識別されている''（{{lang-en-short|over-identified}}）と言う。過剰識別ならば、そのモデルのモーメント条件がデータと適合するかどうかを調べることが出来る。

概念的に、モデルがデータによくフィットしているかは、<math>\hat{m}(\hat\theta)</math> が十分0に近いかどうかで調べることが出来る。一般化モーメント法は方程式 <math>\hat{m}(\theta)=0</math> を解く問題、つまり <math>\theta</math> が制約を確かに満たすように選ぶという問題を最小化計算に置き換えている。この最小化は <math>m(\theta_0)=0</math> を満たすような <math>\theta_0</math> が存在しないとしても、常に実行可能である。J検定はこの制約が成立しているかを確かめることができる。J検定は''過剰識別制約についての検定''とも呼ばれる。

以下の統計的仮説を考えよう。
* <math>H_0:\ m(\theta_0)=0</math> （モデルが妥当であるという[[帰無仮説]]）
* <math>H_1:\ m(\theta)\neq 0,\ \forall \theta\in\Theta</math> （モデルが妥当でないという[[対立仮説]]。データは制約を満たすほど近づかない。）

仮説 <math>H_0</math> の下で以下のJ検定統計量は漸近的に自由度 ''k-l'' の[[カイ2乗分布]]に従う。
:<math>J \equiv T \cdot \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\hat\theta)\bigg)' \hat{W}_T \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\hat\theta)\bigg)\ \xrightarrow{d}\ \chi^2_{k-\ell}</math> &nbsp; under <math>H_0,</math>
ここで <math>\hat\theta</math> はパラメーター <math>\theta_0</math> の一般化モーメント法による推定量、''k'' はモーメント条件の数（ベクトル ''g'' の次元）、''l'' は推定パラメーターの数（ベクトル ''θ'' の次元）である。行列 <math>\hat{W}_T</math> は <math>\Omega^{-1}</math> に確率収束しなくてはならない。<math>\Omega^{-1}</math> は効率的な加重行列である（以前、推定量が効率的であるためには、''W'' は <math>\Omega^{-1}</math> に比例することだけが必要だった。しかし、J検定を行うには、''W'' は <math>\Omega^{-1}</math> と一致せねばならず、単純に比例するだけではいけない）。

対立仮説 <math>H_1</math> の下で、J検定等計量は漸近的に非有界である。
: <math>J\ \xrightarrow{p}\ \infty</math>&nbsp; under <math>H_1</math>

検定を行う為に、データから ''J'' の値を計算しなくてはならない。''J'' は非負である。''J'' を（例えば）<math>\chi^2_{k-\ell}</math> 分布の95%[[分位点]]と比較する。
* もし <math> J > q_{0.95}^{\chi^2_{k-\ell}} </math> ならば、帰無仮説 <math>H_0</math> は[[有意水準]]5%で棄却される。
* もし <math> J < q_{0.95}^{\chi^2_{k-\ell}} </math> ならば、帰無仮説 <math>H_0</math> は[[有意水準]]5%で棄却できない。

== 用例 ==
他の多くの推定法は一般化モーメント法の意味で解釈できる。
* [[最小二乗法]]（{{lang-en-short|Ordinary least squares, OLS}}）は一般化モーメント法と以下のモーメント条件で同値となる。
:<math>\operatorname{E}[\,x_t(y_t - x_t'\beta)\,]=0</math>
* {{仮リンク|一般化最小二乗法|en|Generalized least squares}}（{{lang-en-short|Generalized least squares, GLS}}）
:<math>\operatorname{E}[\,x_t(y_t - x_t'\beta)/\sigma^2(x_t)\,]=0</math>
* [[操作変数法]]（{{lang-en-short|Instrumental variables regression, IV}}）
:<math>\operatorname{E}[\,z_t(y_t - x_t'\beta)\,]=0</math>
* [[非線形最小二乗法]]（{{lang-en-short|Non-linear least squares, NLS}}）
:<math>\operatorname{E}[\,\nabla_{\!\beta}\, g(x_t,\beta)\cdot(y_t - g(x_t,\beta))\,]=0</math>
* [[最尤法]]（{{lang-en-short|Maximum likelihood estimation, MLE}}）
:<math>\operatorname{E}[\,\nabla_{\!\theta} \ln f(x_t,\theta) \,]=0</math>
ただし、NLSやMLEは上のモーメント条件だけではパラメーターを識別できない場合がある。つまり、上のモーメント条件を満たすパラメーターが複数存在する可能性がある。従って、NLSやMLEの場合、GMMによる推定はできる限り避けることを推奨する。例えば、Dominguez and Lobato (2004)のシュレーションでは、非線形回帰のGMM推定が非常に不安定になることが示されている。

== 実装例 ==
*[[Wikibooks:R Programming/Method of Moments|R Programming wikibook, Method of Moments]]
*[http://cran.r-project.org/web/packages/gmm/gmm.pdf R]
*[https://www.stata.com/manuals13/rgmm.pdf Stata]
*[http://www.eviews.com/EViews7/ev7features.html EViews]
*[http://support.sas.com/documentation/cdl/en/etsug/60372/HTML/default/viewer.htm#etsug_model_sect035.htm SAS]

== 参考文献 ==
{{Reflist}}
{{refbegin}}
* {{Citation
 |last = Faciane Jr. | first = Kirby Adam
 |title = Statistics for Empirical and Quantitative Finance
 |year = 2006
 |publisher = H.C. Baird
 |location = Philadelphia
 |isbn = 0-9788208-9-4
}}
* {{Citation
 |last = Hall | first = Alastair R.
 |title = Generalized Method of Moments (Advanced Texts in Econometrics)
 |year = 2005
 |publisher = Oxford University Press
 |isbn = 0-19-877520-2
}}
* {{Citation
 |authorlink = ラース・ハンセン
 |first = Lars Peter |last = Hansen
 |year = 1982
 |title = Large Sample Properties of Generalized Method of Moments Estimators
 |journal = Econometrica
 |volume = 50
 |issue = 4
 |pages = 1029–1054
 |jstor=1912775
}}
* {{Citation
 |authorlink = ラース・ハンセン
 |first = Lars Peter |last = Hansen
 |year = 2002
 |contribution = Method of Moments
 |title = International Encyclopedia of the Social and Behavior Sciences
 |editor1-last = Smelser | editor1-first = N. J
 |editor2-last = Bates   | editor2-first = P. B
 |publisher = Pergamon
 |location = Oxford
}}
* {{Citation
  | last1 = Hansen | first1 = Lars Peter
  | last2 = Heaton | first2 = John
  | last3 = Yaron  | first3 = Amir
  | authorlink1 = ラース・ハンセン
  | year = 1996
  | title = Finite-sample properties of some alternative GMM estimators
  | journal = Journal of Business & Economic Statistics
  | volume = 14 | issue = 3
  | pages = 262–280
  | jstor = 1392442
  | ref = harv | doi=10.1080/07350015.1996.10524656
  }}
* {{Citation
  | last1 = Imbens  | first1 = Guido W.
  | last2 = Spady   | first2 = Richard H.
  | last3 = Johnson | first3 = Phillip
  | year = 1998
  | title = Information theoretic approaches to inference in moment condition models
  | journal = Econometrica
  | volume = 66 | issue = 2
  | pages = 333–357
  | jstor = 2998561
  | ref = harv
  }}
* {{Citation
  | last1 = Newey | first1 = W
  | last2 = McFadden | first2 = D
  | authorlink2 = ダニエル・マクファデン
  | year = 1994
  | contribution = Large sample estimation and hypothesis testing
  | title = Handbook of Econometrics
  | publisher = Elsevier Science
  | ref = harv
}}
* Special issues of Journal of Business and Economic Statistics: [http://www.jstor.org/stable/i260360 vol. 14, no. 3] and [http://pubs.amstat.org/toc/jbes/20/4 vol. 20, no. 4].
{{refend}}

== 関連項目 ==
* [[最尤法]]
* {{仮リンク|一般化経験尤度|en|Generalized empirical likelihood}}

{{デフォルトソート:いつはんかもーめんとほう}}
[[Category:統計学]]
[[Category:計量経済学]]
[[Category:ノーベル経済学賞]]