統計に対する決定論的アプローチは、詳細な説明を提供します。 平方差は、(正当に採用される可能性がある場合はいつでも)考慮する必要のある統計的手順の大幅な簡素化につながる幅広い損失関数のプロキシであると言います。
残念ながら、これが何を意味するのかを説明し、それが本当である理由を示すには、多くの設定が必要です。表記はすぐに理解不能になる可能性があります。ここで私が目指しているのは、少しの手間をかけずに主要なアイデアをスケッチすることです。完全なアカウントについては、参照を参照してください。
データ標準の豊富なモデルは、分布Fが自然の状態である分布のある集合Ωの要素であることが知られている(実際のベクトル値の)ランダム変数Xの実現であると仮定します。統計的手順は、一連の決定D、決定空間で値を取るxの関数tです。xXFΩtxD
たとえば、予測または分類の問題では、は「トレーニングセット」と「データのテストセット」の結合で構成され、tはxをテストセットの予測値のセットにマッピングします。すべての可能な予測値のセットはDになります。xtxD
手順の完全な理論的な議論は無作為化された手順を収容しなければなりません。ランダム化された手順は、いくつかの確率分布(データ依存)に従って2つ以上の可能な決定から選択します。データが2つの選択肢を区別しないように見える場合、その後、明確な選択肢を決定するために「コインを裏返す」という直感的な考え方を一般化します。多くの人々は、ランダム化された手順を嫌い、そのような予測不可能な方法で決定を下すことに反対しています。x
決定理論の際立った特徴は、損失関数 使用です。W 任意の自然状態のためのと判断D ∈ D、損失F∈Ωd∈D
W(F,d)
自然の真の状態がFである場合に決定を行うことがどの程度「悪い」かを表す数値です。小さな損失は良好、大きな損失は不良です。たとえば、仮説検定の状況では、Dには2つの要素 "accept"および "reject"(帰無仮説)があります。損失関数は正しい決定を行うことを強調しています。決定が正しい場合はゼロに設定され、そうでない場合は定数wに設定されます。(これは「0 − 1損失関数:」と呼ばれます。すべての悪い決定は等しく悪く、すべての良い決定は等しく良いです。)特に、W (F 、 accept )= 0のときdFDw0−1W(F, accept)=0帰無仮説とである W (F 、 リジェクト)= 0のとき Fは、対立仮説です。FW(F, reject)=0F
手順を使用する場合、自然の真の状態がFであるときのデータxの損失は、W (F 、t (x ))と書くことができます。これにより、損失W (F 、t (X ))は、分布が(未知の)Fによって決定される確率変数になります。txFW(F,t(x))W(F,t(X))F
手続きの予想損失その呼び出されたリスクを、Rのトン。期待値は真の自然状態Fを使用するため、期待値演算子の添え字として明示的に表示されます。リスクをFの関数と見なし、表記法でそれを強調します。trtFF
rt(F)=EF(W(F,t(X))).
より良い手順はより低いリスクを持ちます。 したがって、リスク関数の比較は、優れた統計的手順を選択するための基礎となります。すべてのリスク関数を共通の(正の)定数で再スケーリングしても比較は変わらないため、のスケールに違いはありません。好きな正の値を掛けることができます。具体的には、乗算の際にWをすることにより1 / wは、我々は常にとることができるwは= 1のために0 - 1損失関数(その名前を正当化します)。WW1/ww=10−1
説明仮説検定の例、続行するには損失関数を、これらの定義は、いずれかの危険性を暗示Fの任意のリスク一方で「拒否」決定がされる可能性を帰無仮説でをされてFの代替ではあります決定が「受け入れられる」チャンス。(全てにわたる最大値F帰無仮説では)試験であるサイズ対立仮説上に定義されたリスク関数の一部が試験の補数であるが、電源(パワーT(F )= 1 - のR T(F )0−1FFFpowert(F)=1−rt(F))。これでは、古典的(頻度論的)仮説検定理論の全体が、特別な種類の損失のリスク関数を比較する特定の方法にどのようになるかがわかります。
ところで、これまでに紹介したものはすべて、ベイジアンパラダイムを含むすべての主流の統計と完全に互換性があります。さらに、ベイジアン分析では、に対する「事前」確率分布を導入し、これを使用してリスク関数の比較を簡素化します。潜在的に複雑な関数r tは、事前分布に関する期待値で置き換えることができます。したがって、すべての手順tは、単一の数値r tによって特徴付けられます。ベイズ手順(通常は一意)はr tを最小化します。損失関数は、r tの計算において依然として重要な役割を果たします。Ωrttrtrtrt
損失関数の使用をめぐるいくつかの(避けられない)論争があります。 どのように選ぶのですか?仮説検定では本質的にユニークですが、他のほとんどの統計設定では多くの選択肢が可能です。それらは意思決定者の価値を反映しています。たとえば、データが医療患者の生理学的測定値であり、決定が「治療する」または「治療しない」である場合、医師はいずれかのアクションの結果を考慮し、バランスをとる必要があります。結果がどのように評価されるかは、患者自身の希望、年齢、生活の質、および他の多くのものに依存する可能性があります。損失関数の選択は複雑で、深く個人的なものです。通常、統計学者に任せるべきではありません!W
私たちが知りたいことの1つは、損失が変わったときに最良の手順の選択がどのように変わるかということです。 多くの一般的で実際的な状況では、最適な手順を変更せずに一定量の変動を許容できることがわかります。これらの状況は、次の条件によって特徴付けられます。
決定空間は凸集合(多くの場合、数字の間隔)です。これは、2つの決定の間にある値も有効な決定であることを意味します。
可能な限り最良の決定が行われた場合の損失はゼロであり、そうでない場合は増加します(行われた決定と真の(ただし未知の)自然状態に対して行われる可能性のある最良の決定との矛盾を反映するため)。
損失は、決定の微分可能な関数です(少なくともローカルで最良の決定に近い)。それは仕方ジャンプしない-これは、それが連続している意味の損失はありませんが-それはまた、決定が近い最高のものにする場合、それは比較的小さな変化することを意味します。0−1
これらの条件が満たされると、リスク関数の比較に伴ういくつかの合併症がなくなります。の微分可能性と凸性により、ジェンセンの不等式を適用して、W
(1)ランダム化された手順を考慮する必要はありません[Lehmann、corollary 6.2]。
(2)1つの手続きがそのようなWに対して最高のリスクを持っていると考えられる場合、十分な統計量のみに依存し、少なくともすべてのそのようなWに対してリスク関数と同等の手続きt ∗に改善できます。、p。151]。tWt∗ W
例として、が平均μ(および単位分散)を持つ正規分布のセットであると仮定します。これにより、すべての実数のセットでΩが識別されます。したがって、表記法を使用して、「μ」を使用して、平均μでΩの分布を識別します。ましょXは、サイズのIIDサンプルでnはこれらのディストリビューションの一つから。μを推定することが目的だとします。これは、μ(任意の実数)のすべての可能な値で決定空間Dを識別します。まかせμ独断を指定し、損失関数でありますΩμΩμΩμXnμDμμ^
W(μ,μ^)≥0
の場合に限りμ = μ。上記の仮定は、(テイラーの定理を介して)W(μ,μ^)=0μ=μ^
W(μ,μ^)=w2(μ^−μ)2+o(μ^−μ)2
一定の正の数。( "表記少し-O 、O (Y )P手段任意の関数" Fの制限値F (Y )/ Y pがある0としてY → 0。)先に述べたように、我々が再スケールに自由であるWせるwは2 = 1。この家族のためにΩ、平均X書かれた、ˉ Xは、十分統計量です。前の結果(キーファーから引用)は、w2o(y)pff(y)/yp0y→0Ww2=1ΩXX¯のいくつかの任意の関数とすることができ、 n個の変数(X 1、... 、xはn個)そのような適して Wは、のみに依存推定に変換することができる ˉ X、少なくともそのようなすべてのための良好なようであり、 W。μn(x1,…,xn)Wx¯W
この例で成し遂げられたのは典型的なものです:もともと変数のおそらくランダム化された関数で構成されていた非常に複雑な可能性のある手続きのセットは、単一変数の非ランダム化関数または、十分な統計が多変量である場合、少なくとも少数の変数)。そしてこれは、意思決定者の損失関数が何であるかを正確に心配することなく行うことができます。ただし、凸関数で微分可能である場合に限ります。n
最も単純なそのような損失関数は何ですか? もちろん、残りの項を無視するもので、純粋に2次関数にします。この同じクラスの他の損失関数には、べき乗= | μ - μ | 2よりも大きい(質問で言及されている2.1 、e 、およびπなど)、exp (z )− 1 − zなど。z=|μ^−μ|22.1,e,πexp(z)−1−z
青(上の)曲線はをプロットし、赤(下)の曲線はz 2をプロットします。青い曲線も最小値が0で微分可能で凸であるため、2次損失(赤い曲線)が享受する統計的手法の優れた特性の多くは、2(exp(|z|)−1−|z|)z20(全体的に指数関数でも二次関数とは異なる動作をします)。
これらの結果は、(明らかに課せられた条件によって制限されるが)二次損失は、統計理論と実践に遍在である理由助けを説明:限られた範囲で、それがために解析的に便利なプロキシである任意の凸微分損失関数。
二次損失は、考慮すべき唯一の損失でも最良の損失でもありません。 確かに、リーマンはそれを書いています
凸損失関数は、推定問題のいくつかの単純化につながると見られています。しかし、そのような損失関数が現実的である可能性が高いかどうか疑問に思うかもしれません。場合はあなたがすべて失われたら、あなたはそれ以上失うことができない。ただ不正確性の尺度が、実際の(例えば、金融)の損失をない表し、1はそのようなすべての損失が制限されていることを主張することがあります。...W(F,d)
... [F]急成長する損失関数は、[仮定された分布の]テールの振る舞いについて行われた仮定に敏感になる傾向がある推定量を導きます。これらの仮定は通常、ほとんど情報に基づいていないため、あまり重要ではありません信頼性のある。
二乗誤差損失によって生成される推定量は、この点で不愉快なほど敏感であることが多いことがわかります。
[リーマン、セクション1.6。表記を少し変更しました。]
代替損失を考慮すると、豊富な可能性のセットが開かれます。分位点回帰、M推定量、ロバストな統計など、すべてをこの決定理論的な方法でフレーム化し、代替損失関数を使用して正当化できます。簡単な例については、パーセンタイル損失関数を参照してください。
参照資料
Jack Carl Kiefer、統計的推論の紹介。 Springer-Verlag 1987。
EL Lehmann、ポイント推定の理論。ワイリー1983。