Geometric nested sampling: sampling from distributions defined on non-trivial geometries

Abstract

Metropolis Hastings nested sampling evolves a Markov chain, accepting new points along the chain according to a version of the Metropolis Hastings acceptance ratio, which has been modified to satisfy the nested sampling likelihood constraint. The geometric nested sampling algorithm I present here is based on the Metropolis Hastings method, but treats parameters as though they represent points on certain geometric objects, namely circles, tori and spheres. For parameters which represent points on a circle or torus, the trial distribution is "wrapped" around the domain of the posterior distribution such that samples cannot be rejected automatically when evaluating the Metropolis ratio due to being outside the sampling domain. Furthermore, this enhances the mobility of the sampler. For parameters which represent coordinates on the surface of a sphere, the algorithm transforms the parameters into a Cartesian coordinate system before sampling which again makes sure no samples are automatically rejected, and provides a physically intuitive way of the sampling the parameter space.

Peer reviewed and published in JOSS. arXiv admin note: substantial text overlap with arXiv:1905.09110

References in corpus (9)

• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• PolyChord: next-generation nested sampling
• A Digital Correlator Upgrade for the Arcminute MicroKelvin Imager
• Bayesian modelling of clusters of galaxies from multi-frequency pointed Sunyaev--Zel'dovich observations
• Sunyaev--Zel'dovich profile fitting with joint AMI-Planck analysis
• Physical modelling of galaxy clusters detected by Planck
• Comparison of physical and observational galaxy cluster modelling
• Physical modelling of galaxy clusters and Bayesian inference in astrophysics