Section 3 —
ProvedDantzigSelector.Oracle.gaussian_max_tailconcentrationgaussianp2o-batch-pfp1ap2o-gran-per-chapterp2o-plan-paperp2o-v1probability
Let have unit-normed columns , , and let be a vector of independent standard normal random variables on a probability space. Put , which is , and . Then for every ,
With this bounds the probability that the noise violates the orthogonality condition (3.1) for all ; it is the only probabilistic input to Theorems 1.1 and 1.2.
Formalization Note Section 3 normalizes , so the noise coordinates have law and are mutually independent. The event is written "some has ", which is the same event since there are finitely many ; the probability is the measure of that set.
Preamble
import Mathlib import Definitions.Def_CandesTao_Decoding_Norms import Definitions.Def_CandesTao_Decoding_RestrictedIsometry import Definitions.Def_DantzigSelector_Sparse_Model import Definitions.Def_DantzigSelector_Oracle_Model open MeasureTheory ProbabilityTheory CandesTao.Decoding DantzigSelector.Sparse
Formal statement
namespace DantzigSelector.Oracle
/-- Candès–Tao (2007), Section 3, p. 15 (σ = 1): if the columns of `X` are unit-normed and
`z_1, …, z_n` are independent `N(0,1)` variables, then `Z_j := ⟨z, X_j⟩` obeys, for every
`u > 0`, `P(sup_j |Z_j| > u) ≤ 2p · φ(u)/u` with `φ(u) = (2π)^{-1/2} e^{-u²/2}`. -/
theorem gaussian_max_tail {n p : ℕ} {Ω : Type*} [MeasurableSpace Ω] (P : Measure Ω)
[IsProbabilityMeasure P] (X : Matrix (Fin n) (Fin p) ℝ) (hX : UnitNormColumns X)
(z : Fin n → Ω → ℝ) (hz : ∀ i, HasLaw (z i) (gaussianReal 0 1) P) (hind : iIndepFun z P)
(u : ℝ) (hu : 0 < u) :
P {ω | ∃ j : Fin p, u < |∑ i, X i j * z i ω|} ≤
ENNReal.ofReal
(2 * (p : ℝ) * ((Real.sqrt (2 * Real.pi))⁻¹ * Real.exp (-u ^ 2 / 2)) / u) := by sorry
end DantzigSelector.Oracle
Source
Candès & Tao, The Dantzig Selector: Statistical Estimation When p Is Much Larger than n, arXiv:math/0506081v3, p. 15, Section 3 (sentence after Eq. (3.1))
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.