ボンフェローニ補正

提供: testwiki
ナビゲーションに移動 検索に移動

統計学において、ボンフェローニ補正(ボンフェローニほせい、テンプレート:Lang-en-short)は、多重比較問題に対抗するために使われるいくつかの手法のうちの1つである。

背景

本手法の名称はボンフェローニの不等式[1]を使用することにちなむ。本手法の信頼区間への拡張はテンプレート:仮リンクによって提唱された[2]

統計的仮説検定は、観察されたデータの帰無仮説の下での尤度が低ければ帰無仮説を棄却することに基づく。複数の仮説が検定されるとすると、稀な事象を観察する可能性が高まり、その結果として、帰無仮説を誤って棄却する(すなわち第一種過誤を犯す)可能性が高まる[3]

ボンフェローニ補正は、有意水準α/mαは望ましい全体としてのα水準、mは仮説の数)で個々の仮説を検証することによって第一種過誤を犯す可能性の高まりを補償する[4]。例えば、1回の試行がm=20個の仮説を望むα=0.05の水準で検定しているとすると、ボンフェローニ補正は個別の仮説をα=0.05/20=0.0025の水準で検定する。同じように、複数の信頼区間を構築する時、同じ現象が表われる。

定義

H1,,Hmを仮説の族(ファミリー)、p1,,pmをそれらの対応するp値とする。mを帰無仮説の総数、m0を真である帰無仮説の数とする。ファミリーワイズエラー率(FWER)は少くとも1つの真であるHiを棄却する確率、すなわち少くとも1つの第一種過誤を犯す確率である。ボンフェローニ補正はpiαmで帰無仮説を棄却することで、FWERを水準αで制御する。この制御の証明は以下のようにブールの不等式から得られる。

FWER=P{i=1m0(piαm)}i=1m0{P(piαm)}=m0αmmαm=α

この制御はp値間の依存性またはいくつの帰無仮説が真であるかに関していかなる仮定も必要としない[5]

拡張

一般化

個々の検定の水準がデータを見るよりも前に決定されるという条件で、α/m水準で個々の仮説を検証するよりむしろ、仮説は合計αとなる水準のどの組合せでも検定してもよい[6]。例えば、2つの仮説検定について、1つの検定を0.04、もう一方の検定を0.01の水準で実行することによって全体として0.05のαを維持することができる。

信頼区間

ダンによって提唱された手順[2](順位ベースの分散力についてのダンの手順[7]と混同してはならない)は信頼区間を調整するために使うことができる。m個の信頼区間を定め、全体の信頼水準を1αにしたいと望むとすると、個々の信頼区間は1αmの水準に調整することができる[2]

連続問題

連続パラメータ空間における信号を探索する時も、多重比較の問題(どこでも効果)が存在しうる。例えば、ある物理学者が幅広い範囲の質量を考慮することによって未知の質量の粒子を発見したいと見ているとする。これはノーベル賞を受賞したヒッグス粒子の検出の際に当て嵌る。こういった場合、試行の有効数mと事前–事後体積比を関連付けるベイズロジックを利用することによって連続パラメータに対して一般化されたボンフェローニ補正を適用することができる[8]

代替手法

テンプレート:Main ファミリーワイズエラー率を制御するためには複数の代替手法が存在する。例えば、ホルム=ボンフェローニ法シダック補正はボンフェローニ補正よりも普遍的に強力な手順である。これは常に少くとも強力であることを意味する。ボンフェローニの手順とは異なり、これらの手法は族毎の第一種過誤の期待数(族毎の第一種過誤の確率)を制御しない[9]

批判

FWER制御に関して、数多くの検定が存在する時と検定統計量が正に相関している時の両方またはどちらか一方の場合、ボンフェローニ補正は保守的であるかもしれない[10]

ボンフェローニ補正は偽陰性を生む確率を増大させる、すなわちテンプレート:仮リンクを低下させる犠牲を払う[11][10]。全ての場合において仮説族をどのように定義するかについて決定的な意見の一致は存在しないが、調整された検定結果は仮説の族に含められた検定数に依存して変動するかもしれないテンプレート:Citation needed。こういった批判は一般にFWER制御に向けられ、ボンフェローニ補正の特有のものではない。

出典

テンプレート:Reflist

参考文献

テンプレート:Refbegin

テンプレート:Refend

関連項目

外部リンク