ホルム＝ボンフェローニ法

統計学において、**ホルム＝ボンフェローニ法**（ホルム＝ボンフェローニほう、）は、多重比較問題に対抗するために使われる手法である。**ホルムの方法**または**ボンフェローニ＝ホルム法**とも呼ばれる。ファミリーワイズエラー率（FWER）を制御することが意図されており、ボンフェローニ補正よりも（検出力の高い）単純な検定を与える。名称は本手法を体系化したスウェーデンの統計学者スチューレ・ホルム（Sture Holm）とボンフェローニの不等式ひいてはカルロ・エミリオ・ボンフェローニにちなむ。

## 動機
いくつかの仮説を考慮する時、多重性の問題が生じる。つまり、より多くの仮説を調べる程、第一種過誤（偽陽性）を得る確率がより高くなる。ホルム＝ボンフェローニ法は個別の仮説のそれぞれに対する棄却基準を調整することによってファミリーワイズエラー率（1つ以上第一種過誤を犯す確率）を制御するための多くの手法の1つである。

## 定式化
本手法は以下の通りである。
* 小さい順P_1,\ldots,P_mに並べられたm個のp値とそれらに対応する仮説H_1,\ldots,H_mを持っていることを考える。ファミリーワイズエラー率は事前に設定した特定の有意水準 \alpha 以下にしたい。
* P_1 であるならば、H_1を棄却し、次の段階に進む。さもなければここで検定を止める。
* P_2 であるならば、H_2も棄却し、次の段階に進む。さもなければここで検定を止める。
* これを繰り返す。それぞれのP値について、P_k かどうかを検定し、条件を満たせばH_kを棄却し、次のより大きなP値について調べる。さもなければここで検定を止める。

この手法はファミリーワイズエラー率が \leq\alpha であることを保証する。

### 理論的根拠
単純なボンフェローニ補正は、1つ以上の真である帰無仮説を棄却する（すなわち、1つ以上の第一種過誤を犯す）危険が最大でも\alphaであることを保証するために、*p*値が\frac{\alpha}{m}未満の帰無仮説のみを棄却する。この第一種過誤に対する保護の代償は、1つ以上の偽である帰無仮説を棄却し損う（すなわち、1つ以上の第二種過誤を犯す）危険が増大することである。

ホルム＝ボンフェローニ法は水準\alphaで最大ファミリーワイズエラー率も制御するが、古典的なボンフェローニ法よりも第二種過誤の危険の増大がより小さい。ホルム＝ボンフェローニ法は*p*値を小さい順番に並べ、それぞれ\frac{\alpha}{m}から\alphaの名目α水準（すなわち\frac{\alpha}{m},  \frac{\alpha}{m-1}, \ldots ,  \frac{\alpha}{2},  \frac{\alpha}{1}）と比較する。
* 指数kは、棄却を有効とするのに十分低くない最初の*p*値を特定する。結果として、帰無仮説H_{(1)}, \ldots , H_{(k-1)}が既約されるのに対して、帰無仮説 H_{(k)}, ... , H_{(m)}は判断が保留される（棄却されない）。
* もしk = 1ならば、棄却のために十分小さな*p*値はなく、そのため棄却される帰無仮説はない（すなわち全ての帰無仮説について判断が留保される）。
* こういった指数kを見つけることができなかったならば、全ての*p*値が棄却のために十分小さく、したがって全ての帰無仮説が棄却される。

### 証明
ホルム＝ボンフェローニ法はFWERを以下のように制御する。H_{(1)}\ldots H_{(m)}を仮説の族、P_{(1)}\leq P_{(2)}\leq\cdots\leq P_{(m)}を並び換えされたp値とする。I_0をm_0個の成員を有する（未知の）真である仮説に対応する指数の組とする。

真である仮説を誤って棄却することを仮定する。この事象の確率が最大でも\alphaであることを証明しなければならない。hを最初に棄却された真である仮説（ボンフェローニ＝ホルム検定によって与えられる順序での最初）とする。すると、H_{(1)}, \ldots, H_{(h-1)}は全ての棄却された偽である仮説であり、h-1 \leq m -m_0である。そこから、\frac{1}{m-h+1}\leq \frac{1}{m_0} (1) を得る。hは棄却されているため、本検定の定義によりP_{(h)} \leq \frac{\alpha}{m-h+1}を得る。(1) 式を使うと、右辺は最大でも\frac{\alpha}{m_0}である。したがって、もし真である仮説を誤って棄却するならば、最大でも\frac{\alpha}{m_0}のP値を持つ真である仮説が存在しなければならない。

そこで、確率変数A=\left\{ P_i \leq \frac{\alpha}{m_0} \text{ for } i\in I_0\right\}を定義する。真の仮説の（未知の）集合I_0が何であれ、（ボンフェローニの不等式により）\Pr(A)\leq \alphaとなる。その結果、真である仮説を棄却する確率は最大でも\alphaである。

