20 INTRODUCTION
for Y makes the joint density of (X ,Y ) uniform over X ∈{0,1}
V
and Y ∈{[0, ∞) :
(∀{i, j})(y ({i, j}) ≤ min(exp(
β
),1))}.
The Markov chain then is simple: given the state x, choose a new value of y
conditioned on x, and then choose a new x conditioned on y. (This can be viewed as
a deterministic scan Gibbs sampler applied to the new joint distribution.)
Auxiliary chain Input: old state (x,y) Output: new state (x ,y)
1) Draw y conditioned on x
2) Draw x conditioned on y
From the construction, drawing y conditioned on x is straightforward. Drawing x
conditioned on y can be more complex. Consider again the Ising model with
β
> 0.
Then exp(
β
) > 1. When y({i, j}) ∈ [0, 1], then it could be true that x(i)=x( j) or
x(i) = x( j).Butwheny({i, j}) > 1, it must be true that x(i)=x( j). The edges with
y({i, j}) > 1 break the graph into clusters of connected components. Each component
must have the same x value across the component, but otherwise is uniform.
Therefore, break the graph into connected components u sing the edges with
y({i, j}) > 1, then for each connected component, choose a label uniformly from
{−1,1}and apply that label to every node in the component.
This example was one of the first examples of an auxiliary variable chain, and
is due to Swendsen and Wang [121]. For this reason, auxiliary variable chains are
sometimes known as Swendsen-Wang chains.
Swendsen Wang chain Input: old state ( x,y) Output: new state (x,y)
1) For each edge {i, j}∈E
2) Draw y({i, j}) ← Unif([0,exp(
β
·1(x(i)=x ( j)))])
3) Let A ←{{i, j} : y({i, j}) > 1}
4) For each connected component C in the graph (V,A)
5) Draw B ← Unif({−1,1})
6) For each v ∈ B,setx(v) ← B
1.8.4 Slice sampler
The slice sampler (see [95] for m ore details) is a special case of auxiliary vari-
ables that uses just a single component Y with conditional distribution [Y |X = x] ∼
Unif([0, f
X
(X)]).Then[X |Y = y] ∼ Unif(Ω ∩{x : f
X
(x) ≥ y}).
This last set Ω ∩{x : f
X
(x) ≥ y} is a slice through the density of f
X
at height y,
hence the name slice sampler.
1.8.5 Shift chains
Although Gibbs, Metropolis-Hastings, and auxiliary random variables are the most
widely used types of Markov chains, they are not the only types o f chains.
An important type of chain that is useful for repulsive processes is the shift move.
The canonical example is permutations: they are repulsive in the sense that two com-
ponents of the permutation cannot both be assigned the same value. The solution to