example_logisticregression_cost.cpp

Programa para el testeo de las funciones.Para compilar o código example_logisticregression_cost.cpp:

g++ -static -o example_logisticregression_cost example_logisticregression_cost.cpp -lpdsmlmm -lpdsramm -lpdsspmm

Para executar o programa:

./example_logisticregression_cost

Retornando por consola:

W: -0.19761437373493    0.69316252037071        0.69316252037071
IG1: 4.5720683890471e-05
    
W: -0.42161372586748    -1      0.63252353657496
IG2: 0.0092228013606607

Código example_logisticregression_cost.cpp:

 
#include <iostream>
#include <Pds/Ra>
#include <Pds/Ml>
        
int main(void)
{
    // Generating data
    Pds::Matrix X(Pds::Ra::TextFormat,"../test/hard-data-x.txt");
    Pds::Vector Y(Pds::Ra::TextFormat,"../test/hard-data-y.txt");
    Pds::Vector W;
    
    double IG1,IG2;
    
    // Create W using Pds::LogisticModel::GetW0MeanMethod()
    W=Pds::LogisticModel::GetW0MeanMethod(X);
    W.T().Print("W: ");
    
    IG1=Pds::LogisticModel::CostInformationGain(W,X,Y);
    std::cout<<"IG1: "<<IG1<<std::endl;
    
    // Create W using Pds::LogisticModel::GetW0MeanSquareMethod()
    W=Pds::LogisticModel::GetW0MeanSquareMethod(X);
    W.T().Print("W: ");
    
    IG2=Pds::LogisticModel::CostInformationGain(W,X,Y);
    std::cout<<"IG2: "<<IG2<<std::endl;
    
    return 0;
}