Abstract:Error bounds are derived for sampling and estimation using a discretization of an intrinsically defined Langevin diffusion with invariant measure $\text{d}\mu_\phi \propto e^{-\phi} \mathrm{dvol}_g $ on a compact Riemannian manifold. Two estimators of linear functionals of $\mu_\phi $ based on the discretized Markov process are considered: a time-averaging estimator based on a single trajectory and an ensemble-averaging estimator based on multiple independent trajectories. Imposing no restrictions beyond a nominal level of smoothness on $\phi$, first-order error bounds, in discretization step size, on the bias and variance/mean-square error of both estimators are derived. The order of error matches the optimal rate in Euclidean and flat spaces, and leads to a first-order bound on distance between the invariant measure $\mu_\phi$ and a stationary measure of the discretized Markov process. This order is preserved even upon using retractions when exponential maps are unavailable in closed form, thus enhancing practicality of the proposed algorithms. Generality of the proof techniques, which exploit links between two partial differential equations and the semigroup of operators corresponding to the Langevin diffusion, renders them amenable for the study of a more general class of sampling algorithms related to the Langevin diffusion. Conditions for extending analysis to the case of non-compact manifolds are discussed. Numerical illustrations with distributions, log-concave and otherwise, on the manifolds of positive and negative curvature elucidate on the derived bounds and demonstrate practical utility of the sampling algorithm.

Estimating sample paths of Gauss-Markov processes from noisy data

Pathwise Derivatives for Multivariate Distributions

Path properties of a d-dimensional Gaussian process

A variational approach to sampling in diffusion processes

Increments and Sample Path Properties of Gaussian Processes

Conditioning of Gaussian processes and a zero area Brownian bridge

Gaussian Approximation and Moderate Deviations of Poisson Shot Noises with Application to Compound Generalized Hawkes Processes

Deviation frequencies of Brownian path property approximations

Sample Path Regularity of Gaussian Processes from the Covariance Kernel

Sample path properties of reflected Gaussian processes

Estimating the Mixing Coefficients of Geometrically Ergodic Markov Processes

Sampling and estimation on manifolds using the Langevin diffusion

Estimating the Hawkes process from a discretely observed sample path

The saddlepoint approximation factors over sample paths of recursively compounded processes

On the sample path properties of mixed Poisson processes.

Gaussian Processes with Noisy Regression Inputs for Dynamical Systems

Sequential sampling of Gaussian process latent variable models

Conditioning of Banach Space Valued Gaussian Random Variables: An Approximation Approach Based on Martingales

Generalised shot-noise representations of stochastic systems driven by non-Gaussian Lévy processes

Generalised shot noise representations of stochastic systems driven by non-Gaussian Lévy processes

Estimation Method and Its Properties of Moments of First Passage Time for Markov Chains in Reverse Problems