zoukankan      html  css  js  c++  java
  • pandas

        异常值概念:是指那些远离正常值的观测,即“不合群”观测。异常值的出现一般是人为的记录错误或者是设备的故障等,异常值的出现会对模型的创建和预测产生

    严重的后果。当然异常值也不一定是坏事,有些情况下,通过寻找异常值就能够给业务带来良好的发展,如销毁“钓鱼”网站,关闭“薅羊毛”用户的权限等。

    异常值的判定方法:

    1.n个标准差法

    2.箱线图法

    标准差法,就是用以样本均值+样本标准差为基准,如果样本离平均值相差2个标准差以上的就是异常值

    箱线图法:以上下四分位作为参考, x > Q3+nIQR 或者 x < Q1 - nIQR 简单地理解,就是如果样本值不在上下四分位+标准差范围内,就是异常值

    两种异常值判定是,如数据近似服从正态分布是,优先选择n个标准差法,因为数据的分布相对比较对称:否则优先选择箱线图法,因为分位数并不会受极端值的影响。

    异常数据处理方法:

    1.删除法(前提是异常观测的比例不能太大)

    2.替换法(可以考虑使用低于判别上下限的最大值或最小值,均值或中位数替换等)

    python处理异常值实例:

    原数据(部分):

     代码:

    import pandas as pd
    sunspots = pd.read_table(r'D:sunspots.csv',sep = ',')
    xbar = sunspots.counts.mean()
    xstd = sunspots.counts.std()
    print('标准差法异常值上限检测: ',any(sunspots.counts > xbar + 2 * xstd))
    print('标准差法异常值下限检测: ',any(sunspots.counts < xbar - 2 * xstd))

    #异常值 箱线图法
    Q1 = sunspots.counts.quantile(q = 0.25)
    Q3 = sunspots.counts.quantile(q = 0.75)
    IQR = Q3 -Q1
    print('箱线图法异常值上限检测: ',any(sunspots.counts > Q3 + 1.5*IQR))
    print('箱线图法异常值下限检测: ',any(sunspots.counts < Q1 - 1.5*IQR))

    out:

    如上结果所示,不管是标准差检验法还是箱线图法,都发现太阳黑子数据中存在异常值,而且异常值都是超过上限临界值的。接下来,通过绘制太阳黑子数量的直方图和核密度曲线图,用于检测数据是否近似服从正态分布,进而选择一个最终的异常值判别方法:

    import matplotlib.pyplot as plt
    plt.style.use('ggplot')
    sunspots.counts.plot(kind='hist',bins = 30,density = True)
    sunspots.counts.plot(kind='kde')
    plt.show()

    如上图所示,不管是直方图还是核密度曲线,所呈现的数据分布形状都是有偏的,并且属于右偏。基于此,这里选择箱线图法来判定太阳黑子数据中的那些异常值。接下来要做的就是选用删除法或替换法来处理这些异常值,删除法就跟上篇我们讲过的一样处理方式,下面介绍一下替换法,即使用上限下限的最大最小值来替换,代码如下:(接异常值检测部分代码)

    print('异常值替换前的数据统计特征: ',sunspots.counts.describe())
    UL = Q3 + 1.5 * IQR
    print('判别异常值的上限临界值: ',UL)
    replace_vaule = sunspots.counts[sunspots.counts < UL].max()
    print('用以替换异常值的数据: ',replace_vaule)

    sunspots.counts[sunspots.counts > UL] = replace_vaule
    print('异常值替换后的数据统计特征: ',sunspots.counts.describe())

    out:

    经过判别异常值,得知,如果一年内太阳黑子超过 148.85时即为异常值年份,对于这些年份的异常值使用141.7替换。

    本篇知识导图:

         

  • 相关阅读:
    Windows 8.1 应用开发文章汇总
    演绎事件委托
    Azure China
    Orchard 刨析:导航篇
    为ASP.NET MVC应用程序使用高级功能
    兼容Mono的下一代云环境Web开发框架ASP.NET vNext
    商品导航菜单
    纳斯达克对经济泡沫的定义是:当交易价格远大于内在价值的时候,就称为泡沫(转)
    当一个低级问题,第一次解决时,你会感受到成就;第二次解决时,你感受到责任,第三次解决时,你可能更多的感受到无力(转)
    互联网组织的未来:剖析GitHub员工的任性之源(转)
  • 原文地址:https://www.cnblogs.com/tinglele527/p/11955103.html
Copyright © 2011-2022 走看看