[論文レビュー] Interval estimation in the presence of nuisance parameters. 1. Bayesian approach

J. G. Heinrich, C. Blocker+5 more|ArXiv.org|Sep 27, 2004

Statistical Methods and Inference参考文献 8被引用数 14

ひとこと要約

本論文は、素粒子物理学における不確実な検出器効率やバックグラウンド率といった余分母を伴う状況下で区間推定を行うベイズ的手法を提案する。階層ベイズモデルと共役事前分布を用い、余分母を周辺化することで、信号率の正確な後部分布を導出し、平均的頻度的カバレッジを保証するとともに、系統的不確実性に対して頑健な信用区間を提供する。

ABSTRACT

We address the common problem of calculating intervals in the presence of systematic uncertainties. We aim to investigate several approaches, but here describe just a Bayesian technique for setting upper limits. The particular example we study is that of inferring the rate of a Poisson process when there are uncertainties on the acceptance and the background. Limit calculating software associated with this work is available in the form of C functions.

研究の動機と目的

• 系統的不確実性（例：効率とバックグラウンド）が存在する状況での信頼区間を信頼性高く構築する課題に対処すること。
• 補助的測定や点推定を必要としない、余分母を考慮するベイズ枠組みを構築すること。
• 得られた区間が、特に事前分布の上での平均カバレッジを満たす、適切な頻度的カバレッジ特性を維持することを保証すること。
• ポisson計数実験における系統的不確実性を伴う上限推定を計算する実用的で計算可能な手法を提供すること。
• 離散化されたデータが存在する状況において、頻度的手法と比較して、カバレッジと区間長の両面で本手法の性能を評価すること。

提案手法

• 信号率 $ s $、効率 $ \epsilon $、バックグラウンド $ b $ を適切な事前分布を備えた確率変数として扱う階層ベイズモデルを用いる。
• 解析的周辺化を可能にするために、$ s $、$ \epsilon $、$ b $ に対してガンマ分布を共役事前分布として適用する。
• 余分母 $ \epsilon $ と $ b $ を統合することで、$ s $ の周辺後部分布を導出し、超幾何関数を含む複雑な式を得る。
• 平均カバレッジ定理を用いて、事前分布の上での平均をとった場合に、ベイズ信用区間が平均的頻度的カバレッジを満たすことを証明する。
• 数値積分と特殊関数（例：超幾何関数）を用いて、観測度数 $ n $ に対する信用区間を計算し、明示的な公式を提示する。
• 実験的物理学における上限推定の再現可能な計算を可能にするために、C言語による実装コードを提供する。

実験結果

リサーチクエスチョン

• RQ1 系統的不確実性（効率とバックグラウンド）が存在する状況で、ベイズ的手法をどのように信号率の信頼区間を生成するように適応できるか？
• RQ2 余分母が存在する状況におけるベイズ信用区間のカバレッジ特性は何か？頻度的手法と比較するとどうなるか？
• RQ3 正確なカバレッジが真のパラメータ値に依存して変動する場合でも、ベイズ的手法が平均的頻度的カバレッジを維持できるか？
• RQ4 不適切な事前分布（例：一様事前分布）を用いることで、ベイズ区間のカバレッジ特性にどのような影響が生じるか？
• RQ5 特に低カウントのポアソン実験において、得られる区間の長さと妥当性は、頻度的手法の結果とどのように比較できるか？

主な発見

• 平均カバレッジ定理により、事前分布の上での平均をとると、ベイズ的手法が平均的頻度的カバレッジを持つことが証明された。
• 特定の真の $ s $ 値に対して個々のベイズ区間がカバレッジ不足となる場合でも、事前分布の上での平均をとることで、名目水準でのカバレッジが保証される。
• $ n $ が非負の整数で、形状母数が正である場合、$ s $ の周辺後部分布に対して正確な解析的表現（超幾何関数を含む）が得られる。
• 特に統計量が少ない状況で生じる可能性のある、空または退化した区間の問題を回避する。
• $ \omega \to \infty $ の極限では、$ b $ が固定された場合の標準的結果に一致し、既知の解と整合性が確認される。
• C言語による実装コードが提供されており、高エネルギー物理学実験における系統的不確実性を伴う上限推定の直接的計算が可能である。