Chan–Geyer and polynomial-ergodicity CLTs (Jones Cor 2)
ProvedMarkovChainCLT.clt_of_geometric_or_polynomialLet 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 . Assume one of the following three conditions: (1) the chain is geometrically ergodic and for some ; (2) the chain is polynomially ergodic of order with for the rate constant , and for some with ; (3) the chain is polynomially ergodic of order with , and -almost surely 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,
Case (1) is the Chan–Geyer CLT, the most frequently cited sufficient condition in the MCMC literature; cases (2)–(3) are its polynomial analogues.
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 2** (case 1: Chan–Geyer 1994): a Harris ergodic chain satisfying one
of: (1) geometric ergodicity with `E_π |f|^{2+δ} < ∞` for some `δ > 0`;
(2) polynomial ergodicity of order `m` with integrable constant and
`E_π |f|^{2+δ} < ∞` with `mδ > 2+δ`; (3) polynomial ergodicity of order `m > 1`
with integrable constant and `f` bounded `π`-a.s. — satisfies the CLT for every
initial distribution. -/
theorem MarkovChainCLT.clt_of_geometric_or_polynomial {X : Type*} [MeasurableSpace X]
(P : Kernel X X) [IsMarkovKernel P] (π : Measure X) [IsProbabilityMeasure π]
(hP : HarrisErgodic P π) (f : X → ℝ) (hf : Measurable f)
(hcase :
(GeometricallyErgodic P π ∧
∃ δ : ℝ, 0 < δ ∧ Integrable (fun x => |f x| ^ (2 + δ)) π) ∨
(∃ m δ : ℝ, 0 < δ ∧ 2 + δ < m * δ ∧ PolynomiallyErgodicL1 P π m ∧
Integrable (fun x => |f x| ^ (2 + δ)) π) ∨
(∃ m : ℝ, 1 < m ∧ PolynomiallyErgodicL1 P π m ∧
∃ B : ℝ, ∀ᵐ x ∂π, |f x| < B)) :
SatisfiesCLT P π f := by sorry