Unbiased time-average estimators for Markov chains

We consider a time-average estimator of a functional of a Markov chain. Under a coupling assumption, we show that the expectation of has a limit as the number of time-steps goes to infinity. We describe a modification of that yields an unbiased estimator of . It is shown that is square-integrable and has finite expected running time. Under certain conditions, can be built without any precomputations, and is asymptotically at least as efficient as , up to a multiplicative constant arbitrarily close to . Our approach provides an unbiased estimator for the bias of . We study applications to volatility forecasting, queues, and the simulation of high-dimensional Gaussian vectors. Our numerical experiments are consistent with our theoretical findings.
View on arXiv