Abstract:Traditional binned statistics such as $\chi^2$ suffer from information loss and arbitrariness of the binning procedure. We point out that the underlying statistical quantity (the log likelihood $L$) does not require any binning beyond the one implied by instrumental readout channels, and we propose to use it for low-count data. The performance of $L$ in the model classification and point estimation problems is explored by Monte-Carlo simulations of Chandra and XMM X-ray spectra, and is compared to the performances of the binned Poisson statistic ($C$), Pearson's $\chi^2$ and Neyman's $\chi^2_N$, the Kolmogorov- Smirnov, and Kuiper' statistics. It is found that the unbinned log likelihood $L$ performs best with regard to the expected chi-square distance between true and estimated spectra, the chance of a successful identification among discrete candidate models, the area under the receiver-operator curve of reduced (two-model) binary classification problems, and generally also with regard to the mean square errors of individual spectrum parameters. The $\chi^2$ ($\chi^2_{\rm N}$) statistics should only be used if more than 10 (15) predicted counts per bin are available. From the practical point of view, the computational cost of evaluating $L$ is smaller than for any of the alternative methods if the forward model is specified in terms of a Poisson intensity and normalization is a free parameter. The maximum-$L$ method is applied to 14 observations from the Taurus Molecular Cloud, and the unbinned results are compared to binned XSPEC results, and found to generally agree, with exceptions explained by instability under re-binning and by background fine structures. The maximum-$L$ method has no lower limit on the available counts, and allows to treat weak sources which are beyond the means of binned methods.

LOO-PIT: A sensitive posterior test

Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC

A new set of tools for goodness-of-fit validation

Matching Bayesian and frequentist coverage probabilities when using an approximate data covariance matrix

On the Application of Bayesian Leave-One-Out Cross-Validation to Exoplanet Atmospheric Analysis

Are you using test log-likelihood correctly?

Informed Total-Error-Minimizing Priors: Interpretable cosmological parameter constraints despite complex nuisance effects

Statistical significance testing for mixed priors: a combined Bayesian and frequentist analysis

New statistical goodness of fit techniques in noisy inhomogeneous inverse problems - With application to the recovering of the luminosity distribution of the Milky Way

A New Goodness-Of-Fit Test Based On The Laplace Statistic For A Large Class Of Nhpp Models

Unbinned maximum-likelihood estimators for low-count data: Applications to faint X-ray spectra in the Taurus Molecular Cloud

PINT: Maximum-likelihood estimation of pulsar timing noise parameters

Asymptotics for Outlier Hypothesis Testing

On the choice of the splitting ratio for the split likelihood ratio test

Improving Pearson's chi-squared test: hypothesis testing of distributions -- optimally

Integrative conformal p -values for out-of-distribution testing with labelled outliers

On Hypothesis Testing for Comparing Image Quality Assessment Metrics [tips & Tricks]

POPE: post optimization posterior evaluation of likelihood free models

Calibration diagnostics for point process models via the probability integral transform

Study on Outliers in the Big Stellar Spectral Dataset of the Fifth Data Release (DR5) of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST)

Testing the approximations described in "Asymptotic formulae for likelihood-based tests of new physics"