statistics.h

Go to the documentation of this file.

// -*- Mode: C++; tab-width: 2; -*-
// vi: set ts=2:
//
// 

#ifndef STATISTICS
#define STATISTICS

#include <vector>
#include <iostream>
#include <cmath>
#include <BALL/MATHS/parsedFunction.h>

#include <set>
#include <Eigen/Core>

namespace BALL 
{
  namespace QSAR
  {
    class BALL_EXPORT Statistics
    {
      public:
        Statistics(){};
        
        ~Statistics(){};
        
        
        
        static void scaling(std::vector<double>& v);
  
        static void scaling(std::vector<std::vector<double> >& m);
  
        static void centering(std::vector<std::vector<double> >& m);
        
        static void centering(std::vector<double>& v);
        
        static void centering(vector<double>& v, double& mean, double& std);
        
        static double getVariance(const std::vector<double>& v, double mean=-1);
        
        static double getStddev(const vector<double>& v, double mean=-1);
        
        static double getCovariance(const vector<double>& v1, const vector<double>& v2, double mean1, double mean2);
  
        static double getMean(const std::vector<double>& v);
        
        static double calculateRankCorrelation(vector<double>& observed_values, vector<double>& expected_values);
        //--- methods for calculating mean, covar, var of matrix-ROWS ---
        
        static double getRowCovariance(const vector<vector<double> >& v, int row1, int row2, double mean1=-1, double mean2=-1, std::multiset<int>* features_to_use=0);
        
        static double getRowMean(const vector<vector<double> >& v, int row, std::multiset<int>* features_to_use=0);
        
        static double getRowVariance(const vector<vector<double> >& v, int row, double mean=-1, std::multiset<int>* features_to_use=0);
        
        static double getRowStddev(const vector<vector<double> >& v, int row, double mean=-1, std::multiset<int>* features_to_use=0);
                
        //  ------------------------
        
        
        static void centering(Eigen::MatrixXd& m);
  
        static void centering(Eigen::MatrixXd& m, int col);
  
        static double getMean(const Eigen::MatrixXd& m, int col);
  
        static double getVariance(const Eigen::MatrixXd& m, int col, double mean=-1);
        
        static double getStddev(const Eigen::MatrixXd& m, int col, double mean=-1);
        
        static double getCovariance(const Eigen::MatrixXd& m, int col1, int col2, double mean1=-1, double mean2=-1);
        
        static double sq(const Eigen::MatrixXd& m, int col, double mean=-1);
        
        static double distance(const Eigen::MatrixXd& m, int& row1, int& row2, double& p);
      
        static double distance(const Eigen::MatrixXd& m1, const Eigen::MatrixXd& m2, int& row1, int& row2, double& p);

        
        static double distance(const Eigen::MatrixXd& m1, const Eigen::MatrixXd& m2, int& row1, int& row2, String& f, String& g);
        
        
        static double euclDistance(const Eigen::MatrixXd& m1, const Eigen::MatrixXd& m2, int row1, int row2);
        
        static double euclNorm(const Eigen::VectorXd& cv);
    
        static double euclDistance(const Eigen::VectorXd& c1, const Eigen::VectorXd& c2);
        
        static double scalarProduct(const Eigen::VectorXd& cv);
    
    };
  }
}


#endif // STATISTICS