Affiliation:
1. School of Statistics, Southwestern University of Finance and Economics, Chengdu 611130, China
Abstract
Sampling from constrained distributions has posed significant challenges in terms of algorithmic design and non-asymptotic analysis, which are frequently encountered in statistical and machine-learning models. In this study, we propose three sampling algorithms based on Langevin Monte Carlo with the Metropolis–Hastings steps to handle the distribution constrained within some convex body. We present a rigorous analysis of the corresponding Markov chains and derive non-asymptotic upper bounds on the convergence rates of these algorithms in total variation distance. Our results demonstrate that the sampling algorithm, enhanced with the Metropolis–Hastings steps, offers an effective solution for tackling some constrained sampling problems. The numerical experiments are conducted to compare our methods with several competing algorithms without the Metropolis–Hastings steps, and the results further support our theoretical findings.
Subject
General Physics and Astronomy
Reference43 articles.
1. Bayesian analysis of constrained parameter and truncated data problems using Gibbs sampling;Gelfand;J. Am. Stat. Assoc.,1992
2. Latent dirichlet allocation;Blei;J. Mach. Learn. Res.,2003
3. Klein, J.P., and Moeschberger, M.L. (2005). Survival Analysis: Techniques for Censored and Truncated Data, Springer.
4. Johnson, V.E., and Albert, J.H. (2006). Ordinal Data Modeling, Springer.
5. Regularization in regression: Comparing Bayesian and frequentist methods in a poorly informative situation;Celeux;Bayesian Anal.,2012