Stereographic Markov Chain Monte Carlo
- URL: http://arxiv.org/abs/2205.12112v2
- Date: Wed, 21 Feb 2024 04:46:20 GMT
- Title: Stereographic Markov Chain Monte Carlo
- Authors: Jun Yang, Krzysztof {\L}atuszy\'nski, Gareth O. Roberts
- Abstract summary: High-dimensional distributions are notoriously difficult for off-the-shelf MCMC samplers.
We introduce a new class of MCMC samplers that map the original high-dimensional problem in Euclidean space onto a sphere.
In the best scenario, the proposed samplers can enjoy the blessings of dimensionality'' that the convergence is faster in higher dimensions.
- Score: 2.9304381683255945
- License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/
- Abstract: High-dimensional distributions, especially those with heavy tails, are
notoriously difficult for off-the-shelf MCMC samplers: the combination of
unbounded state spaces, diminishing gradient information, and local moves
results in empirically observed ``stickiness'' and poor theoretical mixing
properties -- lack of geometric ergodicity. In this paper, we introduce a new
class of MCMC samplers that map the original high-dimensional problem in
Euclidean space onto a sphere and remedy these notorious mixing problems. In
particular, we develop random-walk Metropolis type algorithms as well as
versions of the Bouncy Particle Sampler that are uniformly ergodic for a large
class of light and heavy-tailed distributions and also empirically exhibit
rapid convergence in high dimensions. In the best scenario, the proposed
samplers can enjoy the ``blessings of dimensionality'' that the convergence is
faster in higher dimensions.
Related papers
- ACMamba: Fast Unsupervised Anomaly Detection via An Asymmetrical Consensus State Space Model [51.83639270669481]
Unsupervised anomaly detection in hyperspectral images (HSI) aims to detect unknown targets from backgrounds.
HSI studies are hindered by steep computational costs due to the high-dimensional property of HSI and dense sampling-based training paradigm.
We propose an Asymmetrical Consensus State Space Model (ACMamba) to significantly reduce computational costs without compromising accuracy.
arXiv Detail & Related papers (2025-04-16T05:33:42Z) - Multilevel Generative Samplers for Investigating Critical Phenomena [3.8160065878097797]
Long-range correlations cause critical slowing down in Markov chain Monte Carlo.
We propose a novel sampler specialized for near-critical systems.
We show that the effective sample size RiGCS is few orders of magnitude higher than state-of-the-art generative model baselines.
arXiv Detail & Related papers (2025-03-11T22:03:54Z) - Symmetry-driven embedding of networks in hyperbolic space [0.4779196219827508]
Hyperbolic models can reproduce the heavy-tailed degree distribution, high clustering, and hierarchical structure of empirical networks.
Current algorithms for finding the hyperbolic coordinates of networks, however, do not quantify uncertainty in the inferred coordinates.
We present BIGUE, a Markov chain Monte Carlo algorithm that samples the posterior distribution of a Bayesian hyperbolic random graph model.
arXiv Detail & Related papers (2024-06-15T18:44:02Z) - Dimension-free Relaxation Times of Informed MCMC Samplers on Discrete Spaces [5.075066314996696]
We develop general mixing time bounds for Metropolis-Hastings algorithms on discrete spaces.
We establish sufficient conditions for a class of informed Metropolis-Hastings algorithms to attain relaxation times independent of the problem dimension.
arXiv Detail & Related papers (2024-04-05T02:40:45Z) - On Cyclical MCMC Sampling [13.002470980542707]
We show that cyclical MCMC converges to the desired probability distribution in settings where the Markov kernels used are fast mixing.
We also show that cyclical MCMC estimates well the local shape of the target distribution around each mode, even when we do not have convergence to the target.
arXiv Detail & Related papers (2024-03-01T02:20:44Z) - Weighted Riesz Particles [0.0]
We consider the target distribution as a mapping where the infinite-dimensional space of the parameters consists of a number of deterministic submanifolds.
We study the properties of the point, called Riesz, and embed it into sequential MCMC.
We find that there will be higher acceptance rates with fewer evaluations.
arXiv Detail & Related papers (2023-12-01T14:36:46Z) - Scaling Riemannian Diffusion Models [68.52820280448991]
We show that our method enables us to scale to high dimensional tasks on nontrivial manifold.
We model QCD densities on $SU(n)$ lattices and contrastively learned embeddings on high dimensional hyperspheres.
arXiv Detail & Related papers (2023-10-30T21:27:53Z) - Chebyshev Particles [0.0]
We are first to consider the posterior distribution of the objective as a mapping of samples in an infinite-dimensional Euclidean space.
We propose a new criterion by maximizing the weighted Riesz polarization quantity, to discretize rectifiable submanifolds via pairwise interaction.
We have achieved high performance from the experiments for parameter inference in a linear state-space model with synthetic data and a non-linear volatility model with real-world data.
arXiv Detail & Related papers (2023-09-10T16:40:30Z) - SIGMA: Scale-Invariant Global Sparse Shape Matching [50.385414715675076]
We propose a novel mixed-integer programming (MIP) formulation for generating precise sparse correspondences for non-rigid shapes.
We show state-of-the-art results for sparse non-rigid matching on several challenging 3D datasets.
arXiv Detail & Related papers (2023-08-16T14:25:30Z) - Orthogonal Matrix Retrieval with Spatial Consensus for 3D Unknown-View
Tomography [58.60249163402822]
Unknown-view tomography (UVT) reconstructs a 3D density map from its 2D projections at unknown, random orientations.
The proposed OMR is more robust and performs significantly better than the previous state-of-the-art OMR approach.
arXiv Detail & Related papers (2022-07-06T21:40:59Z) - Simultaneous Transport Evolution for Minimax Equilibria on Measures [48.82838283786807]
Min-max optimization problems arise in several key machine learning setups, including adversarial learning and generative modeling.
In this work we focus instead in finding mixed equilibria, and consider the associated lifted problem in the space of probability measures.
By adding entropic regularization, our main result establishes global convergence towards the global equilibrium.
arXiv Detail & Related papers (2022-02-14T02:23:16Z) - A Note on Optimizing Distributions using Kernel Mean Embeddings [94.96262888797257]
Kernel mean embeddings represent probability measures by their infinite-dimensional mean embeddings in a reproducing kernel Hilbert space.
We show that when the kernel is characteristic, distributions with a kernel sum-of-squares density are dense.
We provide algorithms to optimize such distributions in the finite-sample setting.
arXiv Detail & Related papers (2021-06-18T08:33:45Z) - What Are Bayesian Neural Network Posteriors Really Like? [63.950151520585024]
We show that Hamiltonian Monte Carlo can achieve significant performance gains over standard and deep ensembles.
We also show that deep distributions are similarly close to HMC as standard SGLD, and closer than standard variational inference.
arXiv Detail & Related papers (2021-04-29T15:38:46Z)
This list is automatically generated from the titles and abstracts of the papers in this site.
This site does not guarantee the quality of this site (including all information) and is not responsible for any consequences.