Stationarity of the two-coordinate Gibbs sampler (source code)

= Stationarity of the two-coordinate Gibbs sampler

The joint target density is stationary for a full Gibbs sweep, since
$$
\int f_X(x)f_{Y\mid X}(y'\mid x)
f_{X\mid Y}(x'\mid y')\,dx
=f_Y(y')f_{X\mid Y}(x'\mid y')
=f_{XY}(x',y').
$$