Binary classification - 聊聊评价指标的那些事儿【实战篇】
分類問題就像披著羊皮的狼,看起來天真無害用起來天雷滾滾。比如在建模前你思考過下面的問題么?
- 你的分類模型輸出的概率只是用來做樣本間的相對排序,還是概率本身?
- 你的訓(xùn)練數(shù)據(jù)本身分布如何是否存在Imbalanced Sample?
要是您都想到了拜拜?。要是有1各您感興趣的問題,那就接著往下看吧。本來是想先回顧一下各個分類問題中可能用到的metric,但是又覺得讀的人可能覺得無聊,就分成了兩章。要是有的指標(biāo)斷片了就來這里回憶一下: 回憶篇
問題1 Rank or Probability?
分類問題可以根據(jù)對輸出形式的要求分成兩類
有人會問,上述兩種需求究竟對解決一個二分類問題有什么影響? 答案是損失函數(shù)/評價指標(biāo)
讓我們來看一個直觀的例子,下圖我們嘗試用LightGBM解決一個二分類問題,我們選擇的擬合指標(biāo)是最大化AUC。
X軸是預(yù)測概率,Y軸是真實(shí)概率,藍(lán)線是LGB的預(yù)測結(jié)果,綠線對應(yīng)真實(shí)概率=預(yù)測概率。為什么模型的AUC高達(dá)98.93%(這里還有ImbalancedSample的影響,讓我們先忽略這一點(diǎn)),但是預(yù)測概率和真實(shí)概率卻差到了姥姥家。
讓我們對預(yù)測概率再做一層處理,黃線可以簡單理解為我們對LGB的預(yù)測結(jié)果做了一層映射 \(\hat{p} \to f(\hat{p})\),這時校準(zhǔn)后的預(yù)測概率和真實(shí)概率基本一致了。但是有趣的是校準(zhǔn)后的預(yù)測概率AUC = 98.94%和原始預(yù)測基本沒差別?!
Duang Duang Duang!敲黑板!AUC是對相對概率排序的檢驗!其實(shí)只要用心(我學(xué)AUC的時候一定沒用心>_<)看一下AUC的計算方式就會發(fā)現(xiàn),AUC只關(guān)心給定各個閾值,把樣本按照預(yù)測概率分成0/1,并計算正負(fù)樣本預(yù)測的準(zhǔn)確率。
舉個最夏天的例子,兩個瓜一個甜一個不甜,我們訓(xùn)練一個西瓜模型來預(yù)測它們甜(1)/不甜(0)。
模型1: 甜的瓜預(yù)測概率是0.8,不甜的瓜預(yù)測概率是0.1,
模型2: 甜的瓜預(yù)測概率是0.51,不甜的瓜預(yù)測概率是0.49
兩個模型的AUC是相同的,因為它們都完美對兩個瓜進(jìn)行了分類。
所以當(dāng)使用最大化AUC作為損失函數(shù)時,當(dāng)正負(fù)樣本的預(yù)測準(zhǔn)確率不再提高,模型就會停止學(xué)習(xí)。這時模型的預(yù)測概率并不是對真實(shí)概率的擬合。那如何才能得到對真實(shí)概率的預(yù)測? 答案是logloss/cros-entropy
\[ \begin{align} L &= \sum_{i=1}^N y_i * log(p_i) + (1-y_i) *log(1-p_i)\\ \end{align} \]
我們可以從兩個角度來理解為什么logloss是對真實(shí)概率的估計
從極大似然估計的角度
logloss可以由極大似然函數(shù)取對數(shù)得到,最小化logloss對應(yīng)的最大化似然函數(shù)。\(p_i\)是對\(p(y_i=1)\)的估計
\[ argmax_p \prod_{i=1}^N {p_i}^{y_i} * {(1-p_i)}^{1-y_i} \]
從信息論的角度
不熟悉信息論的同學(xué)看這里 Intro to Information Theory
logloss也叫cross-entropy(交叉熵),用來衡量兩個分布的相似程度。
交叉熵本身可以分解為P本身的信息熵+分布P和分布q之間的距離。這里P是樣本的真實(shí)分布信息,信息熵一定。所以最小化交叉熵就變成了最小化分布p和q之間的距離,也就是樣本分布和模型估計間的距離,如下
\[ \begin{align} crossentropy &= H(p,q)\\ &= -\sum_{c=1}^C p(c) * log(q(c))\\ & = - \sum_{c=1}^C p(c) * log(p(c)) + \sum_{c=1}^C p(c)[log(p(c) - log(q(c)))] \\ &= H(p) + KL(p||q)\\ \end{align} \]
乍一看會覺得交叉熵和logloss長的不像一家人。因為在訓(xùn)練模型時分布p是從訓(xùn)練樣本的分布中抽象得到的。二分類問題中C=2, 讓我們把上述交叉熵再推一步
\[ \begin{align} H(p,q) &= p *log(q) + (1-p) *log(1-q) \\ p& = \sum_{i=1}^N I(y_i=1)/N \\ H(p,q) &= \frac{1}{N} \sum_i I(y_i=1) *log(q)+ I(y_i=0) *log(1-q) \\ \end{align} \]
所以下次解決分類問題,如果你的目標(biāo)是計算對真實(shí)概率的估計的話,別選錯指標(biāo)喲�
興趣卡片- 預(yù)測概率校準(zhǔn)
其實(shí)黃線用了Isotonic Regression來校準(zhǔn)預(yù)測概率。是一種事后將預(yù)測概率根據(jù)真實(shí)概率進(jìn)行校準(zhǔn)的方法。感興趣的可以看一下Reference里面的材料1,2。原理并不復(fù)雜,但在分析特定算法,尤其是boosting,bagging類的集合算法為什么使用loggloss對概率估計依舊會有偏的部分蠻有趣的
問題2 Imbalanced Sample ?
正負(fù)樣本分布不均大概是分類問題中最常遇到的問題。正確解決Imbalane問題需要注意的并不只是評價指標(biāo),往往還要注意采樣和訓(xùn)練集測試集的劃分。但這里我們只討論在解決樣本分布不均的問題時,我們應(yīng)該選擇什么指標(biāo)來評價模型表現(xiàn)。讓我們挨個來剔除不好用的指標(biāo)。
舉個極端的例子,100個樣本里只有1個正樣本
Accuracy
這種情況下即便我們?nèi)款A(yù)測為負(fù),我們的準(zhǔn)確率依舊高達(dá)99%。所以Accuracy只適用于正負(fù)樣本均勻分布的情況,因為它把正負(fù)樣本的預(yù)測準(zhǔn)確率柔和在一起看了。
AUC
AUC是fpr和tpr(recall)組成的ROC的曲線下面積。還記得我們在【回憶篇】里面說過fpr,tpr是分別衡量在正負(fù)樣本上的準(zhǔn)確率的。
而fpr和tpr之間的trade-off,在正樣本占比很小的情況下,這種trad-off會被樣本量更大的一方主導(dǎo)。所以當(dāng)正樣本占比很小的時候,AUC往往會看起來過于優(yōu)秀。
但就像硬幣的正反面一樣,從另一個角度看這也是AUC的優(yōu)點(diǎn),就是AUC本身不會很大的受到樣本實(shí)際分布的影響,相同的模型相同的樣本,你把正樣本downsample /upsample 1倍,AUC不會有很大的改變。
下圖來自An introduction to ROC analysis, 上面的AUC和PR是正負(fù)樣本1:1的預(yù)測表現(xiàn),下面是1:10的表現(xiàn)。我們會發(fā)現(xiàn)AUC基本沒有變化,但是precision-recall發(fā)生了劇烈變化。
AP/AUCPR
AP是recall和precision組成的PR的曲線下面積。這里recall和precision分別從真實(shí)分布和預(yù)測分布兩個角度衡量了對正樣本的預(yù)測準(zhǔn)確率。說到這里已經(jīng)有人反應(yīng)過來了。是的這一對trade-off指標(biāo)都是針對正樣本的,在計算中沒有用到True negative.所以當(dāng)你的數(shù)據(jù)集存在Imbalance的時候,AP一般會是更好的選擇。
...你還遇到過啥問題嘞?歡迎留言
Reference
轉(zhuǎn)載于:https://www.cnblogs.com/gogoSandy/p/11123688.html
總結(jié)
以上是生活随笔為你收集整理的Binary classification - 聊聊评价指标的那些事儿【实战篇】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: oracle函数 trunc(x[,y
- 下一篇: oracle自动选择索引