[論文レビュー] KLTS: A rigorous method to compute the confidence intervals for the Three-Cornered Hat and for Groslambert Covariance

Éric Lantz, Claudio Calosso+5 more|arXiv (Cornell University)|Apr 11, 2019

Scientific Measurement and Uncertainty Evaluation参考文献 10被引用数 6

ひとこと要約

本論文は、カルフネン＝ローベ変換と十分統計量を用いた厳密なベイズ的手法であるKLTSを紹介する。この手法により、3クロック系におけるアラン分散推定値の正確な信頼区間が、Three-Cornered Hat法およびGroslambert共分散法を用いて計算可能となる。信頼区間は1自由度でも厳密に正であり、多数のモンテカルロシミュレーションによる検証と、冷凍サファイア発振器への応用が行われた。

ABSTRACT

The three-cornered hat / Groslambert Covariance methods are widely used to estimate the stability of each individual clock in a set of three, but no method gives reliable confidence intervals for large integration times. We propose a new KLTS (Karhunen-Lo\\`eve Tansform using Sufficient statistics) method which uses these estimators to take into account the statistics of all the measurements between the pairs of clocks in a Bayesian way. The resulting Cumulative Density Function (CDF) yields confidence intervals for each clock AVAR. This CDF provides also a stability estimator which is always positive. Checked by massive Monte-Carlo simulations, KLTS proves to be perfectly reliable even for one degree of freedom. An example of experimental measurement is given.

研究の動機と目的

• 長時間積分における等価自由度（EDF）が低い場合のThree-Cornered Hat法およびGroslambert共分散法推定値に対する信頼区間の信頼性の欠如に対処すること。
• 従来のThree-Cornered Hat法およびGCov法で一般的に見られる負の分散推定値の問題を解消すること。
• 推定値のみに依存するのではなく、元の測定値の全統計量を用いることで、頑健な信頼区間を導出するベイズ的手法の開発。
• GCov推定値が負になる場合でも有効な、正の中央値に基づく安定性推定値を提供すること。
• 低EDF、高測定ノイズ、クロック安定性の不均衡といった多様な条件下での手法の妥当性を検証すること。

提案手法

• KLTSは、2つのクロック間の比較から得られる十分統計量に基づくベイズ的逆問題を用い、推定分散に依存しない。
• すべての測定された周波数偏差の同時尤度を用いて、ノイズ分散の累積分布関数（CDF）を構築する。
• データの共分散構造を効率的に表現し、十分統計量を抽出するために、カルフネン＝ローベ変換を用いる。
• 信頼区間は、導出されたCDFの分位数から得られる：95％信頼区間は2.5百分位数と97.5百分位数から、中央値推定値は50百分位数から得られる。
• さまざまなEDF、測定ノイズレベル、クロック安定性の不均衡を想定した広範なモンテカルロシミュレーションにより、手法の妥当性を検証した。
• 高EDF（>100）の場合、以前に提案されたKLTG手法が信頼できる代替手段であることが示され、KLTSと併用可能な相補的関係にある。

実験結果

リサーチクエスチョン

• RQ1 十分統計量に基づくベイズ的手法が、特に1自由度を含む極めて低いEDFでも、クロック安定性推定値の信頼区間を信頼できるものにできるか？
• RQ2 提案手法が、Three-Cornered Hat法およびGCov法で一般的に見られる負の分散推定値の発生を排除できるか？
• RQ3 測定ノイズのレベルやクロック安定性の不均衡の違いに対して、KLTSの性能が理論的期待値およびモンテカルロシミュレーションとどの程度一致するか？
• RQ4 導出されたCDFの中央値が、GCov推定値が負であっても、物理的に意味のある、常に正の安定性推定値として機能するか？
• RQ5 測定ノイズが、特に短い積分時間における信頼区間の幅に及ぼす実用的影響は何か？

主な発見

• KLTSは、1自由度でも近似や簡略化を一切用いずに、アラン分散推定値の正確な95％信頼区間を提供する。
• 導出されたCDFの中央値は常に正であり、GCov推定値が負になる場合でも物理的に意味のある安定性推定値を提供する。
• モンテカルロシミュレーションにより、KLTSの信頼区間が、すべてのテストされたEDFレベル、測定ノイズ条件、クロック安定性の不均衡に対して理論的期待値と一致することが確認された。
• 積分時間が5000秒以上（EDF < 100）になると、KLTSの信頼区間が著しく広がり、安定性が低いクロックでは下限が0に近づくため、上界のみが特定可能になる。
• 1秒の積分時間では、測定ノイズが支配的となり、信頼区間が非常に広がる（例：CSO Cでは真の安定性の5倍以上に広がる）ため、高EDFであっても高ノイズが精度を制限することが示された。
• この手法は、GCov推定値が信頼できない（特に負値のとき）場合にそれを特定し、信頼区間による意味のある上界に置き換えることに成功した。