【多重比較補正】Bonferroni・Holm・FDRの選択と論文記述
複数の仮説検定を同時に実施すると、偶然による偽陽性(Type I error)が蓄積します。α = 0.05で20回検定すれば、平均1回は偽陽性が生じます。多重比較補正はこの問題を制御する重要な手法であり、特にオミクス研究・サブグループ解析・探索的研究で必須です。
1. 多重比較問題とは
α = 0.05で独立した検定をk回実施した場合の偽陽性確率:
- k=1: 5%
- k=5: 1-(0.95)⁵ = 23%
- k=10: 1-(0.95)¹⁰ = 40%
- k=20: 1-(0.95)²⁰ = 64%
- k=100(ゲノム研究): ≈ 99%
補正なしでは、多くの「有意な」発見が偶然による偽陽性です。
2. Family-wise Error Rate(FWER)制御
| 手法 | 補正方法 | 特徴 | 適用場面 |
|---|---|---|---|
| Bonferroni | α’ = α/k | 最も保守的・単純 | 少数の独立した検定(k≤10) |
| Holm(逐次Bonferroni) | 順位付けして逐次補正 | Bonferroniより検出力高い | Bonferroniと同じ場面、より多い検定 |
| Šídák | α’ = 1-(1-α)^(1/k) | 独立検定に正確 | 完全独立な検定のみ |
5つのアウトカムを検定する場合: α’ = 0.05/5 = 0.01
各検定でP < 0.01を有意と判断します。
報告例(英語): “To control the family-wise error rate at 0.05 across five primary outcomes, a Bonferroni correction was applied, with statistical significance set at P < 0.010 for individual tests.”
3. False Discovery Rate(FDR)制御
| 手法 | 制御する指標 | 特徴 | 適用場面 |
|---|---|---|---|
| Benjamini-Hochberg (BH) | FDR = E[V/R] | FWERより検出力高い | ゲノム研究、多数のアウトカム |
| Benjamini-Yekutieli | FDR(相関ある検定) | 検定間の相関を考慮 | 相関する多数の検定 |
| q値 | 各検定のFDR | 個々の発見のFDP推定 | マイクロアレイ・プロテオミクス |
Benjamini-Hochberg法の手順:
- k個のP値を昇順に並べ: P₍₁₎ ≤ P₍₂₎ ≤ … ≤ P₍k₎
- 各P₍ᵢ₎について: i/k × α を計算
- P₍ᵢ₎ ≤ i/k × α を満たす最大のiを見つける
- そのiまでのすべての仮説を有意と判断
報告例(英語): “Multiple comparisons were controlled using the Benjamini-Hochberg procedure to maintain a false discovery rate of 5% (q < 0.05).”
4. 手法の比較と選択ガイド
| 状況 | 推奨手法 | 理由 |
|---|---|---|
| 主要アウトカムが3〜5個(確証的研究) | Bonferroni / Holm | FWER制御が必須 |
| 主要アウトカムが10個以上 | Holm(Bonferroniより検出力高い) | 保守的すぎない |
| ゲノム研究(1000個以上の検定) | BH(FDR制御) | FWERは保守的すぎる |
| 探索的研究・サブグループ解析 | BH(FDR制御) | 発見の機会を増やす |
| ANOVAの事後検定 | Tukey HSD / Dunnett | 特化した事後検定を使用 |
多重比較補正のジャーナル要件確認
ジャーナルによって多重比較補正の要件が異なります。投稿規程の読み方ガイドで確認してください。また統計報告完全ガイドで補正後P値の記述方法を確認しましょう。
5. ANOVAの事後検定(Post-hoc Tests)
| 事後検定 | 特徴 | 適用場面 |
|---|---|---|
| Tukey HSD | 全ペア比較、等分散仮定 | 標準的な全対全比較 |
| Dunnett | 対照群との比較のみ | 介入vs対照の多数比較 |
| Games-Howell | 等分散不要 | 等分散性が満たされない場合 |
| Bonferroni(事後) | 最も保守的 | 少数の計画された比較 |
| Šídák | Bonferroniより若干緩い | 独立した比較 |
6. 論文での多重比較補正の記述
“For the analysis of five pre-specified primary outcomes, a Bonferroni correction was applied to control the family-wise error rate at 5%, with statistical significance defined as P < 0.010. For exploratory analyses involving 47 secondary outcomes, the Benjamini-Hochberg procedure was applied to maintain a false discovery rate of 5% (q < 0.05). Post-hoc pairwise comparisons following significant ANOVA results were conducted using Tukey’s Honestly Significant Difference test.”
多重比較補正の選択・適用で悩んでいませんか?
どの補正法を使うべきか、結果の報告方法はどうすべきか——Editverse Japanの統計専門家が、あなたの研究デザインに最適な多重比較補正の実装から論文記述まで完全サポートします。
参考PMID
| PMID | 著者 | 内容 |
|---|---|---|
| 27225100 | Wasserstein & Lazar (2016) | ASA声明 — 多重比較補正の必要性 |
| 26186117 | Ioannidis (2015) | 多重比較による偽陽性率の膨張問題 |
| 25892050 | Open Science Collaboration (2015) | 多重比較が再現性危機に与える影響 |
| 28817453 | Benjamin et al. (2018) | 有意水準変更と多重比較補正の関係 |
| 30581843 | Amrhein et al. (2019) | 多重比較補正後のCI報告の重要性 |
多重比較の見落としは査読で必ず指摘されます
サブグループ解析や複数アウトカムの研究では、多重比較補正の欠如が査読却下の原因になります。Editverse Japan統計サポートで投稿前に確認しましょう。
関連サービス:医学統計解析サポート
多重比較の補正漏れは、査読で指摘されます
「多重比較の補正が行われていません」は頻出の指摘です。Bonferroni、Holm、FDR など研究デザインに応じた補正手法の選択と実施、結果の再提示までを定額¥79,800で。再現可能なRスクリプトと出力ログを納品するため、追加の指摘にも対応できます。執筆まで含む「データから論文へ」パッケージは定額¥238,000です。