testpca.cpp

Go to the documentation of this file.

#include "PCA.h"
#include "statutils.h"
#include "UniformRNG.h"
#include "VecMat.h"
#include <iostream>
#include <fstream>
#include <vector>
#include <itpp/base/fastica.h>
#include <itpp/base/mat.h>
#include <itpp/base/sigfun.h>

int main()
{
        UniformRNG Random;
                std::ifstream infile ("in.pca");
        std::vector<double> invector;
        copy(std::istream_iterator<double>(infile),std::istream_iterator<double>(),back_inserter(invector));
        const size_t nobs = invector.size()/2;
        const size_t nchan = 2;
        rmat input(nchan,nobs);
        cmat evec(nchan,nchan);
        cvec eval(nchan);
        itpp::mat ittpInput(nchan,nobs);
        
        for (int i = 0; i < invector.size()/2; i+=2)
        {
                input(0,i/2) = invector.at(i);
                input(1,i/2) = invector.at(i+1);
        }
    for (size_t i = 0; i < input.size2(); ++i)
        for (size_t j = 0; j < input.size1(); ++j)
        {
                //input(j,i) = Random.GetNumber();
                ittpInput(j,i) = input(j,i);
        }
   
        itpp::Fast_ICA fastica(ittpInput);
        //fastica.set_nrof_independent_components(nchan);
        //fastica.set_non_linearity(FICA_NONLIN_TANH);
        //fastica.set_approach(FICA_APPROACH_DEFL);
        fastica.separate();
        PCA(input,evec,eval);
        std::cout << "Input: " << input << std::endl;
        std::cout << "pca evec: " << evec << std::endl;
        std::cout << "pca eval: " << eval << std::endl;
        std::cout << "pca WhM: " << WhiteMat(evec, eval) << std::endl;
        std::cout << "pca DeWhM: " << DeWhiteMat(evec, eval) << std::endl;
        std::cout << "Itpp PEV: " << fastica.get_principal_eigenvectors() << std::endl;
        std::cout << "Itpp WhM: " << fastica.get_whitening_matrix() << std::endl;
        
        
}

---