[論文レビュー] Rare Event Statistics Applied to Fast Radio Bursts

Scott Vander Wiel, Sarah Burke-Spolaor+3 more|arXiv (Cornell University)|Dec 2, 2016

Pulsars and Gravitational Waves Research参考文献 2被引用数 9

ひとこと要約

本稿は、確率的レートが不確実なポアソン度数に対して、特に負の二項分布を用いたレアイベント統計を、高速電波バースト（FRB）検出調査に適用する。ベイズ的手法と頻度主義的手法を用いて、FRB発生レートの推定、検出確率の評価、複数の調査間でのレート比較を実施し、高銀緯度におけるFRBレートが顕著に高いことの強力な証拠を示し、高い信頼性で検出を実現するための必要な観測時間を定量的に評価している。

ABSTRACT

Statistical interpretation of sparsely sampled event rates has become vital for new transient surveys, particularly those aimed at detecting fast radio bursts (FRBs). We provide an accessible reference for a number of simple, but critical, statistical questions relevant for current transient and FRB research and utilizing the negative binomial model for counts in which the count rate parameter is uncertain or randomly biased from one study to the next. We apply these methods to re-assess and update results from previous FRB surveys, finding as follows. 1) Thirteen FRBs detected across five high-Galactic-latitude (> 30$^\circ$) surveys are highly significant $(p = 5 imes 10^{-5})$ evidence of a higher rate relative to the single FRB detected across four low-latitude (< 5$^\circ$) surveys, even after accounting for effects that dampen Galactic plane sensitivity. High- vs. mid-latitude (5 to 15$^\circ$) is marginally significant $(p = 0.03)$. 2) A meta analysis of twelve heterogeneous surveys gives an FRB rate of 2866 sky$^{-1}$ day$^{-1}$ above 1 Jy at high Galactic latitude (95% confidence 1121 to 7328) and 285 sky$^{-1}$ day$^{-1}$ at low/mid latitudes (95% from 48 to 1701). 3) Using the Parkes HTRU high-latitude setup requires 193 observing hours to achieve 50% probability of detecting an FRB and 937 hours to achieve 95% probability, based on the ten detections of (Champion et al. 2016) and appropriately accounting for uncertainty in the unknown Poisson rate. 4) Two quick detections at Parkes from a small number of high-latitude fields (Ravi et al. 2015; Petroff et al. 2015) tentatively favor a look long survey style relative to the scan wide HTRU survey, but only at $p = 0.07$ significance.

研究の動機と目的

• 複数の調査において度数が少なく、あるいは0件の検出にとどまるような、希少事象のレート推定における統計的課題に取り組む。
• 研究者がFRBレートを推定し、検出までの時間を予測し、異なる調査設定間でレートを比較できる、アクセス可能で実用的な統計ツールを提供する。
• 未知または変動するポアソンレートに起因するFRBレート推定の不確実性を定量的に評価する。特に、調査感度や天の川銀河平面周辺のスカイカバレッジに差がある場合に有効である。
• 天の川銀河平面の感度低下効果を考慮した上で、高銀緯度（>30°）と低・中銀緯度（<5°）の調査間でFRB発生レートが有意に異なるかどうかを評価する。
• 異種のFRB調査をメタアナリシスにより統合し、信頼区間を基準としたFRBレート関数の統合的推定を可能にする。

提案手法

• ポアソンレートパラメータが不確実または調査間で変動する場合に、過分散を許容できるように負の二項分布を度数データのモデルとして用いる。
• 中立的なガンマ事前分布（α=1/3, β=0）を用いたベイズ推論により、FRBレートの事後分布を推定し、検出までの時間に関する直接的な確率的記述を可能にする。
• F分布の分位点を用いた頻度主義的手法により、将来の度数の保守的95％予測区間を計算し、特に少なくとも1件のパルスを検出する場合に有効である。
• 複数の調査が共通の潜在的レートを持つかどうかを検証するため、尤度比検定を用い、高銀緯度および低銀緯度サンプル間の整合性を評価する。
• 負の二項回帰を用いて、FRBレートをフラックス感度、銀河平面からの距離、および余分なポアソン分散の関数としてモデル化し、源数インデックスを推定する。
• 検出時間推定におけるベイズ的手法と頻度主義的手法を比較し、解釈性とより良いキャリブレーションの観点からベイズ手法を優先する。

実験結果

リサーチクエスチョン

• RQ1 既存の検出度数とレートの不確実性を踏まえて、パーキーズ多ビーム受信機を用いた観測で、50％または95％の確率でFRBを検出するには、何時間の観測が必要か？
• RQ2 天の川銀河平面の感度低下を補正した上で、高銀緯度（>30°）におけるFRB検出レートが、低銀緯度調査（<5°）と比べ顕著に高いとされるか？
• RQ3 個々の点向の長時間観測（「長時間観測」）による調査と、広範囲をカバーする広域スキャン（「広域スキャン」）による調査とでは、利用可能な検出データに基づき、FRBの検出レートが顕著に異なるか？
• RQ4 複数の異種調査の集団において、共通の潜在的FRBレートを仮定できるか、それとも顕著な変動があり、より柔軟なモデルが必要か？
• RQ5 統合された調査データに基づき、フラックス感度と銀河平面からの距離に応じたFRBレートの推定値は何か？また、その変動はどのように変化するか？

主な発見

• 高銀緯度調査（|b| > 30°）5件で13件のFRBが検出されたが、低銀緯度調査（|b| < 5°）4件では1件の検出にとどまり、天の川銀河平面の感度低下を補正した後でも、FRBレートが顕著に高いことの強い証拠（p = 5 × 10⁻⁵）が得られた。
• 12件の異種調査のメタアナリシスにより、高銀緯度では2,866 sky⁻¹ day⁻¹（95％信頼区間：1,121～7,328）、低・中銀緯度では285 sky⁻¹ day⁻¹（95％信頼区間：48～1,701）のFRBレートが推定された。
• パーキーズHTRU調査では、50％の確率でFRBを検出するには193時間の観測が必要であり、95％の確率では937時間の観測が必要である。これは10件の検出とレートの不確実性モデルに基づく。
• 小さな高銀緯度領域からの2件の迅速検出（Ravi et al., 2015; Petroff et al., 2015）は、「長時間観測」戦略が「広域スキャン」HTRUに勝る可能性を示唆しているが、これは僅かに有意（p = 0.07）にとどまる。
• 尤度比検定の結果、高銀緯度調査と統合された低・中銀緯度調査の間では共通のレートが妥当であるが、ユークリッドスケーリング仮定下では全調査セットに対しては共通レートは不適切である。
• 負の二項回帰により、他の既存の研究結果と整合するレート関数が得られ、フラックス感度および銀河平面からの距離依存性を反映した源数インデックスが推定され、余分なポアソン分散を考慮した。