DOMINATED COUPLING FROM THE PAST FOR POINT PROCESSES 129
the dominating process Y
t
contains every point in X
t
and possibly some extra ones
because it allows every birth into the state.
Locally stable dominating update Input: x, y Output: x, y
1) Draw T
B
← Exp(M ·m(S))
2) Draw T
D
← Exp(#y)
3) If T
D
< T
B
(then there is a death)
4) Draw v uniformly from y, y ← y −v, x ← x −v
5) Else
6) Draw v uniformly over S
7) Draw U ← Unif([0,1])
8) y ← y + v
8) If U < f (x + v)/[Mf(x)]
9) x ← x + v
This update function is monotonic. To see why, suppose that the input satisfies
x ⊆y.WhenT
D
< T
B
, that means that a point in y is about to die. Since x cannot have
a birth without a birth in y, the death of the point in y must also come before any birth
in x. If this point in y that is about to die is also in x, then it also dies in x. (Note that
if v /∈ x then x −v = x). Either way, after the death event it still holds that x ⊆ y .
On the other hand, if T
B
≤ T
D
, then there is a point v that is added to y. It might
or might not (depending on the value of U) be added to x, but it is always added to y .
Hence, at the end of a birth event, x ⊆ y.
So Y
t
is a dominating process for X
t
. Moreover, Y
t
is a very simple birth-death
process. The stationary distribution for the Y
t
process is a PPP with rate
λ
(x)=M
for all x ∈ S, so it is easy to generate a sample from Y
t
.
Formally Y
t
is a marked process, because for each point in Y
t
, there is a choice of
uniform U
t
called a mark. The marks are sufficient to determine the X
t
,soitisthe
marked Y
t
and X
t
processes that are coupled together.
Now in order to used dominated CFTP, it is also necessary to be able to run the
dominating pr ocess Y
t
backwards in time. This is straightforward from a computa-
tional point of view, since the dominating process is reversible. However, it can get a
bit confusing, as when time is running backward, a “birth” in the reverse process is
a “death” in the forward process. And a “death” in the backward process is a “birth”
in the forward process.
The dominating process can be thought of as maintaining a list of birth and death
events that happen backwards in time. Call this list EL for event list. Let D(0) denote
the configuration of the dominating process at time 0. Then sometime before time 0,
a birth or death event occurred. Let D(−1) denote the state of the dominating process
immediately before this event occurred.
Similarly, let D(−i) denote the state of the dominating process just before time
t
i
, where exactly i birth o r death events o ccurred in interval [t
i
,0). Suppose that the
current event list holds the first n events that happened b efore time 0. This is neces-
sary to be able to extend this list backwards in time to n
> n events. This is done by
the following pseudocode: