User:Atastorino/sandbox

Statistical learning theory is a framework for machine learning drawing from the fields of statistics and functional analysis. Statistical learning theory deals with the problem of finding a predictive function based on data. Statistical learning theory lead to successful applications in fields such as computer vision, speech recognition, and bioinformatics. It is the theoretical framework underlying support vector machines.

=Introduction=

The goal of learning is prediction. Learning falls into many categories, including supervised learning, unsupervised learning, online learning, and reinforcement learning. From the perspective of statistical learning theory, supervised learning is best understood. Supervised learning involves learning from a training set of data. Every point in the training is an input-output pair, where the input maps to an output. The learning problem consists of inferring the function that maps between the input and the output in a predictive fashion, such that the learned function can be used to predict output from future input.

Depending of the type of output, supervised learning problems are either problems of regression or problems of classification. If the output takes a continuous range of values, it is a regression problem. Using Ohm's Law as an example, a regression could be performed with voltage as input and current as output. The regression would find the functional relationship between voltage and current to be $\frac{1}{R}$, such that

I = \frac{1}{R} V $$ Classification problems are those for which the output will be an element from a discrete set of labels. Classification is very common for machine learning applications. In facial recognition, for instance, a picture of a person's face would be the input, and the output label would be that person's name. The input would be represented by a large multidimensional vector, in which each dimension represents the value of one of the pixels.

After learning a function based on the training set data, that function is validated on a test set of data, data that did not appear in the training set. Classification functions can use the percentage of inputs that are correctly classified as a metric for how predictive the learned function is, while regression functions must use some distance metric, called a loss function, for how accurate the predicted value is. A familiar example of a loss function is the square of the difference between the actual value and the predicted value; this is the loss function used in ordinary least squares regression.

=Formal Description=

Take $$X$$ to be the vector space of all possible inputs, and $$Y$$ to be the vector space of all possible outputs. Statistical learning theory takes the perspective that there is some unknown probability distribution over the product space $$Z = X \otimes Y$$, i.e. there exists some unknown $$p(z) = p(\vec{x},y)$$. The training set is made up of $$n$$ samples from this probability distribution, and is notated
 * $$S = \{(\vec{x}_1,y_1), \dots ,(\vec{x}_n,y_n)\} = \{\vec{z}_1, \dots ,\vec{z}_n\}$$

Every $$\vec{x}_i$$ is an input vector from the training data, and $$y_i$$ is the output that corresponds to it.

In this formalism, the inference problem consists of finding a function $$f: X \mapsto Y$$ such that $$f(\vec{x}) \sim y$$. Let $$\mathcal{H}$$ be a space of functions $$f: X \mapsto Y$$ called the hypothesis space. The hypothesis space is the space of functions the algorithm will search through. Let $$V(f(\vec{x}),y)$$ be the loss function, a metric for the difference between the predicted value $$f(\vec{x})$$ and the actual value $$y$$. The expected risk is defined to be
 * $$I[f] = \displaystyle \int_{X \otimes Y} V(f(\vec{x}),y) p(\vec{x},y) d\vec{x} dy$$

The target function, the best possible function $$f$$ that can be chosen, is given by the $$f$$ that satisfies
 * $$\inf_{f \in \mathcal{H}} I[f]$$

Because the probability distribution $$p(\vec{x},y)$$ is unknown, a proxy measure for the expected risk must be used. This measure is based on the training set, a sample from this unknown probability distribution. It is called the empirical risk
 * $$I_S[f] = \frac{1}{n} \displaystyle \sum_{i=1}^n V(\vec{x}_i,y_i)$$

A learning algorithm that chooses the function $$f_S$$ that minimizes the empirical risk is called empirical risk minimization.

=Loss Functions= The choice of loss function is a determining factor on the function $$f_S$$ that will be chosen by the learning algorithm. The loss function also affects the convergence rate for an algorithm. It is important for the loss function to be convex.

Different loss functions are used depending on whether the problem is one of regression or one of classification.

Regression
The most common loss function for regression is the square loss function. This familiar loss function is used in ordinary least squares regression. The form is:
 * $$V(f(\vec{x}),y) = (y - f(\vec{x}))^2$$

The absolute value loss is also sometimes used:
 * $$V(f(\vec{x}),y) = |y - f(\vec{x})|$$

Classification
In some sense the 0-1 loss is the most natural loss function for classification. It takes the value 0 if the predicted output is the same as the actual output, and it takes the value 1 if the predicted output is different from the actual output. For binary classification, this is:
 * $$V(f(\vec{x},y)) = \theta(- y f(\vec{x}))$$

where $$\Theta$$ is the Heavyside step function.

The 0-1 loss function, however, is not convex. The hinge loss is thus often used:
 * $$V(f(\vec{x},y)) = (- y f(\vec{x}))_+$$

=Regularization= In machine learning problems, a major problem that arises is that of overfitting. Because learning is a prediction problem, the goal is not to find a function that most closely fits the data, but to find one that will most accurately will predict output from future input. Empirical risk minimization runs this risk of overfitting: finding a function that matches the data exactly but does not predict future output well.

Overfitting is symptomatic of unstable solutions; a small pertubation in the training set data would cause a large variation in the learned function. It can be shown that if the stability for the solution can be guaranteed, generalization and consistency are guaranteed as well. Regularization can solve the overfitting problem and give the problem stability.

Regularization can be accomplished by restricting the hypothesis space $$\mathcal{H}$$. A common example would be restricting $$\mathcal{H}$$ to linear functions: this can be seen as a reduction to the standard problem of linear regression. $$\mathcal{H}$$ could also be restricted to polynomial of degree $$p$$, exponentials, or bounded functions on L1. Restriction of the hypothesis space avoids overfitting because the form of the potential functions are limited, and so does not allow for the choice of a function that gives empirical risk arbitrarily close to zero.

Regularization can also be accomplished through Tikhonov regularization. This consists of minimizing
 * $$\frac{1}{n} \displaystyle \sum_{i=1}^n V(f(\vec{x}_i,y_i)) + \gamma

where $$\gamma$$ is a fixed and positive parameter, the regularization parameter. Tikhonov regularization ensures existence, uniqueness, and stability of the solution.
 * f||_{\mathcal{H}}^2$$

=Related Articles=

Reproducing kernel Hilbert spaces are a useful choice for $$\mathcal{H}$$.

Regularized least squares is an example of of Tikhonov regularization in use.

=References=