Python筛选异常值

# -*- coding:utf-8 -*-
"""
作者：sunli
日期：2022年01月04日
"""
# 导入第三方模块
import matplotlib
import pandas as pd
import matplotlib.pyplot as plt
import numpy as np
from scipy import stats
# % matplotlib inline
# 创建数据
from scipy.interpolate import lagrange

fdata1 = "D21.txt"
data = pd.Series(np.loadtxt(fdata1, dtype=np.float64))
u = np.mean(data)  # 计算均值
std = np.std(data)  # 计算标准差
print(u,std)
stats.kstest(data, 'norm', (u, std))  # 正态性检验
plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签

# 绘制数据密度曲线
# 绘制数据密度曲线
fig = plt.figure(figsize=(10, 6))

ax1 = fig.add_subplot(2, 1, 1)

data.plot(kind='kde', style='--k', grid=True, title='密度曲线')

plt.axvline(3*std, linestyle='--', color='r')

plt.axvline(-3*std, linestyle='--', color='r')

# 筛选出异常值和正常值
error = data[np.abs(data - u) > 3*std]
d