十分な統計

統計学において十分性とは、標本データセットを用いて計算される統計量が、そのデータセットのパラメトリックモデルと関連して持つ特性です。十分統計量には、データセットが提供するモデルパラメータに関するすべての情報が含まれています。これは、モデルパラメータに関する情報を含まない補助統計量、およびパラメータに関する情報のみを含み補助情報を含まない完全統計量の概念と密接に関連しています。

関連する概念として線形十分性がある。これは十分性よりも弱い概念だが、十分な統計量がない場合にも適用できる。ただし、線形推定量に限定される。[1]コルモゴロフ構造関数は個々の有限データを扱い、それに関連する概念はアルゴリズム的十分統計量である。

この概念は1920年にロナルド・フィッシャー卿によって提唱されました。[2] スティーブン・スティグラーは1973年に、十分性の概念は分布形式の仮定に強く依存しているため記述統計では好まれなくなっていると指摘しました(下記のピットマン・クープマン・ダルモワの定理を参照)。しかし、理論的研究においては依然として非常に重要な概念です。[3]

背景

大まかに言えば、未知のパラメータ を条件とする、独立かつ同一に分布するデータの集合が与えられた場合、十分な統計量とは、その値がパラメータの推定値(例えば、最尤推定値)を計算するために必要なすべての情報を含む関数です。因数分解定理(下記参照)により、十分な統計量 の場合、確率密度は と表すことができます。この因数分解から、 の最大尤度推定値はを通じてのみ と相互作用することが容易にわかります。通常、十分な統計量はデータの単純な関数、例えばすべてのデータポイントの合計です。

より一般的には、「未知のパラメータ」は未知の量のベクトルを表す場合もあれば、モデルに関する未知または完全には指定されていないすべてのものを表す場合もあります。このような場合、十分な統計量は、共同十分統計量と呼ばれる関数の集合である場合があります。通常、関数の数はパラメータの数と同じです。例えば、平均分散が未知のガウス分布の場合、両方のパラメータの最大尤度推定値を推定できる共同十分統計量は、すべてのデータポイントの合計とすべてのデータポイントの2乗の合計(または、標本平均標本分散)という2つの関数で構成されます。

言い換えれば、データの結合確率分布は、パラメータの十分統計量の値が与えられた場合、パラメータから条件付きで独立である統計量と基礎パラメータはどちらもベクトルである可能性がある。

数学的な定義

統計量t  =  T ( X )が基礎パラメータθに対して十分であるとは、統計量t  =  T ( X )が与えられたデータXの条件付き確率分布がパラメータθに依存しない場合に限る[4]

あるいは、 θ上のすべての事前分布に対して、θT(X)間の相互情報量がθX間の相互情報量に等しい場合、統計量T ( X ) はθに対して十分であると言うこと もできます[5]言い換えれば、データ処理不等式は等式になります。

例えば、分散が既知の正規分布の平均μ (未知)を求めるには、標本平均値だけで十分です。標本平均値が分かれば、標本自体からはμに関する更なる情報を得ることはできません。一方、任意の分布の場合、中央値は平均値を求めるのに十分ではありません。標本の中央値が分かっていても、標本自体を知ることで母平均に関する更なる情報が得られます。例えば、中央値を下回る観測値はわずかに小さいだけであるのに対し、中央値を大きく上回る観測値は母平均に関する推論に影響を与える可能性があります。

フィッシャー・ネイマン因数分解定理

フィッシャーの因数分解定理あるいは因数分解基準は、十分統計量の便利な特徴付けを提供する。確率密度関数が ƒ θ ( x ) である場合、 T がθ に対して十分であることは、非負関数ghが次式を満たす場合に限ります。

すなわち、密度 ƒ は、一方の因子hがθに依存せず、もう一方の因子(θに依存する)がT ( x )を通じてのみxに依存するような積に因数分解できる。この一般的な証明はHalmosとSavage [6]によって与えられており、この定理はHalmos–Savage因数分解定理と呼ばれることもある。[7]以下の証明は特殊なケースを扱っているが、同様の方法で別の一般的な証明を行うこともできる。[8]多くの単純なケースでは、確率密度関数は、および、およびによって完全に指定される(例を参照)。

