+ //1%
+ lowCutOff = sqrt(((quantiles[i][j] - average - 3) * (quantiles[i][j] - average + 3)) / (float)(quantiles[i].size() - 1));
+//cout << "high = " << highCutOff << " low = " << lowCutOff << " de = " << quantiles[i][j] << endl;
+ //if this is below the highcutff and above the lowcutoff
+ if ((quantiles[i][j] < highCutOff) && (quantiles[i][j] > lowCutOff)) {
+
+ newQuanI.push_back(quantiles[i][j]);
+
+ }else { cout << "removed outlier: high = " << highCutOff << " low = " << lowCutOff << " de = " << quantiles[i][j] << endl; }