piaip 的 (lib)SVM 簡易入門

最新推荐文章于 2017-06-04 13:01:40 发布

chvalrous

最新推荐文章于 2017-06-04 13:01:40 发布

阅读量813

点赞数

分类专栏： Data mining 文章标签： svm

Data mining 专栏收录该内容

21 篇文章 0 订阅

订阅专栏

本文转载自http://ntu.csie.org/~piaip/svm/svm_tutorial.html#

piaip at csie dot ntu dot edu dot tw,
Hung-Te Lin
Fri Apr 18 15:04:53 CST 2003
$Id: svm_tutorial.html,v 1.13 2007/10/02 05:51:55 piaip Exp piaip $原作：林弘德，轉載請保留原出處

Why this tutorial is here

我一直覺得 SVM 是個很有趣的東西，不過也一直沒辦法 (mostly 衝堂) 去聽林智仁老師的 Data mining 跟 SVM 的課；後來看了一些網路上的文件跟聽 kcwu 講了一下 libsvm 的用法後，就想整理一下，算是對於並不需要知道完整 SVM 理論的人提供使用 libsvm 的入門。原始 libsvm 的 README 跟 FAQ 也是很好的文件，不過你可能要先對 svm 跟流程有點了解才看得懂 (我在看時有這樣的感覺)；這篇入門就是為了從零開始的人而寫的。

後來還有一些人提供意見，所以在此要感謝：

kcwu, biboshen, puffer, somi

不過請記得底下可能有些說法不一定對，但是對於只是想用 SVM 的人來說我覺得這樣說明會比較易懂。

這篇入門原則上是給會寫基本程式的人看的，也是給我自己一個備忘, 不用太多數學底子，也不用對 SVM 有任何先備知識。

還看不懂的話有三個情形, 一是我講的不夠清楚, 二是你的常識不足, 三是你是小白 ^^;

我自己是以完全不懂的角度開始的，這篇入門也有不少一樣不懂 SVM 的人看過、而且看完多半都有一定程度的理解，所以假設情況一不會發生，那如果不懂一定是後兩個情況 :P 也所以, 有問題別問我。

SVM: What is it and what can it do for me?

SVM, Support Vector Machine , 簡而言之它是個起源跟類神經網路有點像的東西，不過現今最常拿來就是做分類 (classification)。也就是說，如果我有一堆已經分好類的東西（可是分類的依據是未知的！），那當收到新的東西時， SVM 可以預測 (predict) 新的資料要分到哪一堆去。

聽起來是很神奇的事（如果你覺得不神奇，請重想一想這句話代表什麼：分類的依據是未知的！，還是不神奇的話就請你寫個程式解解看這個問題），也很像要 AI 之類的高等技巧... 不過 SVM 基於統計學習理論可以在合理的時間內漂亮的解決這個問題。

以圖形化的例子來說明(by SVMToy), 像假定我在空間中標了一堆用顏色分類的點, 點的顏色就是他的類別, 位置就是他的資料, 那 SVM 就可以找出區隔這些點的方程式, 依此就可以分出一區區的區域; 拿到新的點(資料) 時, 只要對照該位置在哪一區就可以(predict) 找出他應該是哪一顏色(類別)了:

原始資料分佈	SVM找出來的區域

當然 SVM 不是真的只有畫圖分區那麼簡單, 不過看上面的例子應該可以了解 SVM 大概在作什麼.

要對 SVM 再多懂一點點，可以參考 cjlin 在 data mining 課的 slides: pdf or ps 。
底下我試著在不用看那個 slide 的情況解釋及使用 libsvm。

所以, 我們可以把 SVM 當個黑盒子, 資料丟進去讓他處理然後我們再來用就好了.

How do I get SVM?

林智仁(cjlin)老師的 libsvm 當然是最完美的工具.

Download libsvm

下載處:

libsvm.zip or libsvm.tar.gz

.zip 跟 .tar.gz 基本上是一樣的, 只是看你的 OS; 習慣上 Windows 用 .zip 比較方便 (因為有 WinZIP, 不過我都用 WinRAR), UNIX 則是用 .tar.gz

Build libsvm

解開來後, 假定是 UNIX 系統, 直接打 make 就可以了; 編不出來的話請詳讀說明和運用常識. 因為這是 tutorial, 所以我不花時間細談, 而且會編不出來的情形真是少之又少, 通常一定是你的系統有問題或你太笨了. 其他的子目錄可以不管, 只要 svm-train, svm-scale, svm-predict 三個執行檔有編出來就可以了.

Windows 的用戶要自己重編當然也是可以, 不過已經有編好的 binary 在裡面了: 請檢查 windows 子目錄, 應該會有 svmtrain.exe, svmscale.exe, svmpredict.exe, svmtoy.exe .

Using SVM

libsvm 有很多種用法, 這篇 tutorial 只打算講簡單的部分.

The programs

解釋一下幾個主要執行檔的作用: (UNIX/Windows 下檔名稍有不同, 請用常識理解我在講哪個)

svmtrain

