2.9 清理异常值

import numpy as np
import matplotlib.pyplot as plt

def is_outlier(points,threshold=3.5):
    '''Return a boolen array with True if points are out liers and Fa;se otherwise
    data point with a modified zscore fjdsklfadlkfa dfjalskerwwerwerwerwe
    rw
    :parameterwe
    :raisewer
    :keyworder
    we
    r
    we
    r
    etr
xiaxiede
    :returne
    2wwe
    :returnff
    fwr
    wr
    '''
    if len(points.shape) == 1:
        points = points[:, None]

    median = np.median(points, axis=0)

    diff = np.sum((points-median)**2, axis=-1)
    diff = np.sqrt(diff)

    med_abs_deviation = np.median(diff)
    modified_z_score=0.6745*diff/med_abs_deviation
    return modified_z_score>threshold

x=np.random.random(100)

buckets=50

x=np.r_[x, -49, 95, 100, -100]
filtered=x[~is_outlier(x)]
plt.figure()
plt.subplot(211)
plt.hist(x,buckets)
plt.xlabel('Raw')

plt.subplot(212)
plt.hist(filtered,buckets)
plt.xlabel('Cleaned')
plt.show()

 

转载于:https://www.cnblogs.com/boooobao/p/6445354.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值