[論文レビュー] Expected Coverage of Bayesian Confidence Intervals for the Mean of a Poisson Statistic in Measurements with Background

Ilya V. Narsky|ArXiv.org|May 11, 2000

Scientific Measurement and Uncertainty Evaluation参考文献 4被引用数 3

ひとこと要約

本稿では、バックグラウンドを伴うポアソン平均のベイジアン信頼区間の期待被覆率と区間長を評価し、平坦、$1/(s+b)$、$1/\sqrt{s}$、$1/\sqrt{s+b}$ の4つの目的的事前分布と、Feldman-Cousins統合的手法および古典的手法の2つの頻度主義的手法を比較している。$1/\sqrt{s+b}$ 事前分布が、信号タイプにわたって最もバランスの取れた被覆率を示し、名目被覆率レベルを維持する点で優れていることが判明した。

ABSTRACT

Expected coverage and expected length of 90% upper and lower limit and 68.27% central intervals are plotted as functions of the true signal for various values of expected background. Results for several objective priors are shown, and formulas for calculation of confidence intervals are obtained. For comparison, expected coverage and length of frequentist intervals constructed with the unified approach of Feldman and Cousins and a simple classical method are also shown. It is assumed that the expected background is accurately known prior to performing an experiment. The plots of expected coverage and length are provided for values of signal and background typical for particle experiments with small numbers of events.

研究の動機と目的

• 既知のバックグラウンドを伴うポアソン分布に従う信号のベイジアン信頼区間の期待被覆率と区間長を評価すること。
• 平坦、$1/(s+b)$、$1/\sqrt{s}$、$1/\sqrt{s+b}$ の4つの目的的事前分布の被覆率と区間長という観点から、それらの性能を比較すること。
• 頻度主義的手法（Feldman-Cousins統合的手法および古典的手法）がベイジアン手法に対してどのように性能を発揮するかを評価すること。
• 異なる信号およびバックグラウンドレベルにおいて、最も信頼性が高く一貫性のある被覆率を提供する事前分布を特定すること。
• 片側区間および中央区間において、区間長と被覆率の正確性のトレードオフを定量化すること。

提案手法

• 既知のバックグラウンド $b$ を伴うポアソン尤度 $f(n|s) = e^{-(s+b)}(s+b)^n/n!$ を使用する。
• ガンマ関数および不完全ガンマ関数を用いて、平坦、$1/(s+b)$、$1/\sqrt{s}$、$1/\sqrt{s+b}$ の4つの目的的事前分布の事後分布を導出する。
• 事後中央分位数を用いて信頼区間を計算する：$s_1(n)$ および $s_2(n)$ が $\int_0^{s_1} \pi(s|n)ds = \alpha_1$、$\int_{s_2}^\infty \pi(s|n)ds = \alpha_2$ を満たすようにする。
• 期待被覆率 $C(s) = \sum_n I(n,s) f(n|s)$ および期待区間長 $L(s) = \sum_n (s_2(n) - s_1(n)) f(n|s)$ を計算し、$L_1(s)$ は下限用とする。
• $n=100$ まで数値的合計を用いて、$b = 0, 1, 3, 6$ の場合に $0 \leq s \leq 10$ の範囲で期待被覆率と区間長をプロットする。
• Feldman-Cousins統合的手法および古典的頻度主義手法と比較する。ネイマンベルト構築法を用いる。

実験結果

リサーチクエスチョン

• RQ1 異なる目的的事前分布は、バックグラウンドを伴うポアソン測定におけるベイジアン信頼区間の期待被覆率にどのように影響するか？
• RQ2 どの事前分布が信号値や区間タイプ（上側、下側、中央）にわたって最も一貫性のある被覆率を提供するか？
• RQ3 区間長は、事前分布および信号レベルに応じてどのように変化するか？
• RQ4 ベイジアン手法と頻度主義的手法（Feldman-Cousinsおよび古典的手法）は、被覆率と区間長の観点でどのように比較できるか？
• RQ5 名目被覆率を維持しながら過剰被覆または不十分被覆を最小限に抑える観点で、$1/\sqrt{s+b}$ 事前分布は他の事前分布を上回るか？

主な発見

• $1/\sqrt{s+b}$ 事前分布は、すべての区間タイプおよび信号値において最もバランスの取れた期待被覆率を示し、顕著な過剰被覆や不十分被覆を回避している。
• 平坦事前分布は下限に対して不十分被覆を示し、上限に対して過剰被覆を示すが、ジェフレーズの事前分布（$1/(s+b)$）は逆の傾向を示す。
• $1/\sqrt{s}$ 事前分布は、特に低信号値において、上限に対して不十分被覆を示し、下限に対して過剰被覆を示す傾向がある。
• 中央区間に関しては、すべてのベイジアン手法が名目信頼水準の周りを振動的に被覆するが、$1/\sqrt{s+b}$ が最も安定した挙動を示している。
• すべてのベイジアン手法の期待区間長は、信号が増加するに従い漸近的に収束する。平坦事前分布は最も長い区間長を生じ、ジェフレーズの事前分布は最も短い区間長を生じる。
• Feldman-Cousins統合的手法は、高信号値において他の手法の片側区間と比較してより短い両側区間を生成するが、両側性のため長さの直接比較はできない。