### libsdca

Multiclass classification based on stochastic dual coordinate ascent

view repo

In order to push the performance on realistic computer vision tasks, the number of classes in modern benchmark datasets has significantly increased in recent years. This increase in the number of classes comes along with increased ambiguity between the class labels, raising the question if top-1 error is the right performance measure. In this paper, we provide an extensive comparison and evaluation of established multiclass methods comparing their top-k performance both from a practical as well as from a theoretical perspective. Moreover, we introduce novel top-k loss functions as modifications of the softmax and the multiclass SVM losses and provide efficient optimization schemes for them. In the experiments, we compare on various datasets all of the proposed and established methods for top-k error optimization. An interesting insight of this paper is that the softmax loss yields competitive top-k performance for all k simultaneously. For a specific top-k error, our new top-k losses lead typically to further improvements while being faster to train than the softmax.

READ FULL TEXT VIEW PDFMulticlass classification based on stochastic dual coordinate ascent

view repo

The number of classes
is rapidly growing in modern computer vision benchmarks
[47, 62].
Typically, this also leads to
ambiguity in the labels as classes start to overlap.
Even for humans, the error rates in top- performance are often quite high
( on SUN 397 [60]).
While previous research focuses on minimizing the top- error,
we address top- error optimization in this paper.
We are interested in two cases:
a) achieving small top- error for *all* reasonably small ;
and
b) minimization of a specific top- error.

While it is argued in [2] that the one-versus-all (OVA)
SVM scheme performs on par in top- and top- accuracy
with the other SVM variations based on ranking losses,
we have recently shown in [28]
that minimization of the *top- hinge loss*
leads to improvements in top- performance compared to OVA SVM,
multiclass SVM, and other ranking-based formulations.
In this paper,
we study top- error optimization from a wider perspective.
On the one hand,
we compare OVA schemes and direct multiclass losses
in extensive experiments,
and on the other,
we present theoretical discussion
regarding their calibration for the top- error.
Based on these insights,
we suggest new families of loss functions for the top- error.
Two are smoothed versions of the top- hinge losses
[28],
and the other two are top- versions
of the softmax loss.
We discuss their advantages and disadvantages,
and for the convex losses provide an efficient implementation
based on stochastic dual coordinate ascent (SDCA)
[48].

We evaluate a battery of loss functions on

datasets of different tasks ranging from text classification to large scale vision benchmarks, including fine-grained and scene classification. We systematically optimize and report results separately for each top-

accuracy. One interesting message that we would like to highlight is that the softmax loss is able to optimizeRelated work.
Top- optimization has recently received revived attention
with the advent of large scale problems
[20, 28, 30, 31].
The top- error in multiclass classification,
which promotes good ranking of class *labels* for each example,
is closely related to the
precision@ metric in information retrieval,
which counts the fraction of positive instances among the top-
ranked *examples*.
In essence, both approaches enforce a desirable
*ranking* of items [28].

The classic approaches optimize pairwise ranking with
[25, 53],
[11],
and [7].
An alternative direction was proposed by
Usunier *et al*. [54], who described a general family
of convex loss functions for ranking and classification.
One of the loss functions that we consider
( [28])
also falls into that family.
Weston *et al*. [59] then introduced Wsabie,
which optimizes an approximation of a ranking-based loss from
[54].
A Bayesian approach was suggested by [51].

Recent works focus on the top of the ranked list [1, 9, 39, 46], scalability to large datasets [20, 28, 30], explore transductive learning [31] and prediction of tuples [45].

Method | Name | Loss function | Conjugate | SDCA update | Top- calibrated | Convex |

One-vs-all (OVA) SVM | [48] | [48] | no (Prop. 2.2) | yes | ||

yes (Prop. 2.2) | ||||||

Multiclass SVM | [28, 48] | [28, 48] | no (Prop. 2.2) | |||

Softmax (maximum entropy) | Prop. 2.4 | Prop. 3 | yes (Prop. 2.2) | |||

Top- hinge () | [28] | [28] | open
question for |
|||

Top- hinge () | ||||||

