[논문 리뷰] Rare Event Statistics Applied to Fast Radio Bursts

Scott Vander Wiel, Sarah Burke-Spolaor외 3명|arXiv (Cornell University)|2016. 12. 02.

Pulsars and Gravitational Waves Research참고 문헌 2인용 수 9

한 줄 요약

이 논문은 희귀 사건 통계, 특히 불확실한 비율을 가진 포isson 카운트에 대한 음이이항 모형을 활용하여 빠른 라디오 불확실성(_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 발생률을 추정하고, 탐지 시간을 예측하며, 다양한 조사 구성 간 발생률을 비교할 수 있도록 접근하기 쉬우며 실용적인 통계 도구를 제공한다.
• 특히 조사 감도와 천구 영역 커버리지가 다를 경우, 알려지지 않은 또는 변동성이 있는 포isson 발생률로 인한 FRB 발생률 추정의 불확실성을 정량화한다.
• 은하수 평면 감도 영향을 고려한 후에도 고은하위 경도(>30°) 조사와 저은하위/중간은하위 경도 조사 간 FRB 발생률이 유의미하게 다름을 평가한다.
• 이질적인 FRB 조사를 메타분석하여 통합된, 신뢰구간 기반의 FRB 발생률 함수 추정치를 도출할 수 있도록 한다.

제안 방법

• 포isson 발생률 매개변수의 불확실성 또는 조사 간 변동성으로 인해 과분산이 발생할 수 있는 상황에서 카운트 데이터를 모델링하기 위해 음이이항 분포를 사용한다.
• 중립적인 감마 prior(α=1/3, β=0)를 적용하여 베이지안 추론을 수행하여 FRB 발생률에 대한 사후 분포를 추정하고, 탐지 시간에 대한 직접적인 확률적 진술이 가능하게 한다.
• F-분포 분위수를 사용하여 흔한 방법의 예측 구간을 계산하여 미래의 카운트에 대해 보수적인 95% 예측 구간을 산출하며, 특히 최소한 한 개의 펄스를 탐지하는 데 초점을 맞춘다.
• 다양한 조사가 동일한 기초 발생률을 공유하는지 평가하기 위해 우도 비율 검정을 적용하며, 고은하위 및 저은하위 샘플 간 일관성을 검증한다.
• 음이이항 회귀를 사용하여 FRB 발생률을 감도, 은하위 경도, 그리고 초포isson 변동성의 함수로 모델링하고, 소스 수 계수를 추정한다.
• 탐지 시간 추정에 대해 베이지안 및 흔한 방법을 비교하며, 이해 가능성과 더 나은 校정성 덕분에 베이지안 방법을 우선시한다.

실험 결과

연구 질문

• RQ1 기존 탐지 카운트와 발생률에 대한 불확실성을 고려할 때, 파크스 다중빔 수신기로 FRB를 탐지할 확률이 50% 또는 95%가 되기 위해 필요한 관측 시간은 얼마나 되는가?
• RQ2 은하수 평면 감도 저하 영향을 고려한 후에도 고은하위 경도(>30°)에서의 FRB 탐지 발생률이 저은하위 경도(<5°)보다 유의미하게 높은가?
• RQ3 더 긴 개별 포인팅 시간을 가진 조사(‘보내기 길게’)가 넓은 영역을 스캔하는 조사(‘넓게 스캔하기’)보다 실제 탐지 발생률이 더 높은가, 가용 탐지 데이터를 바탕으로 판단할 수 있는가?
• RQ4 다양한 조사로 구성된 집합에서 공통된 기초 발생률를 가정할 수 있는가, 아니면 유의미한 변동성이 존재하여 더 민첩한 모형이 필요한가?
• RQ5 합쳐진 조사 데이터를 기반으로, 감도에 따른 FRB 발생률 추정치는 어떻게 되며, 은하위 경도에 따라 어떻게 변화하는가?

주요 결과

• 13개의 FRB가 5개의 고은하위 경도 조사(|b| > 30°)에서 탐지되어, 은하수 평면 감도 저하 보정 후에도 1개의 FRB가 4개의 저은하위 경도 조사(|b| < 5°)에서 탐지된 것과 비교해 강력한 증거(p = 5 × 10⁻⁵)를 제공한다.
• 12개의 이질적인 조사를 메타분석한 결과, 고은하위 경도에서의 FRB 발생률은 2,866 sky⁻¹ day⁻¹(95% CI: 1,121 ~ 7,328)로 추정되었고, 저은하위/중간은하위 경도에서는 285 sky⁻¹ day⁻¹(95% CI: 48 ~ 1,701)로 추정되었다.
• 파크스 HTRU 조사의 경우, 50% 확률로 FRB를 탐지하기 위해 193시간의 관측이 필요하고, 95% 확률을 확보하기 위해서는 937시간이 필요하다. 이는 10회의 탐지와 발생률 불확실성 모델링을 기반으로 한다.
• 작은 고은하위 경도 영역에서의 두 가지 빠른 탐지(Ravi et al., 2015; Petroff et al., 2015)는 ‘보내기 길게’ 전략이 ‘넓게 스캔하기’ HTRU보다 더 나은 성능을 낼 수 있음을 시사하지만, 이는 오직 약간의 유의미성(p = 0.07)에 그친다.
• 우도 비율 검정 결과, 고은하위 경도와 통합된 저은하위/중간은하위 경도 조사 간 공통 발생률가 타당한 것으로 나타났지만, 유클리드 스케일링 가정 하에 전체 집합에서는 그렇지 않다.
• 음이이항 회귀 분석은 다른 논문에서 보고된 결과와 일치하는 발생률 함수를 추정하였으며, 감도 및 은하위 경도 의존성을 반영한 소스 수 계수를 추정하면서 초포isson 변동성도 고려했다.