確率・統計 · 検出力

検出力

検出力とは、ある特定の対立仮説が真のときに、検定が帰無仮説H0を棄却する確率です。ここでは機械が平均500 gずつ袋に詰めるはずだとします。真の平均をドラッグし、n袋の調査がどのくらいの頻度で差を検出するかを見てみましょう。

この図で見る検出力

検出力とは、ある特定の対立仮説が真のときに、検定が帰無仮説H0H_0を棄却する確率です。ここでH0H_0は、機械が平均μ=500\mu = 500 gずつ袋に詰めるという仮説です。各袋の重量は独立で、σ=10\sigma = 10 gの正規分布に従います。1回の調査でnn袋を量り、平均xˉ\bar xが臨界値を越えたら、有意水準α\alphaでH0H_0を棄却します。真の平均がμ\muのときの検出力は、真の曲線のうち臨界値を越えた部分の面積です。残りのβ\betaは差を見逃す確率、つまり第2種の過誤の確率です。μ=500\mu = 500では見逃す差がなく、検出力はα\alphaそのものになります。

検出力は何で決まるのでしょうか。H1 ⁣:μ<500H_1\!: \mu < 500の検定で不足量がδ=500−μ\delta = 500 - \muなら、検出力はΦ(δn/σ−z1−α)\Phi(\delta\sqrt{n}/\sigma - z_{1-\alpha})です。差が大きいほど、袋数が多いほど高くなります。α\alphaを厳しくすると臨界値が500から遠ざかり、検出力は下がります。nnはn\sqrt{n}の形でしか効かないので、同じ検出力に必要な袋数は1/δ21/\delta^2に比例して増えます。方向も同じくらい重要です。実際には機械が多く詰めている場合、μ<500\mu < 500の検定の検出力は0に近づき、袋数を増やすほどさらに下がります。両側検定はα\alphaを両方の裾に分けるので、どちら側の差も検出できますが、片側検定が見る側では、その片側検定より検出力が低くなります。

有意でない結果が「差がない」ことを示さないのはなぜでしょうか。検出力が低い検定は、実際の差を見逃すことが多いからです。25袋、α=0.05\alpha = 0.05の片側検定で2 gの不足を検出する確率は約4分の1にすぎず、ほとんどの調査は有意になりません。有意でない結果をどう受け止めるかは、見つける価値のある差に対する検出力によります。それはデータを見る前に計算します。重要な最小の差を決め、十分な検出力、よく使われる目安なら80%を得られるnnを選びます。

検出力が低いと、有意な推定値は平均して差を過大評価するのはなぜでしょうか。実際に充填量が不足している機械について、不足を調べる検定を考えます。平均がたまたま500よりかなり低くなった調査だけが臨界値を越えるため、有意な調査は選ばれた一群です。その平均推定値は真の曲線の臨界値を越えた部分の平均で、実際より大きな不足量に対応します。臨界値自体が下限を設けています。25袋、α=0.05\alpha = 0.05では、実際の不足量によらず、約3.29 gを超える不足と推定しないと有意になりません。nnとα\alphaを固定すると、真の不足量が小さいほど検出力が低くなり、有意な推定値が平均して何倍に過大評価するかも大きくなります。これがM型(大きさ)の誤りです。両側検定では有意なのに方向まで逆の結果もあり、S型(符号)の誤りと呼ばれます。有意な結果だけを報告する研究の蓄積には、この偏りが含まれます。

注目するところ

  • 上段の2つの曲線は、どのnnでも同じ高さで描かれます。色の付いた面積は、それぞれの曲線全体に対する割合として読んでください。ばらつきの減少を示すのは高さではなく幅の縮まりです。
  • 検出力は、機械の充填量が不足している確率ではありません。真の平均がページで設定したμ\muだと仮定して計算します:P(棄却H0∣μ)P(\text{棄却} H_0 \mid \mu)。

操作方法

  • 青い曲線やハンドルをドラッグして真の平均を動かすか、下段を押してその位置へ移します。臨界値の四角いハンドルをドラッグするとα\alphaを設定できます。
  • 「検出力80%」の行は、現在の平均、α\alpha、方向で検出力が初めて80%に達する標本サイズにnnを設定します。

関連項目