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

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+1 \[ Sn = Cov(X1,...,Xn) + ε I, \] that is, the sample covariance matrix of the history of the chain plus a (small) constant ε>0 multiple of the identity matrix I. The lower bound on the eigenvalues of Sn induced by the factor ε I is theoretically convenient, but practically cumbersome, as a good value for the parameter ε may not always be easy to choose. This article considers variants of the AM algorithm that do not explicitly bound the eigenvalues of Sn away from zero. The behaviour of Sn is studied in detail, indicating that the eigenvalues of Sn do not tend to collapse to zero in general.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…