Sampling from manifold-restricted distributions using tangent bundle projections

Abstract

A common problem in Bayesian inference is the sampling of target probability distributions at sufficient resolution and accuracy to estimate the probability density, and to compute credible regions. Often by construction, many target distributions can be expressed as some higher-dimensional closed-form distribution with parametrically constrained variables, i.e., one that is restricted to a smooth submanifold of Euclidean space. I propose a derivative-based importance sampling framework for such distributions. A base set of $n$ samples from the target distribution is used to map out the tangent bundle of the manifold, and to seed $nm$ additional points that are projected onto the tangent bundle and weighted appropriately. The method essentially acts as an upsampling complement to any standard algorithm. It is designed for the efficient production of approximate high-resolution histograms from manifold-restricted Gaussian distributions, and can provide large computational savings when sampling directly from the target distribution is expensive.

Published version; Python implementation available at https://github.com/alvincjk/sampling-manifold-restricted-gaussians

References in corpus (5)

• Observation of Gravitational Waves from a Binary Black Hole Merger
• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Properties of the Binary Black Hole Merger GW150914
• Reduced-order modeling with artificial neurons for gravitational-wave inference
• An investigation into the Multiple Optimised Parameter Estimation and Data compression algorithm