example_funcgmm_probability.cpp

 
#include <iostream>
#include <Pds/Ra>
#include <Pds/Ml>
    
int main(void)
{
    Pds::Octave::XLabel="x1";
    Pds::Octave::YLabel="x2";
    Pds::Octave::DAspect=true;
    
    Pds::IterationConf Conf; Conf.Show=true; Conf.SetMinError(1.0e-6);
    unsigned int L=1000;
    Pds::Matrix X;
    Pds::Vector Y;
    Pds::Vector P;
    
    // Generating data X=[x1 x2]
    Pds::DataSet::Classification::LoadKGaussian(L,3,X,Y);
    
    Pds::Octave::MarkerSize=3;
    Pds::Octave::Plot::PointsX2D(X,"testando.m","example_funcgmm_probability_data.png");
    Pds::Octave::MarkerSize=12;
    
    // Mu, Sigma estimation
    Pds::DataGMM Dat=Pds::Clustering::GMM(Conf,X,3,1.0);
    Dat.Print("Dat:\n");
    /*
    Dat.ExportFile("Dat.txt");
    Pds::DataGMM Dat2("Dat.txt");
    Dat2.Print("\nDat2:\n");
    std::ofstream myfile;
    myfile.open("Dat2.txt");
    Dat2.SaveInStream(myfile);
    myfile.close();
    Pds::DataGMM Dat3;
    std::ifstream myinfile("Dat2.txt", std::ifstream::in);
    Dat3.LoadFromStream(myinfile);
    myinfile.close();
    Dat3.Print("\nDat3:\n");
    */
    
    // Gaussian Mixture Model density
    P=Pds::Probability::GMMDensity(X,Dat.Pi,Dat.Mu,Dat.Sinv);
    Pds::Octave::Plot::ScatterX2DY(X,P,"testando.m","example_funcgmm_probability_prob.png");
    
    P=Pds::Probability::GMMDensity(X,Dat.Pi,Dat.Mu,Dat.Sinv,Dat.Det);
    Pds::Octave::Plot::ScatterX2DY(X,P,"testando.m","example_funcgmm_probability_prob2.png");
    
    P=Pds::Probability::GMMDensity(X,Dat);
    Pds::Octave::Plot::ScatterX2DY(X,P,"testando.m","example_funcgmm_probability_prob3.png");
    
    // Gaussian Mixture Model partial density
    Pds::Matrix A;
    Pds::Vector ID;
    
    A=Pds::Probability::GMMPartial(X,Dat.Pi,Dat.Mu,Dat.Sinv);
    ID=A.OperateRows(Pds::ArgMax);
    Pds::Octave::Plot::ScatterX2DY(X,ID,"testando.m","example_funcgmm_probability_clusters.png");
    
    A=Pds::Probability::GMMPartial(X,Dat.Pi,Dat.Mu,Dat.Sinv,Dat.Det);
    ID=A.OperateRows(Pds::ArgMax);
    Pds::Octave::Plot::ScatterX2DY(X,ID,"testando.m","example_funcgmm_probability_clusters2.png");
    
    A=Pds::Probability::GMMPartial(X,Dat);
    ID=A.OperateRows(Pds::ArgMax);
    Pds::Octave::Plot::ScatterX2DY(X,ID,"testando.m","example_funcgmm_probability_clusters3.png");
 
    return 0;
}