CLT under geometric drift: , (Jones Thm 1(i))
ProvedMarkovChainCLT.clt_of_geometric_driftLet be a Markov chain with transition kernel on a state space , Harris ergodic with invariant probability distribution , and let be measurable. Write for the sample average and . Suppose is measurable, is a measurable small set, and are constants, the geometric drift condition
holds with integrable under every , and pointwise.
Then the chain satisfies the central limit theorem for : there is a single asymptotic variance such that for every initial distribution of the chain,
This is the workhorse CLT of applied Markov chain Monte Carlo: drift towards a small set is the standard checkable route to a CLT for a specific sampler (Meyn–Tweedie, Theorem 17.0.1).
Formalization Note "Harris ergodic" is encoded by its total-variation characterization: is invariant for and for every starting point (equivalent to the classical aperiodic, -irreducible, positive Harris recurrent definition; the "every " quantifier is exactly the Harris property). Convergence in distribution is weak convergence of laws, and is read as the point mass at , which absorbs the source's "" caveat. The -algebra of the state space is additionally assumed countably generated, the standard general-state-space setting of Meyn and Tweedie.
import Definitions.Def_MarkovErgodicity import Definitions.Def_MarkovDriftMinorization import Definitions.Def_MarkovChainPathMeasure open MeasureTheory ProbabilityTheory Filter open scoped ENNReal NNReal Topology ProbabilityTheory /-- **Theorem 1, condition 1** (Meyn–Tweedie 1993, Theorem 17.0.1): a Harris ergodic chain satisfying the geometric drift condition towards a small set, with `f² ≤ V`, satisfies the CLT for every initial distribution. -/
theorem MarkovChainCLT.clt_of_geometric_drift {X : Type*} [MeasurableSpace X]
[MeasurableSpace.CountablyGenerated X]
(P : Kernel X X) [IsMarkovKernel P] (π : Measure X) [IsProbabilityMeasure π]
(hP : HarrisErgodic P π) (f : X → ℝ) (hf : Measurable f)
(V : X → ℝ) (hV : Measurable V) (hV1 : ∀ x, 1 ≤ V x)
(C : Set X) (hC : MeasurableSet C) (hsmall : IsSmallSet P C)
(d b : ℝ) (hd : 0 < d) (hdrift : GeoDriftCondition P V d b C)
(hfV : ∀ x, f x ^ 2 ≤ V x) :
SatisfiesCLT P π f := by sorry