50 COUPLING FROM THE PAST
occurs, and t is doubled to 2. These two steps can be viewed as the step from X
−3
to
X
−2
and the step from X
−2
to X
−1
If coalescence occurs, then one more step in the Markov chain is taken, and recall
this was the step from X
−1
up to X
0
. The two blocks together have moved the chain
from step X
−3
up to step X
0
.
If coalescence does not occur. then now t = 4, and steps X
−6
→ X
−5
→ X
−4
→
X
−3
are taken. Then as the recursion unpacks, the two steps from X
−3
→X
−2
→X
−1
are taken, finally followed by the step X
−1
→ X
0
.
No matter how far back the recursion goes, the final output of the procedure (from
this point of view) is X
0
∼
π
. This way of viewing CFTP becomes very important for
designing perfect simulation algorithms for spatial point processes (see Chapter 7.)
3.3 Monotonic CFTP
To use CFTP requires an update function that is stationary with respect to
π
,anda
set A such that for U ∈ A,
φ
(Ω,U)={x}. The stationary requirement is is usually
easy, as it can be done by constructing a Markov chain whose stationary distribution
is
π
. The hard part is creating the sets A
t(i)
together with a method for determining if
U ∈ A
t(i)
.
It is this difficulty that prevents CFTP from being used for every Monte Carlo
application. In the remainder of this chapter, different methods for finding such A
t(i)
will be discussed, together with applications. The first method considered is the use
of a monotonic update function.
Monotonicity was used by Johnson [74] to detect the speed of convergence for
the Gibbs chain for the Ising model, although the term monotonicity was not used
there. Propp and Wilson [112] saw that this same idea could be used as the basis of a
CFTP-type algorithm for perfect simulation.
The idea is as follows. Suppose that the state space Ω admits a partial order.
Recall Definition 1.24, which states that the binary relation is a partial order if it is
reflexive (so s s), antisymmetric (so s t and t s gives s = t), and transitive (so
s t and t r gives s r.) But unlike a total order (such as ≤ over R ), there could
be elements s and t for which neither s t or t s holds.
As an example, consider two Ising configurations x and y. Say that x y if
for all nodes v, x(v) ≤ y(v). So if there are four nodes, (−1, −1, 1,1) precedes
(
−1,−1,−1,1). However, the states (−1 , −1, 1, 1) and (1,1,−1,−1) are incompa-
rable; neither one precedes the other.
Suppose further that Ω has a largest element x
max
, and a smallest element x
min
.
Continuing the Ising model example, the configuration of all −1 values is the small-
est element, and the configuration of all 1 values is the largest element.
Definition 3.3. Let
φ
be an update function for a single step of a Markov chain. If
φ
has the property that for all x y and u ∈ [0,1], it holds that
φ
(x,u)
φ
(y,u),then
call
φ
monotonic.
Consider running two Markov chains {X
t
} and {Y
t
} wh ere X
0
= x
0
and Y
0
= y
0
with x
0
y
0
. Suppose at each step, X
t+1
=
φ
(X
t
,U
t
) and Y
t+1
=
φ
(Y
t
,U
t
),where
U
1
,U
2
,... is an iid stream of Unif([0,1]) random variables.