Delta method for Markov chain statistics (Lemma EC.5), measurable estimator
ProvedMarkovChainCLT.delta_method_of_uniformly_ergodic_of_measurableLemma EC.5 of arXiv:2407.19618 (linearization / delta method for Markov chain statistics), with the estimator required to be measurable.
Let be a uniformly ergodic Harris chain on with stationary distribution , and let be a vector of square-integrable statistics over a finite index set . Write
Let be Fréchet-differentiable at with derivative . Then
where is the Markov chain asymptotic variance — the variance of the sample average of the linearized observable , i.e. in the notation of the paper.
Why the extra hypothesis. The conclusion asserts convergence in distribution, and in Mathlib that carries an almost-everywhere measurability obligation on each . Differentiability of at the single point does not supply it: a function of the form with non-measurable is Fréchet-differentiable at with derivative and is not measurable anywhere else. This statement therefore carries Measurable g explicitly. The paper is on the right side of this already — Definition 1 speaks of differentiable estimators and Appendix EC.3.2 uses " is differentiable (and thus continuous)" — so the hypothesis costs nothing in the intended application, where the estimator is a genuinely differentiable, hence continuous, hence measurable, function of the statistics.
What drives the proof. Because is finite, is a finite linear combination of coordinates, so is the sample average of the scalar observable and inherits the Markov chain central limit theorem directly. The difference between the two sequences is times the first-order remainder of at , and it vanishes in probability because the scaled deviations are tight — uniform ergodicity gives a second-moment bound uniform in — so a Slutsky-type argument transfers the limit law. This is Step 2 of Appendix EC.3.1 combined with Lemma EC.5.
import Definitions.Def_MarkovAsymptoticVariance import Definitions.Def_MarkovErgodicity import Definitions.Def_MarkovChainPathMeasure open MeasureTheory ProbabilityTheory Filter open scoped NNReal ENNReal Topology
theorem MarkovChainCLT.delta_method_of_uniformly_ergodic_of_measurable {X : Type*}
[MeasurableSpace X]
(P : Kernel X X) [IsMarkovKernel P] (π : Measure X) [IsProbabilityMeasure π]
(hP : MarkovChainCLT.HarrisErgodic P π)
(huni : MarkovChainCLT.UniformlyErgodic P π)
{ι : Type*} [Fintype ι]
(u : X → ι → ℝ) (hu : ∀ p, Measurable (fun x => u x p))
(hL2 : ∀ p, MemLp (fun x => u x p) 2 π)
(g : (ι → ℝ) → ℝ) (hgm : Measurable g) (g' : (ι → ℝ) →L[ℝ] ℝ)
(hg : HasFDerivAt g g' (fun p => ∫ x, u x p ∂π)) :
TendstoInDistribution
(fun (n : ℕ) (ω : ℕ → X) =>
Real.sqrt n * (g (fun p => MarkovChainCLT.sampleAvg (fun x => u x p) n ω)
- g (fun p => ∫ x, u x p ∂π)))
atTop (id : ℝ → ℝ) (fun _ => MarkovChainCLT.chainMeasure P π)
(gaussianReal 0
(MarkovChainCLT.asymptoticVariance P π (fun x => g' (u x))).toNNReal) := by sorry