English

Can the Adaptive Metropolis Algorithm Collapse Without the Covariance Lower Bound?

Probability 2011-02-09 v1 Statistics Theory Computation Statistics Theory

Abstract

The Adaptive Metropolis (AM) algorithm is based on the symmetric random-walk Metropolis algorithm. The proposal distribution has the following time-dependent covariance matrix at step n+1n+1 Sn=Cov(X1,...,Xn)+ϵI, S_n = Cov(X_1,...,X_n) + \epsilon I, that is, the sample covariance matrix of the history of the chain plus a (small) constant ϵ>0\epsilon>0 multiple of the identity matrix II. The lower bound on the eigenvalues of SnS_n induced by the factor ϵI\epsilon I is theoretically convenient, but practically cumbersome, as a good value for the parameter ϵ\epsilon may not always be easy to choose. This article considers variants of the AM algorithm that do not explicitly bound the eigenvalues of SnS_n away from zero. The behaviour of SnS_n is studied in detail, indicating that the eigenvalues of SnS_n do not tend to collapse to zero in general.

Keywords

Cite

@article{arxiv.0911.0522,
  title  = {Can the Adaptive Metropolis Algorithm Collapse Without the Covariance Lower Bound?},
  author = {Matti Vihola},
  journal= {arXiv preprint arXiv:0911.0522},
  year   = {2011}
}

Comments

31 pages, 1 figure

R2 v1 2026-06-21T14:06:48.675Z