快速找到離群值的三種方法
本文將介紹3個在數(shù)據(jù)集中查找離群值的Python方法。
離群值(Outliers)是指在數(shù)據(jù)集中與其他數(shù)據(jù)點明顯不同或者異常的數(shù)據(jù)點。這些數(shù)據(jù)點可能比其他數(shù)據(jù)點要遠離數(shù)據(jù)集的中心,或者具有異常的數(shù)值。離群值可能是由于數(shù)據(jù)采集錯誤、異常事件、測量誤差或者其他未知因素引起的。
離群值的存在可以對數(shù)據(jù)分析和統(tǒng)計建模產(chǎn)生重要影響,因為它們可能導(dǎo)致模型不準(zhǔn)確或者產(chǎn)生誤導(dǎo)性的結(jié)果。
我們先創(chuàng)建一個演示的數(shù)據(jù)
import pandas as pd
import matplotlib.pyplot as plt
name = ['John', 'Victor', 'Carlos', 'Leo', 'Kevin', 'Silva', 'Johnson', 'Lewis', 'George', 'Daniel', 'Harry', 'Jordan', 'James']
salary = [4000, 1000, 2000, 100000, 3500, 6000, 1500, 3000, 2500, 3600, 2100, 1700, 1600]
df = pd.DataFrame({'Name': name, 'Salary': salary})
plt.boxplot(df['Salary'])
plt.show()
可以看到上面的點就是離群值,下面我們將介紹快速找到它的方法。
四分位極差法
首先找到第一和第三個四分位數(shù)值,通常記為Q1和Q3。然后用Q3減去Q1計算四分位差(IQR)。
通過減去/增加1.5倍IQR來計算下界和上界。上下邊界外的值就是離群值
q1 = df['Salary'].quantile(0.25)
q3 = df['Salary'].quantile(0.75)
iqr = q3 - q1
lower_bound = q1 - 1.5 * iqr
upper_bound = q3 + 1.5 * iqr
outliers = df[(df['Salary'] < lower_bound) | (df['Salary'] > upper_bound)]
1.5倍IQR的閾值通常是一種常用的標(biāo)準(zhǔn),閾值的選擇可以根據(jù)具體情況進行調(diào)整。有時候,也可以選擇更嚴(yán)格或更寬松的閾值,以適應(yīng)特定的數(shù)據(jù)分析需求。
標(biāo)準(zhǔn)偏差法
標(biāo)準(zhǔn)偏差法(Standard Deviation Method)使用數(shù)據(jù)的標(biāo)準(zhǔn)差來判斷數(shù)據(jù)點是否偏離了數(shù)據(jù)。上界和下界是均值和3倍標(biāo)準(zhǔn)差的加減。
他的方法如下:
計算平均值和標(biāo)準(zhǔn)偏差: 首先,計算數(shù)據(jù)的平均值(Mean)和標(biāo)準(zhǔn)偏差(Standard Deviation)。平均值代表了數(shù)據(jù)的中心位置,標(biāo)準(zhǔn)偏差衡量了數(shù)據(jù)的分散程度。
確定閾值: 定義一個閾值,通常是標(biāo)準(zhǔn)偏差的倍數(shù)(通常為2或3倍標(biāo)準(zhǔn)偏差)。這個閾值決定了什么樣的數(shù)據(jù)點被認(rèn)為是離群值。
識別離群值: 計算每個數(shù)據(jù)點與平均值之間的差值,然后將這個差值與閾值比較。如果差值超過了閾值,數(shù)據(jù)點被認(rèn)為是離群值。
mean = df.Salary.mean()
std = df.Salary.std()
upper_bound = mean + 3 * std
lower_bound = mean - 3 * std
outliers = df[(df['Salary'] < lower_bound) | (df['Salary'] > upper_bound)]
標(biāo)準(zhǔn)偏差法的優(yōu)點在于簡單易懂,而且不需要假設(shè)數(shù)據(jù)分布的形狀。但需要注意以下幾點:
- 通常情況下,閾值使用2或3倍標(biāo)準(zhǔn)偏差作為閾值,但這個值可能需要根據(jù)具體情況進行調(diào)整。
- 這種方法對于正態(tài)分布的數(shù)據(jù)集效果較好,但對于偏斜分布的數(shù)據(jù),可能會導(dǎo)致誤判。
- 標(biāo)準(zhǔn)偏差法可能不適用于小樣本,因為標(biāo)準(zhǔn)偏差在小樣本中可能不夠穩(wěn)定。
Z-分?jǐn)?shù)法
Z-分?jǐn)?shù)(Z-Score)法測量了數(shù)據(jù)點與數(shù)據(jù)集平均值之間的偏差,以標(biāo)準(zhǔn)化方式表示這個偏差。對于每個數(shù)據(jù)點,計算它與平均值之間的差值,然后將這個差值除以標(biāo)準(zhǔn)偏差,得到Z-分?jǐn)?shù)。如果z-score大于3.0或小于-3.0,則該值可歸類為離群值。
我們可以直接使用scipy提供的函數(shù)來進行計算
from scipy import stats
df['Salary_zscore'] = stats.zscore(df['Salary'])
filtered_df = df[(df['Salary_zscore'] <= 3) & (df['Salary_zscore'] >= -3)]
Z-分?jǐn)?shù)法適用于各種類型的數(shù)據(jù)分布,不需要假設(shè)數(shù)據(jù)分布的形狀。并且提供了標(biāo)準(zhǔn)化的度量,使得不同數(shù)據(jù)集之間的離群值比較更加容易。
總結(jié)
以上是可以快速找到離群值的統(tǒng)計學(xué)方法,除此以外,還有一些機器學(xué)習(xí)的方法例如:
DBSCAN(Density-Based Spatial Clustering of Applications with Noise): DBSCAN是一種密度聚類算法,也可用于檢測離群值。它根據(jù)數(shù)據(jù)點的密度來識別離群值,將密度較低的點視為離群值。
LOF(Local Outlier Factor): LOF是一種局部離群值因子方法,用于檢測局部區(qū)域內(nèi)的離群值。它考慮了每個數(shù)據(jù)點周圍的局部密度與相鄰點的密度之間的比率,從而識別離群值。
Isolation Forest: Isolation Forest是一種基于隨機森林的離群值檢測方法,它通過構(gòu)建樹結(jié)構(gòu)來識別離群值。由于使用了隨機性,它對高維數(shù)據(jù)和大數(shù)據(jù)集非常有效。
但是這些方法執(zhí)行的速度會很慢,如果對于速度要求比較嚴(yán)格還是需要謹(jǐn)慎選擇。