Theorem 3 (Davis–Kahan, Yu–Wang–Samworth), p. 46, case s = 1 — ∃ orthogonal Ô, ‖U − ÛÔ‖_F ≤ 2^{3/2}‖Σ̂ − Σ‖_F / g
ProvedRobustSpecClust.Round.davis_kahanThis is the Davis–Kahan theorem in the Frobenius-norm form of Yu, Wang and Samworth, for the leading eigenspace.
Let be symmetric. Let have orthonormal columns that are eigenvectors of , , and suppose there is an eigengap below them: for every orthogonal to the columns of and every ,
i.e. (in particular the are the largest eigenvalues). Let consist of orthonormal top- eigenvectors of . Then there is an orthogonal matrix with
In the proof of Theorem 2 this bounds the distance between the top eigenvectors of the estimated clustering matrix and those of the true clustering matrix , step (i) of (G.6).
Formalization Note The paper states the theorem for eigenvectors with gap ; this item is the case (where ), the only case the paper uses. Any with is allowed; taking equal to the gap gives (G.5). When the complement is and the gap hypothesis is vacuous, as on the page.
import Mathlib import Definitions.Def_RobustSpecClust_Round_Setting
namespace RobustSpecClust.Round
open Matrix
theorem davis_kahan {n r : ℕ} (S Sh : Matrix (Fin n) (Fin n) ℝ) (hS : Sᵀ = S) (hSh : Shᵀ = Sh)
(U Uh : Matrix (Fin n) (Fin r) ℝ) (lam : Fin r → ℝ) (hU : IsTopEigvecsWith S U lam)
(g : ℝ) (hg : 0 < g)
(hgap : ∀ v : Fin n → ℝ, Uᵀ *ᵥ v = 0 → ∀ k, v ⬝ᵥ (S *ᵥ v) ≤ (lam k - g) * (v ⬝ᵥ v))
(hUh : IsTopEigvecs Sh Uh) :
∃ O ∈ Matrix.orthogonalGroup (Fin r) ℝ,
Real.sqrt (frobSq (U - Uh * O)) ≤ 2 ^ (3 / 2 : ℝ) * Real.sqrt (frobSq (Sh - S)) / g := by sorry
end RobustSpecClust.Round
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.