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;
}