[論文レビュー] On the Combining Significances

Sergey Bityukov, Nikolai Krasnikov+1 more|ArXiv.org|Dec 18, 2006

Random Matrices and Applications被引用数 6

ひとこと要約

本稿では、独立した実験からの複数の信号の有意水準を、観測された有意水準を単位分散の正規分布に従う確率変数として扱うことで、統計的枠組みを提案する。シグナルの有意水準を実験数の平方根で割った和を用いることで、妥当な統合有意水準が得られることを示し、モンテカルロシミュレーションおよび信頼度密度解析によって、計数およびp値ベースの有意水準推定器の両方で検証されている。

ABSTRACT

We present the statistical approach to the combining of signal significances.

研究の動機と目的

• 独立した実験からの複数の観測有意水準を統合する厳密な統計的枠組みを確立すること。
• 初期（期待）有意水準と観測（実現）有意水準（確率変数としての）の違いを明確にすること。
• ポアソン分布に従う事象数からの観測有意水準が、単位分散の正規分布に近似されることを検証すること。
• 正規分布の性質を用いて、複数の部分的有意水準を一つの統合有意水準に統合する手法を開発すること。
• 多数のモンテカルロシミュレーションおよび信頼度密度解析を通じて、手法の妥当性を示すこと。

提案手法

• 著者らは、ポアソン分布に従う事象数から得られる観測有意水準を、2つの有意水準定義（$\hat{S}_{c12} = 2(\sqrt{\hat{N}_{\text{obs}}} - \sqrt{N_b})$ および $\hat{S}_{cP}$、Poissonのp値をガウスのシグマに変換）に基づく確率変数の実現値としてモデル化する。
• 初期有意水準の均一なスキャンを実施し、パラメータごとに30,000回のシミュレーションを実行して、初期有意水準の関数としての観測有意水準の信頼度密度を構築する。
• この手法は、正規分布の統計的自己双対性に依存しており、初期有意水準の信頼度密度が観測有意水準の分布と一致することを可能にする。
• 有意水準の統合に関して、著者らは $\hat{S}_{\text{comb}} = \hat{S}_{\text{sum}} / \sqrt{\text{var}(\hat{S}_{\text{sum}})}$ を定義し、$\text{var}(\hat{S}_{\text{sum}}) = \sum_{i=1}^n \text{var}(\hat{S}_i)$ となるようにする。また、ポアソン過程では分散が概ね1に近いことを示している。
• 最大300万回の試行を用いたモンテカルロシミュレーションにより、複数の実験における有意水準の和と正規化された分布を比較することで、手法の妥当性を検証している。
• この手法は $\hat{S}_{c12}$ および $\hat{S}_{cP}$ の両方のケースに適用され、観測有意水準の分布が単位分散の正規分布でよく近似されることを確認している。

実験結果

リサーチクエスチョン

• RQ1 独立した実験からの観測有意水準を、意味のある方法で一つの統合有意水準に統合する方法は何か？
• RQ2 ポアソン分布に従う事象数からの観測有意水準は、どの程度単位分散の正規分布に従うか？
• RQ3 初期（期待）有意水準と観測（実現）有意水準（確率変数としての）との統計的関係は何か？
• RQ4 観測有意水準の実現値から初期有意水準の信頼度密度をどのように構築できるか？
• RQ5 ポアソン統計のもとで、$\bar{S} = \Sigma S_i / \sqrt{n}$ によって有意水準を統合する標準的手法は、有効な推定器を提供するか？

主な発見

• 観測有意水準 $\hat{S}_{c12}$ および $\hat{S}_{cP}$ は、分散が1に近い正規分布でよく近似される確率変数の実現値である。
• 初期有意水準 $S_{c12}$ の信頼度密度は、観測有意水準 $\hat{S}_{c12}$ の分布と一致しており、$\hat{S}_{c12}$ が $S_{c12}$ の有効な推定器であることが確認された。
• モンテカルロシミュレーションにより、独立した実験における観測有意水準の和の分布が概ね正規分布に従うことが確認され、分散は実験数に等しい。
• 正規化された和 $\hat{S}_{\text{comb}} = \hat{S}_{\text{sum}} / \sqrt{n}$ は、真の統合有意水準を中心として分布し、有意水準統合手法の妥当性が裏付けられた。
• 純粋なバックグラウンド（信号なし）の場合、観測有意水準 $\hat{S}_{c12}$ および $\hat{S}_{cP}$ は標準正規分布に従うことが確認され、近似の妥当性が裏付けられた。
• 異なるバックグラウンドおよび信号レベルにおいても、複数のスキャンパラメータおよびシミュレーション試行で一貫した結果が得られたことから、手法のロバスト性が示された。