Globally-centered autocovariances in MCMC

09/03/2020
by   Medha Agarwal, et al.
0

Autocovariances are a fundamental quantity of interest in Markov chain Monte Carlo (MCMC) simulations with autocorrelation function (ACF) plots being an integral visualization tool for performance assessment. Unfortunately, for slow mixing Markov chains, the empirical autocovariance can highly underestimate the truth. For multiple-chain MCMC sampling, we propose a globally-centered estimator of the autocovariance function (G-ACvF) that exhibits significant theoretical and empirical improvements. We show that the bias of the G-ACvF estimator is smaller than the bias of the current state-of-the-art. The impact of this improved estimator is evident in three critical output analysis applications: (1) ACF plots, (2) estimates of the Monte Carlo asymptotic covariance matrix, and (3) estimates of the effective sample size. Under weak conditions, we establish strong consistency of our improved asymptotic covariance estimator, and obtain its large-sample bias and variance. The performance of the new estimators is demonstrated through various examples.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
07/08/2020

Estimating Monte Carlo variance from multiple Markov chains

The ever-increasing power of the personal computer has led to easy paral...
research
05/21/2018

Weighted batch means estimators in Markov chain Monte Carlo

This paper proposes a family of weighted batch means variance estimators...
research
09/12/2018

Lugsail lag windows and their application to MCMC

Lag windows are commonly used in the time series, steady state simulatio...
research
06/09/2015

Measuring Sample Quality with Stein's Method

To improve the efficiency of Monte Carlo estimation, practitioners are t...
research
03/19/2023

Provable Convergence of Variational Monte Carlo Methods

The Variational Monte Carlo (VMC) is a promising approach for computing ...
research
10/11/2018

Efficient estimation of autocorrelation spectra

The performance of Markov chain Monte Carlo calculations is determined b...
research
11/03/2019

Estimating accuracy of the MCMC variance estimator: a central limit theorem for batch means estimators

The batch means estimator of the MCMC variance is a simple and effective...

Please sign up or login with your details

Forgot password? Click here to reset