クラメール・ラオの限界

推定理論英語版統計学におけるクラメール・ラオの限界(CRB)(クラメールラオのげんかい、: Cramér–Rao bound)(クラメール・ラオの下限(CRLB)クラメール・ラオの不等式Frechet–Darmois–Cramér–Rao 不等式情報不等式とも)とは、ある確率分布の未知母数を推定する不偏推定量には、その分散についてある下限値が存在することを示すものである。名称は、1940年代にそれぞれ独立に推定精度に関する限界を見出した、ハラルド・クラメールC. R.ラオ英語版モーリス・ルネ・フレシェジョルジュ・ダルモア英語版フランス語版にちなむ[1][2][3][4][5][6][7]

最も単純に述べると、『任意の不偏推定量の分散は、 そのフィッシャー情報量の逆数以上になる』というものである。不偏な推定量がこの下限を達成するとき、その推定量は(完全な)有効推定量英語版であるという。この場合、その推定量はあらゆる不偏推定量の中で平均二乗誤差が最小のものとなるため、必然的に最小分散不偏推定量英語版(MVU推定量)にもなる。

しかしながら、どんな不偏推定量を考えても分散が決してクラメール・ラオの下限に到達できないようなケースもある(MVU推定量が存在するときでもこれは起こりえる)。

クラメール・ラオの限界には、不偏でない推定量に対するバージョンもある。不偏性の条件を取り除くことで、推定量の分散・平均二乗誤差が、不偏の場合のクラメール・ラオの下限を「下回る」ようなケースも存在する。推定量の偏り英語版も参照。

主張 編集

ここでは、母数が1つ・推定量が不偏である場合から始めて、いくつかのかなり一般的な場合へと拡張していく。どのバージョンでもある種の正規性の仮定をおくが、それはほとんどの「普通のふるまいをする」確率分布については成り立つものである。この条件については後述する。

母数が1つで推定量が不偏の場合 編集

何らかの確率密度関数   に従って分布する量   の観測値から、未知母数   を推定することを考える。このとき、   に対する任意の不偏な推定量   の分散は、フィッシャー情報量   の逆数以上になる:

 

フィッシャー情報量  

 

と定義される。ここで、  尤度自然対数をとったもの(なお  スコア関数英語版という)で、   は平均を表す。

不偏推定量   の有効度は、推定量の分散がこの下限にどの程度接近しているかを測る指標で、次のように定義される。

 

不偏推定量の分散の下限値を、実際の分散で割った値、ともいえる。クラメール・ラオの下限より   となる。

母数が1つで、母数の関数の値を推定する場合 編集

より一般に、確率変数   の関数   を用いて、母数の関数   を推定することを考える。  であるとする。このときの分散の下限は、

 

ここで     による微分、   はフィッシャー情報量である。

母数が1つで推定量が不偏とは限らない場合 編集

母数   の推定量    だけの偏りがあるとする。

  と置いて前項の結果を使うと、

 

不偏のときの不等式は、   とした特別な場合である。

分散を小さくすることだけを考えるなら、定数関数となる「推定量」をとれば、分散はゼロである。しかし上記の式から、推定量の平均二乗誤差には

 

という下限が存在することになる。ここで、平均二乗誤差の標準的な分解式

 

を用いた。

注意:もし   であれば、不偏のときのクラメール・ラオの下限   を下回ることもある。例えば、後述する例では、   となる。

母数が複数(ベクトル値)の場合 編集

クラメール・ラオの限界を、母数が複数の場合にも拡張しよう。母数ベクトル

 

とし(右肩の T は転置を表す(以下同じ))、それによって決まる確率密度関数   を考える。  は後述の正規性の条件をみたすものとする。 フィッシャー情報行列は、   行列で、その成分  

 

で定まる行列のことである。   を、母数ベクトルの任意の推定量としよう:  。ここで、各成分の平均を並べた平均ベクトル    と記す。

このとき、  分散共分散行列に対するクラメール・ラオの限界は、

 

となる。ここで、

  • 行列に対する不等式   は、行列の差  非負定値であるということである。
  •  ヤコビ行列  成分が   )である。


もし    の不偏推定量であれば(つまり   であれば)クラメール・ラオの限界は

 

のようになる。フィッシャー情報行列の逆行列を計算するのが面倒な場合は、単に対応する対角成分の逆数をとることで、(より緩いかもしれないが)1つの下限が得られる[8]

 

正規性の条件 編集

クラメール・ラオの不等式が成り立つための確率密度関数   と推定量   に関する2つの弱い十分条件は、次のとおりである:

  • フィッシャー情報量が常に定義されていること。言い換えると、次式を   で積分した値が有限値として存在すること。
 
  •   の期待値について、   についての積分と、   についての偏微分が交換可能である、つまり
 
が、右辺が存在する限り成り立つこと。
この条件は、以下のいずれかの場合が成り立つことをもって確認されることが多い:
  1. 関数   は、   に依らない有界な関数の台(非ゼロとなる定義域)を持つ。
  2.   に依らない可積分関数   が存在して   を上から抑える。つまり、
 

フィッシャー情報量の変形 編集

   で2階偏微分可能であるとすると、フィッシャー情報量は

 

(3番目の等号の箇所で

 

であることを用いた)

と変形でき、クラメール・ラオの不等式は次のようにも書ける。

 

こちらの公式の方が下限を評価するのにより有用な場合がある。

母数が1つの場合の証明 編集

母数が1つの場合のクラメール・ラオの不等式を一般的に証明する。

  を、確率密度関数が   となる確率分布に従う確率変数とし、   の関数で、母数   の関数である   の不偏推定量であるとする。つまり、  

