Kernel perceptron
# Kernel perceptron

In machine learning, the kernel perceptron is a variant of the popular perceptron learning algorithm that can learn kernel machines, i.e. non-linear classifiers that employ a kernel function to compute the similarity of unseen samples to training samples. The algorithm was invented in 1964, making it the first kernel classification learner.

## Preliminaries

### The perceptron algorithm

The perceptron algorithm is an online learning algorithm that operates by a principle called "error-driven learning". It iteratively improves a model by running it on training samples, then updating the model whenever it finds it has made an incorrect classification with respect to a supervised signal. The model learned by the standard perceptron algorithm is a linear binary classifier: a vector of weights w (and optionally an intercept term b, omitted here for simplicity) that is used to classify a sample vector x as class "one" or class "minus one" according to

    ${\hat {y}}=\operatorname {sgn}(\mathbf {w} ^{\top }\mathbf {x} )$

where a zero is arbitrarily mapped to one or minus one. (The "hat" on ŷ denotes an estimated value.)

In pseudocode, the perceptron algorithm is given by:

    Initialize w to an all-zero vector of length p, the number of predictors (features).
    For some fixed number of iterations, or until some stopping criterion is met:
       For each training example $x_i$ with ground truth label $y_i$ ∈ {-1, 1}:
          Let ŷ = $sgn(w^T$ $x_i).$
          If ŷ ≠ $y_i,$ update w ← w + $y_i$ $x_i.$

### Kernel Methods

By contrast with the linear models learned by the perceptron, a kernel method is a classifier that stores a subset of its training examples $x_i,$ associates with each a weight $α_i,$ and makes decisions for new samples x' by evaluating

    $\operatorname {sgn} \sum _{i}\alpha _{i}y_{i}K(\mathbf {x} _{i},\mathbf {x'} )$.

Here, K is some kernel function. Formally, a kernel function is a non-negative semidefinite kernel (see Mercer's condition), representing an inner product between samples in a high-dimensional space, as if the samples had been expanded to include additional features by a function Φ: K(x, x') = Φ(x) · Φ(x'). Intuitively, it can be thought of as a similarity function between samples, so the kernel machine establishes the class of a new sample by weighted comparison to the training set. Each function x' ↦ $K(x_i,$ x') serves as a basis function in the classification.

## Algorithm

To derive a kernelized version of the perceptron algorithm, we must first formulate it in dual form, starting from the observation that the weight vector w can be expressed as a linear combination of the n training samples. The equation for the weight vector is

    $\mathbf {w} =\sum _{i}^{n}\alpha _{i}y_{i}\mathbf {x} _{i}$

where $α_i$ is the number of times $x_i$ was misclassified, forcing an update w ← w + $y_i$ $x_i.$ Using this result, we can formulate the dual perceptron algorithm, which loops through the samples as before, making predictions, but instead of storing and updating a weight vector w, it updates a "mistake counter" vector α. We must also rewrite the prediction formula to get rid of w:

    ${\begin{aligned}{\hat {y}}&=\operatorname {sgn}(\mathbf {w} ^{\mathsf {T}}\mathbf {x} )\\&=\operatorname {sgn} \left(\sum _{i}^{n}\alpha _{i}y_{i}\mathbf {x} _{i}\right)^{\mathsf {T}}\mathbf {x} \\&=\operatorname {sgn} \sum _{i}^{n}\alpha _{i}y_{i}(\mathbf {x} _{i}\cdot \mathbf {x} )\end{aligned}}$

Plugging these two equations into the training loop turn it into the dual perceptron algorithm.

Finally, we can replace the dot product in the dual perceptron by an arbitrary kernel function, to get the effect of a feature map Φ without computing Φ(x) explicitly for any samples. Doing this yields the kernel perceptron algorithm:

    Initialize α to an all-zeros vector of length n, the number of training samples.
    For some fixed number of iterations, or until some stopping criterion is met:
       For each training example $x_j,$ $y_j:$
          Let ${\hat {y}}=\operatorname {sgn} \sum _{i}^{n}\alpha _{i}y_{i}K(\mathbf {x} _{i},\mathbf {x} _{j})$
          If ŷ ≠ $y_j,$ perform an update by incrementing the mistake counter:
             $α_j$ ← $α_j$ + 1

## Variants and extensions

One problem with the kernel perceptron, as presented above, is that it does not learn sparse kernel machines. Initially, all the $α_i$ are zero so that evaluating the decision function to get ŷ requires no kernel evaluations at all, but each update increments a single $α_i,$ making the evaluation increasingly more costly. Moreover, when the kernel perceptron is used in an online setting, the number of non-zero $α_i$ and thus the evaluation cost grow linearly in the number of examples presented to the algorithm.

The forgetron variant of the kernel perceptron was suggested to deal with this problem. It maintains an active set of examples with non-zero $α_i,$ removing ("forgetting") examples from the active set when it exceeds a pre-determined budget and "shrinking" (lowering the weight of) old examples as new ones are promoted to non-zero $α_i.$

Another problem with the kernel perceptron is that it does not regularize, making it vulnerable to overfitting. The NORMA online kernel learning algorithm can be regarded as a generalization of the kernel perceptron algorithm with regularization. The sequential minimal optimization (SMO) algorithm used to learn support vector machines can also be regarded as a generalization of the kernel perceptron.

The voted perceptron algorithm of Freund and Schapire also extends to the kernelized case, giving generalization bounds comparable to the kernel SVM.

