Chen characterization: CLT bounded in probability (Jones Thm 4)
OpenMarkovChainCLT.clt_iff_boundedInProbabilityLet be a Markov chain with transition kernel , Harris ergodic with invariant probability , and let be measurable with and . Then, for the stationary chain (initial distribution ), the following are equivalent:
Chen's characterization shows the Markov chain CLT is equivalent to mere tightness of the normalized averages — the sharpest available dividing line for square-integrable functionals. By the source's Remark 2, the CLT side is then automatic for every initial distribution.
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_MarkovChainPathMeasure import Definitions.Def_MixingCoefficients open MeasureTheory ProbabilityTheory Filter open scoped ENNReal NNReal Topology ProbabilityTheory /-- **Theorem 4** (Chen 1999): for a Harris ergodic chain and a centered square-integrable functional, the stationary chain satisfies a CLT `√n f̄_n →d N(0, σ²)` for some `σ² ≥ 0` **iff** the sequence `√n f̄_n` is bounded in probability. -/
theorem MarkovChainCLT.clt_iff_boundedInProbability {X : Type*} [MeasurableSpace X]
[MeasurableSpace.CountablyGenerated X]
(P : Kernel X X) [IsMarkovKernel P] (π : Measure X) [IsProbabilityMeasure π]
(hP : HarrisErgodic P π) (f : X → ℝ) (hf : Measurable f)
(hcent : ∫ x, f x ∂π = 0) (hL2 : MemLp f 2 π) :
(∃ v : ℝ≥0, TendstoInDistribution
(fun (n : ℕ) (ω : ℕ → X) => Real.sqrt n * sampleAvg f n ω)
atTop (id : ℝ → ℝ) (fun _ => chainMeasure P π) (gaussianReal 0 v))
↔ BoundedInProbability
(fun (n : ℕ) (ω : ℕ → X) => Real.sqrt n * sampleAvg f n ω)
(chainMeasure P π) := by sorry
Read-back
What the Lean code literally says, in plain math · claude-fable-5
Setting. is an arbitrary type with a measurable-space structure, together with the typeclass assumption that its -algebra is countably generated (generated by some countable family of sets); is a kernel from to assumed (typeclass) to be a Markov kernel, i.e. is a probability measure for every ; is a measure on assumed (typeclass) to be a probability measure. Hypotheses. (1) Harris ergodicity of , which in this bundle means the conjunction of: is invariant for (the measure equals ), and for every starting point , as ; here is the -fold iterate of defined by the identity kernel and , and is this bundle's total-variation distance — with no factor , measure values converted to reals by a convention sending to (immaterial for probability measures), and the supremum taken as a real-number supremum whose junk value is if the set were unbounded (the set always contains , via ). (2) is measurable. (3) . (4) in the MemLp sense: is -almost-everywhere strongly measurable and . Conclusion — an exact if-and-only-if. Left side: there exists a nonnegative real such that the functions converge in distribution to as , under the path measure . Here is the average of over the chain states at times — the time- state is never used, and at Lean's convention together with gives . No centering term is subtracted inside — centering is supplied only by hypothesis (3). (chainMeasure P \pi) is the probability law on trajectory space of the time-homogeneous Markov chain with one-step transition kernel whose initial state is drawn from itself (the stationary start; no other initial distribution appears anywhere in this statement), each being drawn from via the Ionescu–Tulcea trajectory construction. TendstoInDistribution (against the identity map on the probability space ) means for every bounded continuous ; gaussianReal 0 v is the normal law of mean and variance , which is the Dirac point mass at when — so a degenerate limit satisfies the left side. Right side: the same sequence is bounded in probability under , which in this bundle means: for every real there exists a real (not required to be positive) such that for every simultaneously, (the measure value converted to a real, with — immaterial since is a probability measure); the case is included, where .
Confirmed by the mission captain (proposal self-audit).