Physical Sciences and Mathematics Commons^™

Fixed-Width Output Analysis For Markov Chain Monte Carlo, Galin L. Jones, Murali Haran, Brian S. Caffo, Ronald Neath

Johns Hopkins University, Dept. of Biostatistics Working Papers

Markov chain Monte Carlo is a method of producing a correlated sample in order to estimate features of a complicated target distribution via simple ergodic averages. A fundamental question in MCMC applications is when should the sampling stop? That is, when are the ergodic averages good estimates of the desired quantities? We consider a method that stops the MCMC sampling the first time the width of a confidence interval based on the ergodic averages is less than a user-specified value. Hence calculating Monte Carlo standard errors is a critical step in assessing the output of the simulation. In particular, we …