[論文レビュー] Bayesian Statistics at Work: the Troublesome Extraction of the CKM Phase alpha

Jérôme Charles, A. Höcker+3 more|ArXiv.org|Jul 22, 2006

Particle physics theoretical and experimental studies参考文献 9被引用数 12

ひとこと要約

この論文は、CKM位相𝛼のベイズ推論が、特にミラー解と不変でないヤコビアンのため、パrameterizationと事前分布の選択に極めて敏感であることを示している。特定のparameterizationにおける一様事前分布は、𝛼=0°で不適切で発散する事後分布をもたらし、物理的対称性に反し、頻度主義的解と整合しない結果をもたらす。

ABSTRACT

In Bayesian statistics, one's prior beliefs about underlying model parameters are revised with the information content of observed data from which, using Bayes' rule, a posterior belief is obtained. A non-trivial example taken from the isospin analysis of B-->PP (P = pi or rho) decays in heavy-flavor physics is chosen to illustrate the effect of the naive "objective" choice of flat priors in a multi-dimensional parameter space in presence of mirror solutions. It is demonstrated that the posterior distribution for the parameter of interest, the phase alpha, strongly depends on the choice of the parameterization in which the priors are uniform, and on the validity range in which the (un-normalizable) priors are truncated. We prove that the most probable values found by the Bayesian treatment do not coincide with the explicit analytical solution, in contrast to the frequentist approach. It is also shown in the appendix that the alpha-->0 limit cannot be consistently treated in the Bayesian paradigm, because the latter violates the physical symmetries of the problem.

研究の動機と目的

• B→PP崩壊からのCKM位相𝛼を抽出する際のベイズ推論の頑健性を、アイソスピン解析を用いて調査すること。
• 一様で情報のない事前分布を用いる際、事後分布がparameterizationの選択に依存することを明らかにすること。
• 数学的アーティファクトに起因して、ベイズ的手法が𝛼=0°で物理的に不適切な結果——具体的には発散する事後分布——を生じることを示すこと。
• ミラー解が存在する場合、ベイズ的結果と頻度主義的解とが一致しないことを対比して示すこと。
• 振幅モデルの物理的対称性に反するため、ベイズ枠組みが𝛼→0極限を一貫して取り扱えないことの失敗を強調すること。

提案手法

• ベイズの定理を用いて、B→PP崩壊振幅に関する事前信念を、CKM位相𝛼の事後分布に更新する。
• 複数のparameterization（例：RIおよびMA）で一様で情報のない事前分布を適用し、事前分布への感受性を評価する。
• 振幅parameterization間の変換におけるヤコビアンを導出し、事後分布の形状に影響を与える度量依存性を明らかにする。
• RI parameterizationにおける𝛼の事後確率密度関数を分析し、|sin𝛼|に反比例するスケーリングを示し、𝛼=0°および𝛼=180°で発散することを示す。
• 非正規化された事前分布を有限範囲に切断し、範囲を広げると事後分布のピークが端点に近づくことを観察する。
• 複数のparameterization間で結果を比較し、MA parameterizationは発散を回避するが、RI parameterizationは不適切な事後分布を生じることを示す。

実験結果

リサーチクエスチョン

• RQ1 ベイズ的解析におけるB→PP崩壊のCKM位相𝛼の事後分布は、parameterizationの選択にどのように影響を受けるか？
• RQ2 なぜRI parameterizationにおいて𝛼=0°で事後分布が発散するのか？この特異性の原因は何か？
• RQ3 ミラー解が存在する状況下で、多次元parameter spaceにおける一様事前分布は、物理的に不適切な結果をどれほど引き起こすか？
• RQ4 ベイズ的事後分布における𝛼は、頻度主義的手法で得られる明示的な解析解とどのように比較できるか？
• RQ5 振幅モデルの基本的物理的対称性を考慮すると、ベイズ枠組み内で𝛼→0極限を一貫して取り扱えるか？

主な発見

• RI parameterizationにおける𝛼の事後分布は|sin𝛼|に反比例し、𝛼=0°および𝛼=180°で発散するため、不適切な確率密度関数となる。
• 発散は、parameterizationに依存するヤコビアンに起因し、再parameterizationに対して不変でない。
• 振幅の事前分布の範囲を広げると、𝛼の事後分布のピークが0°および180°に近づくことが示され、切断限界への強い感受性が示された。
• ベイズ的取り扱いによって得られた𝛼の最も確実な値は、頻度主義的手法で得られた解析解と一致しない。
• RI parameterizationは正規化されていない事後分布を生じるが、MA parameterizationは有限で良好に振る舞う事後分布を生じ、𝛼=0°で消えるため、物理的期待と整合する。
• 𝛼=0のとき、ベイズ枠組みは𝛼→−𝛼の物理的対称性を尊重できず、標準模型の振幅構造と矛盾することを示している。