15 #ifndef __MLPACK_METHODS_KERNEL_PCA_NAIVE_METHOD_HPP
16 #define __MLPACK_METHODS_KERNEL_PCA_NAIVE_METHOD_HPP
23 template<
typename KernelType>
38 arma::mat& transformedData,
42 KernelType kernel = KernelType())
45 arma::mat kernelMatrix;
47 kernelMatrix.set_size(data.n_cols, data.n_cols);
52 for (
size_t i = 0; i < data.n_cols; ++i)
54 for (
size_t j = i; j < data.n_cols; ++j)
57 kernelMatrix(i, j) = kernel.Evaluate(data.unsafe_col(i),
63 for (
size_t i = 1; i < data.n_cols; ++i)
64 for (
size_t j = 0; j < i; ++j)
65 kernelMatrix(i, j) = kernelMatrix(j, i);
72 arma::rowvec rowMean = arma::sum(kernelMatrix, 0) / kernelMatrix.n_cols;
73 kernelMatrix.each_col() -= arma::sum(kernelMatrix, 1) / kernelMatrix.n_cols;
74 kernelMatrix.each_row() -= rowMean;
75 kernelMatrix += arma::sum(rowMean) / kernelMatrix.n_cols;
78 arma::eig_sym(eigval, eigvec, kernelMatrix);
82 for (
size_t i = 0; i < floor(eigval.n_elem / 2.0); ++i)
83 eigval.swap_rows(i, (eigval.n_elem - 1) - i);
86 eigvec = arma::fliplr(eigvec);
88 transformedData = eigvec.t() * kernelMatrix;
89 transformedData.each_col() /= arma::sqrt(eigval);
Linear algebra utility functions, generally performed on matrices or vectors.
static void ApplyKernelMatrix(const arma::mat &data, arma::mat &transformedData, arma::vec &eigval, arma::mat &eigvec, const size_t, KernelType kernel=KernelType())
Construct the exact kernel matrix.
Include all of the base components required to write MLPACK methods, and the main MLPACK Doxygen docu...