example_logisticregression_meansquare_hard.cpp

 
#include <iostream>
#include <Pds/Ra>
#include <Pds/Ml>
 
        
int main(void)
{
    // Generating data
    Pds::Matrix X(Pds::Ra::TextFormat,"../test/hard2-data-x.txt");
    Pds::Vector Y(Pds::Ra::TextFormat,"../test/hard2-data-y.txt");
    
    
    Pds::Vector W;
    
    Pds::IterationConf Conf;
    Pds::Vector d=Y.BalancedBinaryWeight(0.5);
    
    
    // Create W using training data
    Conf.Show=true; Conf.SetMinError(1e-07);    Conf.SetAlpha(0.1);
    
    W=Pds::LogisticModel::FittingLogitMeanSquare(Conf,X,Y,0.0001);
    W.T().Print("W_ms_conf:\n");
    
    Pds::Octave::XLabel="x_1";
    Pds::Octave::YLabel="x_2";
    Pds::Octave::Plot::PointsX2DYW(X,Y,W,"testando.m","example_logisticregression_meansquare_hard_ms-conf.png");
    
    // Create W using training data
    W=Pds::LogisticModel::FittingLogitMeanSquare(X,Y,0.0001);
    W.T().Print("W_ms:\n");
    
    Pds::Octave::Plot::PointsX2DYW(X,Y,W,"testando.m","example_logisticregression_meansquare_hard_ms.png");
    
    
    // Create W using training data
    Conf.Show=true; Conf.SetMinError(1e-07);    Conf.SetAlpha(0.1);
    
    W=Pds::LogisticModel::FittingLogitWeightedMeanSquare(Conf,X,Y,d,0.0001);
    W.T().Print("W_ms_conf:\n");
    
    Pds::Octave::Plot::PointsX2DYW(X,Y,W,"testando.m","example_logisticregression_meansquare_hard_wms-conf.png");
    
    // Create W using training data
    W=Pds::LogisticModel::FittingLogitWeightedMeanSquare(X,Y,d,0.0001);
    W.T().Print("W_ms_w:\n");
    
    Pds::Octave::Plot::PointsX2DYW(X,Y,W,"testando.m","example_logisticregression_meansquare_hard_wms.png");
 
    
    return 0;
}