Chain CLT from a TV rate: , (Jones Cor 1)
ProvedMarkovChainCLT.clt_of_tv_rateLet 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 the total-variation rate bound holds for all and all , with integrable with respect to and nonincreasing, and that for some ,
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 source's master corollary (its eq. (11)): any total-variation rate plus a matching moment yields the CLT, uniformly over initial distributions; all remaining chain CLTs of the mission are specializations.
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.
import Definitions.Def_MarkovErgodicity
import Definitions.Def_MarkovChainPathMeasure
open MeasureTheory ProbabilityTheory Filter
open scoped ENNReal NNReal Topology ProbabilityTheory
/-- **Corollary 1**: a Harris ergodic chain with total-variation rate `γ`
(nonnegative, nonincreasing) and integrable constant `M`, and a functional with
`E_π |f|^{2+δ} < ∞` for a `δ > 0` such that `∑_n γ(n)^{δ/(2+δ)} < ∞`, satisfies the
CLT for every initial distribution. -/
theorem MarkovChainCLT.clt_of_tv_rate {X : Type*} [MeasurableSpace X]
(P : Kernel X X) [IsMarkovKernel P] (π : Measure X) [IsProbabilityMeasure π]
(hP : HarrisErgodic P π) (f : X → ℝ) (hf : Measurable f)
(M : X → ℝ) (hM0 : ∀ x, 0 ≤ M x) (hM : Integrable M π)
(γ : ℕ → ℝ) (hγ0 : ∀ n, 0 ≤ γ n) (hγa : Antitone γ)
(hrate : ErgodicWithRate P π M γ)
(δ : ℝ) (hδ : 0 < δ) (hmom : Integrable (fun x => |f x| ^ (2 + δ)) π)
(hsum : Summable (fun n => γ n ^ (δ / (2 + δ)))) :
SatisfiesCLT P π f := by sorry