Eqs. (2)–(3): the Metropolis–Hastings kernel satisfies detailed balance iff its off-diagonal part does
ProvedTierneyMH.Reversibility.mhKernel_reversible_iff_offDiagonalLet be a probability measure on a measurable space , let be a Markov transition kernel on , let be measurable, and let be the Metropolis–Hastings kernel . Then satisfies the detailed balance relation
if and only if
both sides being measures on . That is, the diagonal (rejection) component of does not matter.
This reduces reversibility of the Metropolis–Hastings kernel to an identity between two measures built from , and alone, which is the starting point of Theorem 2.
Formalization Note Detailed balance (2) is Mathlib's Kernel.IsReversible P π: for all measurable , which is (2) evaluated on measurable rectangles. The measure is (π ⊗ₘ Q).withDensity α, and the right side of (3) is its image under the swap . is ℝ≥0∞-valued with .
import Mathlib import Definitions.Def_TierneyMH_Shared_mhKernel open TierneyMH.Shared open MeasureTheory ProbabilityTheory open scoped ENNReal
namespace TierneyMH.Reversibility
/-- **§2, Eqs. (2)–(3)** (Tierney 1998, p. 2). The Metropolis–Hastings kernel `P = mhKernel Q α`
of Eq. (1) satisfies detailed balance (2) `π(dx) P(x, dy) = π(dy) P(y, dx)` if and only if
(3) `π(dx) Q(x, dy) α(x, y) = π(dy) Q(y, dx) α(y, x)` holds as an identity of measures on
`E × E`: the diagonal (rejection) component does not matter.
Detailed balance (2) is Mathlib's `Kernel.IsReversible P π`
(`∫_A P(x, B) π(dx) = ∫_B P(x, A) π(dx)` for measurable `A, B`). The measure
`π(dx) Q(x, dy) α(x, y)` is `(π ⊗ₘ Q).withDensity α`, and `π(dy) Q(y, dx) α(y, x)` is its image
under `Prod.swap`. -/
theorem mhKernel_reversible_iff_offDiagonal {E : Type*} [MeasurableSpace E]
(π : Measure E) [IsProbabilityMeasure π] (Q : Kernel E E) [IsMarkovKernel Q]
(α : E × E → ℝ≥0∞) (hα_meas : Measurable α) (hα_le : ∀ p, α p ≤ 1) :
Kernel.IsReversible (mhKernel Q α) π ↔
((π ⊗ₘ Q).withDensity α).map Prod.swap = (π ⊗ₘ Q).withDensity α := by sorry
end TierneyMH.Reversibility
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.