example_pca.cpp

 
#include <iostream>
#include <Pds/Ra>
#include <Pds/Ml>
    
int main(void)
{
    Pds::Matrix X(Pds::Ra::TextFormat,"../test/pca-data-x.txt");
    Pds::Matrix Y(Pds::Ra::TextFormat,"../test/pca-data-y.txt");
 
    Pds::Octave::Plot::ScatterX3DY(X,Y,"testando3d.m","test_pca_3d.png");    
    
    std::cout<<"X.Size():"<<X.Size()<<"\tY.Size():"<<Y.Size()<<"\n\n";
    
    Pds::PCA Pca(X);
    
    Pca.Print("----- PCA :\n");
    
    Pds::Matrix X1=Pca.Transform(X,2);
    
    Pds::Octave::Plot::ScatterX2DY(X1,Y,"testando2d.m","test_pca_2d.png");
 
    return 0;
}