Prove2Me
Navigate
DiscoverFormalpediaBlogsUsersMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Variance asymptotics: n−1E[Sn2]→σ2n^{-1}E[S_n^2]\to\sigma^2n−1E[Sn2​]→σ2

Proved
MarkovChainCLT.tendsto_inv_mul_integral_sq_partialSum

by Gabewhigham · Sep 5, 2026 · Mathlib c5ea003 (Lean v4.30.0)

central-limit-theoremprobabilitystationary-processes

Let Y=(Yi)i≥0Y=(Y_i)_{i\ge0}Y=(Yi​)i≥0​ be a measurable, strictly stationary real sequence on a probability space with Y0∈L2Y_0\in L^2Y0​∈L2, write Sn=∑i<nYiS_n=\sum_{i<n}Y_iSn​=∑i<n​Yi​ and γk=E[Y0Yk]\gamma_k=E[Y_0Y_k]γk​=E[Y0​Yk​], and assume that the positive-lag autocovariance series ∑k≥1γk\sum_{k\ge1}\gamma_k∑k≥1​γk​ converges. Then the variance of the partial sums grows linearly, with slope the asymptotic variance:

1n E[Sn2] ⟶ σ2=γ0+2∑k≥1γk.\frac1n\,E\bigl[S_n^2\bigr]\ \longrightarrow\ \sigma^2 = \gamma_0 + 2\sum_{k\ge1}\gamma_k .n1​E[Sn2​] ⟶ σ2=γ0​+2k≥1∑​γk​.

This is the standard variance asymptotics underlying every central limit theorem for stationary sequences: it is what identifies the limiting Gaussian variance and, when σ2>0\sigma^2>0σ2>0, what makes the two natural normalizations n\sqrt nn​ and E[Sn2]\sqrt{E[S_n^2]}E[Sn2​]​ interchangeable.

Stationarity gives E[YiYi+k]=γkE[Y_iY_{i+k}]=\gamma_kE[Yi​Yi+k​]=γk​, so expanding the square yields E[Sn2]=nγ0+2∑k=1n−1(n−k)γkE[S_n^2]=n\gamma_0+2\sum_{k=1}^{n-1}(n-k)\gamma_kE[Sn2​]=nγ0​+2∑k=1n−1​(n−k)γk​, and after dividing by nnn the weights 1−k/n1-k/n1−k/n increase to 111 while being dominated by the summable sequence ∣γk∣|\gamma_k|∣γk​∣.

Formalization Note Sequences are indexed from 000, so Sn=Y0+⋯+Yn−1S_n=Y_0+\dots+Y_{n-1}Sn​=Y0​+⋯+Yn−1​; the asymptotic variance is the platform's seqAsymptoticVariance, namely E[Y02]+2∑k≥0E[Y0Yk+1]E[Y_0^2]+2\sum_{k\ge0}E[Y_0Y_{k+1}]E[Y02​]+2∑k≥0​E[Y0​Yk+1​].

Preamble
import Definitions.Def_MixingCoefficients

open MeasureTheory ProbabilityTheory Filter MarkovChainCLT
open scoped ENNReal NNReal Topology ProbabilityTheory
Formal statement
theorem MarkovChainCLT.tendsto_inv_mul_integral_sq_partialSum
    {Ω : Type*} [MeasurableSpace Ω]
    (P : Measure Ω) [IsProbabilityMeasure P] (Y : ℕ → Ω → ℝ)
    (hY : ∀ n, Measurable (Y n)) (hstat : IsStrictlyStationary P Y)
    (hL2 : MemLp (Y 0) 2 P)
    (hsum : Summable fun k : ℕ => ∫ ω, Y 0 ω * Y (k + 1) ω ∂P) :
    Tendsto (fun n : ℕ => (n : ℝ)⁻¹ * ∫ ω, (∑ i ∈ Finset.range n, Y i ω) ^ 2 ∂P) atTop
      (𝓝 (seqAsymptoticVariance P Y)) := by sorry
Source
G. L. Jones, "On the Markov Chain Central Limit Theorem", Probability Surveys 1 (2004) 299-320, arXiv math/0409112v2, Section 4 (the variance sigma^2 of Theorems 5-8, eq. (12)); the identity E[S_n^2] = n gamma_0 + 2 sum_{k<n} (n-k) gamma_k with Cesaro passage to the limit is classical, see P. Billingsley, Probability and Measure, 3rd ed. (1995), Section 27.

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactJoin Slack© 2026 Prove2Me