Low-rank matrix approximations are essential tools in the application of kernel methods to large-scale learning problems. Kernel methods (for instance, support vector machines or Gaussian processes) project data points into a high-dimensional or infinite-dimensional feature space and find the optimal splitting hyperplane. In the kernel method the data is represented in a kernel matrix (or, Gram matrix). Many algorithms can solve machine learning problems using the kernel matrix. The main problem of kernel method is its high computational cost associated with kernel matrices. The cost is at least quadratic in the number of training data points, but most kernel methods include computation of matrix inversion or eigenvalue decomposition and the cost becomes cubic in the number of training data. Large training sets cause large storage and computational costs. While low rank decomposition methods (Cholesky decomposition) reduce this cost, they still require computing the kernel matrix. One of the approaches to deal with this problem is low-rank matrix approximations. The most popular examples of them are the Nyström approximation and randomized feature maps approximation methods. Both of them have been successfully applied to efficient kernel learning.
Nyström approximation Kernel methods become computationally unfeasible when the number of points D {\displaystyle D} is so large such that the kernel matrix K {\displaystyle K} cannot be stored in memory. If D {\displaystyle D} is the number of training examples, the storage and computational cost required to find the solution of the problem using general kernel method is O ( D 2 ) {\displaystyle O(D^{2})} and O ( D 3 ) {\displaystyle O(D^{3})} respectively. The Nyström approximation can allow a significant speed-up of the computations. This speed-up is achieved by using, instead of the kernel matrix, its approximation K ~ {\displaystyle {\tilde {K}}} of rank d {\displaystyle d} . An advantage of the method is that it is not necessary to compute or store the whole kernel matrix, but only a submatrix of size d × D {\displaystyle d\times D} . It reduces the storage and complexity requirements to O ( D d ) {\displaystyle O(Dd)} and O ( D d 2 ) {\displaystyle O(Dd^{2})} respectively. The method is named "Nyström approximation" because it can be interpreted as a case of the Nyström method from integral equation theory.
Kernel approximation
… excerpt ends here. Continue reading the full article.
