for any symmetric dominating
ProvedTierneyMH.Mixture.alphaMH_density_formmarkov-chain-monte-carlomarkov-kernelmetropolis-hastingsp2o-batch-p200bp2o-gran-per-chapterp2o-plan-paperp2o-v1
Let be a probability measure and a Markov proposal kernel on , and let . Let be any -finite measure on that is symmetric, , and dominates , and let . Then, as measures on ,
This is the first step of the proof of Proposition 5, where it is applied to the mixture proposal: the part of the joint law of a move that is accepted has density with respect to every symmetric dominating measure, not only the one used to define .
Formalization Note The page writes the first equality as for the mixture ; the statement gives the resulting min form, for an arbitrary Markov proposal kernel. is defined through the specific measure ; the theorem relates it to an arbitrary symmetric dominating .
Preamble
import Mathlib import Definitions.Def_TierneyMH_Mixture_alphaMH open MeasureTheory ProbabilityTheory open scoped ENNReal
Formal statement
namespace TierneyMH.Mixture
/-- **Density form of `μ α_MH`** (Tierney 1998, proof of Proposition 5, p. 8, first two lines
of the display). Let `ν` be any σ-finite symmetric measure on `E × E` (`ν(dx, dy) = ν(dy, dx)`)
dominating `μ(dx, dy) = π(dx) Q(x, dy)`, and let `h = dμ/dν`. Then
`π(dx) Q(x, dy) α_MH(x, y) = min{h(y, x), h(x, y)} ν(dx, dy)`
as measures on `E × E`. -/
theorem alphaMH_density_form {E : Type*} [MeasurableSpace E]
(π : Measure E) [IsProbabilityMeasure π] (Q : Kernel E E) [IsMarkovKernel Q]
(ν : Measure (E × E)) [SigmaFinite ν] (hν_symm : ν.map Prod.swap = ν)
(hμν : π ⊗ₘ Q ≪ ν) :
(π ⊗ₘ Q).withDensity (alphaMH π Q) =
ν.withDensity (fun p => min ((π ⊗ₘ Q).rnDeriv ν p.swap) ((π ⊗ₘ Q).rnDeriv ν p)) := by sorry
end TierneyMH.Mixture
Source
L. Tierney, A Note on Metropolis–Hastings Kernels for General State Spaces, Ann. Appl. Probab. 8(1) (1998) 1–9, DOI 10.1214/aoap/1027961031, p. 8, proof of Proposition 5 (first two lines of the display)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.