Strongly mixing stationary sequences: CLT uniformly integrable (Jones Thm 3)
ProvedMarkovChainCLT.clt_iff_uniformlyIntegrable_of_alpha_mixingLet be a centered, strictly stationary sequence of real random variables on a probability space, with partial sums and . Suppose , the sequence is strongly mixing (), and . Then the following are equivalent:
This characterization (Cogburn; Denker; Mori–Yoshihara) explains exactly what can fail for dependent sequences with second moments: the CLT is equivalent to uniform integrability of the normalized squares, not implied by moments alone.
Formalization Note For the (finitely many) indices with the normalized quantities are interpreted as . Sequences are indexed from , so and the past -algebras used by the mixing coefficients start at ; under strict stationarity this agrees with the source, which indexes from . Absolute convergence of the covariance series is expressed as unconditional summability, and the limit statement is weak convergence of the laws of .
import Definitions.Def_MixingCoefficients
import Mathlib.MeasureTheory.Function.ConvergenceInDistribution
import Mathlib.Probability.Distributions.Gaussian.Real
import Mathlib.MeasureTheory.Function.UniformIntegrable
open MeasureTheory ProbabilityTheory Filter
open scoped ENNReal NNReal Topology ProbabilityTheory
/-- **Theorem 3** (Cogburn 1960; Denker 1986; Mori–Yoshihara 1986): for a centered
strictly stationary strongly mixing square-integrable sequence with
`σ_n² = E[S_n²] → ∞`, the normalized sums `S_n / σ_n` converge in distribution to
`N(0,1)` **iff** the family `{S_n² / σ_n²}` is uniformly integrable. -/
theorem MarkovChainCLT.clt_iff_uniformlyIntegrable_of_alpha_mixing {Ω : Type*} [MeasurableSpace Ω]
(P : Measure Ω) [IsProbabilityMeasure P] (Y : ℕ → Ω → ℝ)
(hY : ∀ n, Measurable (Y n)) (hstat : IsStrictlyStationary P Y)
(hcent : ∫ ω, Y 0 ω ∂P = 0) (hL2 : MemLp (Y 0) 2 P)
(hmix : Tendsto (fun n => alphaMixingCoef P Y n) atTop (𝓝 0))
(hvar : Tendsto (fun n => ∫ ω, (∑ i ∈ Finset.range n, Y i ω) ^ 2 ∂P) atTop atTop) :
TendstoInDistribution
(fun (n : ℕ) ω => (∑ i ∈ Finset.range n, Y i ω)
/ Real.sqrt (∫ ω', (∑ i ∈ Finset.range n, Y i ω') ^ 2 ∂P))
atTop (id : ℝ → ℝ) (fun _ => P) (gaussianReal 0 1)
↔ UniformIntegrable
(fun (n : ℕ) ω => (∑ i ∈ Finset.range n, Y i ω) ^ 2
/ ∫ ω', (∑ i ∈ Finset.range n, Y i ω') ^ 2 ∂P) 1 P := by sorry