軒田機器學習基石課程學習筆記11 — Linear Models for Classification

红色石头發表於2018-07-24

上一節課,我們介紹了Logistic Regression問題,建立cross-entropy error,並提出使用梯度下降演算法gradient descnt來獲得最好的logistic hypothesis。本節課繼續介紹使用線性模型來解決分類問題。

一、Linear Models for Binary Classification

之前介紹幾種線性模型都有一個共同點,就是都有樣本特徵x的加權運算,我們引入一個線性得分函式s:

s=w^Tx

三種線性模型,第一種是linear classification。線性分類模型的hypothesis為h(x)=sign(s),取值範圍為{-1,+1}兩個值,它的err是0/1的,所以對應的E_{in}(w)是離散的,並不好解,這是個NP-hard問題。第二種是linear regression。線性迴歸模型的hypothesis為h(x)=s,取值範圍為整個實數空間,它的err是squared的,所以對應的E_{in}(w)是開口向上的二次曲線,其解是closed-form的,直接用線性最小二乘法求解即可。第三種是logistic regression。邏輯迴歸模型的hypothesis為h(x)=\theta(s),取值範圍為(-1,1)之間,它的err是cross-entropy的,所有對應的E_{in}(w)是平滑的凸函式,可以使用梯度下降演算法求最小值。

從上圖中,我們發現,linear regression和logistic regression的error function都有最小解。那麼可不可以用這兩種方法來求解linear classification問題呢?下面,我們來對這三種模型的error function進行分析,看看它們之間有什麼聯絡。

對於linear classification,它的error function可以寫成:

err_{0/1}(s,y)=|sign(s)\neq y|=|sign(ys)\neq 1|

對於linear regression,它的error function可以寫成:

err_{SQR}(s,y)=(s-y)^2=(ys-1)^2

對於logistic regression,它的error function可以寫成:

err_{CE}(s,y)=ln(1+exp(-ys))

上述三種模型的error function都引入了ys變數,那麼ys的物理意義是什麼?ys就是指分類的正確率得分,其值越大越好,得分越高。

下面,我們用圖形化的方式來解釋三種模型的error function到底有什麼關係:

從上圖中可以看出,ys是橫座標軸,err_{0/1}是呈階梯狀的,在ys>0時,err_{0/1}恆取最小值0。err_{SQR}呈拋物線形式,在ys=1時,取得最小值,且在ys=1左右很小區域內,err_{0/1}err_{SQR}近似。err_{CE}是呈指數下降的單調函式,ys越大,其值越小。同樣在ys=1左右很小區域內,err_{0/1}err_{CE}近似。但是我們發現err_{CE}並不是始終在err_{0/1}之上,所以為了計算討論方便,我們把err_{CE}做幅值上的調整,引入err_{SCE}=log_2(1+exp(-ys))=\frac1{ln2}err_{CE},這樣能保證err_{SCE}始終在err_{0/1}上面,如下圖所示:

由上圖可以看出:

err_{0/1}(s,y)\leq err_{SCE}(s,y)=\frac1{ln2}err_{CE}(s,y)

E_{in}^{0/1}(w)\leq E_{in}^{SCE}(w)=\frac1{ln2}E_{in}^{CE}(w)

E_{out}^{0/1}(w)\leq E_{out}^{SCE}(w)=\frac1{ln2}E_{out}^{CE}(w)

那麼由VC理論可以知道:
從0/1出發:

E_{out}^{0/1}(w)\leq E_{in}^{0/1}(w)+\Omega^{0/1}\leq \frac1{ln2}E_{in}^{CE}(w)+\Omega^{0/1}

從CE出發:

E_{out}^{0/1}(w)\leq \frac1{ln2}E_{out}^{CE}(w)\leq \frac1{ln2}E_{in}^{CE}(w)+\frac1{ln2}\Omega^{CE}

通過上面的分析,我們看到err 0/1是被限定在一個上界中。這個上界是由logistic regression模型的error function決定的。而linear regression其實也是linear classification的一個upper bound,只是隨著sy偏離1的位置越來越遠,linear regression的error function偏差越來越大。綜上所述,linear regression和logistic regression都可以用來解決linear classification的問題。

下圖列舉了PLA、linear regression、logistic regression模型用來解linear classification問題的優點和缺點。通常,我們使用linear regression來獲得初始化的w_0,再用logistic regression模型進行最優化解。

二、Stochastic Gradient Descent

之前介紹的PLA演算法和logistic regression演算法,都是用到了迭代操作。PLA每次迭代只會更新一個點,它每次迭代的時間複雜度是O(1);而logistic regression每次迭代要對所有N個點都進行計算,它的每時間複雜度是O(N)。為了提高logistic regression中gradient descent演算法的速度,可以使用另一種演算法:隨機梯度下降演算法(Stochastic Gradient Descent)。

隨機梯度下降演算法每次迭代只找到一個點,計算該點的梯度,作為我們下一步更新w的依據。這樣就保證了每次迭代的計算量大大減小,我們可以把整體的梯度看成這個隨機過程的一個期望值。

