partial-sum variance for square-integrable observables
ProvedMarkovChainCLT.integral_sq_sum_coord_le_of_sq_integrableThe variance bound for partial sums, for square-integrable observables. If the -step kernel satisfies with , then for every with — no boundedness assumed —
Removing boundedness is what makes this usable. The underlying covariance estimate is proved through the conditional-expectation form of the Markov property, which is stated for bounded observables. But the intended application is precisely to an unbounded remainder: in the central limit theorem for a square-integrable one truncates, , proves the theorem for the bounded part, and needs
for the unbounded remainder , whose sup norm does not tend to zero even though its norm does.
Proof. Truncate: let and , so each is bounded and centred and pointwise (dominated convergence gives ). The bounded case applies to each :
On the right, by dominated convergence, with dominating function . On the left, the integrands converge pointwise, so Fatou's lemma gives
Fatou also shows the left-hand side is finite, so the partial sum really is square-integrable. Note that no domination is needed on the left — which is essential, since a dominating function there would require exactly the kind of bound being proved.
import Definitions.Def_MarkovChainPathMeasure import Definitions.Def_MarkovErgodicity import Definitions.Def_MarkovIterKernel import Definitions.Def_TotalVariationDist import Mathlib.Probability.Kernel.Invariance import Mathlib.MeasureTheory.Integral.Bochner.Set open Filter Finset Function MeasurableSpace MeasureTheory ProbabilityTheory open MarkovChainCLT open scoped ENNReal NNReal Topology
theorem MarkovChainCLT.integral_sq_sum_coord_le_of_sq_integrable {X : Type*}
[MeasurableSpace X] (P : Kernel X X) [IsMarkovKernel P] (π : Measure X)
[IsProbabilityMeasure π] (hinv : Kernel.Invariant P π) (N : ℕ) (hN : 1 ≤ N) (ρ : ℝ)
(hρ0 : 0 ≤ ρ) (hρ : 4 * ρ ≤ 1 / 4) (hrate : ∀ x, tvDist (iterKernel P N x) π ≤ ρ)
(r : X → ℝ) (hr : Measurable r) (hL2 : Integrable (fun x => (r x) ^ 2) π)
(hmean : ∫ x, r x ∂π = 0) (n : ℕ) :
∫ ω, (∑ k ∈ Finset.range n, r (ω (k + 1))) ^ 2 ∂(chainMeasure P π)
≤ 4 * N * n * ∫ x, (r x) ^ 2 ∂π := by sorry