Proposition 5: the maximal kernel of a mixture proposal dominates the mixture of maximal kernels,
ProvedTierneyMH.Mixture.proposition5_mixtureLet be a measurable space with measurable singletons and a countably generated -algebra, and let be a probability measure on (the target). Let be a finite or countable family of Markov proposal kernels on and let with . Let be the maximal Metropolis–Hastings kernel for , and let be the maximal Metropolis–Hastings kernel for the mixture proposal . Then
that is, for -almost every ,
There are two ways of combining Metropolis–Hastings samplers built on proposals : mix the samplers, or run a single sampler with the mixture proposal. Combined with Peskun's theorem for general state spaces (Theorem 4 of the paper), the proposition shows that the second has asymptotic variances of sample-path averages no larger than the first, whatever the weights.
Formalization Note Two hypotheses implicit in the paper are made explicit: singletons are measurable, as the sets and point masses require; and the -algebra is countably generated. The second is an addition: the proof in the paper shows an inequality between measures on , and passing to the kernel-level statement, with one -null set of exceptional for all simultaneously, uses a countable generating family. The weights are in ℝ≥0 with HasSum β 1, and the index type is countable.
import Mathlib import Definitions.Def_TierneyMH_Shared_OffDiagDominates import Definitions.Def_TierneyMH_Mixture_maxMHKernel import Definitions.Def_TierneyMH_Mixture_mixKernel open TierneyMH.Shared open MeasureTheory ProbabilityTheory open scoped ENNReal NNReal
namespace TierneyMH.Mixture
/-- **Proposition 5** (Tierney 1998, p. 7). Let `Qᵢ` be a (countable) sequence of Markov
proposal kernels and `βᵢ ≥ 0` with `∑ βᵢ = 1`. Let `Pᵢ` be the maximal Metropolis–Hastings
kernel for `Qᵢ` and `P` the maximal Metropolis–Hastings kernel for the mixture proposal
`Q = ∑ βᵢ Qᵢ`, all for the target `π`. Then `P ⪰ ∑ βᵢ Pᵢ`: for `π`-almost every `x`,
`P(x, A \ {x}) ≥ ∑ βᵢ Pᵢ(x, A \ {x})` for all measurable `A`.
Added hypotheses: measurable singletons (implicit in the paper's `A \ {x}`), and a countably
generated σ-algebra on `E`, which makes the exceptional `π`-null set uniform over `A`. -/
theorem proposition5_mixture {E : Type*} [MeasurableSpace E] [MeasurableSingletonClass E]
[MeasurableSpace.CountablyGenerated E]
(π : Measure E) [IsProbabilityMeasure π]
{ι : Type*} [Countable ι] (Q : ι → Kernel E E) [∀ i, IsMarkovKernel (Q i)]
(β : ι → ℝ≥0) (hβ : HasSum β 1) :
OffDiagDominates π (maxMHKernel π (mixKernel β Q))
(mixKernel β fun i => maxMHKernel π (Q i)) := by sorry
end TierneyMH.Mixture
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.