Mean shift cluster recognition method implementation in the nested sampling algorithm
arXiv:2002.01431 · doi:10.3390/e22020185
Abstract
Nested sampling is an efficient algorithm for the calculation of the Bayesian evidence and posterior parameter probability distributions. It is based on the step-by-step exploration of the parameter space by Monte Carlo sampling with a series of values sets called live points that evolve towards the region of interest, i.e. where the likelihood function is maximal. In presence of several local likelihood maxima, the algorithm converges with difficulty. Some systematic errors can also be introduced by unexplored parameter volume regions. In order to avoid this, different methods are proposed in the literature for an efficient search of new live points, even in presence of local maxima. Here we present a new solution based on the mean shift cluster recognition method implemented in a random walk search algorithm. The clustering recognition is integrated within the Bayesian analysis program NestedFit. It is tested with the analysis of some difficult cases. Compared to the analysis results without cluster recognition, the computation time is considerably reduced. At the same time, the entire parameter space is efficiently explored, which translates into a smaller uncertainty of the extracted value of the Bayesian evidence.
References in corpus (27)
- Cosmological parameters from CMB and other data: a Monte-Carlo approach
- MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
- dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
- Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
- GW190425: Observation of a Compact Binary Coalescence with Total Mass
- Bayes in the sky: Bayesian inference and model selection in cosmology
- PolyChord: next-generation nested sampling
- A Nested Sampling Algorithm for Cosmological Model Selection
- Bayesian coherent analysis of in-spiral gravitational wave signals with a detector network
- A statistical test for Nested Sampling algorithms
- DIAMONDS: a new Bayesian Nested Sampling tool. Application to Peak Bagging of solar-like oscillations
- Properties of Nested Sampling
- A super-Earth and two sub-Neptunes transiting the nearby and quiet M dwarf TOI-270
- Direct comparison between Bayesian and frequentist uncertainty quantification for nuclear reactions
- Bayesian Calibrated Significance Levels Applied to the Spectral Tilt and Hemispherical Asymmetry
- Resolution of Single Spin-Flips of a Single Proton
- Measurement of the charged pion mass using X-ray spectroscopy of exotic atoms
- Constant-pressure nested sampling with atomistic dynamics
- Observation of the 2p3/2 -> 2s1/2 intra-shell transition in He-like uranium
- Bayesian data analysis tools for atomic physics
- Incorporating Uncertainties in Atomic Data Into the Analysis of Solar and Stellar Observations: A Case Study in Fe XIII
- New Test of Modulated Electron Capture Decay of Hydrogen-Like Pm Ions: Precision Measurement of Purely Exponential Decay
- Bayesian inference of a negative quantity from positive measurement results
- The Nested_fit data analysis program
- Line shape analysis of the K transition in muonic hydrogen
- Measurement of the charged pion mass using a low-density target of light atoms
- Bayesian Analysis of Femtosecond Pump-Probe Photoelectron-Photoion Coincidence Spectra with Fluctuating Laser Intensities