PCA.h

Go to the documentation of this file.

/*
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or
 * (at your option) any later version.
 *
 * Written (W) 1999-2008 Gunnar Raetsch
 * Written (W) 1999-2008,2011 Soeren Sonnenburg
 * Copyright (C) 1999-2009 Fraunhofer Institute FIRST and Max-Planck-Society
 * Copyright (C) 2011 Berlin Institute of Technology
 */

#ifndef PCA_H_
#define PCA_H_
#include <shogun/lib/config.h>
#ifdef HAVE_LAPACK
#include <shogun/mathematics/lapack.h>
#include <stdio.h>
#include <shogun/preprocessor/DimensionReductionPreprocessor.h>
#include <shogun/features/Features.h>
#include <shogun/lib/common.h>

namespace shogun
{
enum EPCAMode
{
    THRESHOLD,
    VARIANCE_EXPLAINED,
    FIXED_NUMBER
};

class CPCA: public CDimensionReductionPreprocessor
{
    public:

        CPCA(bool do_whitening=false, EPCAMode mode=FIXED_NUMBER, float64_t thresh=1e-6);

        virtual ~CPCA();

        virtual bool init(CFeatures* features);

        virtual void cleanup();

        virtual SGMatrix<float64_t> apply_to_feature_matrix(CFeatures* features);

        virtual SGVector<float64_t> apply_to_feature_vector(SGVector<float64_t> vector);

        SGMatrix<float64_t> get_transformation_matrix();

        SGVector<float64_t> get_eigenvalues();

        SGVector<float64_t> get_mean();

        virtual const char* get_name() const { return "PCA"; }

        virtual EPreprocessorType get_type() const { return P_PCA; }

    protected:

        void init();

    protected:

        SGMatrix<float64_t> m_transformation_matrix;
        int32_t num_dim;
        int32_t num_old_dim;
        SGVector<float64_t> m_mean_vector;
        SGVector<float64_t> m_eigenvalues_vector;
        bool m_initialized;
        bool m_whitening;
        EPCAMode m_mode;
        float64_t thresh;
};
}
#endif
#endif