F ( t ) が1対1関数であり、 T が十分統計量である場合、F ( T ) が十分統計量であることは容易に理解できます。特に、十分統計量に非ゼロの定数を乗じることで、別の十分統計量を得ることができます。

尤度原理の解釈

この定理の含意は、尤度に基づく推論を用いる場合、十分統計量T ( X ) の値が同一となる2つのデータセットは、常にθに関する推論結果が同じになるという点である。因数分解基準によれば、尤度のθへの依存性はT ( X )のみに依存する。どちらの場合もこれが同じであるため、 θへの依存性も同じとなり、同一の推論結果が導かれる。

証拠

HoggとCraigによる。[9]は、 ι  <  θ  <  δに対して確率密度関数f ( xθ )を持つ分布からのランダムサンプルを表す。Y 1  =  u 1 ( X 1X 2 , ...,  X n )を、確率密度関数g 1 ( y 1θ )である統計量とする。証明したいのは、Y 1  =  u 1 ( X 1 , X 2 , ...,  X n )がθに対する十分統計量であるための必要十分条件であり、かつ、ある関数Hに対して、

まず、

i  = 1, ...,  nに対してy i  =  u i ( x 1x 2 , ...,  x n )という変換を行い、その逆関数x i  =  w i ( y 1y 2 ,  ...  y n )ヤコビ行列求める。したがって、

左側の要素は、Y 1 = u 1 ( X 1 , ..., X n ) 、... 、Y n = u n ( X 1 , ..., X n ) の結合確率密度関数g ( y 1 , y 2 , ..., y n ; θ ) です。右側の要素は の確率密度関数なので、は と の商です。つまり、が与えられたの条件付き確率密度関数です

しかし、 であり、したがって は に依存しないことが示されていますは変換に導入されておらず、したがってヤコビ行列 にも導入されていないため、は に依存しないことが示されは の十分な統計量となります

逆は次のように証明されます。

ここで はに依存しません。なぜならは のみに依存し、 はによって条件付けられる場合から独立しているからです。これは仮説により十分統計量となります。ここで、両方の要素を非零ヤコビアン の絶対値で割りの関数 を に置き換えます。これは次式を得ます。

ここで、は の項における値に置き換えられたヤコビアンである。左辺は必然的にの結合確率密度関数となる、そして はに依存しないので、

は に依存しない関数です

もう一つの証拠

より単純でわかりやすい証明は次のとおりです。ただし、これは離散的な場合にのみ適用されます。

の結合確率密度を で表すために、略記法を用いる。 は の決定論的関数であるためある限り が成り立ち、 でなければ0 となる。したがって、

最後の等式は十分統計量の定義により真となる。したがって、およびとなる

逆に、の場合には、

最初の等式は多重変数の pdf の定義により、2 番目は上記のコメントにより、3 番目は仮説により、4 番目は合計が を超えていないためです

与えられたの条件付き確率密度をとします。すると、これを明示的に表す式を導出できます。

最初の等式は条件付き確率密度の定義により、2番目の等式は上記の注釈により、3番目の等式は上で証明された等式により、そして4番目の等式は簡略化により成立する。この式は に依存しないため、十分な統計量となる。[10]

最低限の充足

十分統計量は、他の十分統計量の関数として表せる場合、最小十分統計量である。言い換えれば、 S ( X ) が最小十分統計量となるのは、[11]

  1. S ( X ) で十分であり、
  2. T ( X )が十分であれば、 S ( X )= f ( T ( X ))となる関数fが存在する。

直感的に言えば、最小限の十分統計量は、パラメータθに関するすべての可能な情報を最も効率的に捕捉します。

最小十分性の有用な特徴付けは、密度f θが存在するとき、S ( X )が最小十分性を満たす場合、

θに依存しない :S ( x )= S ( y )

これは、上で述べたフィッシャーの因数分解定理の結果として生じます。

