test_matrix_statistics.cpp

 
#include <Pds/Ra>
    
int main(void)
{
    unsigned int L=4;
    Pds::Matrix A(L,L);
    A.FillRandC(0.5);
    A.Print("A:\n");
    
    Pds::Matrix B=A+2.0*Pds::RandU(L,L);
    
    std::cout<<"A.MultiplyValues(): "<<A.MultiplyValues() <<std::endl;
    std::cout<<"A.Max()  : "<<A.Max() <<std::endl;
    std::cout<<"A.Min()  : "<<A.Min() <<std::endl;
    std::cout<<"A.Sum()  : "<<A.Sum() <<std::endl;
    std::cout<<"A.Mean() : "<<A.Mean()<<std::endl;
    std::cout<<"A.MeanInCols() :\n"<<A.MeanInCols();
    std::cout<<"A.Std()  : "<<A.Std() <<std::endl;
    std::cout<<"A.Corr(B): "<<A.Corr(B)<<std::endl;
    
    Pds::Matrix C("-2 2 2 \n-1 2 -1\n 1 1 3");
    std::cout<<"C.MultiplyValues(): "<<C.MultiplyValues() <<std::endl;
 
    Pds::Matrix D("1 1\n1 -1");
    
    C.Print("C:\n");
    D.Print("D:\n");
    
    C.XCorr(D).Print("C.XCorr(D):\n");
    C.XCorr(D,true).Print("C.XCorr(D,true):\n");
 
    
    A.Print("A:\n");
    std::vector<Pds::Matrix> G(2);
    G[0]=A.GetRows(0,2);    G[0].Print("G[0]:\n");
    G[1]=A.GetRows(3,L-1);  G[1].Print("G[1]:\n");
    double IG=A.InformationGain(G,0.5);
    std::cout<<"IG: "<<IG<<"\n";
    std::cout<<"IG: "<< A.BinaryEntropy(0.5)
                        -(G[0].Nel()*1.0/A.Nel())*G[0].BinaryEntropy(0.5)
                        -(G[1].Nel()*1.0/A.Nel())*G[1].BinaryEntropy(0.5)<<"\n";
    
    C.Print("\nC:\n");
    C.MinInCols().Print("C.MinInCols():\n");
    C.MaxInCols().Print("C.MaxInCols():\n");
    
    std::vector<unsigned int> stdlin;
    
    C.Print("\nC:\n");
    C.MaxInCols(stdlin).Print("C.MaxInCols(stdlin):\n");
    std::cout<<"stdlin: "<<Pds::Ra::IndicesToString(stdlin,"\t")<<"\n";
    C.MinInCols(stdlin).Print("C.MinInCols(stdlin):\n");
    std::cout<<"stdlin: "<<Pds::Ra::IndicesToString(stdlin,"\t")<<"\n";
    
    Pds::Array<unsigned int> Lin;
    
    C.Print("\nC:\n");
    C.MaxInCols(Lin).Print("C.MaxInCols(Lin):\n");
    Lin.Print("Lin:\n");
    C.MinInCols(Lin).Print("C.MinInCols(Lin):\n");
    Lin.Print("Lin:\n");
    
    return 0;
}