Algorithms for Sparse Support Vector Machines

Abstract

Many problems in classification involve huge numbers of irrelevant features. Variable selection reveals the crucial features, reduces the dimensionality of feature space, and improves model interpretation. In the support vector machine literature, variable selection is achieved by ${\ell }_1$ penalties. These convex relaxations seriously bias parameter estimates toward 0 and tend to admit too many irrelevant features. The current paper presents an alternative that replaces penalties by sparse-set constraints. Penalties still appear, but serve a different purpose. The proximal distance principle takes a loss function $L\left(\beta \right)$ and adds the penalty $\frac{\rho }{2}\text{dist}{\left(\beta ,S_k\right)}^2$ capturing the squared Euclidean distance of the parameter vector $\beta$ to the sparsity set $S_k$ where at most k components of $\beta$ are nonzero. If ${\beta }_{\rho }$ represents the minimum of the objective $f_{\rho }\left(\beta \right)=L\left(\beta \right)+\frac{\rho }{2}\text{dist}{\left(\beta ,S_k\right)}^2$, then ${\beta }_{\rho }$ tends to the constrained minimum of $L\left(\beta \right)$ over $S_k$ as $\rho$ tends to $\infty$. We derive two closely related algorithms to carry out this strategy. Our simulated and real examples vividly demonstrate how the algorithms achieve better sparsity without loss of classification power.

Keywords: Julia; discriminant analysis; sparsity; unsupervised learning.