最小十分統計量が存在しないケースは、Bahadur (1954) によって示されました。[12]しかし、軽度の条件下では、最小十分統計量は常に存在します。特にユークリッド空間では、( に関連付けられた)確率変数がすべて離散的であるか、すべて連続的である場合、これらの条件は常に成立します。

最小十分統計量が存在する場合(通常はそうである)、すべての完全十分統計量は必然的に最小十分統計量となる[13](ただし、この記述は、完全十分統計量が存在する一方で最小十分統計量が存在しないという病的なケースを排除するものではないことに注意)。最小十分統計量が存在しないケースを見つけるのは難しいが、完全十分統計量が存在しないケースを見つけるのはそれほど難しくない。

パラメータ空間が離散的である場合、尤度比の集合は最小十分統計量です

ベルヌーイ分布

X 1、...、  X nが期待値pを持つ独立したベルヌーイ分布の確率変数である場合、合計T ( X ) =  X 1  + ... +  X nはpに対する十分な統計量です(ここで、「成功」はX i  = 1 に対応し、「失敗」は X i  = 0 に対応します。したがって、Tは成功の総数です)。

これは結合確率分布を考慮するとわかります。

観測は独立しているので、これは次のように書ける。

そして、pと1 −  pの累乗を集めると、

これは因数分解の基準を満たしており、h ( x ) = 1 は単なる定数です。

重要な特徴に注意してください。未知のパラメータpは、統計T ( x )=Σxi 介してのみデータxと相互作用します。

具体的な応用として、これは公平なコインと偏ったコインを区別する手順を示します

均一分布

X 1、...、X nが独立しており、区間 [0、θ ]上で均一に分布している場合、 T ( X ) = max( X 1、...、X n ) は θ に対して十分です。つまり、標本の最大値は母集団の最大値に対する十分な統計量です。

これを理解するには、 X   ( X 1 ,..., X n )の結合確率密度関数を考えてみましょう。観測値は独立しているので、確率密度関数は個々の密度の積として表すことができます。

ここで、1 { ... }は指示関数です。したがって、密度はフィッシャー・ネイマン分解定理によって要求される形をとります。ここで、h ( x ) =  1 {min{ x i }≥0}となり、式の残りの部分はθT ( x ) = max{ x i } のみの関数となります。

実際、θの最小分散不偏推定量(MVUE)

これは、バイアスを補正するためにスケーリングされた標本最大値でありレーマン・シェッフェの定理によりMVUEとなる。スケーリングされていない標本最大値T ( X )は、θ最大尤度推定値である。

均一分布(2つのパラメータを持つ)

が独立しており、区間およびは未知のパラメータ)上で均一に分布している場合、 は の2 次元の十分な統計量です

これを理解するには、の結合確率密度関数を考えてみましょう。観測値は独立しているので、pdfは個々の密度の積として表すことができます。つまり、

サンプルの結合密度は、フィッシャー・ネイマン因数分解定理によって求められる形をとり、

はパラメータに依存せず関数を通じてのみ依存するため

フィッシャー・ネイマン因数分解定理によれば、は の十分な統計量である

ポアソン分布

X 1、...、  X nが独立しており、パラメータλを持つポアソン分布に従う場合、合計T ( X ) =  X 1  + ... +  X nはλに対する十分な統計量です 

これを理解するには、結合確率分布を考えてみましょう。

観測は独立しているので、これは次のように書ける。

これは次のように書ける。

これは因数分解基準が満たされていることを示しています。ここで、h ( x ) は階乗の積の逆数です。パラメータ λ は、その和T ( X ) を通じてのみデータと相互作用することに注意してください。

正規分布

が独立で、期待値(パラメータ)と既知の有限分散を持つ正規分布である場合

十分な統計量である

これを理解するには、の結合確率密度関数を考えてみましょう。観測値は独立しているので、pdfは個々の密度の積として表すことができます。つまり、

サンプルの結合密度は、フィッシャー・ネイマン因数分解定理によって求められる形をとり、

はパラメータに依存せず関数を通じてのみ依存するため

フィッシャー・ネイマン因数分解定理によれば、は の十分な統計量である

