Multiclass Learning with Simplex Coding
Abstract
In this paper we discuss a novel framework for multiclass learning, defined by a suitable coding/decoding strategy, namely the simplex coding, that allows to generalize to multiple classes a relaxation approach commonly used in binary classification. In this framework, a relaxation error analysis can be developed avoiding constraints on the considered hypotheses class. Moreover, we show that in this setting it is possible to derive the first provably consistent regularized method with training/tuning complexity which is independent to the number of classes. Tools from convex analysis are introduced that can be used beyond the scope of this paper.
1 Introduction
As bigger and more complex datasets are available, multiclass learning is becoming increasingly important in machine learning.
While theory and algorithms for solving binary classification problems are well established, the problem of multicategory classification is much less understood. Practical multiclass algorithms often reduce the problem to a collection of binary classification problems. Binary classification algorithms are often based
on a relaxation approach: classification is posed as a nonconvex minimization problem and hence relaxed to a convex one, defined by suitable convex loss functions. In this context, results in statistical learning theory quantify the error incurred by relaxation and in particular derive comparison inequalities explicitly relating the excess misclassification risk with the excess expected loss, see for example [2, 27, 14, 29] and [18] Chapter 3 for an exhaustive presentation as well as generalizations.
Generalizing the above approach and results to more than two classes is not straightforward.
Over the years, several computational solutions have been proposed (among others, see [10, 6, 5, 25, 1, 21]. Indeed, most of the above methods can be interpreted as a kind of relaxation. Most proposed
methods have complexity which is more than linear in the number of classes and simple onevs all in practice offers
a good alternative both in terms of performance and speed [15].
Much fewer works have focused on deriving theoretical guarantees. Results in this sense have been pioneered by [28, 20], see also [11, 7, 23].
In these works the error due to relaxation is studied asymptotically and under constraints on the function class to be considered.
More quantitative results in terms of comparison inequalities are given in [4] under similar restrictions
(see also [19]).
Notably, the above results show that seemigly intuitive extensions of binary classification
algorithms might lead to methods which are not consistent.
Further, it is interesting to note that these restrictions on the function class, needed to prove the theoretical guarantees, make the computations in the corresponding algorithms more involved and are in fact often ignored in practice.
In this paper we dicuss a novel framework for multiclass learning,
defined by a suitable coding/decoding strategy, namely the simplex coding,
in which a relaxation error analysis can be developed avoiding constraints on the considered hypotheses class.
Moreover, we show that in this framework it is possible to derive the first provably consistent
regularized method with training/tuning complexity which is independent to the number of classes.
Interestingly, using the simplex coding, we can naturally generalize results, proof techniques and methods from the binary case, which is recovered as a special case of our theory. Due to space restriction in this paper we focus on extensions of least squares, and SVM loss functions, but our analysis can be generalized to large class of simplex loss functions, including extension of logistic and exponential loss functions (used in boosting). Tools from convex analysis are developed in the longer version of the paper and can be useful beyond the scopes of this paper, and in particular in structured prediction.
The rest of the paper is organized as follow. In Section 2 we discuss problem statement and background. In Section 3 we discuss the simplex coding framework that we analyze in Section 4. Algorithmic aspects and numerical experiments are discussed in Section 5 and Section 6, respectively. Proofs and supplementary technical results are given in the longer version of the paper.
2 Problem Statement and Previous Work
Let be two random variables with values in two measurable spaces and , . Denote by , the law of on , and by , the conditional probabilities for . The data is a sample , from identical and independent copies of .We can think of as a set of possible inputs and of , as a set of labels describing a set of semantic categories/classes the input can belong to. A classification rule is a map , and its error is measured by the misclassification risk
The optimal classification rule that minimizes is the Bayes rule,
Computing the Bayes rule by directly minimizing the risk , is not possible since the probability distribution is unknown.
In fact one could think of minimizing the empirical risk (ERM), ,
which is an unbiased estimator of the , but the corresponding
optimization problem is in general not feasible.
In binary classification, one of the most common way to obtain computationally efficient methods
is based on a relaxation approach. We recall this approach in the next section and describe its extension to multiclass in the rest of the paper.
Relaxation Approach to Binary Classification. If , we can set .
Most modern machine learning algorithms for binary classification consider a convex relaxation of the ERM functional . More precisely:
1) the indicator function in is replaced by non negative loss which is convex in the second argument and is sometimes called a surrogate loss;
2) the classification rule replaced by a real valued measurable function .
A classification rule is then obtained by considering the sign of .
It often suffices to consider a special class of loss functions, namely
large margin loss functions of the form .
This last expression is suggested by the observation that the misclassification risk, using the labels , can be written as
where is the heavy side step function. The quantity , sometimes called the margin, is a natural pointwise measure of the classification error. Among other examples of large margin loss functions (such as the logistic and exponential loss),
we recall the hinge loss used in support vector machine, and the square loss used in regularized least squares (note that ).
Using surrogate large margin loss functions it is possible to design effective learning algorithms replacing the empirical risk with regularized empirical risk minimization
(1) 
where is a suitable regularization functional and is the regularization parameter, see Section 5.
2.1 Relaxation Error Analysis
As we replace the misclassification loss with a convex surrogate— loss, we are effectively changing the problem:
the misclassification risk is replaced by the expected loss,
. The expected loss can be seen as a functional on a large space of functions
, which depend on and . Its minimizer, denoted by , replaces the Bayes rule as the target of our algorithm.
The question arises of the price we pay by a considering a relaxation approach:
“What is the relationship between and ?” More generally, “What is the approximation we incur into
by estimating the expected risk rather than the misclassification risk?”
The relaxation error for a given loss function can be quantified by the following two requirements:
1) Fisher Consistency. A loss function is Fisher consistent if almost surely (this property is related to the notion of classificationcalibration [2]).
2) Comparison inequalities. The excess misclassification risk, and the excess expected loss are related by a comparison inequality
for any function , where is a suitable function that depends on , and possibly on the data distribution.
In particular should be such that as , so that if is a (possibly random) sequence of functions, such that (possibly in probability), then the corresponding sequences of classification rules is Bayes consistent, i.e. (possibly in probability). If is explicitly known, then bounds on the excess expected loss yields bounds on the excess misclassification risk.
The relaxation error in the binary case has been thoroughly studied in [2, 14].
In particular, Theorem 2 in [2] shows that if a large margin surrogate
loss is convex, differentiable and decreasing in a neighborhood of , then the loss is Fisher consistent.
Moreover, in this case it is possible to give an explicit expression of the function .
In particular, for the hinge loss the target function is exactly the Bayes rule and . For least squares, , and .
The comparison inequality for the square loss can be improved for a suitable class of probability distribution satisfying the so called Tsybakov noise condition [22],
Under this condition the probability of points such that decreases polynomially.
In this case the comparison inequality for the square loss is given by , see [2, 27].
Previous Works in Multiclass Classification.
From a practical perspective, over the years, several computational solutions to multiclass learning have been proposed.
Among others, we mention for example [10, 6, 5, 25, 1, 21].
Indeed, most of the above methods can be interpreted as a kind of relaxation of the original multiclass problem.
Interestingly, the study in [15] suggests that the simple onevs all schemes should be a practical benchmark for multiclass algorithms
as it seems to experimentally achive performances that are similar or better to more sophisticated methods.
As we previously mentioned from a theoretical perspective a general account of a large class of multiclass methods
has been given in [20], building on results in [2] and [28].
Notably, these results show that seemingly intuitive extensions of binary classification algorithms might lead to inconsistent methods.
These results, see also [11, 23],
are developed in a setting where a classification rule is found by applying a suitable prediction/decoding map to a function
where is found considering a loss function
The considered functions have to satisfy the constraint , for all .
The latter requirement is problematic since it makes the computations in the corresponding algorithms more involved
and is in fact often ignored, so that practical algorithms often come with no consistency guarantees.
In all the above papers relaxation is studied in terms of Fisher and Bayes consistency
and the explicit form of the function is not given. More quantitative results in terms of explicit comparison inequality are given in [4] and (see also [19]), but also need to to impose the ”sum to zero” constraint on the considered function class.
3 A Relaxation Approach to Multicategory Classification
In this section we propose a natural extension of the relaxation approach that avoids constraining the class of functions to be considered, and allows to derive explicit comparison inequalities. See Remark 1 for related approaches.
Simplex Coding.
We start considering a suitable coding/decoding strategy. A coding map turns a label into a code vector. The corresponding decoding map given a vector returns a label in . Note that, this is what we implicitly did while treating binary classification encoding the label space using the coding , so that the naturally decoding strategy is simply . The coding/decoding strategy we study is described by the following definition.
Definition 1 (Simplex Coding).
The simplex coding is a map where the code vectors satisfy: 1) , , 2) for with , and 3) . The corresponding decoding is the map
The simplex coding corresponds to the most separated vectors on the
hypersphere in , that is the vertices of the simplex (see Figure 1).
For binary classification it reduces to the coding and the decoding map is equivalent to taking the sign of .
The decoding map has a natural geometric interpretation: an input point is mapped to a vector by a function ,
and hence assigned to the class having closer code
vector (for and , we have .
Relaxation for Multiclass Learning.
We use the simplex coding to propose an extension of the binary classification approach.
Following the binary case, the relaxation can be described in two steps:

using the simplex coding, the indicator function is upper bounded by a nonnegative loss function , such that for all , and ,

rather than we consider functions with values in , so that , for all and .
In the next section we discuss several loss functions satisfying the above definitions and we study in particular the extension of the least squares and SVM loss functions.
Multiclass Simplex Loss Functions.
Several loss functions for binary classification can be naturally extended to multiple classes using the simplex coding. Due to space restriction, in this paper we focus on extensions of least squares, and SVM loss functions, but our analysis can be generalized to large class of simplex loss functions, including extension of logistic and exponential loss functions( used in boosting). The Simplex Least Square loss (SLS) is given by ,
and reduces to the usual least square approach to binary classification for
One natural extension of the SVM’s hinge loss in this setting would be to consider the Simplex Half space SVM loss (SHSVM)
.
We will see in the following that while this loss function
would induce efficient algorithms in general is not Fisher consistent unless further constraints are assumed.
In turn, this latter constraint would considerably slow down the computations.
Then we consider a second loss function Simplex Cone SVM (SCSVM), related to the hinge loss, which is defined as
The latter loss function is related to the one considered in the multiclass SVM proposed in [10].
We will see that it is possible to quantify the relaxation error of the loss function without requiring further constraints.
Both the above SVM loss functions reduce to the binary SVM hinge loss if .
Remark 1 (Related approaches).
The simplex coding has been considered in [8],[26], and [16]. In particular, a kind of SVM loss is considered in [8] where and with . More recently [26] considered the loss function , and a simplex multiclass boosting loss was introduced in [16], in our notation While all those losses introduce a certain notion of margin that makes use of the geometry of the simplex coding, it is not to clear how to derive explicit comparison theorems and moreover the computational complexity of the resulting algorithms scales linearly with the number of classes in the case of the losses considered in [16, 26] and for losses considered in [8] .
4 Relaxation Error Analysis
If we consider the simplex coding, a function taking values in , and the decoding operator , the misclassification risk can also be written as: . Then, following a relaxation approach we replace the misclassification loss by the expected risk induced by one of the loss functions defined in the previous section. As in the binary case we consider the expected loss Let ,
The following theorem studies the relaxation error for SHSVM, SCSVM, and SLS loss functions.
Theorem 1.
For SHSVM, SCSVM, and SLS loss functions, there exists a such that is convex and continuous. Moreover,

The minimizer of over exists and .

For any , where the expressions of and are given in Table 1.
Loss  

SHSVM  
SCSVM  
SLS 
The proof of this theorem is given in the longer version of the paper.
The above theorem can be improved for Least Squares under certain classes of distribution .
Toward this end we introduce the following notion of misclassification noise that generalizes Tsybakov’s noise condition.
Definition 2.
Fix , we say that the distribution satisfy the multiclass noise condition with parameter , if
(2) 
where .
If a distribution is characterized by a very large , then, for each , is arbitrarily close to one of the coding vectors. For , the above condition reduces to the binary Tsybakov noise. Indeed, let , and , if , , and if , .
The following result improves the exponent of simplexleast square to :
Theorem 2.
For each , if holds, then for SLS we have the following inequality,
(3) 
for a constant .
Remark 2.
Note that the comparison inequalities show a tradeoff between the exponent and the constant , for SLS and SVM losses. While the constant is order for SVM it is order for SLS, on the other hand the exponent is for SVM losses and for SLS. The latter could be enhanced to for close to separable classification problems by virtue of the Tsybakov noise condition.
Remark 3.
Comparison inequalities given in Theorems 1 and 3 can be used to derive generalization bounds on the excess misclassification risk. For least square minmax sharp bound, for vector valued regression are easy to derive.Standard techniques for deriving sample complexity bound in binary classification extended for multiclass SVM losses could be found in [7] and could be adapted to our setting. The obtained bound are not known to be tight, better bounds akin to those in [18], will be subject to future work.
5 Computational Aspects and Regularization Algorithms
In this section we discuss some computational implications of the framework we presented.
Regularized Kernel Methods.
We consider regularized methods of the form (1), induced by simplex loss functions and where
the hypotheses space is a vector valued reproducing kernel Hilbert spaces (VVRKHSs) and the regularizer the corresponding norm.
See Appendix D. for a brief introduction to VVRKHSs.
In the following, we consider a class of kernels such that the corresponding RKHS
is given by the completion of the span , where we note that the coefficients are vectors in . While other choices are possible this is the kernel more directly related to a one vs all approach.
We will discuss in particular the case where the kernel is induced by a finite dimensional feature map,
, and is the inner product in .
In this case we can write each function in as , where .
It is known [12, 3] that the representer theorem [9] can be easily extended to a vector valued setting, so that
that minimizer of a simplex version of Tikhonov regularization is given by , for all , where the explicit expression of the coefficients depends on the considered loss function.
We use the following notations:
Simplex Regularized Least squares (SRLS).
SRLS is obtained considering the simplex least square loss in the Tikhonov functionals.
It is easy to see [15] that in this case the coefficients must satisfy
either or in the linear case, where and .
Interestingly, the classical results from [24] can be extended to show that
the value , obtained computing the solution removing the point
from the training set (the leave one out solution), can be computed in closed form.
Let .
Let and .
Define , such that: ,
One can show similarly to [15], that , where is the Hadamard product.
Then, the leaveoneout error
can be minimized at essentially no extra cost by precomputing the eigen decomposition of (or ).
Simplex Cone Support Vector Machine (SCSVM). Using standard reasoning it is easy to show that (see Appendix C.),
for the SCSVM the coefficients in the representer theorem are given by
where solve the quadratic programming (QP) problem
(4)  
where and , , for and is the Kronecker delta.
Simplex Halfspaces Support Vector Machine (SHSVM). A similar, yet more more complicated procedure, can be derived for the SHSVM.
Here, we omit this derivation and observe instead that if we neglect the convex hull
constraint from Theorem 1, requiring for almost all ,
then the SHSVM has an especially simple formulation at the price of loosing consistency guarantees.
In fact, in this case the coefficients are given by
where , with solve the quadratic programming (QP) problem
where . The latter formulation could be trained at the same complexity of the binary SVM (worst case ) but lacks consistency.
Online/Incremental Optimization
The regularized estimators induced by the simplex loss functions can be computed
by mean of online/incremental first order (sub) gradient methods.
Indeed, when considering finite dimensional feature maps,
these strategies offer computationally feasible solutions to train estimators
for large datasets where neither a by or an by matrix fit in memory.
Following [17] we can alternate a step of stochastic descent on a data point : and a projection on the Frobenius ball
(See Algorithn C. for details.)
The algorithm depends on the used loss function through the computation
of the (pointwise) subgradient . The latter can be easily computed for all the loss functions previously discussed.
For the SLS loss we have
while for the SCSVM loss we have
where .
For the SHSVM loss we have:
5.1 Comparison of Computational Complexity
The cost of solving SRLS for fixed is in the worst case (for example via Choleski decomposition). If we are interested into computing the regularization path for regularization parameter values, then as noted in [15] it might be convenient to perform an eigendecomposition of the kernel matrix rather than solving the systems times. For explicit feature maps the cost is , so that the cost of computing the regularization path for simplex RLS algorithm is and hence independent of . One can contrast this complexity with the one of a näive One Versus all (OVa) approach that would lead to a complexity. Simplex SVMs can be solved using solvers available for binary SVMs that are considered to have complexity with (actually the complexity scales with the number of support vectors) . For SCSVM, though, we have rather than unknowns and the complexity is . SHSVM where we omit the constraint, could be trained at the same complexity of the binary SVM (worst case ) but lacks consistency. Note that unlike for SRLS, there is no straightforward way to compute the regularization path and the leave one out error for any of the above SVMs . The online algorithms induced by the different simplex loss functions are essentially the same, in particular each iteration depends linearly on the number of classes.
6 Numerical Results
We conduct several experiments to evaluate the performance of our batch and online algorithms, on 5 UCI datasets as listed in Table 2, as well as on Caltech101 and Pubfig83. We compare the performance of our algorithms to on versus all svm (libsvm) , as well as the simplex based boosting [16].
For UCI datasets we use the raw features, on Caltech101 we use hierarchical features

Landsat  Optdigit  Pendigit  Letter  Isolet  Ctech  Pubfig83 

SCSVM Online (ho)  
SHSVM Online (ho)  
SLS Online (ho)  
SLS Batch (loo)  
SLS rbf Batch (loo)  
SVM batch ova (ho)  
SVM rbf batch ova (ho)  
Simplex boosting [16] 
As suggested by the theory, the consistent methods SCSVM and SLS have a big advantage over SHSVM (where we omitted the convex hull constraint) . Batch methods are overall superior to online methods, with online SCSVM achieving the best results. More generally, we see that rbf S LS has the best performance among the simplex methods including the simplex boosting [16]. When compared to One Versus All SVMrbf, we see that SLS rbf achieves essentially the same performance.
Footnotes
 The data set will be made available upon acceptance.
References
 Erin L. Allwein, Robert E. Schapire, and Yoram Singer. Reducing multiclass to binary: a unifying approach for margin classifiers. Journal of Machine Learning Research, 1:113–141, 2000.
 Peter L. Bartlett, Michael I. Jordan, and Jon D. McAuliffe. Convexity, classification, and risk bounds. Journal of the American Statistical Association, 101(473):138–156, 2006.
 A. Caponnetto and E. De Vito. Optimal rates for regularized leastsquares algorithm. Foundations of Computational Mathematics, 2006.
 D. Chen and T. Sun. Consistency of multiclass empirical risk minimization methods based in convex loss. Journal of machine learning, X, 2006.
 Crammer.K and Singer.Y. On the algorithmic implementation of multiclass kernelbased vector machines. JMLR, 2001.
 Thomas G. Dietterich and Ghulum Bakiri. Solving multiclass learning problems via errorcorrecting output codes. Journal of Artificial Intelligence Research, 2:263–286, 1995.
 Yann Guermeur. Vc theory of large margin multicategory classiÞers. Journal of Machine Learning Research, 8:2551–2594, 2007.
 Simon I. Hill and Arnaud Doucet. A framework for kernelbased multicategory classification. J. Artif. Int. Res., 30(1):525–564, December 2007.
 G. Kimeldorf and G. Wahba. A correspondence between bayesian estimation of stochastic processes and smoothing by splines. Ann. Math. Stat., 41:495–502, 1970.
 Lee.Y, L.Yin, and Wahba.G. Multicategory support vector machines: Theory and application to the classification of microarray data and satellite radiance data. Journal of the American Statistical Association, 2004.
 Liu.Y. Fisher consistency of multicategory support vector machines. Eleventh International Conference on Artificial Intelligence and Statistics, 289296, 2007.
 C.A. Micchelli and M. Pontil. On learning vector–valued functions. Neural Computation, 17:177–204, 2005.
 N. Pinto, Z. Stone, T. Zickler, and D.D. Cox. Scalingup biologicallyinspired computer vision: A casestudy on facebook. 2011.
 M.D. Reid and R.C. Williamson. Composite binary losses. JMLR, 11, September 2010.
 Rifkin.R and Klautau.A. In defense of one versus all classification. journal of machine learning, 2004.
 Saberian.M and Vasconcelos .N. Multiclass boosting: Theory and algorithms. In NIPS 2011, 2011.
 Shai ShalevShwartz, Yoram Singer, and Nathan Srebro. Pegasos: Primal estimated subgradient solver for svm. In Proceedings of the 24th ICML, ICML ’07, pages 807–814, New York, NY, USA, 2007. ACM.
 I. Steinwart and A. Christmann. Support vector machines. Information Science and Statistics. Springer, New York, 2008.
 Van de Geer.S Tarigan.B. A moment bound for multicategory support vector machines. JMLR 9, 21712185, 2008.
 A. Tewari and P. L. Bartlett. On the consistency of multiclass classification methods. In Proceedings of the 18th Annual Conference on Learning Theory, volume 3559, pages 143–157. Springer, 2005.
 I. Tsochantaridis, T. Joachims, T. Hofmann, and Y. Altun. Large margin methods for structured and interdependent output variables. JMLR, 6(2):1453–1484, 2005.
 Alexandre B. Tsybakov. Optimal aggregation of classifiers in statistical learning. Annals of Statistics, 32:135–166, 2004.
 Elodie Vernet, Robert C. Williamson, and Mark D. Reid. Composite multiclass losses. In Proceedings of Neural Information Processing Systems (NIPS 2011), 2011.
 G. Wahba. Spline models for observational data, volume 59 of CBMSNSF Regional Conference Series in Applied Mathematics. SIAM, Philadelphia, PA, 1990.
 Weston and Watkins. Support vector machine for multi class pattern recognition. Proceedings of the seventh european symposium on artificial neural networks, 1999.
 Tong Tong Wu and Kenneth Lange. Multicategory vertex discriminant analysis for highdimensional data. Ann. Appl. Stat., 4(4):1698–1721, 2010.
 Y. Yao, L. Rosasco, and A. Caponnetto. On early stopping in gradient descent learning. Constructive Approximation, 26(2):289–315, 2007.
 T. Zhang. Statistical analysis of some multicategory large margin classification methods. Journal of Machine Learning Research, 5:1225–1251, 2004.
 Tong Zhang. Statistical behavior and consistency of classification methods based on convex risk minimization. The Annals of Statistics, Vol. 32, No. 1, 56Ð134, 2004.