ネイマン・ピアソンの補題とは何ですか? ネイマン-ピアソンの補題は、統計的仮説検定の分野における基本定理であり、単純な仮説に対する最も強力な検定を決定するための枠組みを提供します。1930 年代にイェジ・ネイマンとエゴン・ピアソンによって導入され、現代の統計的推論の基礎と.. ここでは、そんな仮説検定とネイマンピアソンの基本定理について、 仮説検定とは? 仮説検定で使われる用語として挙げられる「検定統計量・有意水準・P値」について、 棄却域や採択域から判定を誤る「第一種の錯誤と第二種の錯誤」について、
ネイマン・ピアソンの補題(ネイマン・ピアソンのほだい)とは、統計学的仮説検定に関する補題。
2つの仮説
H0: θ=θ0とH1: θ=θ1
の間で仮説検定を行う際に、H1を支持しH0を排除するような、次に示す尤度比による尤度比検定:
Λ
(
x
)
=
L
(
x
∣
θ
0
)
L
(
x
∣
θ
1
)
≤
k
{\displaystyle \Lambda (x)={\frac {L(x\mid \theta _{0})}{L(x\mid \theta _{1})}}\leq k}
(ただしここで
P
r
(
Λ
(
X
)
≤
k
|
H
0
)
=
α
{\displaystyle Pr(\Lambda (X)\leq k|H_{0})=\alpha }
とする)が、サイズ(危険率、第一種過誤)
α
{\displaystyle \alpha }
の仮説検定の中で最もパワー(検出力)
1
−
β
{\displaystyle 1-\beta }
が大きい、というものである(
β
{\displaystyle \beta }
は第二種過誤)。
θ
0
{\displaystyle \theta _{0}}
、
θ
1
{\displaystyle \theta _{1}}
が単純仮説であれば一様最強力検定となる。
「αを決めておき、その中で検出力が最も大きい検定法を選択する」という方針をネイマン・ピアソンの基準という。この補題はその方法を具体的に与えるものである。ただしこの尤度比検定法が直接用いられるよりも、近似が用いられることが多い。
証明
(尤度は確率密度関数
f
{\displaystyle f}
で与えられているものとする。また、この補題において構成される尤度比検定は、危険率が
α
{\displaystyle \alpha }
以下のあらゆる検定法の中で検出力が最大となっていることを示す。)
ネイマン・ピアソンの補題における帰無仮説の棄却域は
R
N
P
=
{
x
:
f
(
x
∣
θ
0
)
f
(
x
∣
θ
1
)
≤
k
}
{\displaystyle R_{NP}=\left\{x:{\frac {f(x\mid \theta _{0})}{f(x\mid \theta _{1})}}\leq k\right\}}
である。ここで
k
{\displaystyle k}
は
P
r
(
X
∈
R
N
P
|
H
0
)
=
α
{\displaystyle Pr(X\in R_{NP}|H_{0})=\alpha }
となるようとった定数。
危険率が
α
{\displaystyle \alpha }
以下の検定法を任意に取ってきて、その棄却域を
R
A
{\displaystyle R_{A}}
とする。
α
=
P
r
(
X
∈
R
N
P
|
H
0
)
≥
P
r
(
X
∈
R
A
|
H
0
)
.
{\displaystyle \alpha =Pr(X\in R_{NP}|H_{0})\geq Pr(X\in R_{A}|H_{0}).}
任意の母数に対し、確率の値は次のような2項の和に分けることができる。
P
r
(
X
∈
R
N
P
|
θ
)
=
P
r
(
X
∈
R
N
P
∩
R
A
|
θ
)
+
P
r
(
X
∈
R
N
P
∩
R
A
c
|
θ
)
P
r
(
X
∈
R
A
|
θ
)
=
P
r
(
X
∈
R
N
P
∩
R
A
|
θ
)
+
P
r
(
X
∈
R
N
P
c
∩
R
A
|
θ
)
{\displaystyle {\begin{aligned}Pr(X\in R_{NP}|\theta )&=Pr(X\in R_{NP}\cap R_{A}|\theta )+Pr(X\in R_{NP}\cap R_{A}^{c}|\theta )\\Pr(X\in R_{A}|\theta )&=Pr(X\in R_{NP}\cap R_{A}|\theta )+Pr(X\in R_{NP}^{c}\cap R_{A}|\theta )\end{aligned}}}
θ
=
θ
0
{\displaystyle \theta =\theta _{0}}
とすれば、
P
r
(
X
∈
R
N
P
∩
R
A
c
|
H
0
)
≥
P
r
(
X
∈
R
N
P
c
∩
R
A
|
H
0
)
{\displaystyle Pr(X\in R_{NP}\cap R_{A}^{c}|H_{0})\geq Pr(X\in R_{NP}^{c}\cap R_{A}|H_{0})}
である。
それぞれの検定法の検出力は
P
r
(
X
∈
R
N
P
|
H
1
)
{\displaystyle Pr(X\in R_{NP}|H_{1})}
,
P
r
(
X
∈
R
A
|
H
1
)
{\displaystyle Pr(X\in R_{A}|H_{1})}
であるので
P
r
(
X
∈
R
N
P
|
H
1
)
≥
P
r
(
X
∈
R
A
|
H
1
)
{\displaystyle Pr(X\in R_{NP}|H_{1})\geq Pr(X\in R_{A}|H_{1})}
を証明すれば良い。上記の等式よりこれは
P
r
(
X
∈
R
N
P
∩
R
A
c
|
H
1
)
≥
P
r
(
X
∈
R
N
P
c
∩
R
A
|
H
1
)
{\displaystyle Pr(X\in R_{NP}\cap R_{A}^{c}|H_{1})\geq Pr(X\in R_{NP}^{c}\cap R_{A}|H_{1})}
と同値である。そこで以下、この不等式を示す。
P
r
(
X
∈
R
N
P
∩
R
A
c
|
H
1
)
=
∫
R
N
P
∩
R
A
c
f
(
x
∣
θ
1
)
d
x
≥
1
k
∫
R
N
P
∩
R
A
c
f
(
x
∣
θ
0
)
d
x
by definition of
R
N
P
=
1
k
P
r
(
X
∈
R
N
P
∩
R
A
c
|
H
0
)
≥
1
k
P
r
(
X
∈
R
N
P
c
∩
R
A
|
H
0
)
=
1
k
∫
R
N
P
c
∩
R
A
f
(
x
∣
θ
0
)
d
x
>
∫
R
N
P
c
∩
R
A
f
(
x
∣
θ
1
)
d
x
by definition of
R
N
P
=
P
r
(
X
∈
R
N
P
c
∩
R
A
|
H
1
)
{\displaystyle {\begin{aligned}Pr(X\in R_{NP}\cap R_{A}^{c}|H_{1})&=\int _{R_{NP}\cap R_{A}^{c}}f(x\mid \theta _{1})\,dx\\[4pt]&\geq {\frac {1}{k}}\int _{R_{NP}\cap R_{A}^{c}}f(x\mid \theta _{0})\,dx&&{\text{by definition of }}R_{NP}\\[4pt]&={\frac {1}{k}}Pr(X\in R_{NP}\cap R_{A}^{c}|H_{0})\\[4pt]&\geq {\frac {1}{k}}Pr(X\in R_{NP}^{c}\cap R_{A}|H_{0})\\[4pt]&={\frac {1}{k}}\int _{R_{NP}^{c}\cap R_{A}}f(x\mid \theta _{0})\,dx\\[4pt]&>\int _{R_{NP}^{c}\cap R_{A}}f(x\mid \theta _{1})\,dx&&{\text{by definition of }}R_{NP}\\[4pt]&=Pr(X\in R_{NP}^{c}\cap R_{A}|H_{1})\end{aligned}}}
関連項目
- イェジ・ネイマン
- エゴン・ピアソン
参考文献
- 「自然科学の統計学」(東京大学出版会)ISBN 4-13-042067-4
わかりやすい仮説検定とネイマンピアソンの定理 Yosshi Labo.
【機械学習】「異常検知と変化検知 1章」ネイマンピアソンの補題の式の行間を埋めてみる MachineLearning Qiita
okuslidestat22 Speaker Deck
【統計検定対策】ネイマンピアソンの基本定理|KOYO 統計検定教材
okuslidestat22 Speaker Deck
自然科学の統計学6章:章末問題回答 Rinahuのブログ
自然科学の統計学6章:章末問題回答 Rinahuのブログ
okuslidestat22 Speaker Deck
ネイマン・ピアソンの基本定理の証明あるノマドの知の旅路~数学・統計学への道
okuslidestat22 Speaker Deck
okuslidestat22 Speaker Deck
自然科学の統計学6章:章末問題回答 Rinahuのブログ
確率・統計(統計検定2級対応)26:検出力関数とネイマン・ピアソンの定理 数式に埋もれる
okuslidestat22 Speaker Deck
一様最強力検定を用いた二項分布に関する導出の流れの確認とネイマン・ピアソンの補題の解釈 あつまれ統計の森
自然科学の統計学6章:章末問題回答 Rinahuのブログ
雑記: ネイマンピアソンの補題とカーリンルビンの定理 クッキーの日記
尤度比検定例題 ネイマンピアソンの補題 GWCA
一様最強力検定を用いた二項分布に関する導出の流れの確認とネイマン・ピアソンの補題の解釈 あつまれ統計の森
わかりやすい仮説検定とネイマンピアソンの定理 Yosshi Labo.
Press ←or →to navigate between chapters Press Sor /to search in the book. まとめ 本記事では、ネイマン・ピアソンの基本補題を導出し、最強力検定の構造を明らかにしました。 ネイマン・ピアソンの基本補題の本質は次のように要約できます。 有意水準 α α の制約の下で検出力を最大化する検定は、尤度比 f 1 (x) / f 0 (x) f 1(x)/f 0(x) が閾値 k k を超える領域を棄却域と.