Conditional sampling in Turing

I’m not entirely sure I understand what you want to do. But one thing you might be able to do is to model the equality of x and y in terms of a Bernoulli distributed RV.

I don’t quite understand what you mean by “condition on the observation that x and y take the same values”.