数理統計 第7章 §7.5 ★★★ Neyman–Pearson の基本定理・単調尤度比と一様最強力検定・不偏性。2025年 統計数理 問4 の題材に直結
§7.3 の尤度比検定は統計量を作る手続きで、それが最良かどうかは別の問題である。この節では「同じ水準の検定の中で検出力が最大」という基準を定め、その最良の検定を構成する。
最強力検定と Neyman–Pearson の基本定理
定義 7.10(最強力検定・一様最強力検定)
単純仮説どうしの検定 H 0 : θ = θ 0 対 H 1 : θ = θ 1 H_0:\theta=\theta_0\ \text{対}\ H_1:\theta=\theta_1 H 0 : θ = θ 0 対 H 1 : θ = θ 1 を考える。水準 α \alpha α の検定 φ ∗ \varphi^* φ ∗ が、同じ水準の任意の検定に対して
E θ 1 [ φ ∗ ] ≥ E θ 1 [ φ ] E_{\theta_1}[\varphi^*]\ge E_{\theta_1}[\varphi] E θ 1 [ φ ∗ ] ≥ E θ 1 [ φ ]
を満たすとき、φ ∗ \varphi^* φ ∗ を最強力検定 (MP 検定)と呼ぶ。
対立仮説が複合仮説 Θ 1 \Theta_1 Θ 1 のときは、その全体で同時に次が成り立つ水準 α \alpha α の検定を一様最強力検定 (UMP 検定)と呼ぶ。
E θ [ φ ∗ ] ≥ E θ [ φ ] ( 水準 α の任意の検定 φ ) . E_\theta[\varphi^*]\ge E_\theta[\varphi]\qquad(\text{水準}\ \alpha\ \text{の任意の検定}\ \varphi). E θ [ φ ∗ ] ≥ E θ [ φ ] ( 水準 α の任意の検定 φ ) .
帰無仮説と対立仮説のもとでのデータの密度を、それぞれ f 0 f_0 f 0 、f 1 f_1 f 1 と書く(離散なら確率関数)。
定理 7.11(Neyman–Pearson の基本定理 ★★★)
0 < α < 1 0<\alpha<1 0 < α < 1 とする。
(存在)定数 k ≥ 0 k\ge0 k ≥ 0 と γ ∈ [ 0 , 1 ] \gamma\in[0,1] γ ∈ [ 0 , 1 ] をうまく選ぶと、
φ ∗ ( x ) = { 1 ( f 1 ( x ) > k f 0 ( x ) ) γ ( f 1 ( x ) = k f 0 ( x ) ) 0 ( f 1 ( x ) < k f 0 ( x ) ) (7.14) \varphi^*(\bm x)=\begin{cases}1&(f_1(\bm x)>kf_0(\bm x))\\ \gamma&(f_1(\bm x)=kf_0(\bm x))\\ 0&(f_1(\bm x)<kf_0(\bm x))\end{cases}\tag{7.14} φ ∗ ( x ) = ⎩ ⎨ ⎧ 1 γ 0 ( f 1 ( x ) > k f 0 ( x )) ( f 1 ( x ) = k f 0 ( x )) ( f 1 ( x ) < k f 0 ( x )) ( 7.14 )
のサイズ E 0 [ φ ∗ ] E_0[\varphi^*] E 0 [ φ ∗ ] をちょうど α \alpha α にできる。2. (十分性)ある k ≥ 0 k\ge0 k ≥ 0 について次の2条件を満たし、サイズがちょうど α \alpha α である検定を考える。これは同じ水準の最強力検定である。
f 1 > k f 0 の点で φ ∗ = 1 , f 1 < k f 0 の点で φ ∗ = 0. f_1>kf_0\ \text{の点で}\ \varphi^*=1,\qquad f_1<kf_0\ \text{の点で}\ \varphi^*=0. f 1 > k f 0 の点で φ ∗ = 1 , f 1 < k f 0 の点で φ ∗ = 0.
3. (必要性)k > 0 k>0 k > 0 とする。水準 α \alpha α の最強力検定は、集合 { f 1 ≠ k f 0 } \{f_1\ne kf_0\} { f 1 = k f 0 } の上でほとんど至るところ (1) の検定と一致する。
証明 (1) 尤度比 Λ 1 ( X ) = f 1 ( X ) / f 0 ( X ) \Lambda_1(\bm X)=f_1(\bm X)/f_0(\bm X) Λ 1 ( X ) = f 1 ( X ) / f 0 ( X ) の、帰無仮説のもとでの分布を考える。f 0 = 0 f_0=0 f 0 = 0 の点は帰無仮説のもとで確率 0 なので無視してよい。次のようにおく。
G ( t ) = P 0 ( Λ 1 > t ) . G(t)=P_0(\Lambda_1>t). G ( t ) = P 0 ( Λ 1 > t ) .
関数 G G G は右連続な減少関数で、t → ∞ t\to\infty t → ∞ で 0 に近づく。また G ( 0 − ) = 1 G(0^-)=1 G ( 0 − ) = 1 である。k = inf { t ≥ 0 : G ( t ) ≤ α } k=\inf\{t\ge0:G(t)\le\alpha\} k = inf { t ≥ 0 : G ( t ) ≤ α } とすると
G ( k ) ≤ α ≤ G ( k − ) = P 0 ( Λ 1 ≥ k ) . G(k)\le\alpha\le G(k^-)=P_0(\Lambda_1\ge k). G ( k ) ≤ α ≤ G ( k − ) = P 0 ( Λ 1 ≥ k ) .
γ \gamma γ は、境目の確率 P 0 ( Λ 1 = k ) = G ( k − ) − G ( k ) P_0(\Lambda_1=k)=G(k^-)-G(k) P 0 ( Λ 1 = k ) = G ( k − ) − G ( k ) が正かどうかで決める。
正なら γ = ( α − G ( k ) ) / P 0 ( Λ 1 = k ) \gamma=(\alpha-G(k))/P_0(\Lambda_1=k) γ = ( α − G ( k )) / P 0 ( Λ 1 = k ) とおく。これは [ 0 , 1 ] [0,1] [ 0 , 1 ] に入る。
0 なら γ = 0 \gamma=0 γ = 0 とおく。
どちらの場合も E 0 [ φ ∗ ] = G ( k ) + γ P 0 ( Λ 1 = k ) = α E_0[\varphi^*]=G(k)+\gamma P_0(\Lambda_1=k)=\alpha E 0 [ φ ∗ ] = G ( k ) + γ P 0 ( Λ 1 = k ) = α となる。
(2) φ \varphi φ を水準 α \alpha α の任意の検定とする。点ごとに次が成り立つ。
f 1 > k f 0 f_1>kf_0 f 1 > k f 0 の点では φ ∗ = 1 ≥ φ \varphi^*=1\ge\varphi φ ∗ = 1 ≥ φ 。
f 1 < k f 0 f_1<kf_0 f 1 < k f 0 の点では φ ∗ = 0 ≤ φ \varphi^*=0\le\varphi φ ∗ = 0 ≤ φ 。
よってどの点でも ( φ ∗ − φ ) ( f 1 − k f 0 ) ≥ 0 (\varphi^*-\varphi)(f_1-kf_0)\ge0 ( φ ∗ − φ ) ( f 1 − k f 0 ) ≥ 0 である。積分して
0 ≤ ∫ ( φ ∗ − φ ) ( f 1 − k f 0 ) d x = ( E 1 [ φ ∗ ] − E 1 [ φ ] ) − k ( E 0 [ φ ∗ ] − E 0 [ φ ] ) . (7.15) 0\le\int(\varphi^*-\varphi)(f_1-kf_0)\,\mathrm{d}\bm x=\bigl(E_1[\varphi^*]-E_1[\varphi]\bigr)-k\bigl(E_0[\varphi^*]-E_0[\varphi]\bigr).\tag{7.15} 0 ≤ ∫ ( φ ∗ − φ ) ( f 1 − k f 0 ) d x = ( E 1 [ φ ∗ ] − E 1 [ φ ] ) − k ( E 0 [ φ ∗ ] − E 0 [ φ ] ) . ( 7.15 )
ここで E 0 [ φ ∗ ] = α ≥ E 0 [ φ ] E_0[\varphi^*]=\alpha\ge E_0[\varphi] E 0 [ φ ∗ ] = α ≥ E 0 [ φ ] かつ k ≥ 0 k\ge0 k ≥ 0 なので、右辺の第2項で引く量は 0 以上である。よって E 1 [ φ ∗ ] ≥ E 1 [ φ ] E_1[\varphi^*]\ge E_1[\varphi] E 1 [ φ ∗ ] ≥ E 1 [ φ ] 。
(3) φ \varphi φ も最強力なら E 1 [ φ ] = E 1 [ φ ∗ ] E_1[\varphi]=E_1[\varphi^*] E 1 [ φ ] = E 1 [ φ ∗ ] である。このとき (7.15) の右辺は − k ( α − E 0 [ φ ] ) ≤ 0 -k(\alpha-E_0[\varphi])\le0 − k ( α − E 0 [ φ ]) ≤ 0 となる。左辺の 0 ≤ 0\le 0 ≤ と合わせて、積分は 0 である。被積分関数は非負なので、ほとんど至るところ
( φ ∗ − φ ) ( f 1 − k f 0 ) = 0 (\varphi^*-\varphi)(f_1-kf_0)=0 ( φ ∗ − φ ) ( f 1 − k f 0 ) = 0
が成り立つ。すなわち f 1 ≠ k f 0 f_1\ne kf_0 f 1 = k f 0 の点では φ = φ ∗ \varphi=\varphi^* φ = φ ∗ 。
定理の使い方は決まっている。
尤度比 f 1 / f 0 f_1/f_0 f 1 / f 0 を計算し、それが十分統計量 T T T の単調増加関数になっていることを確かめる。
すると「f 1 / f 0 > k f_1/f_0>k f 1 / f 0 > k 」は「T > c T>c T > c 」と同値になる。最強力検定は、T > c T>c T > c なら棄却し、境目では確率 γ \gamma γ で棄却する形になる。
帰無仮説のもとでの統計量の分布から、サイズが指定の水準になるように c c c と γ \gamma γ を決める。
T T T が連続分布なら P 0 ( T = c ) = 0 P_0(T=c)=0 P 0 ( T = c ) = 0 なので、確率化は不要である。離散分布なら一般に確率化が必要になる。
正規平均
正規分布 N ( μ , 1 ) \mathrm{N}(\mu,1) N ( μ , 1 ) の標本で、H 0 : μ = 0 対 H 1 : μ = μ 1 ( μ 1 > 0 ) H_0:\mu=0\ \text{対}\ H_1:\mu=\mu_1\ (\mu_1>0) H 0 : μ = 0 対 H 1 : μ = μ 1 ( μ 1 > 0 ) とする。尤度比
f 1 f 0 = exp ( μ 1 ∑ x i − n μ 1 2 2 ) \frac{f_1}{f_0}=\exp\!\left(\mu_1\sum x_i-\frac{n\mu_1^2}{2}\right) f 0 f 1 = exp ( μ 1 ∑ x i − 2 n μ 1 2 )
は標本平均の増加関数なので、最強力検定は X ˉ > z α / n \bar X>z_\alpha/\sqrt n X ˉ > z α / n で棄却する。この棄却域は対立仮説の値 μ 1 \mu_1 μ 1 によらない。
単調尤度比と一様最強力検定
上の例では、最強力検定の棄却域が μ 1 \mu_1 μ 1 によらなかった。これが起きる条件を一般化したものが単調尤度比である。
定義 7.12(単調尤度比)
1次元のパラメータをもつ分布族 { f ( x ; θ ) } \{f(\bm x;\theta)\} { f ( x ; θ )} と、統計量 T ( x ) T(\bm x) T ( x ) を考える。任意の θ 1 < θ 2 \theta_1<\theta_2 θ 1 < θ 2 に対して、比
f ( x ; θ 2 ) f ( x ; θ 1 ) \frac{f(\bm x;\theta_2)}{f(\bm x;\theta_1)} f ( x ; θ 1 ) f ( x ; θ 2 )
が T ( x ) T(\bm x) T ( x ) の非減少関数になるとき、分布族は T T T について単調尤度比 (MLR)をもつという。分母が 0 で分子が正の点では、比を + ∞ +\infty + ∞ とみなす。
1母数の指数型分布族を例にとる。
f ( x ; θ ) = h ( x ) exp { η ( θ ) T ( x ) − A ( θ ) } f(\bm x;\theta)=h(\bm x)\exp\{\eta(\theta)T(\bm x)-A(\theta)\} f ( x ; θ ) = h ( x ) exp { η ( θ ) T ( x ) − A ( θ )}
で、η \eta η が増加関数とする。このとき比は
exp { ( η ( θ 2 ) − η ( θ 1 ) ) T ( x ) − ( A ( θ 2 ) − A ( θ 1 ) ) } \exp\{(\eta(\theta_2)-\eta(\theta_1))T(\bm x)-(A(\theta_2)-A(\theta_1))\} exp {( η ( θ 2 ) − η ( θ 1 )) T ( x ) − ( A ( θ 2 ) − A ( θ 1 ))}
で、T T T の増加関数なので MLR をもつ。正規・二項・ポアソン・指数分布はすべてこの形である。
定理 7.13(Karlin–Rubin の定理 ★★★)
分布族が T T T について MLR をもつとする。仮説 H 0 : θ ≤ θ 0 対 H 1 : θ > θ 0 H_0:\theta\le\theta_0\ \text{対}\ H_1:\theta>\theta_0 H 0 : θ ≤ θ 0 対 H 1 : θ > θ 0 に対し、次の検定を考える。
φ ( x ) = { 1 ( T ( x ) > c ) γ ( T ( x ) = c ) 0 ( T ( x ) < c ) , E θ 0 [ φ ] = α \varphi(\bm x)=\begin{cases}1&(T(\bm x)>c)\\ \gamma&(T(\bm x)=c)\\ 0&(T(\bm x)<c)\end{cases},\qquad E_{\theta_0}[\varphi]=\alpha φ ( x ) = ⎩ ⎨ ⎧ 1 γ 0 ( T ( x ) > c ) ( T ( x ) = c ) ( T ( x ) < c ) , E θ 0 [ φ ] = α
これは水準 α \alpha α の一様最強力検定である。さらに検出力関数 β ( θ ) = E θ [ φ ] \beta(\theta)=E_\theta[\varphi] β ( θ ) = E θ [ φ ] は非減少である。
証明 準備 。任意の θ ′ < θ ′ ′ \theta'<\theta'' θ ′ < θ ′′ をとる。
検定 φ \varphi φ が、単純仮説 θ ′ 対 θ ′ ′ \theta'\ \text{対}\ \theta'' θ ′ 対 θ ′′ の Neyman–Pearson 型の検定であることを示す。比 f ( x ; θ ′ ′ ) / f ( x ; θ ′ ) f(\bm x;\theta'')/f(\bm x;\theta') f ( x ; θ ′′ ) / f ( x ; θ ′ ) は統計量の値だけで決まるので、その関数を r r r と書く。これは非減少関数なので、k = r ( c ) k=r(c) k = r ( c ) とおくと次が成り立つ。
f ( ⋅ ; θ ′ ′ ) > k f ( ⋅ ; θ ′ ) f(\cdot;\theta'')>kf(\cdot;\theta') f ( ⋅ ; θ ′′ ) > k f ( ⋅ ; θ ′ ) なら r ( T ) > r ( c ) r(T)>r(c) r ( T ) > r ( c ) である。よって T > c T>c T > c で、φ = 1 \varphi=1 φ = 1 。
f ( ⋅ ; θ ′ ′ ) < k f ( ⋅ ; θ ′ ) f(\cdot;\theta'')<kf(\cdot;\theta') f ( ⋅ ; θ ′′ ) < k f ( ⋅ ; θ ′ ) なら T < c T<c T < c で、検定関数は 0。
したがって定理 7.11 (2) より、φ \varphi φ は同じサイズの検定の中で、点 θ ′ ′ \theta'' θ ′′ での検出力が最大となる。
残りを3段で示す。
検出力の単調性 。定数の検定 φ 0 ≡ β ( θ ′ ) \varphi_0\equiv\beta(\theta') φ 0 ≡ β ( θ ′ ) もサイズ β ( θ ′ ) \beta(\theta') β ( θ ′ ) をもつ。準備の結果から
β ( θ ′ ′ ) = E θ ′ ′ [ φ ] ≥ E θ ′ ′ [ φ 0 ] = β ( θ ′ ) . \beta(\theta'')=E_{\theta''}[\varphi]\ge E_{\theta''}[\varphi_0]=\beta(\theta'). β ( θ ′′ ) = E θ ′′ [ φ ] ≥ E θ ′′ [ φ 0 ] = β ( θ ′ ) .
2. 水準 。検出力関数が非減少なので、次が成り立つ。よって φ \varphi φ は水準 α \alpha α の検定である。
sup θ ≤ θ 0 β ( θ ) = β ( θ 0 ) = α . \sup_{\theta\le\theta_0}\beta(\theta)=\beta(\theta_0)=\alpha. θ ≤ θ 0 sup β ( θ ) = β ( θ 0 ) = α .
3. 一様最強力性 。ψ \psi ψ を水準 α \alpha α の任意の検定とする。特に E θ 0 [ ψ ] ≤ α E_{\theta_0}[\psi]\le\alpha E θ 0 [ ψ ] ≤ α である。任意の θ 1 > θ 0 \theta_1>\theta_0 θ 1 > θ 0 をとる。準備の結果を θ ′ = θ 0 \theta'=\theta_0 θ ′ = θ 0 、θ ′ ′ = θ 1 \theta''=\theta_1 θ ′′ = θ 1 に使う。すると E θ 1 [ φ ] ≥ E θ 1 [ ψ ] E_{\theta_1}[\varphi]\ge E_{\theta_1}[\psi] E θ 1 [ φ ] ≥ E θ 1 [ ψ ] を得る。θ 1 \theta_1 θ 1 は任意なので φ \varphi φ は一様最強力である。
対立仮説が H 1 : θ < θ 0 H_1:\theta<\theta_0 H 1 : θ < θ 0 の向きでは、不等号を逆にした「T < c T<c T < c で棄却」が一様最強力になる。
両側検定と不偏性
両側の問題 H 0 : μ = μ 0 対 H 1 : μ ≠ μ 0 H_0:\mu=\mu_0\ \text{対}\ H_1:\mu\ne\mu_0 H 0 : μ = μ 0 対 H 1 : μ = μ 0 では、一様最強力検定は存在しない。理由は次のとおりである。
μ 1 > μ 0 \mu_1>\mu_0 μ 1 > μ 0 に対する最強力検定は、右側の棄却域をもつ。
μ 1 < μ 0 \mu_1<\mu_0 μ 1 < μ 0 に対する最強力検定は、左側の棄却域をもつ。
定理 7.11 (3) によりどちらも本質的に一意だから、両方を同時に満たす検定がない(確認テスト 7.5-5)。
図 7.2 の右片側検定は、μ > 0 \mu>0 μ > 0 で最も検出力が高い。しかし μ < 0 \mu<0 μ < 0 では検出力が α \alpha α より小さくなる。これは「対立仮説が正しいときに、帰無仮説が正しいときよりも棄却しにくい」という不自然な性質である。
定義 7.14(不偏検定・一致性)
水準 α \alpha α の検定で、次を満たすものを不偏検定 と呼ぶ。
β ( θ ) ≥ α ( すべての θ ∈ Θ 1 ) . \beta(\theta)\ge\alpha\qquad(\text{すべての}\ \theta\in\Theta_1). β ( θ ) ≥ α ( すべての θ ∈ Θ 1 ) .
標本の大きさ n n n の検定の列 φ n \varphi_n φ n が次を満たすとき、一致性 をもつという。
β n ( θ ) → 1 ( n → ∞ ) ( すべての θ ∈ Θ 1 ) . \beta_n(\theta)\to1\quad(n\to\infty)\qquad(\text{すべての}\ \theta\in\Theta_1). β n ( θ ) → 1 ( n → ∞ ) ( すべての θ ∈ Θ 1 ) .
不偏検定に範囲を絞ると最良のものが存在する場合があり、それを一様最強力不偏検定 (UMPU 検定)と呼ぶ。正規分布の両側 z z z 検定、両側 t t t 検定はUMPU 検定であることが知られている(証明は指数型分布族の条件付き検定の議論を使うので省く)。両側 z z z 検定が不偏であることは確認テスト 7.2-4 で示した。
図 7.2 n = 16 n=16 n = 16 の N ( μ , 1 ) \mathrm{N}(\mu,1) N ( μ , 1 ) 標本で H 0 : μ = 0 H_0:\mu=0 H 0 : μ = 0 を水準 0.05 で検定するときの検出力関数。実線は右片側検定 X ˉ > z 0.05 / 4 \bar X>z_{0.05}/4 X ˉ > z 0.05 /4 、破線は両側検定 ∣ X ˉ ∣ > z 0.025 / 4 \lvert \bar X\rvert >z_{0.025}/4 ∣ X ˉ ∣ > z 0.025 /4 。右片側検定は μ > 0 \mu>0 μ > 0 で両側検定より検出力が高いが、μ < 0 \mu<0 μ < 0 では 0.05 を下回る。両側検定は μ = 0 \mu=0 μ = 0 で最小値 0.05 をとる不偏検定である。
例題 7.8 ポアソン分布の最強力検定と確率化 ★★★
オリジナル X 1 , … , X 4 X_1,\dots,X_4 X 1 , … , X 4 は独立に P o ( λ ) \mathrm{Po}(\lambda) Po ( λ ) に従う。H 0 : λ = 0.5 対 H 1 : λ = 1 H_0:\lambda=0.5\ \text{対}\ H_1:\lambda=1 H 0 : λ = 0.5 対 H 1 : λ = 1 の有意水準 0.05 の最強力検定を求め、その検出力を求めよ。また、この検定が H 0 : λ ≤ 0.5 対 H 1 : λ > 0.5 H_0:\lambda\le0.5\ \text{対}\ H_1:\lambda>0.5 H 0 : λ ≤ 0.5 対 H 1 : λ > 0.5 の一様最強力検定でもあることを説明せよ。
解答 T = ∑ i = 1 4 X i T=\sum_{i=1}^4X_i T = ∑ i = 1 4 X i とおく。尤度比は
f 1 ( x ) f 0 ( x ) = ∏ i = 1 4 e − 1 1 x i / x i ! e − 0.5 0.5 x i / x i ! = e − 2 ⋅ 2 T \frac{f_1(\bm x)}{f_0(\bm x)}=\prod_{i=1}^4\frac{e^{-1}1^{x_i}/x_i!}{e^{-0.5}0.5^{x_i}/x_i!}=e^{-2}\cdot2^{T} f 0 ( x ) f 1 ( x ) = i = 1 ∏ 4 e − 0.5 0. 5 x i / x i ! e − 1 1 x i / x i ! = e − 2 ⋅ 2 T
で、T T T の増加関数である。定理 7.11 より、最強力検定は T > c T>c T > c で棄却し、境目では確率 γ \gamma γ で棄却する形になる。ポアソン分布の再生性から、帰無仮説のもとで T ∼ P o ( 2 ) T\sim\mathrm{Po}(2) T ∼ Po ( 2 ) 。
P 0 ( T ≥ 5 ) = 1 − e − 2 ( 1 + 2 + 2 + 4 3 + 2 3 ) = 1 − 7 e − 2 = 0.05265 , P_0(T\ge5)=1-e^{-2}\left(1+2+2+\tfrac43+\tfrac23\right)=1-7e^{-2}=0.05265, P 0 ( T ≥ 5 ) = 1 − e − 2 ( 1 + 2 + 2 + 3 4 + 3 2 ) = 1 − 7 e − 2 = 0.05265 ,
P 0 ( T ≥ 6 ) = P 0 ( T ≥ 5 ) − P 0 ( T = 5 ) = 0.05265 − e − 2 2 5 5 ! = 0.05265 − 0.03609 = 0.01656. P_0(T\ge6)=P_0(T\ge5)-P_0(T=5)=0.05265-e^{-2}\tfrac{2^5}{5!}=0.05265-0.03609=0.01656 . P 0 ( T ≥ 6 ) = P 0 ( T ≥ 5 ) − P 0 ( T = 5 ) = 0.05265 − e − 2 5 ! 2 5 = 0.05265 − 0.03609 = 0.01656.
P 0 ( T ≥ 6 ) < 0.05 < P 0 ( T ≥ 5 ) P_0(T\ge6)<0.05<P_0(T\ge5) P 0 ( T ≥ 6 ) < 0.05 < P 0 ( T ≥ 5 ) なので、c = 5 c=5 c = 5 とする。次の式を解いて γ \gamma γ を決める。
0.01656 + γ × 0.03609 = 0.05 ⟹ γ = 0.9265. 0.01656+\gamma\times0.03609=0.05\quad\Longrightarrow\quad\gamma=0.9265 . 0.01656 + γ × 0.03609 = 0.05 ⟹ γ = 0.9265.
最強力検定は次のとおりである。
φ ( T ) = { 1 ( T ≥ 6 ) 0.9265 ( T = 5 ) 0 ( T ≤ 4 ) \varphi(T)=\begin{cases}1&(T\ge6)\\ 0.9265&(T=5)\\ 0&(T\le4)\end{cases} φ ( T ) = ⎩ ⎨ ⎧ 1 0.9265 0 ( T ≥ 6 ) ( T = 5 ) ( T ≤ 4 )
検出力は、H 1 H_1 H 1 のもとで T ∼ P o ( 4 ) T\sim\mathrm{Po}(4) T ∼ Po ( 4 ) なので
β = P 1 ( T ≥ 6 ) + γ P 1 ( T = 5 ) = 0.21487 + 0.9265 × 0.15629 = 0.3597. \beta=P_1(T\ge6)+\gamma P_1(T=5)=0.21487+0.9265\times0.15629=0.3597 . β = P 1 ( T ≥ 6 ) + γ P 1 ( T = 5 ) = 0.21487 + 0.9265 × 0.15629 = 0.3597.
後半:ポアソン分布の同時確率関数は
∏ e − λ λ x i x i ! = ( ∏ 1 x i ! ) exp { T log λ − 4 λ } \prod\frac{e^{-\lambda}\lambda^{x_i}}{x_i!}=\Bigl(\prod\frac{1}{x_i!}\Bigr)\exp\{T\log\lambda-4\lambda\} ∏ x i ! e − λ λ x i = ( ∏ x i ! 1 ) exp { T log λ − 4 λ }
である。η ( λ ) = log λ \eta(\lambda)=\log\lambda η ( λ ) = log λ は増加関数だから、T T T について MLR をもつ。上の検定は T T T の右側で棄却し、λ = 0.5 \lambda=0.5 λ = 0.5 でサイズ 0.05 である。よって定理 7.13 より一様最強力である。
例題 7.9 一様分布の上端の片側検定 ★★
オリジナル X 1 , … , X n X_1,\dots,X_n X 1 , … , X n は独立に一様分布 U ( 0 , θ ) \mathrm{U}(0,\theta) U ( 0 , θ ) に従う。H 0 : θ ≤ 1 対 H 1 : θ > 1 H_0:\theta\le1\ \text{対}\ H_1:\theta>1 H 0 : θ ≤ 1 対 H 1 : θ > 1 を考える。
分布族が X ( n ) = max i X i X_{(n)}=\max_iX_i X ( n ) = max i X i について MLR をもつことを示せ。
棄却域が { X ( n ) > c } \{X_{(n)}>c\} { X ( n ) > c } の形の、水準 α \alpha α の一様最強力検定を考える。c c c と検出力関数を求めよ。
標本の大きさ 5、水準 0.05 のとき、θ = 1.2 \theta=1.2 θ = 1.2 での検出力を求めよ。また θ = 1.2 \theta=1.2 θ = 1.2 での検出力を 0.9 以上にする最小の n n n を求めよ。
解答
すべての x i ≥ 0 x_i\ge0 x i ≥ 0 のもとで、同時密度は f ( x ; θ ) = θ − n 1 { x ( n ) ≤ θ } f(\bm x;\theta)=\theta^{-n}\bm 1\{x_{(n)}\le\theta\} f ( x ; θ ) = θ − n 1 { x ( n ) ≤ θ } である。θ 1 < θ 2 \theta_1<\theta_2 θ 1 < θ 2 に対して
f ( x ; θ 2 ) f ( x ; θ 1 ) = { ( θ 1 / θ 2 ) n ( x ( n ) ≤ θ 1 ) + ∞ ( θ 1 < x ( n ) ≤ θ 2 ) \frac{f(\bm x;\theta_2)}{f(\bm x;\theta_1)}=\begin{cases}(\theta_1/\theta_2)^n&(x_{(n)}\le\theta_1)\\ +\infty&(\theta_1<x_{(n)}\le\theta_2)\end{cases} f ( x ; θ 1 ) f ( x ; θ 2 ) = { ( θ 1 / θ 2 ) n + ∞ ( x ( n ) ≤ θ 1 ) ( θ 1 < x ( n ) ≤ θ 2 )
で、x ( n ) x_{(n)} x ( n ) について非減少である(x ( n ) > θ 2 x_{(n)}>\theta_2 x ( n ) > θ 2 はどちらの密度も 0 で考えなくてよい)。2. 最大値の分布関数は P θ ( X ( n ) ≤ x ) = ( x / θ ) n P_\theta(X_{(n)}\le x)=(x/\theta)^n P θ ( X ( n ) ≤ x ) = ( x / θ ) n (0 ≤ x ≤ θ 0\le x\le\theta 0 ≤ x ≤ θ )である。境界 θ = 1 \theta=1 θ = 1 でのサイズは 1 − c n 1-c^n 1 − c n なので
1 − c n = α ⟹ c = ( 1 − α ) 1 / n . 1-c^n=\alpha\quad\Longrightarrow\quad c=(1-\alpha)^{1/n}. 1 − c n = α ⟹ c = ( 1 − α ) 1/ n .
定理 7.13 より一様最強力で、検出力関数は
β ( θ ) = { 0 ( θ ≤ c ) 1 − ( c / θ ) n = 1 − ( 1 − α ) θ − n ( θ > c ) . \beta(\theta)=\begin{cases}0&(\theta\le c)\\ 1-(c/\theta)^n=1-(1-\alpha)\theta^{-n}&(\theta>c)\end{cases}. β ( θ ) = { 0 1 − ( c / θ ) n = 1 − ( 1 − α ) θ − n ( θ ≤ c ) ( θ > c ) .
3. 検出力は β ( 1.2 ) = 1 − 0.95 / 1.2 5 = 1 − 0.95 / 2.48832 = 0.618 \beta(1.2)=1-0.95/1.2^5=1-0.95/2.48832=0.618 β ( 1.2 ) = 1 − 0.95/1. 2 5 = 1 − 0.95/2.48832 = 0.618 。必要な標本の大きさは次のように求める。
1 − 0.95 / 1.2 n ≥ 0.9 ⟺ 1.2 n ≥ 9.5 ⟺ n ≥ log 9.5 / log 1.2 = 12.35. \begin{aligned}1-0.95/1.2^n\ge0.9
&\iff 1.2^n\ge9.5\\
&\iff n\ge\log9.5/\log1.2=12.35 .\end{aligned} 1 − 0.95/1. 2 n ≥ 0.9 ⟺ 1. 2 n ≥ 9.5 ⟺ n ≥ log 9.5/ log 1.2 = 12.35.
よって最小の n n n は 13。
確認テスト
オリジナル
7.5-1 右片側 z z z 検定は X ˉ > μ 0 + z α σ / n \bar X>\mu_0+z_\alpha\sigma/\sqrt n X ˉ > μ 0 + z α σ / n で棄却する。これを両側の問題 H 0 : μ = μ 0 対 H 1 : μ ≠ μ 0 H_0:\mu=\mu_0\ \text{対}\ H_1:\mu\ne\mu_0 H 0 : μ = μ 0 対 H 1 : μ = μ 0 に使うと、不偏検定にならない。その理由を1行で答えよ。
7.5-2 1個の観測 X X X について、次の2つの仮説を考える(どちらも 0 < x < 1 0<x<1 0 < x < 1 上の密度)。
H 0 : f 0 ( x ) = 1 , H 1 : f 1 ( x ) = 2 x . H_0:\ f_0(x)=1,\qquad H_1:\ f_1(x)=2x. H 0 : f 0 ( x ) = 1 , H 1 : f 1 ( x ) = 2 x .
水準 α \alpha α の最強力検定を求め、検出力を α \alpha α の式で表せ。水準 0.05 での値も求めよ。3. 7.5-3 X 1 , … , X 5 X_1,\dots,X_5 X 1 , … , X 5 は独立に指数分布 E x ( 1 / θ ) \mathrm{Ex}(1/\theta) Ex ( 1/ θ ) に従う。平均は θ \theta θ である。仮説 H 0 : θ ≤ 2 対 H 1 : θ > 2 H_0:\theta\le2\ \text{対}\ H_1:\theta>2 H 0 : θ ≤ 2 対 H 1 : θ > 2 の水準 0.05 の一様最強力検定を求めよ。また θ = 4 \theta=4 θ = 4 での検出力を求めよ。4. 7.5-4 両側 z z z 検定(確認テスト 7.2-4 の検定)が一致性をもつことを示せ。5. 7.5-5 N ( μ , 1 ) \mathrm{N}(\mu,1) N ( μ , 1 ) の標本で、仮説 H 0 : μ = 0 対 H 1 : μ ≠ 0 H_0:\mu=0\ \text{対}\ H_1:\mu\ne0 H 0 : μ = 0 対 H 1 : μ = 0 を考える。水準 α \alpha α の一様最強力検定が存在しないことを、定理 7.11 (3) を用いて示せ。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com