隨機梯度下降可以看成是真實的梯度加上均值為零的隨機噪聲方向。單次迭代看,好像會對每一步找到正確梯度方向有影響,但是整體期望值上看,與真實梯度的方向沒有差太多,同樣能找到最小值位置。隨機梯度下降的優點是減少計算量,提高運算速度,而且便於online學習;缺點是不夠穩定,每次迭代並不能保證按照正確的方向前進,而且達到最小值需要迭代的次數比梯度下降演算法一般要多。

對於logistic regression的SGD,它的表示式為:

w_{t+1}\leftarrow w_t+\eta\theta(-y_nw_t^Tx_n)(y_nx_n)

我們發現,SGD與PLA的迭代公式有類似的地方,如下圖所示:

我們把SGD logistic regression稱之為’soft’ PLA,因為PLA只對分類錯誤的點進行修正,而SGD logistic regression每次迭代都會進行或多或少的修正。另外,當\eta=1,且w_t^Tx_n足夠大的時候,PLA近似等於SGD。

除此之外,還有兩點需要說明:1、SGD的終止迭代條件。沒有統一的終止條件,一般讓迭代次數足夠多;2、學習速率\eta\eta的取值是根據實際情況來定的,一般取值0.1就可以了。

三、Multiclass via Logistic Regression

之前我們一直講的都是二分類問題,本節主要介紹多分類問題,通過linear classification來解決。假設平面上有四個類,分別是正方形、菱形、三角形和星形,如何進行分類模型的訓練呢?

首先我們可以想到這樣一個辦法,就是先把正方形作為正類,其他三種形狀都是負類,即把它當成一個二分類問題,通過linear classification模型進行訓練,得出平面上某個圖形是不是正方形,且只有{-1,+1}兩種情況。然後再分別以菱形、三角形、星形為正類,進行二元分類。這樣進行四次二分類之後,就完成了這個多分類問題。

但是,這樣的二分類會帶來一些問題,因為我們只用{-1,+1}兩個值來標記,那麼平面上某些可能某些區域都被上述四次二分類模型判斷為負類,即不屬於四類中的任何一類;也可能會出現某些區域同時被兩個類甚至多個類同時判斷為正類,比如某個區域又判定為正方形又判定為菱形。那麼對於這種情況,我們就無法進行多類別的準確判斷,所以對於多類別,簡單的binary classification不能解決問題。

針對這種問題,我們可以使用另外一種方法來解決:soft軟性分類,即不用{-1,+1}這種binary classification,而是使用logistic regression,計算某點屬於某類的概率、可能性,去概率最大的值為那一類就好。

soft classification的處理過程和之前類似,同樣是分別令某類為正,其他三類為負,不同的是得到的是概率值,而不是{-1,+1}。最後得到某點分別屬於四類的概率,取最大概率對應的哪一個類別就好。效果如下圖所示:

這種多分類的處理方式,我們稱之為One-Versus-All(OVA) Decomposition。這種方法的優點是簡單高效,可以使用logistic regression模型來解決;缺點是如果資料類別很多時,那麼每次二分類問題中,正類和負類的數量差別就很大,資料不平衡unbalanced,這樣會影響分類效果。但是,OVA還是非常常用的一種多分類演算法。

四、Multiclass via Binary Classification

上一節,我們介紹了多分類演算法OVA,但是這種方法存在一個問題,就是當類別k很多的時候,造成正負類資料unbalanced,會影響分類效果,表現不好。現在,我們介紹另一種方法來解決當k很大時,OVA帶來的問題。

這種方法呢,每次只取兩類進行binary classification,取值為{-1,+1}。假如k=4,那麼總共需要進行C_4^2=6次binary classification。那麼,六次分類之後,如果平面有個點,有三個分類器判斷它是正方形,一個分類器判斷是菱形,另外兩個判斷是三角形,那麼取最多的那個,即判斷它屬於正方形,我們的分類就完成了。這種形式就如同k個足球對進行單迴圈的比賽,每場比賽都有一個隊贏,一個隊輸,贏了得1分,輸了得0分。那麼總共進行了C_k^2次的比賽,最終取得分最高的那個隊就可以了。

這種區別於OVA的多分類方法叫做One-Versus-One(OVO)。這種方法的優點是更加高效,因為雖然需要進行的分類次數增加了,但是每次只需要進行兩個類別的比較,也就是說單次分類的數量減少了。而且一般不會出現資料unbalanced的情況。缺點是需要分類的次數多,時間複雜度和空間複雜度可能都比較高。

五、總結

本節課主要介紹了分類問題的三種線性模型:linear classification、linear regression和logistic regression。首先介紹了這三種linear models都可以來做binary classification。然後介紹了比梯度下降演算法更加高效的SGD演算法來進行logistic regression分析。最後講解了兩種多分類方法,一種是OVA,另一種是OVO。這兩種方法各有優缺點,當類別數量k不多的時候,建議選擇OVA,以減少分類次數。

註明:

文章中所有的圖片均來自臺灣大學林軒田《機器學習基石》課程

更多AI資源請關注公眾號:AI有道(ID:redstonewill)

相關文章