Cargando…
Mean Shift Cluster Recognition Method Implementation in the Nested Sampling Algorithm
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...
Autores principales: | , |
---|---|
Formato: | Online Artículo Texto |
Lenguaje: | English |
Publicado: |
MDPI
2020
|
Materias: | |
Acceso en línea: | https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7516612/ https://www.ncbi.nlm.nih.gov/pubmed/33285961 http://dx.doi.org/10.3390/e22020185 |
_version_ | 1783587041384071168 |
---|---|
author | Trassinelli, Martino Ciccodicola, Pierre |
author_facet | Trassinelli, Martino Ciccodicola, Pierre |
author_sort | Trassinelli, Martino |
collection | PubMed |
description | 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. |
format | Online Article Text |
id | pubmed-7516612 |
institution | National Center for Biotechnology Information |
language | English |
publishDate | 2020 |
publisher | MDPI |
record_format | MEDLINE/PubMed |
spelling | pubmed-75166122020-11-09 Mean Shift Cluster Recognition Method Implementation in the Nested Sampling Algorithm Trassinelli, Martino Ciccodicola, Pierre Entropy (Basel) Article 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. MDPI 2020-02-06 /pmc/articles/PMC7516612/ /pubmed/33285961 http://dx.doi.org/10.3390/e22020185 Text en © 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). |
spellingShingle | Article Trassinelli, Martino Ciccodicola, Pierre Mean Shift Cluster Recognition Method Implementation in the Nested Sampling Algorithm |
title | Mean Shift Cluster Recognition Method Implementation in the Nested Sampling Algorithm |
title_full | Mean Shift Cluster Recognition Method Implementation in the Nested Sampling Algorithm |
title_fullStr | Mean Shift Cluster Recognition Method Implementation in the Nested Sampling Algorithm |
title_full_unstemmed | Mean Shift Cluster Recognition Method Implementation in the Nested Sampling Algorithm |
title_short | Mean Shift Cluster Recognition Method Implementation in the Nested Sampling Algorithm |
title_sort | mean shift cluster recognition method implementation in the nested sampling algorithm |
topic | Article |
url | https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7516612/ https://www.ncbi.nlm.nih.gov/pubmed/33285961 http://dx.doi.org/10.3390/e22020185 |
work_keys_str_mv | AT trassinellimartino meanshiftclusterrecognitionmethodimplementationinthenestedsamplingalgorithm AT ciccodicolapierre meanshiftclusterrecognitionmethodimplementationinthenestedsamplingalgorithm |