Sequential minimal optimization

From Infogalactic: the planetary knowledge core
Jump to: navigation, search
Sequential minimal optimization
Class Optimization algorithm for training support vector machines
Worst case performance O(n³)

Sequential minimal optimization (SMO) is an algorithm for solving the quadratic programming (QP) problem that arises during the training of support vector machines. It was invented by John Platt in 1998 at Microsoft Research.[1] SMO is widely used for training support vector machines and is implemented by the popular LIBSVM tool.[2][3] The publication of the SMO algorithm in 1998 has generated a lot of excitement in the SVM community, as previously available methods for SVM training were much more complex and required expensive third-party QP solvers.[4]

Optimization problem

Consider a binary classification problem with a dataset (x1, y1), ..., (xn, yn), where xi is an input vector and yi ∈ {-1, +1} is a binary label corresponding to it. A soft-margin support vector machine is trained by solving a quadratic programming problem, which is expressed in the dual form as follows:

\max_{\alpha} \sum_{i=1}^n \alpha_i - \frac12 \sum_{i=1}^n \sum_{j=1}^n y_i y_j K(x_i, x_j) \alpha_i \alpha_j,
subject to:
0 \leq \alpha_i \leq C, \quad \mbox{ for } i=1, 2, \ldots, n,
\sum_{i=1}^n y_i \alpha_i = 0

where C is an SVM hyperparameter and K(xi, xj) is the kernel function, both supplied by the user; and the variables \alpha_i are Lagrange multipliers.


SMO is an iterative algorithm for solving the optimization problem described above. SMO breaks this problem into a series of smallest possible sub-problems, which are then solved analytically. Because of the linear equality constraint involving the Lagrange multipliers \alpha_i, the smallest possible problem involves two such multipliers. Then, for any two multipliers \alpha_1 and \alpha_2, the constraints are reduced to:

0 \leq \alpha_1, \alpha_2 \leq C,
y_1 \alpha_1 + y_2 \alpha_2 = k,

and this reduced problem can be solved analytically: one needs to find a minimum of a one-dimensional quadratic function. k is the negative of the sum over the rest of terms in the equality constraint, which is fixed in each iteration.

The algorithm proceeds as follows:

  1. Find a Lagrange multiplier \alpha_1 that violates the Karush–Kuhn–Tucker (KKT) conditions for the optimization problem.
  2. Pick a second multiplier \alpha_2 and optimize the pair (\alpha_1,\alpha_2).
  3. Repeat steps 1 and 2 until convergence.

When all the Lagrange multipliers satisfy the KKT conditions (within a user-defined tolerance), the problem has been solved. Although this algorithm is guaranteed to converge, heuristics are used to choose the pair of multipliers so as to accelerate the rate of convergence. This is critical for large data sets since there are n(n − 1) possible choices for \alpha_i and \alpha_j.

Related Work

The first approach to splitting large SVM learning problems into a series of smaller optimization tasks was proposed by Bernhard E Boser, Isabelle M Guyon, Vladimir Vapnik.[5] It is known as the "chunking algorithm". The algorithm starts with a random subset of the data, solves this problem, and iteratively adds examples which violate the optimality conditions. One disadvantage of this algorithm is that it is necessary to solve QP-problems scaling with the number of SVs. On real world sparse data sets, SMO can be more than 1000 times faster than the chunking algorithm.[1]

In 1997, E. Osuna, R. Freund, and F. Girosi proved a theorem which suggests a whole new set of QP algorithms for SVMs.[6] By the virtue of this theorem a large QP problem can be broken down into a series of smaller QP sub-problems. A sequence of QP sub-problems that always add at least one violator of the Karush–Kuhn–Tucker (KKT) conditions is guaranteed to converge. The chunking algorithm obeys the conditions of the theorem, and hence will converge.[1] The SMO algorithm can be considered a special case of the Osuna algorithm, where the size of the optimization is two and both Lagrange multipliers are replaced at every step with new multipliers that are chosen via good heuristics.[1]

The SMO algorithm is closely related to a family of optimization algorithms called Bregman methods or row-action methods. These methods solve convex programming problems with linear constraints. They are iterative methods where each step projects the current primal point onto each constraint.[1]

See also


  1. 1.0 1.1 1.2 1.3 1.4 Platt, John (1998), Sequential Minimal Optimization: A Fast Algorithm for Training Support Vector Machines, CiteSeerX:<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  2. Chang, Chih-Chung; Lin, Chih-Jen (2011). "LIBSVM: A library for support vector machines". ACM Transactions on Intelligent Systems and Technology. 2 (3).<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  3. Luca Zanni (2006). Parallel Software for Training Large Scale Support Vector Machines on Multiprocessor Systems.
  4. Rifkin, Ryan (2002), "Everything Old is New Again: a Fresh Look at Historical Approaches in Machine Learning", Ph.D. thesis: 18<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  5. Boser, B. E.; Guyon, I. M.; Vapnik, V. N. (1992). "A training algorithm for optimal margin classifiers". Proceedings of the fifth annual workshop on Computational learning theory - COLT '92. p. 144. doi:10.1145/130385.130401. ISBN 089791497X.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  6. Osuna, E.; Freund, R.; Girosi, F. (1997). "An improved training algorithm for support vector machines". Neural Networks for Signal Processing [1997] VII. Proceedings of the 1997 IEEE Workshop. pp. 276–285. doi:10.1109/NNSP.1997.622408. ISBN 0-7803-4256-9.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>