目標は、任意の   に対して

 

を示すことである。

  を次のように定義する(これはスコア関数である):

 

ここで連鎖律を使った。  の期待値はゼロである。なぜなら:

 

ここで積分と偏微分の順序が交換可能であること(正規性の条件の2番目)を使った。

  共分散   は、   だから  、よって次式を得る。

 

ここで再び、積分と微分が交換可能であるという条件(正規性の条件の2番目)を使った。

コーシー・シュワルツの不等式から、

 

よって

 

これが示したかったことである。

観測を繰り返し行うときの注意 編集

確率変数列   を使って推定を行う場合について、未知母数が1つ( )のときに絞って概要を述べる。  と書くことにする。

  • 尤度関数は、結合確率密度関数   で与えられる(標本の値   が代入されたとして   の関数とみなしている)。
  • スコア関数は、尤度関数の自然対数をとってから   で偏微分したものである。
 
これらはいずれも実数値関数であるので、
  • フィッシャー情報量も実数値であり、
 
となる。

本記事でここまでに述べた事柄は、次の置き換えをすれば基本的に全て同じ形式で成り立つ。

 

特に、確率変数列  独立同分布で、その確率密度関数が   であるとすると、

  • 尤度関数は  
  • スコア関数は  
  • フィッシャー情報量は
 

となる。

編集

多変量正規分布 編集

平均値ベクトル  、分散共分散行列   が未知母数ベクトル   で定まるような、一般的な d 次元正規分布   の場合、

フィッシャー情報行列の成分は[9]

 

ここで "tr" は行列のトレースを表す。

より簡単な例として、平均   が未知で分散   が既知の正規分布から、独立に   回抽出してえられる標本量ベクトルを   とする。

 

ここで   は 1 を d 個並べたベクトル、 d 次単位行列である。未知母数が1つなのでフィッシャー情報量は

 

とスカラーで与えられ、クラメール・ラオの下限は

 

平均が既知の正規分布の場合 編集

  を、平均   が既知、分散   が未知の正規分布に従う独立な確率変数(列)だとする。次のような統計量を考えよう:

 

このとき   より、    の不偏推定量になる。

  •  分散は、
 
(2番目の等号は分散の定義)。第1項は正規分布の4次の中心モーメントであり、   に等しい。第2項は分散の2乗、つまり   である。よって
 
  • 一方フィッシャー情報量については、まず、観測1回あたりのスコア関数   が尤度関数   から次のように計算できる。
 
最後の等号は簡単な計算でわかる。この情報量は、   をもう一度偏微分してから平均をとり、マイナス1倍したものに等しい。
 
  回の独立な観測の情報量は、これを単純に   倍したものになり、
 

クラメール・ラオの不等式は   だが、この場合は等号が成り立っているため、推定量が有効英語版であることがわかる。

不偏でない推定量を用いれば、分散及び平均二乗誤差をより小さくすることもできる。例えば   とすれば、分散は明らかにより小さくなる。実際

 

ここで偏りは   であり、平均二乗誤差は、『(平均二乗誤差(MSE))=(分散)+(偏りの2乗)』の分解式から

 

となる。こちらも不偏推定量のときの

 

を下回っている。

正規母集団の平均も分散も未知の場合、分散の推定量の平均二乗誤差が最小になるのは、  を平均の推定量として

 

のときである(分母が n − 1 や n + 2 のときではない)。

関連項目 編集

脚注 編集

  1. ^ Cramér, Harald (1946). Mathematical Methods of Statistics. Princeton, NJ: Princeton Univ. Press. ISBN 0-691-08004-6. OCLC 185436716 
  2. ^ Rao, Calyampudi Radakrishna (1945). “Information and the accuracy attainable in the estimation of statistical parameters”. Bulletin of the Calcutta Mathematical Society 37: 81–89. MR0015748. 
  3. ^ Rao, Calyampudi Radakrishna (1994). S. Das Gupta. ed. Selected Papers of C. R. Rao. New York: Wiley. ISBN 978-0-470-22091-7. OCLC 174244259 
  4. ^ Darmois, Georges (1945). “Sur les limites de la dispersion de certaines estimations”. Rev. Int. Inst. Statist. 13: 9–15. 
  5. ^ Fréchet, Maurice (1943). “Sur l'extension de certaines évaluations statistiques au cas de petit échantillons”. Rev. Inst. Int. Statist. 11: 182–205. 
  6. ^ Gart, John J. (1958). “An extension of the Cramér–Rao inequality”. Ann. Math. Stat. 29: 367–380. 
  7. ^ Malécot, Gustave (1999). “Statistical methods and the subjective basis of scientific knowledge [translated from Année X 1947 by Daniel Gianola]”. Genet. Sel. Evol. 31: 269–298. 
  8. ^ For the Bayesian case, see eqn. (11) of Bobrovsky; Mayer-Wolf; Zakai (1987). “Some classes of global Cramer–Rao bounds”. Ann. Stat. 15 (4): 1421–38. 
  9. ^ Kay, S. M. (1993). Fundamentals of Statistical Signal Processing: Estimation Theory. Prentice Hall. p. 47. ISBN 0-13-042268-1 

参考文献 編集

  • Bos, Adriaan van den (2007). Parameter Estimation for Scientists and Engineers. Hoboken: John Wiley & Sons. pp. 45–98. ISBN 0-470-14781-4 
  • Kay, Steven M. (1993). Fundamentals of Statistical Signal Processing, Volume I: Estimation Theory. Prentice Hall. ISBN 0-13-345711-7 . Chapter 3.
  • Shao, Jun (1998). Mathematical Statistics. New York: Springer. ISBN 0-387-98674-X . Section 3.1.3.