未知であり、 なので、上記の尤度は次のように書き直すことができる。

フィッシャー・ネイマン因数分解定理は依然として成り立ち、が の十分な結合統計量であることを意味します

指数分布

が独立しており、期待値θ (未知の実数値の正のパラメータ)を持つ指数分布である場合、 はθ の十分な統計量です。

これを理解するには、の結合確率密度関数を考えてみましょう。観測値は独立しているので、pdfは個々の密度の積として表すことができます。つまり、

サンプルの結合密度は、フィッシャー・ネイマン因数分解定理によって求められる形をとり、

はパラメータに依存せず関数を通じてのみ依存するため

フィッシャー・ネイマン因数分解定理によれば、は の十分な統計量である

ガンマ分布

が独立で として分布している場合およびはガンマ分布の未知のパラメータ)、 は の2 次元の十分な統計量です

これを理解するには、の結合確率密度関数を考えてみましょう。観測値は独立しているので、pdfは個々の密度の積として表すことができます。つまり、

サンプルの結合密度は、フィッシャー・ネイマン因数分解定理によって求められる形をとり、

はパラメータに依存せず関数を通じてのみ依存するため

フィッシャー・ネイマン因数分解定理によれば、

ラオ・ブラックウェルの定理

十分性はラオ・ブラックウェル定理に有用に応用される。この定理は、g ( X ) がθの何らかの推定値である場合、通常、十分な統計量T ( X ) を与えられたg ( X )の条件付き期待値は、 θのより良い推定値(分散が低いという意味で)であり、決して劣ることはない、と述べている。場合によっては、非常に粗い推定値g ( X )を非常に簡単に構築し、その条件付き期待値を評価することで、様々な意味で最適な推定値を得ることができる。

指数族

ピットマン・クープマン・ダルモワの定理によれば、推定されるパラメータによって定義域が変化しない確率分布族の中で、指数分布族においてのみ、標本数が増加しても次元が有界のままである十分な統計量が存在する。直感的には、これは実数直線上の非指数分布族では、データの情報を完全に捉えるためにノンパラメトリック統計量が必要であることを意味している。

もっと簡単に言えば、が独立した同一分布の実数ランダム変数であり、その分布が によってパラメータ化されたある確率分布族に属することがわかっており、特定の技術的な正則性条件を満たしているとすると、その族が指数族であるためには、サンプルサイズnが増加してもスカラー成分の数が増加しない十分な 値統計量が存在する必要がある[14]

この定理は、有限次元の実ベクトル値の十分な統計量の存在が、実数直線上の分布の族の可能な形式を厳しく制限することを示しています。

パラメータやランダム変数が実数値ではなくなると、状況はさらに複雑になります。[15]

その他の充足性の種類

ベイズ十分性

統計量が十分であるという条件をベイズ流の文脈で定式化すると、データセット全体を用いた場合と統計量のみを用いた場合の事後分布が考慮される。したがって、ほぼすべてのxに対して、

より一般的には、パラメトリックモデルを仮定せずに、統計量T十分に予測可能であると言えるのは

この「ベイズ的十分性」は上記の定式化の結果であることが判明しているが[16] 、無限次元の場合、それらは直接的には同等ではない。[17]ベイズ的文脈における十分性については、様々な理論的結果が得られている。[18]

線形充足性

「線形十分性」と呼ばれる概念は、ベイズ理論の文脈[19]、そしてより一般的には[20]で定式化することができる。まず、ベクトルYのXに基づく最良の線形予測子を と定義する。すると、線形統計量T ( x )が線形十分性を持つのは[21]

参照

