Mean shift cluster recognition method implementation in the nested sampling algorithm
This work addresses convergence issues in Bayesian analysis for researchers using nested sampling, offering an incremental improvement in efficiency and accuracy.
The authors tackled the problem of nested sampling struggling with multiple local likelihood maxima by implementing a mean shift cluster recognition method within a random walk search algorithm, which reduced computation time and improved the uncertainty of Bayesian evidence estimates.
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.