The Annals of Statistics
- Ann. Statist.
- Volume 37, Number 4 (2009), 2011-2055.
On the computational complexity of MCMC-based estimators in large samples
Alexandre Belloni and Victor Chernozhukov
Abstract
In this paper we examine the implications of the statistical large sample theory for the computational complexity of Bayesian and quasi-Bayesian estimation carried out using Metropolis random walks. Our analysis is motivated by the Laplace–Bernstein–Von Mises central limit theorem, which states that in large samples the posterior or quasi-posterior approaches a normal density. Using the conditions required for the central limit theorem to hold, we establish polynomial bounds on the computational complexity of general Metropolis random walks methods in large samples. Our analysis covers cases where the underlying log-likelihood or extremum criterion function is possibly nonconcave, discontinuous, and with increasing parameter dimension. However, the central limit theorem restricts the deviations from continuity and log-concavity of the log-likelihood or extremum criterion function in a very specific manner.
Under minimal assumptions required for the central limit theorem to hold under the increasing parameter dimension, we show that the Metropolis algorithm is theoretically efficient even for the canonical Gaussian walk which is studied in detail. Specifically, we show that the running time of the algorithm in large samples is bounded in probability by a polynomial in the parameter dimension d and, in particular, is of stochastic order d2 in the leading cases after the burn-in period. We then give applications to exponential families, curved exponential families and Z-estimation of increasing dimension.
Article information
Source
Ann. Statist., Volume 37, Number 4 (2009), 2011-2055.
Dates
First available in Project Euclid: 18 June 2009
Permanent link to this document
https://projecteuclid.org/euclid.aos/1245332839
Digital Object Identifier
doi:10.1214/08-AOS634
Mathematical Reviews number (MathSciNet)
MR2533478
Zentralblatt MATH identifier
1175.65015
Subjects
Primary: 65C05: Monte Carlo methods
Secondary: 65C60: Computational problems in statistics
Keywords
Markov chain Monte Carlo computational complexity Bayesian increasing dimension
Citation
Belloni, Alexandre; Chernozhukov, Victor. On the computational complexity of MCMC-based estimators in large samples. Ann. Statist. 37 (2009), no. 4, 2011--2055. doi:10.1214/08-AOS634. https://projecteuclid.org/euclid.aos/1245332839