Train (訓練) data. 跑 SVM 被戲稱為 "開火車" 也是由於這個程式名而來. train 會接受特定格式的輸入, 產生一個 "Model" 檔. 這個 model 你可以想像成 SVM 的內部資料, 因為 predict 要 model 才能 predict, 不能直接吃原始資料. 想想也很合理, 假定 train 本身是很耗時的動作, 而 train 好可以以某種形式存起內部資料, 那下次要 predict 時直接把那些內部資料 load 進來就快多了.

svmpredict

依照已經 train 好的 model, 再加上給定的輸入 (新值), 輸出 predict (預測) 新值所對應的類別 (class).

svmscale

Rescale data. 因為原始資料可能範圍過大或過小, svmscale 可以先將資料重新 scale (縮放) 到適當範圍.

File Format

檔案格式要先交代一下. 你可以參考 libsvm 裡面附的 "heart_scale":

[label] [index1]:[value1] [index2]:[value2] ... [label] [index1]:[value1] [index2]:[value2] ... . .

一行一筆資料，如

+1 1:0.708 2:1 3:1 4:-0.320 5:-0.105 6:-1

label

或說是 class, 就是你要分類的種類，通常是一些整數。

index

是有順序的索引，通常是放連續的整數。

value

就是用來 train 的資料，通常是一堆實數。

每一行都是如上的結構, 意思就是: 我有一排資料, 分別是 value1, value2, .... valueN, (而且它們的順序已由 indexN 分別指定)，這排資料的分類結果就是 label。

或許你會不太懂，為什麼會是 value1,value2,.... 這樣一排呢？這牽涉到 SVM 的原理。你可以這樣想（我沒說這是正確的），它的名字就叫 Support "Vector" Machine，所以輸入的 training data 是 "Vector"(向量), 也就是一排的 x1, x2, x3, ... 這些值就是 valueN，而 x[n] 的 n 就是由 indexN 指定。這些東西又稱為 "attribute"。

真實的情況是，大部份時候我們給定的資料可能有很多 "特徵(feature)" 或說 "屬性(attribute)"，所以輸入會是一組的。舉例來說，以前面畫點分區的例子來說，我們不是每個點都有 X 跟 Y 的座標嗎？所以它就有兩種 attribute。假定我有兩個點： (0,3) 跟 (5,8) 分別在 label(class) 1 跟 2 ，那就會寫成1 1:0 2:3 2 1:5 2:8同理，空間中的三維座標就等於有三組 attribute。

這種檔案格式最大的好處就是可以使用 sparse matrix，或說有些 data 的 attribute 可以不存在。

To Run libsvm

來解釋一下 libsvm 的程式怎麼用。你可以先拿 libsvm 附的 heart_scale 來做輸入，底下也以它為例：

看到這裡你應該也了解，使用 SVM 的流程大概就是：

準備資料並做成指定格式 (有必要時需 svmscale)
用 svmtrain 來 train 成 model
對新的輸入，使用 svmpredict 來 predict 新資料的 class

svmtrain

svmtrain 的語法大致就是:

svmtrain [options] training_set_file [model_file]

training_set_file 就是之前的格式，而 model_file 如果不給就會叫 [training_set_file].model。 options 可以先不要給。

下列程式執行結果會產生 heart_scale.model 檔：(螢幕輸出不是很重要，沒有錯誤就好了)

./svm-train heart_scale optimization finished, #iter = 219 nu = 0.431030 obj = -100.877286, rho = 0.424632 nSV = 132, nBSV = 107 Total nSV = 132

svmpredict

svmpredict 的語法是 :

svmpredict test_file model_file output_file

test_file 就是我們要 predict 的資料。它的格式跟 svmtrain 的輸入，也就是 training_set_file 是一樣的！ predict 完會順便拿 predict 出來的值跟 test_file 裡面寫的值去做比對，這代表： test_file 寫的 label 是真正的分類結果，拿來跟我們 predict 的結果比對就可以知道 predict 有沒有猜對了。

也所以，我們可以拿原 training set 當做 test_file再丟給 svmpredict 去 predict (因為格式一樣)，看看正確率有多高，方便後面調參數。

其它參數就很好理解了： model_file 就是 svmtrain 出來的檔案， output_file 是存輸出結果的檔案。

輸出的格式很簡單，每行一個 label，對應到你的 test_file 裡面的各行。

下列程式執行結果會產生 heart_scale.out：

./svm-predict heart_scale heart_scale.model heart_scale.out Accuracy = 86.6667% (234/270) (classification) Mean squared error = 0.533333 (regression) Squared correlation coefficient = 0.532639(regression)

As you can see，我們把原輸入丟回去 predict，第一行的 Accuracy 就是預測的正確率了。如果輸入沒有 label 的話，那就是真的 predict 了。

看到這裡，基本上你應該已經可以利用 svm 來作事了：你只要寫程式輸出正確格式的資料，交給 svm 去 train，後來再 predict 並讀入結果即可。

Advanced Topics

後面可以說是一些稍微進階的部份，我可能不會講的很清楚，因為我的重點是想表達一些觀念和解釋一些你看相關文件時很容易碰到的名詞。

