使用Python進(jìn)行異常檢測的解決方案
異常檢測是異常值分析中的一項統(tǒng)計任務(wù),但是如果我們開發(fā)一個機(jī)器學(xué)習(xí)模型來自動化地進(jìn)行異常檢測,可以節(jié)省很多時間。
異常檢測有很多用例,包括信用卡欺詐檢測、故障機(jī)器檢測、基于異常特征的硬件系統(tǒng)檢測、基于醫(yī)療記錄的疾病檢測都是很好的例子,除此之外也還有很多的用例。在本文中,我們將使用Python從頭開始實(shí)現(xiàn)異常檢測算法。公式和過程與我之前解釋過的其他機(jī)器學(xué)習(xí)算法相比,我們使用的異常檢測算法要簡單得多。該算法使用均值和方差來計算每個訓(xùn)練數(shù)據(jù)的概率。如果一個訓(xùn)練實(shí)例的概率很高,則是正常的;如果某個訓(xùn)練實(shí)例的概率很低,那就可以是一個異常樣本。對于不同的訓(xùn)練集,高概率和低概率的定義是不同的,這個我們以后再討論。接下來我們來看一下異常檢測的工作過程。使用以下公式計算平均值:
這里m是數(shù)據(jù)集的長度或訓(xùn)練數(shù)據(jù)的數(shù)量,而是一個單獨(dú)的訓(xùn)練例子。如果你有多個訓(xùn)練特征,大多數(shù)情況下都需要計算每個特征的平均值。使用以下公式計算方差:
這里,mu是上一步計算的平均值。現(xiàn)在,用這個概率公式來計算每個訓(xùn)練例子的概率。
不要被這個公式中的求和符號弄糊涂了!這實(shí)際上是Sigma方差。稍后我們實(shí)現(xiàn)該算法時,就會理解它了,F(xiàn)在我們需要找到概率的臨界值。正如我前面提到的,如果一個訓(xùn)練例子的概率很低,那這就是一個異常樣本。多大是低概率呢這沒有統(tǒng)一的標(biāo)準(zhǔn),我們需要為我們的訓(xùn)練數(shù)據(jù)集找出這個閾值。我們從步驟3中得到的輸出中獲取一系列概率值,對于每個概率,通過閾值的設(shè)置來判斷數(shù)據(jù)是否異常然后計算一系列概率的精確度、召回率和f1分?jǐn)?shù)。精度可使用以下公式計算
召回率的計算公式如下:
在這里,True positives(真正例)是指是異常的且算法檢測到也是異常的樣本。False Positives(假正例)是指不是異常的但算法檢測到是異常的樣本。False Negative(假反例)是指不是異常的且算法檢測到也不是異常的樣本。從上面的公式你可以看出,更高的精確度和更高的召回率說明算法性能更好,因?yàn)檫@意味著我們有更多的真正的正例,但同時,假正例和假反例也起著至關(guān)重要的作用,這需要一個平衡點(diǎn),根據(jù)你的行業(yè),你需要決定哪一個對你來說是可以忍受的。一個好辦法是取平均數(shù)。計算平均值有一個獨(dú)特的公式,這就是f1分?jǐn)?shù),f1得分公式為:
這里,P和R分別表示精確性和召回率。如果你對該公式感興趣的話,可以查看:https://towardsdatascience.com/a-complete-understanding-of-precision-recall-and-f-score-concepts-23dc44defef6根據(jù)f1分?jǐn)?shù),你需要選擇你的閾值概率。異常檢測算法我們將使用Andrew Ng的機(jī)器學(xué)習(xí)課程的數(shù)據(jù)集,它具有兩個訓(xùn)練特征。我沒有在本文中使用真實(shí)的數(shù)據(jù)集,因?yàn)檫@個數(shù)據(jù)集非常適合學(xué)習(xí),它只有兩個特征。在任何真實(shí)的數(shù)據(jù)集中,都不可能只有兩個特征。有兩個特性的好處是可以可視化數(shù)據(jù),這對學(xué)習(xí)者非常有用。請從該鏈接下載數(shù)據(jù)集:https://github.com/rashida048/Machine-Learning-With-Python/blob/master/ex8data1.xlsx首先,導(dǎo)入必要的包import pandas as pd
import numpy as np
導(dǎo)入數(shù)據(jù)集。這是一個excel數(shù)據(jù)集。在這里,訓(xùn)練數(shù)據(jù)和交叉驗(yàn)證數(shù)據(jù)存儲在單獨(dú)的表中。df = pd.read_excel('ex8data1.xlsx', sheet_name='X', header=None)
df.head()
讓我們將第0列與第1列進(jìn)行比較。plt.figure()
plt.scatter(df[0], df[1])
plt.show()
你可能通過看這張圖知道哪些數(shù)據(jù)是異常的。檢查此數(shù)據(jù)集中有多少個訓(xùn)練示例:m = len(df)

發(fā)表評論
請輸入評論內(nèi)容...
請輸入評論/評論長度6~500個字
最新活動更多
-
7月22-29日立即報名>> 【線下論壇】第三屆安富利汽車生態(tài)圈峰會
-
7.30-8.1火熱報名中>> 全數(shù)會2025(第六屆)機(jī)器人及智能工廠展
-
7月31日免費(fèi)預(yù)約>> OFweek 2025具身智能機(jī)器人產(chǎn)業(yè)技術(shù)創(chuàng)新應(yīng)用論壇
-
免費(fèi)參會立即報名>> 7月30日- 8月1日 2025全數(shù)會工業(yè)芯片與傳感儀表展
-
即日-2025.8.1立即下載>> 《2024智能制造產(chǎn)業(yè)高端化、智能化、綠色化發(fā)展藍(lán)皮書》
-
8月5日立即報名>> 【在線會議】CAE優(yōu)化設(shè)計:醫(yī)療器械設(shè)計的應(yīng)用案例與方案解析
推薦專題