Abstract
This paper defines an approximation scheme for a solution of the Poisson equation of a geometrically ergodic Metropolis-Hastings chain Φ. The scheme is based on the idea of weak approximation and gives rise to a natural sequence of control variates for the ergodic average Sk(F) = (1/k)Pki=1 F(Φi), where F is the force function in the Poisson equation. The main results show that the sequence of the asymptotic variances (in the CLTs for the control-variate estimators) converges to zero and give a rate of this convergence. Numerical examples in the case of a double-well potential are discussed.
Original language | English |
---|---|
Pages (from-to) | 2401-2428 |
Journal | BERNOULLI |
Volume | 24 |
Issue number | 3 |
Early online date | 19 Apr 2017 |
Publication status | E-pub ahead of print - 19 Apr 2017 |