A Robust Spectral Clustering Algorithm for Sub-Gaussian Mixture Models with Outliers 2: Spectral Rounding Without Outliers Gives Σ|S_k|/n_k ≤ 64(2+ε)‖X⁰ − X̂‖₁/n²minResearch Paper
Motivation
Spectral clustering recovers group labels from a similarity matrix: compute its leading eigenvectors, embed each point as a row of the eigenvector matrix, and cluster the rows with k-means. It is a standard last step of convex-relaxation methods for clustering, which first estimate a clustering matrix (with when and belong to the same cluster) and then have to turn the estimate back into labels.
Srivastava, Sarkar and Hanasusanto (arXiv:1912.07546v3; Operations Research 71(1), 2023) estimate by a linear program on a thresholded Gaussian kernel and then round it this way (Algorithm 1, steps 3–4). Their Theorem 2 controls the error of the rounding step in terms of the error of the estimate. The argument follows Lei and Rinaldo (Ann. Statist. 43(1), 2015), who analysed spectral clustering for stochastic block models with a polynomial-time -approximate k-means step, combined with the Davis–Kahan theorem in the form of Yu, Wang and Samworth (Biometrika 102(2), 2015). The same pattern, "matrix estimation error ⇒ eigenvector error ⇒ misclassification count", recurs throughout the analysis of SDP and spectral methods for community detection and mixture models.
Setting
There are points and clusters. A fixed labelling gives clusters of sizes , and . The true clustering matrix is if and otherwise. Its scaled membership matrix is
The estimate is any symmetric matrix with entries in . Spectral rounding takes a matrix of orthonormal eigenvectors of for its largest eigenvalues (ties broken arbitrarily), and runs k-means on the rows of . A membership matrix is a -matrix with exactly one per row; a pair with is a -approximate k-means solution on the rows of if
Let be an optimal rotation, an orthogonal matrix minimising . With and the rows of and , the misclassified set of cluster is
Norms: (entrywise), .
Formalization targets
Goal: Theorem 2, (15), without outliers
For every , , and optimal as above and every ,
Milestones
- Davis–Kahan (Theorem 3, case ): an eigengap of a symmetric below its top eigenvalues gives an orthogonal with .
- Eigenstructure of : consists of top- eigenvectors of with eigenvalues , and vanishes on the orthogonal complement, so the eigengap is .
- (G.6)(iii): for matrices with entries in .
- (G.6): .
- Lemma 5, (G.4) (Lei–Rinaldo): for and a -approximate k-means solution on any , .
- Centre separation: the rows of satisfy .
Companion: permutation recovery
If , the k-means labels agree with up to one permutation of on every point outside .
Significance
Inequality (15) separates the analysis of the algorithm into an estimation part (how close is to in , Theorem 1 of the paper) and a purely deterministic rounding part. Any estimator of the clustering matrix with entries in — the Robust-LP, the SDP relaxations of Peng–Wei or Amini–Levina, a thresholded kernel — inherits a misclassification bound from its error through this one statement, with an explicit constant. The companion statement shows that the sets really are the misclassified points.
The result is proved on paper; none of it is formalized. Mathlib has the spectral theorem for Hermitian matrices but no eigenspace perturbation bound of Davis–Kahan type, and the platform has no k-means approximation lemma for membership matrices. A machine-checked Davis–Kahan theorem in Frobenius norm and a machine-checked Lei–Rinaldo lemma are reusable well beyond this paper: they are the two standard ingredients of every spectral clustering guarantee for block models and mixtures.
Difficulty
The k-means part is elementary but bookkeeping-heavy: the comparison point in the bound must be a feasible k-means solution, which is why the true embedding has to factor as a membership matrix times a centre matrix. The main obstacle is the Davis–Kahan step. The naive route through individual eigenvectors fails when eigenvalues of or of are repeated (here has eigenvalue with multiplicity whenever cluster sizes coincide), and only the subspace is stable; the bound must be stated up to an orthogonal matrix and proved for the whole eigenspace, with no ordering or separation assumption on the eigenvalues of .
Formalization scope
Points are Fin n, clusters Fin r, matrices Matrix ι κ ℝ; the norms are explicit sums, not Mathlib's matrix norm instances. The labelling is fixed; membership matrices are encoded by assignments , and "" by " the cost of every assignment and every centre matrix". "Top- eigenvectors" is variational: orthonormal columns, eigenvector equations, and every Rayleigh quotient on the orthogonal complement at most each selected eigenvalue. is a finite minimum and requires .
Departures from the printed text, each disclosed in the item concerned:
- No outliers (). The printed Theorem 2 allows outliers; it is false for , because the outlier rows of are zero, so the true embedding is not k-means feasible and Lemma 5 does not apply. A counterexample with two clusters of size , one isolated outlier and violates (15).
- No probability, no , no separation condition. (15) is deterministic in ; the printed probabilistic hypotheses only feed Theorem 1 into the permutation claim.
- The printed is replaced by the weaker .
- Davis–Kahan is stated for the leading eigenspace (), the only case used.
- In Lemma 5 the thresholds are required to satisfy , which the page leaves implicit; for a negative the bound (G.4) fails at .
The goal quantifies over every optimal rotation , not over some orthogonal matrix chosen afterwards (which could make empty), and over every top-eigenvector matrix and every approximate k-means solution; a formalization that fixes to good eigenvectors, compares only against the true labels, or drops the eigenvalue-dominance clause (allowing bottom eigenvectors) proves a different and partly false statement.
Contributions welcome: Davis–Kahan in Frobenius norm for symmetric real matrices (general included), the Lei–Rinaldo lemma in full, and the eigenstructure of block-constant matrices.
Selected references
- P. R. Srivastava, P. Sarkar, G. A. Hanasusanto, A Robust Spectral Clustering Algorithm for Sub-Gaussian Mixture Models with Outliers, arXiv:1912.07546v3, 2021; Operations Research 71(1), 2023. https://arxiv.org/abs/1912.07546
- J. Lei, A. Rinaldo, Consistency of spectral clustering in stochastic block models, Annals of Statistics 43(1), 215–237, 2015. https://doi.org/10.1214/14-AOS1274
- Y. Yu, T. Wang, R. J. Samworth, A useful variant of the Davis–Kahan theorem for statisticians, Biometrika 102(2), 315–323, 2015. https://doi.org/10.1093/biomet/asv008
- C. Davis, W. M. Kahan, The rotation of eigenvectors by a perturbation. III, SIAM Journal on Numerical Analysis 7(1), 1–46, 1970. https://doi.org/10.1137/0707001