Open Access
September 2014 Automatic estimation of flux distributions of astrophysical source populations
Raymond K. W. Wong, Paul Baines, Alexander Aue, Thomas C. M. Lee, Vinay L. Kashyap
Ann. Appl. Stat. 8(3): 1690-1712 (September 2014). DOI: 10.1214/14-AOAS750


In astrophysics a common goal is to infer the flux distribution of populations of scientifically interesting objects such as pulsars or supernovae. In practice, inference for the flux distribution is often conducted using the cumulative distribution of the number of sources detected at a given sensitivity. The resulting “$\log(N>S)$–$\log(S)$” relationship can be used to compare and evaluate theoretical models for source populations and their evolution. Under restrictive assumptions the relationship should be linear. In practice, however, when simple theoretical models fail, it is common for astrophysicists to use prespecified piecewise linear models. This paper proposes a methodology for estimating both the number and locations of “breakpoints” in astrophysical source populations that extends beyond existing work in this field.

An important component of the proposed methodology is a new interwoven EM algorithm that computes parameter estimates. It is shown that in simple settings such estimates are asymptotically consistent despite the complex nature of the parameter space. Through simulation studies it is demonstrated that the proposed methodology is capable of accurately detecting structural breaks in a variety of parameter configurations. This paper concludes with an application of our methodology to the Chandra Deep Field North (CDFN) data set.


Download Citation

Raymond K. W. Wong. Paul Baines. Alexander Aue. Thomas C. M. Lee. Vinay L. Kashyap. "Automatic estimation of flux distributions of astrophysical source populations." Ann. Appl. Stat. 8 (3) 1690 - 1712, September 2014.


Published: September 2014
First available in Project Euclid: 23 October 2014

zbMATH: 1304.85001
MathSciNet: MR3271349
Digital Object Identifier: 10.1214/14-AOAS750

Keywords: $\log N$–$\log S$ , Broken power law , CDFN X-ray survey , interwoven EM algorithm , likelihood computations , Pareto distribution

Rights: Copyright © 2014 Institute of Mathematical Statistics

Vol.8 • No. 3 • September 2014
Back to Top