#ifndef __DEFINITION_OF_STAT_UTIL__
#define __DEFINITION_OF_STAT_UTIL__

#include <set>
#include <vector>
#include <algorithm>

#include <math.h>

bool get_mean_variance(const std::vector<double>& x, double *mean, double *var); 
bool get_mean_variance(const std::set<double>& x, double *mean, double *var); 
bool get_covariance(const std::vector<double>& x, const std::vector<double>& y, double *cov);
bool get_correlation(std::vector<double> *x, std::vector<double>* y, double *corr, double *Rsquare, double *_alpha, double *_beta); 
bool get_robust_correlation(double trim_percent, const std::vector<double>& x, const std::vector<double>& y, double *corr, double *Rsquare, double *_alpha, double *_beta); 
bool get_median(std::vector<double> src, double* dest); 
double rank_z_alpha_from_rank(const std::vector<int>& rank1, const std::vector<int>& rank2, double *_u1, double *_u2, double *_rank_sum1, double *_rank_sum2);
double rank_z_alpha_from_value(const std::vector<double>& value1, const std::vector<double>& value2, double *_u1, double *_u2, double *_rank_sum1, double *_rank_sum2); 
bool remove_by_maha_dist(double trim_perc, const std::vector<double>& srcx, const std::vector<double>& srcy, std::vector<double>* destx, std::vector<double>* desty); 

#endif