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.

Estimators for the exponent and upper limit, and goodness-of-fit tests for (truncated) power-law distributions

A new tool to derive simultaneously exponent and extremes of power-law distributions

Segmental estimation and testing method for power-law distributions and some extensions

Self-consistent dynamical models with a finite extent -- III. Truncated power-law spheres

Power-Law Distributions in Empirical Data

A complete data frame work for fitting power law distributions

Is Least-Squares Inaccurate in Fitting Power-Law Distributions? the Criticism is Complete Nonsense

Discrete Truncated Power‐Law Distributions

Statistical Methods for Investigating the Cosmic Ray Energy Spectrum

Detection of two power-law tails in the probability distribution functions of massive GMCs

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

A Multiple Power Law Distribution for Initial Mass Functions

Near-Gaussian distributions for modelling discrete stellar velocity data with heteroskedastic uncertainties

Statistical Tools for Analyzing the Cosmic Ray Energy Spectrum

New probability distributions in astrophysics: II. The generalized and double truncated Lindley

Power Law Stellar Distributions

Lens galaxies in the Illustris simulation: power-law models and the bias of the Hubble constant from time-delays

The Effect of Limited Sample Sizes on the Accuracy of the Estimated Scaling Parameter for Power-Law-Distributed Solar Data

Spectroscopic Constraints on the Form of the Stellar Cluster Mass Function

Flat to nonflat: Calculating nonlinear power spectra of biased tracers for nonflat $Λ$CDM model

Studying the properties of galaxy cluster morphology estimators