Theorem 2: a Metropolis–Hastings kernel satisfies detailed balance iff a.e. off and a.e. on
ProvedTierneyMH.Reversibility.theorem2_detailed_balance_iffLet be a probability measure on a measurable space , let be a Markov transition kernel on (the proposal), and let be measurable (the acceptance probability). Let be the Metropolis–Hastings kernel , and put and . Let be any symmetric set on which and are mutually absolutely continuous and off which they are mutually singular, and let be any version of with and everywhere (Proposition 1). Then satisfies detailed balance, , if and only if
- is -almost everywhere zero on , and
- satisfies
This gives necessary and sufficient conditions on the proposal kernel and the acceptance probability for a Metropolis–Hastings sampler to be reversible with respect to on a general state space, with no density or dominating measure assumed. The standard Metropolis–Hastings acceptance probability and many special cases of the literature (common dominating measure, deterministic proposals, dimension-changing moves) are instances.
Formalization Note The paper states the theorem for "the detailed balance condition (4)", , which the text on pp. 2–3 identifies with detailed balance (2) of the kernel; the statement here uses (2) for the kernel itself, as Mathlib's Kernel.IsReversible (mhKernel Q α) π. The pair is universally quantified, so the theorem holds for every version. is π ⊗ₘ Q, is its image under the swap, and "-a.e. on " is almost everywhere for the restriction of to . and are ℝ≥0∞-valued.
import Mathlib import Definitions.Def_TierneyMH_Shared_mhKernel import Definitions.Def_TierneyMH_Shared_IsSymmetricSplit import Definitions.Def_TierneyMH_Shared_IsRatioVersion open TierneyMH.Shared open MeasureTheory ProbabilityTheory open scoped ENNReal
namespace TierneyMH.Reversibility
/-- **Theorem 2** (Tierney 1998, p. 3). Let `π` be a probability measure and `Q` a Markov kernel
on a measurable space `E`, let `α : E × E → [0, 1]` be measurable, and put
`μ(dx, dy) = π(dx) Q(x, dy)` (`μ = π ⊗ₘ Q`). Let `(R, r)` be **any** pair given by
Proposition 1 for `μ`. Then the Metropolis–Hastings kernel `mhKernel Q α` of Eq. (1) satisfies
detailed balance with respect to `π` if and only if
* (i) `α = 0` `μ`-almost everywhere on `Rᶜ`, and
* (ii) `α(x, y) r(x, y) = α(y, x)` `μ`-almost everywhere on `R`.
The paper names detailed balance in the form (4) `μ(dx, dy) α(x, y) = μᵀ(dx, dy) α(y, x)`,
which the text on pp. 2–3 identifies with (2) for the kernel (1); the statement here uses (2)
itself, `Kernel.IsReversible (mhKernel Q α) π`. -/
theorem theorem2_detailed_balance_iff {E : Type*} [MeasurableSpace E]
(π : Measure E) [IsProbabilityMeasure π] (Q : Kernel E E) [IsMarkovKernel Q]
(α : E × E → ℝ≥0∞) (hα_meas : Measurable α) (hα_le : ∀ p, α p ≤ 1)
(R : Set (E × E)) (r : E × E → ℝ≥0∞)
(hR : IsSymmetricSplit (π ⊗ₘ Q) R) (hr : IsRatioVersion (π ⊗ₘ Q) R r) :
Kernel.IsReversible (mhKernel Q α) π ↔
((∀ᵐ p ∂((π ⊗ₘ Q).restrict Rᶜ), α p = 0) ∧
(∀ᵐ p ∂((π ⊗ₘ Q).restrict R), α p * r p = α p.swap)) := by sorry
end TierneyMH.Reversibility
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.