Abstract
We study the problem of sampling from a probability distribution π on Rd which has a density w.r.t. the Lebesgue measure known up to a normalization factor x→ e−U(x)/f Rd e−U(y) dy. We analyze a sampling method based on the Euler discretization of the Langevin stochastic differential equations under the assumptions that the potential U is continuously differentiable, ∇U is Lipschitz, and U is strongly concave. We focus on the case where the gradient of the log-density cannot be directly computed but unbiased estimates of the gradient from possibly dependent observations are available. This setting can be seen as a combination of a stochastic approximation (here stochastic gradient) type algorithms with discretized Langevin dynamics. We obtain an upper bound of the Wasserstein-2 distance between the law of the iterates of this algorithm and the target distribution π with constants depending explicitly on the Lipschitz and strong convexity constants of the potential and the dimension of the space. Finally, under weaker assumptions on U and its gradient but in the presence of independent observations, we obtain analogous results in Wasserstein-2 distance.
| Original language | English |
|---|---|
| Pages (from-to) | 1-33 |
| Number of pages | 33 |
| Journal | Bernoulli |
| Volume | 27 |
| Issue number | 1 |
| DOIs | |
| Publication status | Published - 1 Feb 2021 |
Keywords
- L-mixing
- Langevin diffusion
- Monte Carlo methods
- Stochastic approximation
Fingerprint
Dive into the research topics of 'On stochastic gradient Langevin dynamics with dependent data streams in the logconcave case'. Together they form a unique fingerprint.Cite this
- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver