Searching for narrow emission lines in X-ray spectra: Computation and methods

Taeyoung Park, David A. Van Dyk, Aneta Siemiginowska

Research output: Contribution to journalArticle

15 Citations (Scopus)

Abstract

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.

Original languageEnglish
Pages (from-to)807-825
Number of pages19
JournalAstrophysical Journal
Volume688
Issue number2
DOIs
Publication statusPublished - 2008 Dec 1

Fingerprint

inference
Markov chains
Markov chain
x rays
quasars
Monte Carlo method
narrowband
standard deviation
method
deviation
photons
approximation
simulation
energy
distribution
test

All Science Journal Classification (ASJC) codes

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

Park, Taeyoung ; Van Dyk, David A. ; Siemiginowska, Aneta. / Searching for narrow emission lines in X-ray spectra : Computation and methods. In: Astrophysical Journal. 2008 ; Vol. 688, No. 2. pp. 807-825.
@article{51ce7b0e0a8341d5804f94f949d9f155,
title = "Searching for narrow emission lines in X-ray spectra: Computation and methods",
abstract = "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.",
author = "Taeyoung Park and {Van Dyk}, {David A.} and Aneta Siemiginowska",
year = "2008",
month = "12",
day = "1",
doi = "10.1086/591631",
language = "English",
volume = "688",
pages = "807--825",
journal = "Astrophysical Journal",
issn = "0004-637X",
publisher = "IOP Publishing Ltd.",
number = "2",

}

Searching for narrow emission lines in X-ray spectra : Computation and methods. / Park, Taeyoung; Van Dyk, David A.; Siemiginowska, Aneta.

In: Astrophysical Journal, Vol. 688, No. 2, 01.12.2008, p. 807-825.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Searching for narrow emission lines in X-ray spectra

T2 - Computation and methods

AU - Park, Taeyoung

AU - Van Dyk, David A.

AU - Siemiginowska, Aneta

PY - 2008/12/1

Y1 - 2008/12/1

N2 - 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.

AB - 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.

UR - http://www.scopus.com/inward/record.url?scp=56849125148&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=56849125148&partnerID=8YFLogxK

U2 - 10.1086/591631

DO - 10.1086/591631

M3 - Article

AN - SCOPUS:56849125148

VL - 688

SP - 807

EP - 825

JO - Astrophysical Journal

JF - Astrophysical Journal

SN - 0004-637X

IS - 2

ER -