Smooth top- hinge () | Eq. (12) w/ | Prop. 2.3 | Prop. 3 | |||

Smooth top- hinge () | Eq. (12) w/ | |||||

Top- entropy | Prop. 2.4 | Eq. (14) | Prop. 3 | |||

Truncated top- entropy | Eq. (22) | - | - | yes (Prop. 2.5) | no | |

Note that and . | ||||||

We let (binary one-vs-all); , (multiclass); . |

Contributions.
We study the problem of top- error optimization
on a diverse range of learning tasks.
We consider existing methods
as well as propose 4 novel loss functions for minimizing the top- error.
A brief overview of the methods is given in Table 1.
For the proposed convex top- losses,
we develop an efficient optimization scheme based on SDCA^{1}^{1}1
Code available at: https://github.com/mlapin/libsdca
,
which can also be used for training with the softmax loss.
All methods are evaluated empirically in terms of the top- error
and, whenever possible, in terms of classification calibration.
We discover that the softmax loss
and the proposed smooth top- SVM are astonishingly
competitive in all top- errors.
Further small improvements can be obtained with the new top- losses.

We consider multiclass problems with classes where the training set consists of examples along with the corresponding labels . We use and to denote a permutation of (indexes) . Unless stated otherwise,

reorders components of a vector

in descending order,While we consider linear classifiers in our experiments, all loss functions below are formulated in the general setting where a function

is learned and prediction at test time is done via , resp. the top- predictions. For the linear case, all predictors have the form . Let be the stacked weight matrix, be a convex loss function, and be a regularization parameter. We consider the following multiclass optimization problemWe use the Iverson bracket notation , defined as if is true, otherwise; and introduce a shorthand . We generalize the standard zero-one error and allow guesses instead of one. Formally, the top- zero-one loss (top- error) is

(1) |

Note that for we recover the standard zero-one error. Top- accuracy is defined as minus the top- error.

In this section, we establish the best achievable top- error, determine when a classifier achieves it, and define a notion of top- calibration.

The Bayes optimal top- error at is

where . A classifier is top- Bayes optimal at if and only if

where .

Let and be a permutation such that . The expected top- error at is

The error is minimal when is maximal, which corresponds to taking the

largest conditional probabilities

and yields the Bayes optimal top- error at .Since the relative order within is irrelevant for the top- error, any classifier , for which the sets and coincide, is Bayes optimal.

Note that we assumed that there is a clear cut between the most likely classes and the rest. In general, ties can be resolved arbitrarily as long as we can guarantee that the largest components of correspond to the classes (indexes) that yield the maximal sum and lead to top- Bayes optimality. ∎

Optimization of the zero-one loss (and, by extension, the top- error)
leads to hard combinatorial problems.
Instead, a standard approach is to use a convex surrogate loss
which upper bounds the zero-one error.
Under mild conditions on the loss function [3, 52],
the optimal classifier w.r.t. the surrogate yields
a Bayes optimal solution for the zero-one loss.
Such loss is called *classification calibrated*

, which is known in statistical learning theory as a necessary condition for a classifier to be universally Bayes consistent

[3]. We introduce now the notion of calibration for the top- error. A loss function (or a reduction scheme) is called top- calibrated if for all possible data generating measures on and allIf a loss is *not* top- calibrated,
it implies that even in the limit of infinite data,
one does not obtain a classifier with
the Bayes optimal top- error from Lemma 2.1.

The standard multiclass problem is often solved using the one-vs-all (OVA) reduction into a set of binary classification problems. Every class is trained versus the rest which yields classifiers .

Typically, the binary classification problems are formulated with a convex margin-based loss function , where and . We consider in this paper:

(2) | ||||

(3) |

The hinge (2) and
logistic (3) losses correspond to
the SVM and logistic regression respectively.
We now show when the OVA schemes are top- calibrated,
not only for (standard multiclass loss) but *for all* simultaneously.

The OVA reduction is top- calibrated for any if the Bayes optimal function of the convex margin-based loss is a strictly monotonically increasing function of .

For every class , the Bayes optimal classifier for the corresponding binary problem has the form

