169 research outputs found
Randomized Hamiltonian Monte Carlo as Scaling Limit of the Bouncy Particle Sampler and Dimension-Free Convergence Rates
The Bouncy Particle Sampler is a Markov chain Monte Carlo method based on a
nonreversible piecewise deterministic Markov process. In this scheme, a
particle explores the state space of interest by evolving according to a linear
dynamics which is altered by bouncing on the hyperplane tangent to the gradient
of the negative log-target density at the arrival times of an inhomogeneous
Poisson Process (PP) and by randomly perturbing its velocity at the arrival
times of an homogeneous PP. Under regularity conditions, we show here that the
process corresponding to the first component of the particle and its
corresponding velocity converges weakly towards a Randomized Hamiltonian Monte
Carlo (RHMC) process as the dimension of the ambient space goes to infinity.
RHMC is another piecewise deterministic non-reversible Markov process where a
Hamiltonian dynamics is altered at the arrival times of a homogeneous PP by
randomly perturbing the momentum component. We then establish dimension-free
convergence rates for RHMC for strongly log-concave targets with bounded
Hessians using coupling ideas and hypocoercivity techniques.Comment: 47 pages, 2 figure
Langevin and Hamiltonian based Sequential MCMC for Efficient Bayesian Filtering in High-dimensional Spaces
Nonlinear non-Gaussian state-space models arise in numerous applications in
statistics and signal processing. In this context, one of the most successful
and popular approximation techniques is the Sequential Monte Carlo (SMC)
algorithm, also known as particle filtering. Nevertheless, this method tends to
be inefficient when applied to high dimensional problems. In this paper, we
focus on another class of sequential inference methods, namely the Sequential
Markov Chain Monte Carlo (SMCMC) techniques, which represent a promising
alternative to SMC methods. After providing a unifying framework for the class
of SMCMC approaches, we propose novel efficient strategies based on the
principle of Langevin diffusion and Hamiltonian dynamics in order to cope with
the increasing number of high-dimensional applications. Simulation results show
that the proposed algorithms achieve significantly better performance compared
to existing algorithms
- …