当前位置:   article > 正文

常用去除离群值的算法!_怎么滤除离群值 c语言

怎么滤除离群值 c语言

1.引言:

       数据预处理的方法主要包括去极值、标准化、中性化。那么这里介绍几个常用的去离群值的方法,所谓离群值就是异常值,这个和极值其实并不一样。常用的剔除离群值的方法有MAD、3\sigma、百分位法。参考资料:数据预处理(上)之离群值处理、标准化数据预处理(下)之中性化

2.MAD算法

        MAD,即median absolute deviation,可译为绝对中位值偏差。其大致思想是通过判断每一个元素与中位值的偏差是否处于合理的范围内来判断该元素是否为离群值。具体方法如下:

(1)计算所有元素的中位值:X_{median}

(2)计算所有元素与中位值的绝对偏差:bias = |X_{i} - X_{median}|

(3)取得绝对偏差的中位值:MAD = bias_{median}

(4)确定参数n,则可以对所有的数据作如下调整:

在这里我没有调整离群值,而是直接将其剔除了。

代码如下:

  1. import numpy as np
  2. # MAD法: media absolute deviation
  3. def MAD(dataset, n):
  4. median = np.median(dataset) # 中位数
  5. deviations = abs(dataset - median)
  6. mad = np.median(deviations)
  7. remove_idx = np.where(abs(dataset - median) > n * mad)
  8. new_data = np.delete(dataset, remove_idx)
  9. return new_data

2.3\sigma

      3\sigma法又称为标准差法。标准差本身可以体现因子的离散程度,和MAD算法类似,只是3\sigma法用到的不是中位值,而是均值,并且n的取值为3,代码如下:

  1. # 3sigma法
  2. def three_sigma(dataset, n= 3):
  3. mean = np.mean(dataset)
  4. sigma = np.std(dataset)
  5. remove_idx = np.where(abs(dataset - mean) > n * sigma)
  6. new_data = np.delete(dataset, remove_idx)
  7. return new_data

3.百分位法

       百分位计算的逻辑是将因子值进行升序的排序,对排位百分位高于97.5%或排位百分位低于2.5%的因子值,类似于比赛中”去掉几个最高分,去掉几个最低分“的做法。代码如下:这里参数采用的是20%和80%,具体取值,还需具体情况具体分析。

  1. # 百分位法:原始参数 min=0.025, max=0.975
  2. def percent_range(dataset, min= 0.20, max= 0.80):
  3. range_max = np.percentile(dataset, max * 100)
  4. range_min = -np.percentile(-dataset, (1 - min) * 100)
  5. # 剔除前20%和后80%的数据
  6. new_data = []
  7. for value in dataset:
  8. if value < range_max and value > range_min:
  9. new_data.append(value)
  10. return new_data

这三个方法思路简单,易于实现,但是只能处理一维数据,接下来,再考虑更复杂但却更精准的可适用于多维数据的离群值处理方法!参考资料:离群点检测---基于kNN的离群点检测、LOF算法和CLOF算法

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/2023面试高手/article/detail/301914
推荐阅读
相关标签
  

闽ICP备14008679号