where is a strictly monotonically increasing function. The ranking of corresponds to the ranking of and hence the OVA reduction is top- calibrated for any . ∎

Next, we check if the one-vs-all schemes employing hinge and logistic regression losses are top- calibrated. OVA SVM is not top- calibrated.

First, we show that the Bayes optimal function for the binary hinge loss is

We decompose the expected loss as

Thus, one can compute the Bayes optimal classifier pointwise by solving

for every , which leads to the following problem

where . It is obvious that the optimal is contained in . We get

The minimum is attained at the boundary and we get

Therefore, the Bayes optimal classifier for the hinge loss is not a strictly monotonically increasing function of .

To show that OVA hinge is not top- calibrated, we construct an example problem with classes and , . Note that for every class , the Bayes optimal binary classifier is , hence the predicted ranking of labels is arbitrary and may not produce the Bayes optimal top- error. ∎

In contrast, logistic regression is top- calibrated.

OVA logistic regression is top- calibrated.

First, we show that the Bayes optimal function for the binary logistic loss is

As above, the pointwise optimization problem is

The logistic loss is known to be convex and differentiable and thus the optimum can be computed via

Re-writing the first fraction we get

which can be solved as and leads to the formula for the Bayes optimal classifier stated above.

We check now that the function defined as is strictly monotonically increasing.

The derivative is strictly positive on , which implies that is strictly monotonically increasing. The logistic loss, therefore, fulfills the conditions of Lemma 2.2 and is top- calibrated for any . ∎

An alternative to the OVA scheme with binary losses is to use a *multiclass* loss
directly.
We consider two generalizations of the hinge and logistic losses below:

(4) | |||

(5) |

Both the multiclass hinge loss (4) of Crammer & Singer [15] and the softmax loss (5) are popular losses for multiclass problems. The latter is also known as the cross-entropy or multiclass logistic loss and is often used as the last layer in deep architectures [6, 26, 50]. The multiclass hinge loss has been shown to be competitive in large-scale image classification [2], however, it is known to be not calibrated [52] for the top- error. Next, we show that it is not top- calibrated for any .

Multiclass SVM is not top- calibrated.

First, we derive the Bayes optimal function.

Let . Given any , a Bayes optimal function for the loss (4) is

Let , then

Suppose that the maximum of is not unique. In this case, we have

as the term is always active. The best possible loss is obtained by setting for all , which yields an expected loss of . On the other hand, if the maximum is unique and is achieved by , then

As the loss only depends on the gap , we can optimize this with .

As only the minimal enters the last term, the optimum is achieved if all are equal for (otherwise it is possible to reduce the first term without affecting the last term). Let for all . The problem becomes

Let . The solution is

and the associated risk is

If , then the Bayes optimal classifier for all and any . Otherwise, and

Moreover, we have that the Bayes risk at is

It follows, that the multiclass hinge loss is not (top-) classification calibrated at any where as its Bayes optimal classifier reduces to a constant. Moreover, even if for some , the loss is not top- calibrated for as the predicted order of the remaining classes need not be optimal. ∎

Again, a contrast between the hinge and logistic losses.

The softmax loss is top- calibrated.

The multiclass logistic loss is (top-) calibrated for the zero-one error in the following sense. If

then for some and all

which implies

We now prove this result and show that it also generalizes to top- calibration for . Using the identity

and the fact that , we write for a

As the loss is convex and differentiable, we get the global optimum by computing a critical point. We have

for . We note that the critical point is not unique as multiplication leaves the equation invariant for any . One can verify that satisfies the equations for any . This yields a solution

for any fixed . We note that is a strictly monotonically increasing function of the conditional class probabilities. Therefore, it preserves the ranking of and implies that is top- calibrated for any . ∎

The implicit reason for top-

calibration of the OVA schemes and the softmax loss is that one can estimate the probabilities

from the Bayes optimal classifier. Loss functions which allow this are calledWe have established that the OVA logistic regression and the softmax loss are top- calibrated for any , so why should we be interested in defining new loss functions for the top- error? The reason is that calibration is an asymptotic property as the Bayes optimal functions are obtained pointwise. The picture changes if we use linear classifiers, since they obviously cannot be minimized independently at each point. Indeed, most of the Bayes optimal classifiers cannot be realized by linear functions.

