[論文レビュー] A shortcut for Hommel's procedure in linearithmic time

Rosa J. Meijer, Thijmen J. P. Krebs+2 more|arXiv (Cornell University)|Oct 23, 2017

Statistical Methods in Clinical Trials参考文献 28被引用数 4

ひとこと要約

本稿では、Hommelの多重仮説検定手順に対する線形対数時間の正確な高速化手法を提案する。この手法は、Hommel法の検出力とHochberg法の計算効率を統合し、O(m log m)時間で調整済みp値を計算する。臨界関数h(α)を効率的に計算し、2ポインタ法を用いて棄却臨界値を特定することで、大規模な多重仮説検定問題においてもFWER制御をスケーラブルに実現しつつ、検出力の損失なしに実現する。

ABSTRACT

Hommel's and Hochberg's procedures for familywise error control are both derived as shortcuts in a closed testing procedure with the Simes local test. Hommel's shortcut is exact but takes quadratic time in the number of hypotheses. Hochberg's shortcut takes only linearithmic time, but is conservative. In this paper we present an exact shortcut in linearithmic time, combining the strengths of both procedures. The novel shortcut also applies to a robust variant of Hommel's procedure that does not require the assumption of the Simes inequality.

研究の動機と目的

• 元の手法の2次時間計算量を回避する、計算的に効率的で正確なHommelの多重仮説検定手順の高速化手法を開発すること。
• Hommel法の統計的検出力とHochberg法の計算効率を統合し、線形対数時間計算量を達成すること。
• Simes不等式を必要としないHommelのロバストな変種への新しい高速化手法の拡張により、適用範囲を広げること。
• 遺伝学や神経画像解析などの大規模な多重仮説検定シナリオにおいて、Simesに基づく閉検定法のスケーラブルな適用を可能にすること。
• 既存の手法を上回る実行時間性能を発揮しながら正確な誤差制御を維持する、実用的でオープンソースのR実装を提供すること。

提案手法

• 臨界関数h(α)を計算するための新規アルゴリズムを導入し、これは各仮説が棄却される最小のαを決定するもので、O(m log m)時間で実行される。
• ソート済みのp値と調和平均重みに基づき、2ポインタ法を用いて各仮説iの臨界インデックスt_iを効率的に計算する。
• 各仮説iの調整済みp値は、min(s_{t_i} p_i, α_{t_i})として計算される。ここでs_jとα_jは、Simes不等式および調和級数から事前に計算された列である。
• p値の単調性とSimes局所検定の構造を活用することで、全閉検定を避けることができ、計算量をO(2^m)からO(m log m)に削減する。
• Simes不等式の代わりに、同時p値の依存性を仮定しない調和平均に基づく棄却基準を用いることで、Hommelのロバスト手順への一般化を実現する。
• 実装は2パス方式を採用する：まず線形対数時間スキャンによりh(α)を計算し、次にt_iインデックスを用いて全調整済みp値を線形時間で計算する。

実験結果

リサーチクエスチョン

• RQ1 Hommel法の正確な高速化手法を、Hochberg法の効率性に匹敵する線形対数時間で開発可能か？ これによりHommel法の検出力を損なわず、かつ実行速度を向上させられるか？
• RQ2 Simes不等式に依存しないHommel法のロバストな変種への新しい高速化手法を拡張可能か？
• RQ3 大規模な多重仮説検定問題において、実際の計算性能がHommel法およびHochberg法の手順と比較してどうなるか？
• RQ4 R言語でこの新しいアルゴリズムを効率的に実装可能か？ これにより遺伝学や神経画像解析における実世界の応用を支援できるか？

主な発見

• 提案された高速化手法は、調整済みp値をO(m log m)時間で計算する。これはHommel法のO(m²)計算量に比べ顕著な改善である。
• この手法は正確である。つまり、完全な閉検定手順と同一の棄却数を達成するが、Hochberg法の保守的高速化手法とは異なり、検出力の損失がない。
• この新しいアルゴリズムは、Hommel法のオリジナル手法を実行時間面で上回り、m > 5×10⁷の仮説に対しても計算時間が実用的である。
• Rパッケージhommelにおける実装は、p.adjustによるHochberg法の実装と同等またはそれ以上の性能を発揮し、検出力の損失なしに実現する。
• この手法により、遺伝学的関連解析や神経画像解析など、FWER制御が不可欠な大規模な設定においてHommel法の実用的応用が可能になる。
• 2ポインタ法により、h(α)が決定されれば、全調整済みp値を線形時間で効率的に計算できるため、全体の手続きは非常にスケーラブルである。