Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

検出力

検定での意思決定の誤り

検定の意思決定の誤り

母数空間Θ\Thetaが互いに排反な2つの部分集合Θ0,Θ1\Theta_0, \Theta_1に分けられているとする。未知の母数θ\thetaについての帰無仮説H0H_0と対立仮説H1H_1を以下のように表す。

H0:θΘ0 vs H1:θΘ1H_0: \theta \in \Theta_0 \text { vs } H_1: \theta \in \Theta_1

帰無仮説を受容する意思決定をd=0d=0、帰無仮説を棄却する決定をd=1d=1と表すことにし、それらの集合(決定空間)はD={0,1}D=\{0, 1\}と表すことにする。

帰無仮説が正しいとき (θΘ0)(\theta \in \Theta_0) に帰無仮説を枼却する (d=1)(d=1) 誤りを 第1種の過誤(error of the first kind)偽陽性 とよぶ。また対立仮説が正しいとき (θΘ1)(\theta \in \Theta_1) に帰無仮説を受容する (d=0)(d=0) 誤りを 第2種の過誤(error of the second kind)偽陰性 とよぶ。

過誤まとめ
  1. 第1種の過誤(type I error)

    • 帰無仮説が正しいのに帰無仮説を棄却してしまう誤り。

    • これを犯す確率が有意水準α\alpha

  2. 第2種の過誤(type II error)

    • 帰無仮説が正しくないのに帰無仮説を受容してしまう誤り。こちらはβ\betaで表される。

これら2種類の誤りを統一的に表現する関数が検出力関数

損失関数

実際には正しくない仮説を採択する意思決定をとったときを「損失」として表し、意思決定の損失関数を以下のような損失関数L(θ,d)L(\theta, d)で表現することにする。

L(θ,0)={0, if θΘ0a, if θΘ1L(θ,1)={0, if θΘ1b, if θΘ0\begin{aligned} L(\theta, 0) &= \begin{cases} 0, & \text { if } \theta \in \Theta_0 \\ a, & \text { if } \theta \in \Theta_1 \end{cases} \\ L(\theta, 1) &= \begin{cases} 0, & \text { if } \theta \in \Theta_1 \\ b, & \text { if } \theta \in \Theta_0 \end{cases} \end{aligned}

(※第1種の過誤を犯したときaaの損失が、第2種の過誤を犯したときbbの損失が得られるということ。)

検出力関数

検定問題における決定関数 (検定関数 test function とよばれる) をδ:XD\delta: \mathcal{X} \to Dとおく。検定関数δ\deltaX=xX=xを観測したとき帰無仮説を棄却するなら1、受容するなら0をとる関数。

損失関数の期待値である リスク関数R(θ,δ)=Eθ[L(θ,δ(X))]R(\theta, \delta) = E_{\theta}[L(\theta, \delta(X))]と定義される。

ところで、

βδ(θ)=Eθ[δ(X)]=Pθ(δ(X)=1)\beta_\delta(\theta)=E_\theta[\delta(X)]=P_\theta(\delta(X)=1)

と定義すると、δ\deltaのリスク関数は以下のように表すことができる。

R(θ,δ)={βδ(θ), if θΘ0 (=第1種の過誤の確率) 1βδ(θ), if θΘ1 (=第2種の過誤の確率) R(\theta, \delta)= \begin{cases} \beta_\delta(\theta), & \text { if } \theta \in \Theta_0 & \text{ (=第1種の過誤の確率) } \\ 1-\beta_\delta(\theta), & \text { if } \theta \in \Theta_1 & \text{ (=第2種の過誤の確率) } \end{cases}

βδ(θ)\beta_\delta(\theta)検出力関数(power function) といい、その値を 検出力(power) という。

例:不良率の検定

10個の製品を検査して、もし1個以上の不良品があれば帰無仮説 H0:pp0H_0: p \leq p_0 を棄却するとする(なお p0=0.01p_0 = 0.01と想定する)

この場合の検定関数は

δ(x)={1, if x10, otherwise \delta(x) = \begin{cases} 1, & \text { if } x \geq 1 \\ 0, & \text { otherwise } \end{cases}

となる。XXは2項分布B(10,p)B(10, p)に従うと考えればよいので、検出力関数は

βδ(p)=P(X1)=1P(X=0)=1(1p)10\beta_\delta(p)=P(X \geq 1)=1-P(X=0)=1-(1-p)^{10}

で与えられる。p0=0.01p_0 = 0.01のときβδ(p)=1(1p)10=0.0956\beta_\delta(p) = 1-(1-p)^{10} = 0.0956