Documentation

TauCeti.Probability.Distributions.Gaussian.Concentration

Gaussian concentration of smooth functions #

Let μ be a centred Gaussian measure on a separable Banach space E, and let f : E → ℝ be a C¹ function whose derivative has variance at most v under μ at almost every point, that is Var[fderiv ℝ f x; μ] ≤ v. Then f concentrates around its mean with Gaussian tails: f - 𝔼 f is sub-Gaussian with variance proxy π² v / 4, and

μ {x | ε ≤ |f x - 𝔼 f|} ≤ 2 exp (-2 ε² / (π² v)).

For the standard Gaussian measure on a finite-dimensional inner product space the variance of a linear functional is its squared norm, so the hypothesis holds with v = K² when ‖f' x‖ ≤ K, that is, when the C¹ function f is K-Lipschitz. Estimates of this kind, transported to the Euclidean sphere, are the concentration of measure behind Dvoretzky's theorem on almost Euclidean sections of normed spaces.

The proof is the interpolation argument of Maurey and Pisier. For (x, y) ∈ E × E and the rotation R θ (x, y) = (cos θ • x + sin θ • y, -sin θ • x + cos θ • y), the first coordinate of R θ (x, y) runs along a quarter circle from x to y with velocity the second coordinate, so (ContDiff.sub_eq_integral_fderiv_rotation)

f y - f x = ∫ θ in 0..π/2, f' (R θ (x, y)).1 (R θ (x, y)).2.

Averaging over θ and applying Jensen's inequality to the exponential bounds exp (t (f y - f x)) by an average of exp (t π / 2 · f' (R θ (x, y)).1 (R θ (x, y)).2). The product μ ⊗ μ is invariant under each rotation R θ (ProbabilityTheory.IsGaussian.map_rotation_eq_self), so after integration every θ contributes the same amount, which is an exponential moment of the Gaussian variable f' x y in y. Jensen's inequality in y then gives the exponential-moment estimate 𝔼 exp (t (f - 𝔼 f)) ≤ exp (π² t² v / 8), and Chernoff's bound (ProbabilityTheory.HasSubgaussianMGF.measure_ge_le) gives the tails.

Main results #

References #

theorem ProbabilityTheory.IsGaussian.lintegral_exp_mul_sub_le {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] [CompleteSpace E] [SecondCountableTopology E] [MeasurableSpace E] [BorelSpace E] {μ : MeasureTheory.Measure E} [IsGaussian μ] {f : E → ℝ} (hμ : ∫ (x : E), id x ∂μ = 0) (hf : ContDiff ℝ 1 f) (t : ℝ) :
∫⁻ (p : E × E), ENNReal.ofReal (Real.exp (t * (f p.2 - f p.1))) ∂μ.prod μ ≤ ∫⁻ (x : E), ENNReal.ofReal (Real.exp (Real.pi ^ 2 * t ^ 2 / 8 * variance (⇑(fderiv ℝ f x)) μ)) ∂μ