In particular, convexity of the softmax and multiclass hinge losses leads to phenomena where , but . This happens if and adds a bias when working with “rigid” function classes such as linear ones. The loss functions which we introduce in the following are modifications of the above losses with the goal of alleviating that phenomenon.

Recently, we introduced two top- versions of the multiclass hinge loss (4) in [28], where the second version is based on the family of ranking losses introduced earlier by [54]. We use our notation from [28] for direct comparison and refer to the first version as and the second one as . Let , where is the all ones vector, is the -th basis vector, and let be defined componentwise as . The two top- hinge losses are

(6) | ||||

(7) |

where is the -th largest component of . It was shown in [28] that (6) is a tighter upper bound on the top- error than (7), however, both losses performed similarly in our experiments. In the following, we simply refer to them as the top- hinge or the top- SVM loss.

Both losses reduce to the multiclass hinge loss
(4) for .
Therefore, they are unlikely to be top- calibrated,
even though we can currently neither
prove nor disprove this for .
The multiclass hinge loss is not calibrated as
it is non-smooth and does not allow to estimate
the class conditional probabilities .
Our new family of *smooth* top- hinge losses
is based on the Moreau-Yosida regularization
[5, 34].
This technique has been used in [48]
to smooth the binary hinge loss (2).
Interestingly, smooth binary hinge loss fulfills
the conditions of Lemma 2.2
and leads to a top- calibrated OVA scheme.
The hope is that the smooth top- hinge loss
becomes top- calibrated as well.

Smoothing works by adding a quadratic term to the conjugate function^{2}^{2}2
The convex conjugate of is
.
,
which then becomes strongly convex.
Smoothness of the loss, among other things,
typically leads to much faster optimization as
we discuss in Section 3.

OVA smooth hinge is top- calibrated.

In order to derive the smooth hinge loss, we first compute the conjugate of the standard binary hinge loss,

(8) |

The smoothed conjugate is

The corresponding primal smooth hinge loss is given by

(9) |

is convex and differentiable with the derivative

We compute the Bayes optimal classifier pointwise.

Let , the optimal is found by solving

Case . Consider the case ,

This case corresponds to , which follows from the constraint . Next, consider ,

unless , which is already captured by the first case. Finally, consider . Then

where we have if . We obtain the Bayes optimal classifier for as follows:

Note that while is not a continuous function of for , it is still a strictly monotonically increasing function of for any .

Case . First, consider ,

From , we get the condition . Next, consider ,

which is in the range if . Finally, consider ,

where we have if . Overall, the Bayes optimal classifier for is

Note that is again a strictly monotonically increasing function of . Therefore, for any , the one-vs-all scheme with the smooth hinge loss (9) is top- calibrated for all by Lemma 2.2. ∎

Next, we introduce the *multiclass* smooth top- hinge losses,
which extend the top- hinge losses (6)
and (7).
We define the top- simplex ( and ) of radius as

(10) | ||||

(11) |

We also let and . [[28]] The convex conjugates of (6) and (7) are respectively , if , otherwise; and , if , otherwise.

Smoothing applied to the top- hinge loss (6) yields the following smooth top- hinge loss (). Smoothing of (7) is done similarly, but the set is replaced with .

Let be the smoothing parameter. The smooth top- hinge loss () and its conjugate are

(12) | ||||

(13) |

where is the Euclidean projection of on . Moreover, is -smooth.

We take the convex conjugate of the top- hinge loss, which was derived in [28, Proposition 2],

and add the regularizer to obtain the -strongly convex conjugate loss as stated in the proposition. As mentioned above [21] (see also [48, Lemma 2]), the primal smooth top- hinge loss , obtained as the convex conjugate of , is -smooth. We now obtain a formula to compute it based on the Euclidean projection onto the top- simplex. By definition,

For the constraint , we have

The final expression follows from the fact that

∎

There is no analytic expression for (12) and evaluation requires computing a projection onto the top- simplex