Comparison of Step Samplers for Nested Sampling

Abstract

Bayesian inference with nested sampling requires a likelihood-restricted prior sampling method, which draws samples from the prior distribution that exceed a likelihood threshold. For high-dimensional problems, Markov Chain Monte Carlo derivatives have been proposed. We numerically study ten algorithms based on slice sampling, hit-and-run and differential evolution algorithms in ellipsoidal, non-ellipsoidal and non-convex problems from 2 to 100 dimensions. Mixing capabilities are evaluated with the nested sampling shrinkage test. This makes our results valid independent of how heavy-tailed the posteriors are. Given the same number of steps, slice sampling is outperformed by hit-and-run and whitened slice sampling, while whitened hit-and-run does not provide as good results. Proposing along differential vectors of live point pairs also leads to the highest efficiencies, and appears promising for multi-modal problems. The tested proposals are implemented in the UltraNest nested sampling package, enabling efficient low and high-dimensional inference of a large class of practical inference problems relevant to astronomy, cosmology, particle physics and astronomy.

accepted MaxEnt 2022 proceeding, published in Physical Sciences Forum. UltraNest nested sampling package https://johannesbuchner.github.io/UltraNest/

References in corpus (2)

• A statistical test for Nested Sampling algorithms
• Nested sampling for physical scientists

Cited by in corpus (4)

• Cross-Comparison of Sampling Algorithms for Pulse Profile Modeling of PSR J0740+6620
• A comparison of Bayesian sampling algorithms for high-dimensional particle physics and cosmology applications
• Einstein@Home all-sky "bucket" search for continuous gravitational waves in LIGO O3 public data
• Determination of neutron star radius from pulse profile modeling using profile likelihood