package rule;

import java.util.HashMap;

import java.util.List;

import java.util.Map;

//格拉布斯准则

public class GLBS

{

static Map map=new HashMap();

List list;

double mean; //均值

double mean2; //平方均值

int listNum;

double variance; //标准差

public GLBS(List list)

{

this.list=list;

if(map.size()==0)

{

initMap();

}

mean_variance();

}

/**

* 均值和方差

*/

private void mean_variance()

{

listNum= list.size();

for(Double d: list)

{

mean+=d;

mean2+=d*d;

}

mean/=listNum; //均值

mean2/=listNum;

variance=Math.sqrt(mean2-mean*mean); // 标准差

}

private void initMap()

{

double[][] crisis= {

{0,0,0,1.148,1.425,1.602,1.729,1.828,1.909,1.977,2.036,2.088,2.134,2.175,2.213,2.247,2.279,2.309,2.335,

2.361,2.385,2.408,2.429,2.448,2.467,2.486,2.502,2.519,2.534,2.549,2.583,2.577,2.591,2.604,

2.616,2.628,2.639,2.65,2.661,2.671,2.682,2.692,2.7,2.71,2.719,2.727,2.736,2.744,2.753,2.76,

2.768,2.775,2.783,2.79,2.798,2.804,2.811,2.818,2.824,2.831,2.837,2.842,2.849,2.854,2.86,2.866,

2.871,2.877,2.883,2.888,2.893,2.897,2.903,2.908,2.912,2.917,2.922,2.927,2.931,2.935,2.94,2.945,

2.949,2.953,2.957,2.961,2.966,2.97,2.973,2.977,2.981,2.984,2.989,2.993,2.996,3,3.003,3.006,

3.011,3.014,3.017},

{0,0,0,1.153,1.463,1.672,1.822,1.938,2.032,2.11,2.176,2.234,2.285,2.331,2.371,2.409,2.443,2.475,2.501,

2.532,2.557,2.58,2.603,2.624,2.644,2.663,2.681,2.698,2.714,2.73,2.745,2.759,2.773,2.786,2.799,

2.811,2.823,2.835,2.846,2.857,2.866,2.877,2.887,2.896,2.905,2.914,2.923,2.931,2.94,2.948,2.956,

2.943,2.971,2.978,2.986,2.992,3,3.006,3.013,3.019,3.025,3.032,3.037,3.044,3.049,3.055,3.061,

3.066,3.071,3.076,3.082,3.087,3.092,3.098,3.102,3.107,3.111,3.117,3.121,3.125,3.13,3.134,3.139,

3.143,3.147,3.151,3.155,3.16,3.163,3.167,3.171,3.174,3.179,3.182,3.186,3.189,3.193,3.196,3.201,

3.204,3.207},

{0,0,0,1.155,1.481,1.715,1.887,2.02,2.126,2.215,2.29,2.355,2.412,2.462,2.507,2.549,2.585,2.62,2.651,

2.681,2.709,2.733,2.758,2.781,2.802,2.822,2.841,2.859,2.876,2.893,2.908,2.924,2.938,2.952,

2.965,2.979,2.991,3.003,3.014,3.025,3.036,3.046,3.057,3.067,3.075,3.085,3.094,3.103,3.111,3.12,

3.128,3.136,3.143,3.151,3.158,3.166,3.172,3.18,3.186,3.193,3.199,3.205,3.212,3.218,3.224,3.23,

3.235,3.241,3.246,3.252,3.257,3.262,3.267,3.272,3.278,3.282,3.287,3.291,3.297,3.301,3.305,

3.309,3.315,3.319,3.323,3.327,3.331,3.335,3.339,3.343,3.347,3.35,3.355,3.358,3.362,3.365,3.369,

3.372,3.377,3.38,3.383},

{0,0,0,1.155,1.492,1.749,1.944,2.097,2.22,2.323,2.41,2.485,2.55,2.607,2.659,2.705,2.747,2.785,2.821,

2.954,2.884,2.912,2.939,2.963,2.987,3.009,3.029,3.049,3.068,3.085,3.103,3.119,3.135,3.15,3.164,

3.178,3.191,3.204,3.216,3.228,3.24,3.251,3.261,3.271,3.282,3.292,3.302,3.31,3.319,3.329,3.336,

3.345,3.353,3.361,3.388,3.376,3.383,3.391,3.397,3.405,3.411,3.418,3.424,3.43,3.437,3.442,3.449,

3.454,3.46,3.466,3.471,3.476,3.482,3.487,3.492,3.496,3.502,3.507,3.511,3.516,3.521,3.525,3.529,

3.534,3.539,3.543,3.547,3.551,3.555,3.559,3.563,3.567,3.57,3.575,3.579,3.582,3.586,3.589,3.593,

3.597,3.6}};

for(int i=3; i<=100; i++)

{

add(i, 0.90, crisis[0][i]);

add(i, 0.95, crisis[1][i]);

add(i, 0.975, crisis[2][i]);

add(i, 0.99, crisis[3][i]);

}

add(101, 0.90, 3.02);

add(101, 0.95, 3.21);

add(101, 0.975, 3.39);

add(101, 0.99, 3.65);

}

private void add(int i, double d, double e)

{

map.put(i+d, e);

}

private double Subtract(int n, double a, double value)

{

if(n>100)

{

n=101;

}

return value-map.get(n+a);

}

/**

* @return 返回剔除了所有异常数据的list

*/

public List clipList()

{

int i=0;

while(clipOne())

{

if(listNum<3)

{

break;

}

i++;

}

System.out.println(i);

return list;

}

/**

* 至多剔除一个异常数据

* @return

*/

private boolean clipOne()

{

double t1=Math.abs(list.get(0)-mean)/variance;

double tn=Math.abs(list.get(listNum-1)-mean)/variance;

double sub1=Subtract(listNum, 0.90, t1);

double subn=Subtract(listNum, 0.90, tn);

if(sub1>0||subn>0)

{

double remove;

if(sub1>subn)

{

remove=list.get(0);

list.remove(0);

}

else

{

remove=list.get(listNum-1);

list.remove(listNum-1);

}

mean=(listNum*mean-remove)/(listNum-1); //均值更新

mean2=(listNum*mean2-remove*remove)/(listNum-1);

variance=mean2-mean*mean; //方差更新

listNum--; //个数更新

return true;

}

return false;

}

public double getMean()

{

return mean;

}

}