## 例
未調整p値p_1=0.01、p_2=0.04、 p_3=0.03、p_4=0.005を持つ4つの帰無仮説が有意水準\alpha=0.05で検定されることを考える。本手順はステップダウン（下降）式であうため、初めに最小のp値p_4=p_{(1)}=0.005を持つH_4=H_{(1)}を検定する。p値は\alpha/4=0.0125と比較され、この帰無仮説は棄却されて、次の仮説に進む。p_1=p_{(2)}=0.01であるため、H_1=H_{(2)}も同様に棄却され、次に進む。次の仮説H_3はp_3=p_{(3)}=0.03 > 0.025=\alpha/2であるため棄却されない。ここで検定を止め、H_1とH_4は棄却され、H_2とH_3は棄却されないと結論付ける。ここでファミリーワイズエラー率は水準\alpha=0.05で制御されている。ここで留意すべきは、p_2=p_{(4)}=0.04 であるにもかかわらず、H_2は棄却されないという点である。これは、棄却ができなかったらそこで検定手順が停止するためである。

## 拡張

### ホルム＝シダックの方法

仮説検定が負の依存関係にある時、\frac{\alpha}{m},\frac{\alpha}{m-1}, \ldots, \frac{\alpha}{1}を
: 1-(1-\alpha)^{1/m},1-(1-\alpha)^{1/(m-1)},\ldots,1-(1-\alpha)^{1}
で置き換えることが可能となる。これによって、わずかにより強力な検定となる。

### 重み付け版
P_{(1)},\ldots,P_{(m)}を並び換えされた未調整p値とする。H_{(i)}について、0\leq w_{(i)}をP_{(i)}に対応させる。
: P_{(j)}\leq \frac{w_{(j)}}{\sum^m_{k=j} w_{(k)}} \alpha,\quad j=1,\ldots,i
であるならば、H_{(i)}を棄却する。

### 調整''p''値
ホルム＝ボンフェローニ法に対する調整''p''値は
: \widetilde{p}_{(i)}=\max_{j\leq i}\left\{ (m-j+1)p_{(j)}\right\}_1, \text{ where } \{x\}_1 \equiv \min(x,1)
である。

前の例では、調整*p*値は\widetilde{p}_1 = 0.03、\widetilde{p}_2 = 0.06、\widetilde{p}_3 = 0.06、\widetilde{p}_4 = 0.02となる。仮説H_1およびH_4のみが水準\alpha=0.05で棄却される。

重み付けされた*p*値は、
:\widetilde{p}_{(i)}=\max_{j\leq i} \left\{\frac{\sum^m_{k=j}{w_{(k)}}}{w_{(j)}} p_{(j)}\right\}_1, \text{ where } \{x\}_1 \equiv \min(x,1)
である。調整*p*値がα未満である時かつその時に限り仮説は水準αで棄却される。等しい重みを用いた前の例では、調整*p*値は0.03、0.06、0.06、0.02である。これは、α = 0.05を使って、この手順によって仮説1および4のみが棄却されることを見るための別のやり方である。

## 代替法と用法

ホルム＝ボンフェローニ法は古典的なボンフェローニ補正よりも「一様に」より検出力が高い。これは、常に少なくとも同等に検出力が高いことを意味する。

ファミリーワイズエラー率を制御するためのホルム＝ボンフェローニ法よりも強力なその他の手法が存在する。例えば、ホッホベルクのステップアップ手順では、H_{(1)} \ldots H_{(k)}の棄却はP_{(k)} \leq \frac{\alpha}{m+1-k}であるような「最大の」指数kを見つけた後に成される。したがって、ホッホベルクの手順はホルムの手順よりも一様により強力である。しかしながら、ホッホベルクの手順は仮説が独立である、または正の依存性を持つ特定の形式の下にあることを必要とするが、ホルム＝ボンフェローニ法はそういった仮定なしに適用することができる。同様のステップアップ手順にホンメル（Hommel）の手順がある。これはホッホベルクの手順よりも一様により強力である。

## 名称
カルロ・エミリオ・ボンフェローニは本項で記述したホルム＝ボンフェローニ法の考案には関与していない。ホルムは元々本手法を「逐次棄却型ボンフェローニ検定」と呼び、しばらくした後でホルム＝ボンフェローニ法と呼ばれるようになった。自身の手法をボンフェローニに因んで命名したホルムの動機は原論文において以下のように説明されている: 『多重推測理論内でのブールの不等式の使用は大抵ボンフェローニ・テクニックと呼ばれ、この理由からこの検定を逐次棄却型ボンフェローニ検出と呼ぶことにする』。

## 出典