SEARCHING FOR NARROW EMISSION LINES IN X-RAY SPECTRA: COMPUTATION AND METHODS

被引:21
作者
Park, Taeyoung [1 ]
van Dyk, David A. [2 ]
Siemiginowska, Aneta [3 ]
机构
[1] Univ Pittsburgh, Dept Stat, Pittsburgh, PA 15120 USA
[2] Univ Calif Irvine, Dept Stat, Irvine, CA 92697 USA
[3] Smithsonian Astrophys Observ, Cambridge, MA 02138 USA
关键词
methods: statistical; quasars: emission lines;
D O I
10.1086/591631
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
The detection and quantification of narrow emission lines in X-ray spectra is a challenging statistical task. The Poisson nature of the photon counts leads to local random fluctuations in the observed spectrum that often result in excess emission in a narrow band of energy resembling a weak narrow line. From a formal statistical perspective, this leads to a ( sometimes highly) multimodal likelihood. Many standard statistical procedures are based on ( asymptotic) Gaussian approximations to the likelihood and simply cannot be used in such settings. Bayesian methods offer a more direct paradigm for accounting for such complicated likelihood functions, but even here multimodal likelihoods pose significant computational challenges. The new Markov chain Monte Carlo (MCMC) methods developed in 2008 by van Dyk and Park, however, are able to fully explore the complex posterior distribution of the location of a narrow line, and thus provide valid statistical inference. Even with these computational tools, standard statistical quantities such as means and standard deviations cannot adequately summarize inference and standard testing procedures cannot be used to test for emission lines. In this paper, we use new efficient MCMC algorithms to fit the location of narrow emission lines, we develop new statistical strategies for summarizing highly multimodal distributions and quantifying valid statistical inference, and we extend the method of posterior predictive p-values proposed by Protassov and coworkers to test for the presence of narrow emission lines in X-ray spectra. We illustrate and validate our methods using simulation studies and apply them to the Chandra observations of the high-redshift quasar PG 1634+706.
引用
收藏
页码:807 / 825
页数:19
相关论文
共 51 条
[1]  
CASH W, 1979, ASTROPHYS J, V228, P939, DOI 10.1086/156922
[2]   Chandra detects relativistic broad absorption lines from APM 08279+5255 [J].
Chartas, G ;
Brandt, WN ;
Gallagher, SC ;
Garmire, GP .
ASTROPHYSICAL JOURNAL, 2002, 579 (01) :169-175
[3]  
CRITCHLEY F, 1988, BIOMETRIKA, V75, P21
[4]   MAXIMUM LIKELIHOOD FROM INCOMPLETE DATA VIA EM ALGORITHM [J].
DEMPSTER, AP ;
LAIRD, NM ;
RUBIN, DB .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-METHODOLOGICAL, 1977, 39 (01) :1-38
[5]  
Fabian AC, 2006, ESA SPEC PUBL, V604, P463
[6]   Sherpa: a mission-independent data analysis application [J].
Freeman, PE ;
Doe, S ;
Siemiginowska, A .
ASTRONOMICAL DATA ANALYSIS, 2001, 4477 :76-87
[7]   X-ray spectroscopy of quasi-stellar objects with broad ultraviolet absorption lines [J].
Gallagher, SC ;
Brandt, WN ;
Chartas, G ;
Garmire, GP .
ASTROPHYSICAL JOURNAL, 2002, 567 (01) :37-41
[8]  
Gelman A, 1996, STAT SINICA, V6, P733
[9]  
GELMAN A, 1996, MCMC PRACTICE
[10]  
Gelman A, 2003, BAYESIAN DATA ANAL