U-statistics
DefinitionU-statistic#✓fully proved
The U-statistic with kernel h of degree m on a sample x of size n: the average of
h (x ∘ f) over the n.descFactorial m injective maps f : Fin m ↪ Fin n. For m > n there
are no such maps and the value is 0 (Lean's a / 0 = 0).
noncomputable def uStatistic {m n : ℕ} (h : (Fin m → α) → ℝ) (x : Fin n → α) : ℝ :=
(∑ f : Fin m ↪ Fin n, h (x ∘ f)) / n.descFactorial m
Used by 3
DefinitionHoeffding projection of a kernel#✓fully proved
Hoeffding's c-th projection of a kernel h of degree m under μ:
h_c(y) = E h(y₁, …, y_c, X_{c+1}, …, X_m) with X_{c+1}, …, X_m i.i.d. μ.
noncomputable def hoeffdingProj (μ : Measure α) {m : ℕ} (h : (Fin m → α) → ℝ) (c : ℕ)
(y : Fin c → α) : ℝ :=
∫ z : Fin (m - c) → α,
h (fun j => if hj : (j : ℕ) < c then y ⟨j, hj⟩ else z ⟨j - c, by have := j.isLt; omega⟩)
∂(Measure.pi fun _ => μ)
Used by 2
TheoremUnbiasedness#✓fully proved
For an i.i.d. sample of size n ≥ m and an integrable kernel, the
U-statistic has mean E h(X₁, …, X_m).
theorem integral_uStatistic {m n : ℕ} (hmn : m ≤ n) (μ : Measure α) [IsProbabilityMeasure μ]
{h : (Fin m → α) → ℝ} (hh : Integrable h (Measure.pi fun _ => μ)) :
∫ x, uStatistic h x ∂(Measure.pi fun _ : Fin n => μ) =
∫ y, h y ∂(Measure.pi fun _ => μ)
noncomputable def uStatistic {m n : ℕ} (h : (Fin m → α) → ℝ) (x : Fin n → α) : ℝ :=
(∑ f : Fin m ↪ Fin n, h (x ∘ f)) / n.descFactorial m
MathlibExt/Probability/Statistics/UStatistic.lean:70
Uses 1
TheoremHoeffding's variance formula#✓fully proved
Hoeffding's variance formula (1948, §5): for an i.i.d. sample of size n ≥ m
and a symmetric square-integrable kernel h of degree m,
Var U = C(n, m)⁻¹ ∑_{c=1}^{m} C(m, c) C(n - m, m - c) Var h_c.
theorem variance_uStatistic {m n : ℕ} (hmn : m ≤ n) (μ : Measure α) [IsProbabilityMeasure μ]
{h : (Fin m → α) → ℝ} (hsymm : ∀ (σ : Equiv.Perm (Fin m)) y, h (y ∘ σ) = h y)
(hh : MemLp h 2 (Measure.pi fun _ => μ)) :
Var[uStatistic h; Measure.pi fun _ : Fin n => μ] =
(∑ c ∈ Finset.Icc 1 m, (m.choose c : ℝ) * ((n - m).choose (m - c) : ℝ) *
Var[hoeffdingProj μ h c; Measure.pi fun _ : Fin c => μ]) / n.choose m
noncomputable def uStatistic {m n : ℕ} (h : (Fin m → α) → ℝ) (x : Fin n → α) : ℝ :=
(∑ f : Fin m ↪ Fin n, h (x ∘ f)) / n.descFactorial m
MathlibExt/Probability/Statistics/UStatistic.lean:70
noncomputable def hoeffdingProj (μ : Measure α) {m : ℕ} (h : (Fin m → α) → ℝ) (c : ℕ)
(y : Fin c → α) : ℝ :=
∫ z : Fin (m - c) → α,
h (fun j => if hj : (j : ℕ) < c then y ⟨j, hj⟩ else z ⟨j - c, by have := j.isLt; omega⟩)
∂(Measure.pi fun _ => μ)
MathlibExt/Probability/Statistics/UStatistic.lean:75
Uses 2
Loading dependency explorer…
The dependency explorer needs JavaScript and its graph data.