Perfect sampling using bounding chains. (Q1879888)
From MaRDI portal
scientific article
Language | Label | Description | Also known as |
---|---|---|---|
English | Perfect sampling using bounding chains. |
scientific article |
Statements
Perfect sampling using bounding chains. (English)
0 references
15 September 2004
0 references
The basic framework of bounding chains and their application to several interesting problems drawn from statistical mechanics, graph theory and the approximation of NP-complete problems is considered. Successful application of Monte Carlo Markov chain techniques requires to find the mixing time of these chains which is, in general, extremely difficult. Bounding chains give a theoretical and experimental bound on the mixing time. Moreover, the perfect sampling algorithms generate variates exactly from the target distribution without the need to know the mixing time at all. A technique is applied to a finite Markov chain, denoted by \(M\), with a state space \(\Omega \subseteq C^V\) where \(V\) is the set of dimensions and \(C\) is the set of colors. The colorings \(c(v) \in C\) for all \(v \in V\), which satisfy some preassigned restrictions, are considered. The goal is to generate random variates from the stationary distribution \(\pi\) on the set of colorings. The bounding chain \(M'\) with state space \((2^C)^V\), where \(2^C\) is the set of subsets of \(C\), is defined by the requirement that there exists a coupling \((X_t,Y_t)\), \(t=0,1,\dots\), between \(M\) and \(M'\) such that \[ X_t(v) \in Y_t(v) \quad \forall v \in V \Rightarrow X_{t+1}(v) \in Y_{t+1}(v) \quad \forall v \in V, \quad t=0,1,\dots. \] Here \(X_t\) is a stochastic process evolving according to \(M\), so that each \(v \in V\) is given a single \(c \in C\) in \(X_t\), and \(Y_t\) is a stochastic process evolving according to \(M'\), so that each \(v \in V\) is given a subset from \(C\) in \(Y_t\). If \(Y_0\) bounds every state in \(\Omega\), then when \(Y_t\) bounds just one state \(x\) it can be accepted as the variate from the stationary distribution \(\pi\). The bounding chains are presented for transposition chain on permutations, the hard core gas model, proper colorings of a graph, the antiferromagnetic Potts model and sink free orientations of a graph. Estimations of running time are given.
0 references
Monte Carlo Markov chains
0 references
perfect simulation
0 references
coupling from the past
0 references
mixing times
0 references
proper colorings
0 references
Potts model
0 references
sink free orientation
0 references