Scaling

svm-scale 目前不太好用，不過它有其必要性。因為適當的scale有助於參數的選擇(後述)還有解svm的速度。
svmscale 會對每個 attribute 做scale。範圍用 -l, -u 指定，通常是[0,1]或是[-1,1]。輸出在 stdout。
另外要注意的(常常會忘記)是 testing data 和 training data要一起scale。
而 svm-scale 最難用的地方就是沒辦法指定 testing data/training data(不同檔案) 然後一起scale。

Arguments

前面提到，在 train 的時候可以下一些參數。(直接執行 svm-train 不指定輸入檔與參數會列出所有參數及語法說明) 這些參數對應到原始 SVM 公式的一些參數，所以會影響 predict 的正確與否。
舉例來說，改個 c=10:
./svm-train -c 10 heart_scale
再來 predict ，正確率馬上變成 92.2% (249/270)。

Cross Validation

一般而言， SVM 使用的方式(在決定參數時)常是這樣：

先有已分好類的一堆資料
亂數拆成好幾組 training set
用某組參數去 train 並 predict 別組看正確率
正確率不夠的話，換參數再重複 train/predict

等找到一組不錯的參數後，就拿這組參數來建 model 並用來做最後對未知資料的 predict。這整個過程叫 cross validation ，也就是交叉比對。

在我們找參數的過程中，可以利用 svmtrain 的內建 cross validation 功能來幫忙：
-v n: n-fold cross validation
n 就是要拆成幾組，像 n=3 就會拆成三組，然後先拿 1跟2來 train model 並 predict 3 以得到正確率；再來拿 2跟 3 train 並 predict 1，最後 1,3 train 並 predict 2。其它以此類推。

如果沒有交叉比對的話，很容易找到只在特定輸入時好的參數。像前面我們 c=10 得到 92.2%，不過拿 -v 5 來看看：./svm-train -v 5 -c 10 heart_scale ... Cross Validation Accuracy = 80.3704%平均之後才只有 80.37%，比一開始的 86 還差。

What arguments rules?

通常而言，比較重要的參數是 gamma (-g) 跟 cost (-c) 。而 cross validation (-v) 的參數常用 5。

cost 預設值是 1, gamma 預設值是 1/k ，k 等於輸入資料筆數。那我們怎麼知道要用多少來當參數呢？

用試的
是的，別懷疑，就是 Try 參數找比較好的值。

Try 參數的過程常用 exponential 指數成長的方式來增加與減少參數的數值，也就是 2^n (2 的 n 次方)。

因為有兩組參數，所以等於要 try n*n=n^2 次。這個過程是不連續的成長，所以可以想成我們在一個 X-Y 平面上指定的範圍內找一群格子點 (grid，如果你不太明白，想成方格紙或我們把平面上所有整數交點都打個點，就是那樣)，每個格子點的 X 跟 Y 經過換算 (如 2^x, 2^y) 就拿去當 cost 跟 gamma 的值來 cross validation。

所以現在你應該懂得 libsvm 的 python 子目錄下面有個 grid.py 是做啥的了：它把上面的過程自動化，在你給定的範圍內呼叫 svm-train 去 try 所有的參數值。 python 是一種語言，在這裡我不做介紹，因為我會了 :P (just a joke，真正原因是 -- 這是 libsvm 的 tutorial)。 grid.py 還會把結果 plot 出來，方便你尋找參數。 libsvm 有很多跟 python 結合的部份，由此可見 python 是強大方便的工具。很多神奇的功能，像自動登入多台機器去平行跑 grid等等都是 python 幫忙的。不過 SVM 本身可以完全不需要 python，只是會比較方便。

跑 grid (基本上用 grid.py 跑當然是最方便，不過如果你不懂 python 而且覺得很難搞，那你要自己產生參數來跑也是可以的) 通常好的範圍是[c,g]=[2^-10,2^10]*[2^-10,2^10]另外其實 grid 用 [-8,8] 也很夠了。

Regression

另一個值得一提的是 regression。

簡單來說，前面都是拿 SVM 來做分類 (classification), 所以 label 的值都是 discrete data、或說已知的固定值。而 regression 則是求 continuous 的值、或說未知的值。你也可以說，一般是 binary classification, 而 regression是可以預測一個實數。

比如說我知道股市指數受到某些因素影響, 然後我想預測股市.. 股市的指數就是我們的 label, 那些因素量化以後變成 attributes。以後蒐集那些 attributes 給 SVM 它就會預測出指數(可能是沒出現過的數字)，這就要用 regression。那樂透開獎的號碼呢？因為都是固定已知的數字，很明顯我們應該用一般 SVM 的 classification 來 predict。 (註：這是真實的例子 -- llwang 就寫過這樣的東西)

所以說 label 也要 scale, 用svm-scale -y lower upper

但是比較糟糕的情況是 grid.py 不支援 regression ，而且cross validation 對 regression 也常常不是很有效。

總而言之，regression 是非常有趣的東西，不過也是比較進階的用法。在這裡我們不細談了，有興趣的人請再參考 SVM 與 libsvm 的其它文件。