[論文レビュー] Polarization measurements analysis II. Best estimators of polarization fraction and angle

L. Montier, S. Plaszczynski+6 more|arXiv (Cornell University)|Jul 1, 2014

Stellar, planetary, and galactic studies参考文献 19被引用数 28

ひとこと要約

本稿は、ノイズが高く相関するストークスパラメータを伴う状況下で、頻度的推定、ベイズ推定、漸近的手法による偏光度および偏光角の推定を評価・比較する。漸近的手法はバイアスを最小化し、最も正規分布に近い出力分布を生成する一方、低S/N比におけるベイズ推定は強い非対称性を示す。最適化されたベイズ事前分布を用いることで、特に非対称なノイズ共分散を伴う低S/N比領域でのバイアスがさらに低減される。

ABSTRACT

With the forthcoming release of high precision polarization measurements, such as from the Planck satellite, it becomes critical to evaluate the performance of estimators for the polarization fraction and angle. These two physical quantities suffer from a well-known bias in the presence of measurement noise, as has been described in part I of this series. In this paper, part II of the series, we explore the extent to which various estimators may correct the bias. Traditional frequentist estimators of the polarization fraction are compared with two recent estimators: one inspired by a Bayesian analysis and a second following an asymptotic method. We investigate the sensitivity of these estimators to the asymmetry of the covariance matrix which may vary over large datasets. We present for the first time a comparison among polarization angle estimators, and evaluate the statistical bias on the angle that appears when the covariance matrix exhibits effective ellipticity. We also address the question of the accuracy of the polarization fraction and angle uncertainty estimators. The methods linked to the credible intervals and to the variance estimates are tested against the robust confidence interval method. From this pool of estimators, we build recipes adapted to different use-cases: build a mask, compute large maps, and deal with low S/N data. More generally, we show that the traditional estimators suffer from discontinuous distributions at low S/N, while the asymptotic and Bayesian methods do not. Attention is given to the shape of the output distribution of the estimators, and is compared with a Gaussian. In this regard, the new asymptotic method presents the best performance, while the Bayesian output distribution is shown to be strongly asymmetric with a sharp cut at low S/N.Finally, we present an optimization of the estimator derived from the Bayesian analysis using adapted priors.

研究の動機と目的

• 相関的かつ非対称なQ-Uノイズを含む現実的なノイズ条件下での偏光度および偏光角推定器の性能を評価すること。
• 特に低信号対ノイズ比（S/N）におけるノイズによって引き起こされる長年のバイアス問題に取り組むこと。
• 従来の頻度的推定器と最近のベイズおよび漸近的手法を、バイアス、不確実性推定、出力分布の形状の観点から比較すること。
• マスク作成、大規模マップ生成、低S/N比データ解析といった用途に最適化された、ユースケース固有のレシピを提供すること。

提案手法

• 相関的かつ非対称なQおよびUストークスパラメータのノイズを考慮した、尤度に基づくベイズフレームワークを用いて、偏光度および偏光角の最適推定器を導出する。
• 修正ベッセル関数および逐次超幾何関数を用いて、ベイズの最後尤度推定（MB）推定器およびその分散の解析的近似を導出する。
• 漸近的手法を用いて偏光度および偏光角を推定し、従来の推定器と比較して滑らかでより正規分布に近い出力分布を生成することを示す。
• 信頼性を評価するために、信用区間、分散推定値、および頑健な信頼区間を用いて不確実性推定器をテスト・比較する。
• バイアスを低減するために、事前分布を調整することで最適化されたベイズ推定器を導入し、特に低S/N比および高楕円度ノイズ領域で有効である。
• 共分散行列における有効な楕円度が、特に偏光角および偏光度の推定器性能に与える影響を評価する。

実験結果

リサーチクエスチョン

• RQ1 従来の頻度的推定器と、最近のベイズおよび漸近的手法との間で、バイアスおよび不確実性推定の観点から、偏光度および偏光角の推定性能にどのような差があるか？
• RQ2 QおよびUにおける非対称的かつ相関的なノイズが、偏光度および偏光角推定器のバイアスに与える影響は何か？
• RQ3 偏光度および偏光角推定器の出力分布は正規分布からどれほどずれているか？ どの手法が最も対称的かつ良好に制御された分布を生成するか？
• RQ4 ベイズ推定器の解析的近似は、計算速度の向上に信頼性を持って利用できるか？ どのような条件下で精度が保たれるか？
• RQ5 マスク作成、大規模マップ生成、低S/N比データ解析といった異なる天体物理学的ユースケースにおける最適推定器は何か？

主な発見

• 漸近的手法は、低S/N比においても、偏光度の出力分布が最も正規分布に近く、従来の推定器が低S/N比で不連続性を示すのとは対照的に顕著に優れている。
• 従来の頻度的推定器は、特にノイズ共分散が非対称な場合、低S/N比において深刻なバイアスと不連続な分布を示す。
• ベイズ推定器は出力分布に強い非対称性を示し、低S/N比で明確な下限カットオフが見られるが、適応した事前分布を用いることでバイアスを低減できる。
• ベイズ推定器の解析的近似（式55および56）は、S/N < 1で相対誤差が0.15%未満、S/N > 4では0.05%未満であり、中間S/N比の低楕円度領域でも許容誤差（最大8%）以内に収まる。
• ベイズ不確実性推定器は、S/N < 1で正確な値と0.1%未満の一致を示し、低楕円度ケースではS/N > 1でも誤差が8%未満に抑えられ、実用的用途に適している。
• 有効な楕円度ε_eff > 1.1の領域では、ベイズ推定器の解析的近似の誤差が増加し（極端な場合20%に達する）、高精度応用ではノイズの非対称性を慎重にモデル化する必要がある。