注記

  1. ^ Dodge, Y. (2003) — 線形充足性の項目
  2. ^ フィッシャー, RA (1922). 「理論統計の数学的基礎について」. Philosophical Transactions of the Royal Society A . 222 ( 594– 604): 309– 368. Bibcode :1922RSPTA.222..309F. doi : 10.1098/rsta.1922.0009 . hdl : 2440/15172 . JFM  48.1280.02. JSTOR  91208.
  3. ^ スティグラー、スティーブン(1973年12月)「確率統計史研究 第32章 ラプラス、フィッシャー、そして十分性の概念の発見」Biometrika 60 ( 3): 439– 445. doi :10.1093/biomet/60.3.439. JSTOR  2334992. MR  0326872.
  4. ^ Casella, George; Berger, Roger L. (2002).統計的推論, 第2版. Duxbury Press.
  5. ^ Cover, Thomas M. (2006). 『情報理論の要素』 Joy A. Thomas (第2版). ホーボーケン, ニュージャージー: Wiley-Interscience. p. 36. ISBN 0-471-24195-4. OCLC  59879802.
  6. ^ Halmos, PR; Savage, LJ (1949). 「ラドン・ニコディムの定理の十分統計理論への応用」 .数理統計年報. 20 (2): 225– 241. doi : 10.1214/aoms/1177730032 . ISSN  0003-4851.
  7. ^ 「因数分解定理 - 数学百科事典」. encyclopediaofmath.org . 2022年9月7日閲覧
  8. ^ Taraldsen, G. (2022). 「十分性のための因数分解定理」.プレプリント. doi :10.13140/RG.2.2.15068.87687.
  9. ^ ホッグ, ロバート V.; クレイグ, アレン T. (1995).数理統計学入門. プレンティス・ホール. ISBN 978-0-02-355722-4
  10. ^ 「フィッシャー・ネイマン因数分解定理」.Connexions のウェブページ (cnx.org)
  11. ^ Dodge (2003) — 最小限の十分統計量のエントリ
  12. ^ レーマンとカセラ(1998年)、点推定理論、第2版、シュプリンガー、p 37
  13. ^ Lehmann and Casella (1998)、「点推定理論」第2版、Springer、42ページ
  14. ^ Tikochinsky, Y.; Tishby, NZ; Levine, RD (1984-11-01). 「最大エントロピー推論への代替アプローチ」 . Physical Review A. 30 ( 5): 2638– 2644. Bibcode :1984PhRvA..30.2638T. doi :10.1103/physreva.30.2638. ISSN  0556-2791.
  15. ^ アンダーセン、アーリング・ベルンハルト(1970年9月)離散標本空間における十分性と指数族」アメリカ統計学会誌65(331):1248-1255。doi:10.1080/01621459.1970.10481160。ISSN 0162-1459  。
  16. ^ Bernardo, JM ; Smith, AFM (1994). 「セクション5.1.4」.ベイズ理論. Wiley. ISBN 0-471-92416-4
  17. ^ Blackwell, D. ; Ramamoorthi, RV (1982). 「ベイズ統計だが古典的には十分ではない統計量」Annals of Statistics . 10 (3): 1025– 1026. doi : 10.1214/aos/1176345895 . MR  0663456. Zbl  0485.62004.
  18. ^ Nogales, AG; Oyola, JA; Perez, P. (2000). 「ベイズ的観点から見た条件付き独立性と十分性と不変性の関係について」 . Statistics & Probability Letters . 46 (1): 75– 84. doi :10.1016/S0167-7152(99)00089-9. MR  1731351. Zbl  0964.62003.
  19. ^ Goldstein, M.; O'Hagan, A. (1996). 「ベイズ線形十分性と専門家事後評価システム」.英国王立統計学会誌. シリーズB. 58 (2): 301– 316. doi :10.1111/j.2517-6161.1996.tb02083.x. JSTOR  2345978.
  20. ^ Godambe, VP (1966). 「有限母集団からのサンプリングへの新しいアプローチ II 分布に依存しない十分性」.英国王立統計学会誌. シリーズB. 28 (2): 320– 328. doi :10.1111/j.2517-6161.1966.tb00645.x. JSTOR  2984375.
  21. ^ Witting, T. (1987). 「信頼性理論における線形マルコフ性」. ASTIN Bulletin . 17 (1): 71– 84. doi : 10.2143/ast.17.1.2014984 . hdl : 20.500.11850/422507 .

参考文献

Retrieved from "https://en.wikipedia.org/w/index.php?title=Sufficient_statistic&oldid=1312788813"