pcom/statistics.cpp at master · tchronis/pcom · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
/*
 * statistics.cpp
 *
 *  Created on: Jul 19, 2014
 *      Author: tchr
 */

#include<vector>
#include<cmath>

using namespace std;

long double min(vector <long double> r){

	long double minimum=r[0];

	for(long i=1;i<(long)r.size();i++)
		if(r[i]<minimum) minimum=r[i];

	return minimum;
};

long double max(vector <long double> r){

	long double maximum=r[0];

	for(long i=1;i<(long)r.size();i++)
		if(r[i]>maximum) maximum=r[i];

	return maximum;
};

long double mean(vector <long double> r){

	long double sum=0;

	for(long i=0;i<(long)r.size();i++) sum+=r[i];

	return sum/r.size();
};

long double st_dev(vector <long double> r){

	long double sum=0,m=mean(r);

	for(long i=0;i<(long)r.size();i++) sum+=pow(r[i]-m,2);

	return sqrt(sum/(r.size()-1));
};

long double adev(vector <long double> r){

	long double sum=0,m=mean(r);

	for(long i=0;i<(long)r.size();i++) sum+=abs(r[i]-m);

	return sqrt(sum/(r.size()));
};

long double skew(vector <long double> r){

	long double sum=0,m=mean(r);

	for(long i=0;i<(long)r.size();i++) sum+=pow(r[i]-m,3);

	return sum/(r.size()*pow(st_dev(r),3));
};