example_perceptron_kernel_gaussian_yinyang.cpp

 
#include <iostream>
#include <Pds/Ra>
#include <Pds/Ml>
 
    
int main(void)
{
    Pds::ClassificationMetrics Metrics;
    Pds::Vector Yp;
    Pds::Matrix F;
    Pds::IterationConf Conf;
    double Sigma=0.6;
    
    // Generating data
    unsigned int L=1000;
    Pds::Matrix X;
    Pds::Vector Y;
    Pds::Matrix G;
    Pds::DataSet::Classification::LoadDataYinYang(L,X,Y);
    
    // Create Perceptron
    G=Pds::MergeVer( {Pds::Clustering::Kmeans(Conf,X.GetRows(Pds::Find(Y)),60),
                      Pds::Clustering::Kmeans(Conf,X.GetRows(Pds::Find(1-Y)), 60)}
                   ); //X
    F=Pds::Kernel::Gaussian(X,G,Sigma);
    Conf.SetMaxIter(10000);
    Pds::Perceptron Neurona(Conf,F,Y);
    Neurona.Print("\nNeurona:\n");
    
    // Evaluate training data
    Yp=Neurona.Evaluate(F);
    
    // Metrics of training
    Metrics = Pds::ClassificationMetrics::Calculate(0.5,Yp,Y);
    Metrics.Print("\n");
        
    // Evaluate testing data
    Pds::DataSet::Classification::LoadDataYinYang(L,X,Y); // Load new data
    F=Pds::Kernel::Gaussian(X,G,Sigma);
    Yp=Neurona.Evaluate(F);
    
    // Metrics testing
    Metrics = Pds::ClassificationMetrics::Calculate(Metrics.Threshold,Yp,Y);
    Metrics.Print("\n");
    return 0;
}