The Maurey–Pisier inequality. For a centred Gaussian measure μ and a C¹ function f, ∫∫ exp (t (f y - f x)) dμ(x) dμ(y) ≤ ∫ exp (π² t² / 8 · Var[f' x; μ]) dμ(x).

theorem ProbabilityTheory.IsGaussian.lintegral_exp_mul_sub_le_of_variance_fderiv_le {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] [CompleteSpace E] [SecondCountableTopology E] [MeasurableSpace E] [BorelSpace E] {μ : MeasureTheory.Measure E} [IsGaussian μ] {f : E → ℝ} {v : NNReal} (hμ : ∫ (x : E), id x ∂μ = 0) (hf : ContDiff ℝ 1 f) (hv : ∀ᵐ (x : E) ∂μ, variance (⇑(fderiv ℝ f x)) μ ≤ ↑v) (t : ℝ) :
∫⁻ (p : E × E), ENNReal.ofReal (Real.exp (t * (f p.2 - f p.1))) ∂μ.prod μ ≤ ENNReal.ofReal (Real.exp (Real.pi ^ 2 * t ^ 2 / 8 * ↑v))

The Maurey–Pisier inequality when the variances of the derivative are at most v almost everywhere.

theorem ProbabilityTheory.IsGaussian.integrable_exp_mul_of_variance_fderiv_le {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] [CompleteSpace E] [SecondCountableTopology E] [MeasurableSpace E] [BorelSpace E] {μ : MeasureTheory.Measure E} [IsGaussian μ] {f : E → ℝ} {v : NNReal} (hμ : ∫ (x : E), id x ∂μ = 0) (hf : ContDiff ℝ 1 f) (hv : ∀ᵐ (x : E) ∂μ, variance (⇑(fderiv ℝ f x)) μ ≤ ↑v) (t : ℝ) :
MeasureTheory.Integrable (fun (x : E) => Real.exp (t * f x)) μ

A C¹ function whose derivative has variance at most v almost everywhere under a centred Gaussian measure has finite exponential moments of every order.

A C¹ function whose derivative has variance at most v almost everywhere under a centred Gaussian measure is integrable.

theorem ProbabilityTheory.IsGaussian.hasSubgaussianMGF_sub_integral_of_variance_fderiv_le {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] [CompleteSpace E] [SecondCountableTopology E] [MeasurableSpace E] [BorelSpace E] {μ : MeasureTheory.Measure E} [IsGaussian μ] {f : E → ℝ} {v : NNReal} (hμ : ∫ (x : E), id x ∂μ = 0) (hf : ContDiff ℝ 1 f) (hv : ∀ᵐ (x : E) ∂μ, variance (⇑(fderiv ℝ f x)) μ ≤ ↑v) :
HasSubgaussianMGF (fun (x : E) => f x - ∫ (y : E), f y ∂μ) (NNReal.pi ^ 2 / 4 * v) μ

Gaussian concentration, exponential-moment form. If f is C¹ and its derivative has variance at most v under a centred Gaussian measure μ at almost every point, then f - 𝔼 f is sub-Gaussian with variance proxy π² v / 4: 𝔼 exp (t (f - 𝔼 f)) ≤ exp (π² t² v / 8).

theorem ProbabilityTheory.IsGaussian.measureReal_setOf_le_abs_sub_integral_le_of_variance_fderiv_le {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] [CompleteSpace E] [SecondCountableTopology E] [MeasurableSpace E] [BorelSpace E] {μ : MeasureTheory.Measure E} [IsGaussian μ] {f : E → ℝ} {v : NNReal} (hμ : ∫ (x : E), id x ∂μ = 0) (hf : ContDiff ℝ 1 f) (hv : ∀ᵐ (x : E) ∂μ, variance (⇑(fderiv ℝ f x)) μ ≤ ↑v) {ε : ℝ} (hε : 0 ≤ ε) :
μ.real {x : E | ε ≤ |f x - ∫ (y : E), f y ∂μ|} ≤ 2 * Real.exp (-2 * ε ^ 2 / (Real.pi ^ 2 * ↑v))

Gaussian concentration. If f is C¹ and its derivative has variance at most v under a centred Gaussian measure μ at almost every point, then for ε ≥ 0, μ {x | ε ≤ |f x - 𝔼 f|} ≤ 2 exp (-2 ε² / (π² v)).

Under the standard Gaussian measure, a C¹ function f with ‖f' x‖ ≤ K everywhere (that is, a C¹ function which is K-Lipschitz) is sub-Gaussian about its mean with variance proxy π² K² / 4.

Gaussian concentration for the standard Gaussian measure. If f is C¹ with ‖f' x‖ ≤ K everywhere, then for ε ≥ 0, γ {x | ε ≤ |f x - 𝔼 f|} ≤ 2 exp (-2 ε² / (π² K²)) for the standard Gaussian measure γ.