Statistics – Computation
Scientific paper
Jul 2011
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2011aspc..442..439s&link_type=abstract
Astronomical Data Analysis Software and Systems XX. ASP Conference Proceedings, Vol. 442, proceedings of a Conference held at Se
Statistics
Computation
Scientific paper
Typical X-ray spectra have low counts and should be modeled using the Poisson distribution. However, χ2 statistic is often applied as an alternative and the data are assumed to follow the Gaussian distribution. A variety of weights to the statistic or a binning of the data is performed to overcome the low counts issues. However, such modifications introduce biases or/and a loss of information. Standard modeling packages such as XSPEC and Sherpa provide the Poisson likelihood and allow computation of rudimentary MCMC chains, but so far do not allow for setting a full Bayesian model. We have implemented a sophisticated Bayesian MCMC-based algorithm to carry out spectral fitting of low counts sources in the Sherpa environment. The code is a Python extension to Sherpa and allows to fit a predefined Sherpa model to high-energy X-ray spectral data and other generic data. We present the algorithm and discuss several issues related to the implementation, including flexible definition of priors and allowing for variations in the calibration information.
Connors Alanna
Kashyap Vinay
Park Tuson
Refsdal Brian
Siemiginowska Aneta
No associations
LandOfFree
pyblocxs: Bayesian Low-Counts X-ray Spectral Analysis in Sherpa does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with pyblocxs: Bayesian Low-Counts X-ray Spectral Analysis in Sherpa, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and pyblocxs: Bayesian Low-Counts X-ray Spectral Analysis in Sherpa will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-989360