CLT for exp-mixing sequences given variance convergence
ProvedMarkovChainCLT.clt_of_var_limit_of_exp_alphaLet be a centered strictly stationary real-valued sequence on a probability space with exponentially decaying strong mixing coefficients, for some , and with . Assume moreover that the normalized variances converge,
where , and that . Then
This is the blocking half of the Doukhan-Massart-Rio central limit theorem (Jones, Theorem 6): Bernstein big-block/small-block decomposition makes distant blocks asymptotically independent at an exponential rate, so the normalized sum inherits the Gaussian limit from the independent-block approximation once the variance is known to stabilize. It takes the variance limit as a hypothesis, complementing the separately proved variance-convergence lemma. Formalization Note Convergence is weak convergence of the laws under the common probability measure ; the Gaussian variance is the nonnegative-real coercion of , which equals under the positivity hypothesis.
import Definitions.Def_MixingCoefficients import Mathlib.MeasureTheory.Function.ConvergenceInDistribution import Mathlib.Probability.Distributions.Gaussian.Real import Mathlib.Analysis.SpecialFunctions.Log.PosLog open MeasureTheory ProbabilityTheory Filter MarkovChainCLT open scoped ENNReal NNReal Topology ProbabilityTheory /-- The blocking/limit component of the Doukhan-Massart-Rio CLT, given variance convergence. -/
theorem MarkovChainCLT.clt_of_var_limit_of_exp_alpha
{Ω : Type*} [MeasurableSpace Ω]
(P : Measure Ω) [IsProbabilityMeasure P] (Y : ℕ → Ω → ℝ)
(hY : ∀ n, Measurable (Y n)) (hstat : IsStrictlyStationary P Y)
(hcent : ∫ ω, Y 0 ω ∂P = 0)
(c a : ℝ) (ha0 : 0 ≤ a) (ha1 : a < 1)
(hα : ∀ n, alphaMixingCoef P Y n ≤ c * a ^ n)
(hmom : Integrable (fun ω => (Y 0 ω) ^ 2 * Real.posLog |Y 0 ω|) P)
(hvarlim : Tendsto (fun n : ℕ => Var[∑ i ∈ Finset.range n, Y i; P] / (n : ℝ))
atTop (𝓝 (seqAsymptoticVariance P Y)))
(hvar : 0 < seqAsymptoticVariance P Y) :
TendstoInDistribution
(fun (n : ℕ) ω => (Real.sqrt n)⁻¹ * ∑ i ∈ Finset.range n, Y i ω)
atTop (id : ℝ → ℝ) (fun _ => P)
(gaussianReal 0 (seqAsymptoticVariance P Y).toNNReal) := by sorry