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 Yk+1 = Tγ (Yk)+ √ γσ2Zk+1 and Yk+1 = Tγ ( Yk)+ √ γσ2 Zk+1. More precisely, we give nonasymptotic bounds on ρ(L(Yk),L( Yk)), 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∈N 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 Howitt (Ph.D. thesis (2007)) and Eberle and Zimmer (Ann. Inst. Henri Poincaré Probab. Stat. 55 (2019) 2370- 2394). Finally, we apply our result to Bayesian inference of ODE parameters and numerically illustrate them on two particular problems.