Discrete sticky couplings of functional autoregressive processes
In this paper, we provide bounds in Wasserstein and total variation distances between the distributions of the successive iterates of two functional autoregressive processes with isotropic Gaussian noise of the form Y_k+1 = T_γ(Y_k) + √(γσ^2) Z_k+1 and Ỹ_k+1 = T̃_γ(Ỹ_k) + √(γσ^2)Z̃_k+1. More precisely, we give non-asymptotic bounds on ρ(ℒ(Y_k),ℒ(Ỹ_k)), where ρ is an appropriate weighted Wasserstein distance or a V-distance, uniformly in the parameter γ, and on ρ(π_γ,π̃_γ), where π_γ and π̃_γ are the respective stationary measures of the two processes. The class of considered processes encompasses the Euler-Maruyama discretization of Langevin diffusions and its variants. The bounds we derive are of order γ as γ→ 0. To obtain our results, we rely on the construction of a discrete sticky Markov chain (W_k^(γ))_k ∈ℕ which bounds the distance between an appropriate coupling of the two processes. We then establish stability and quantitative convergence results for this process uniformly on γ. In addition, we show that it converges in distribution to the continuous sticky process studied in previous work. Finally, we apply our result to Bayesian inference of ODE parameters and numerically illustrate them on two particular problems.
READ FULL TEXT