10 種機器學習演算法的要點

前言

谷歌董事長施密特曾說過：雖然谷歌的無人駕駛汽車和機器人受到了許多媒體關注，但是這家公司真正的未來在於機器學習，一種讓計算機更聰明、更個性化的技術。

也許我們生活在人類歷史上最關鍵的時期：從使用大型計算機，到個人電腦，再到現在的雲計算。關鍵的不是過去發生了什麼，而是將來會有什麼發生。

工具和技術的民主化，讓像我這樣的人對這個時期興奮不已。計算的蓬勃發展也是一樣。如今，作為一名數據科學家，用複雜的演算法建立數據處理機器一小時能賺到好幾美金。但能做到這個程度可並不簡單！我也曾有過無數黑暗的日日夜夜。

誰能從這篇指南里受益最多？

我今天所給出的，也許是我這輩子寫下的最有價值的指南。

這篇指南的目的，是為那些有追求的數據科學家和機器學習狂熱者們，簡化學習旅途。這篇指南會讓你動手解決機器學習的問題，並從實踐中獲得真知。我提供的是幾個機器學習演算法的高水平理解，以及運行這些演算法的 R 和 Python 代碼。這些應該足以讓你親自試一試了。

我特地跳過了這些技術背後的數據，因為一開始你並不需要理解這些。如果你想從數據層面上理解這些演算法，你應該去別處找找。但如果你想要在開始一個機器學習項目之前做些準備，你會喜歡這篇文章的。

廣義來說，有三種機器學習演算法

1、 監督式學習

工作機制：這個演算法由一個目標變數或結果變數（或因變數）組成。這些變數由已知的一系列預示變數（自變數）預測而來。利用這一系列變數，我們生成一個將輸入值映射到期望輸出值的函數。這個訓練過程會一直持續，直到模型在訓練數據上獲得期望的精確度。監督式學習的例子有：回歸、決策樹、隨機森林、K – 近鄰演算法、邏輯回歸等。

2、非監督式學習

工作機制：在這個演算法中，沒有任何目標變數或結果變數要預測或估計。這個演算法用在不同的組內聚類分析。這種分析方式被廣泛地用來細分客戶，根據干預的方式分為不同的用戶組。非監督式學習的例子有：關聯演算法和 K – 均值演算法。

3、強化學習

工作機制：這個演算法訓練機器進行決策。它是這樣工作的：機器被放在一個能讓它通過反覆試錯來訓練自己的環境中。機器從過去的經驗中進行學習，並且嘗試利用了解最透徹的知識作出精確的商業判斷。 強化學習的例子有馬爾可夫決策過程。

常見機器學習演算法名單

這裡是一個常用的機器學習演算法名單。這些演算法幾乎可以用在所有的數據問題上：

線性回歸

邏輯回歸

決策樹

SVM

樸素貝葉斯

K最近鄰演算法

K均值演算法

隨機森林演算法

降維演算法

Gradient Boost 和 Adaboost 演算法

1、線性回歸

線性回歸通常用於根據連續變數估計實際數值（房價、呼叫次數、總銷售額等）。我們通過擬合最佳直線來建立自變數和因變數的關係。這條最佳直線叫做回歸線，並且用 Y= a *X + b 這條線性等式來表示。

理解線性回歸的最好辦法是回顧一下童年。假設在不問對方體重的情況下，讓一個五年級的孩子按體重從輕到重的順序對班上的同學排序，你覺得這個孩子會怎麼做？他（她）很可能會目測人們的身高和體型，綜合這些可見的參數來排列他們。這是現實生活中使用線性回歸的例子。實際上，這個孩子發現了身高和體型與體重有一定的關係，這個關係看起來很像上面的等式。

在這個等式中：

Y：因變數

a：斜率

x：自變數

b ：截距

係數 a 和 b 可以通過最小二乘法獲得。

參見下例。我們找出最佳擬合直線y=0.2811x+13.9。已知人的身高，我們可以通過這條等式求出體重。

線性回歸的兩種主要類型是一元線性回歸和多元線性回歸。一元線性回歸的特點是只有一個自變數。多元線性回歸的特點正如其名，存在多個自變數。找最佳擬合直線的時候，你可以擬合到多項或者曲線回歸。這些就被叫做多項或曲線回歸。

Python 代碼

#Import Library

#Import other necessary libraries like pandas, numpy...

fromsklearnimportlinear_model

#Load Train and Test datasets

#Identify feature and response variable(s) and values must be numeric and numpy arrays

x_train=input_variables_values_training_datasets

y_train=target_variables_values_training_datasets

x_test=input_variables_values_test_datasets

# Create linear regression object

linear=linear_model.LinearRegression()

# Train the model using the training sets and check score

linear.fit(x_train,y_train)

linear.score(x_train,y_train)

#Equation coefficient and Intercept

print("Coefficient: n",linear.coef_)

print("Intercept: n",linear.intercept_)

#Predict Output

predicted=linear.predict(x_test)

R代碼

#Load Train and Test datasets

#Identify feature and response variable(s) and values must be numeric and numpy arrays

x_train

y_train

x_test

x

# Train the model using the training sets and check score

linear

summary(linear)

#Predict Output

predicted=predict(linear,x_test)

、

2、邏輯回歸

別被它的名字迷惑了！這是一個分類演算法而不是一個回歸演算法。該演算法可根據已知的一系列因變數估計離散數值（比方說二進位數值 0 或 1 ，是或否，真或假）。簡單來說，它通過將數據擬合進一個邏輯函數來預估一個事件出現的概率。因此，它也被叫做邏輯回歸。因為它預估的是概率，所以它的輸出值大小在 0 和 1 之間（正如所預計的一樣）。

讓我們再次通過一個簡單的例子來理解這個演算法。

假設你的朋友讓你解開一個謎題。這隻會有兩個結果：你解開了或是你沒有解開。想像你要解答很多道題來找出你所擅長的主題。這個研究的結果就會像是這樣：假設題目是一道十年級的三角函數題，你有 70%的可能會解開這道題。然而，若題目是個五年級的歷史題，你只有30%的可能性回答正確。這就是邏輯回歸能提供給你的信息。

從數學上看，在結果中，幾率的對數使用的是預測變數的線性組合模型。

odds=p/(1-p)=probability of eventoccurrence/probability ofnotevent occurrence

ln(odds)=ln(p/(1-p))

logit(p)=ln(p/(1-p))=b0+b1X1+b2X2+b3X3....+bkXk

在上面的式子里，p 是我們感興趣的特徵出現的概率。它選用使觀察樣本值的可能性最大化的值作為參數，而不是通過計算誤差平方和的最小值（就如一般的回歸分析用到的一樣）。

現在你也許要問了，為什麼我們要求出對數呢？簡而言之，這種方法是複製一個階梯函數的最佳方法之一。我本可以更詳細地講述，但那就違背本篇指南的主旨了。

Python代碼

#Import Library

fromsklearn.linear_modelimportLogisticRegression

#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset

# Create logistic regression object

model=LogisticRegression()

# Train the model using the training sets and check score

model.fit(X,y)

model.score(X,y)

#Equation coefficient and Intercept

print("Coefficient: n",model.coef_)

print("Intercept: n",model.intercept_)

#Predict Output

predicted=model.predict(x_test)

R代碼

x

# Train the model using the training sets and check score

logistic

summary(logistic)

#Predict Output

predicted=predict(logistic,x_test)

更進一步：

你可以嘗試更多的方法來改進這個模型：

加入交互項

精簡模型特性

使用正則化方法

使用非線性模型

3、決策樹

這是我最喜愛也是最頻繁使用的演算法之一。這個監督式學習演算法通常被用於分類問題。令人驚奇的是，它同時適用於分類變數和連續因變數。在這個演算法中，我們將總體分成兩個或更多的同類群。這是根據最重要的屬性或者自變數來分成儘可能不同的組別。想要知道更多，可以閱讀：簡化決策樹。

（來源： statsexchange）

在上圖中你可以看到，根據多種屬性，人群被分成了不同的四個小組，來判斷 「他們會不會去玩」。為了把總體分成不同組別，需要用到許多技術，比如說 Gini、Information Gain、Chi-square、entropy。

理解決策樹工作機制的最好方式是玩Jezzball，一個微軟的經典遊戲（見下圖）。這個遊戲的最終目的，是在一個可以移動牆壁的房間里，通過造牆來分割出沒有小球的、盡量大的空間。

因此，每一次你用牆壁來分隔房間時，都是在嘗試著在同一間房裡創建兩個不同的總體。相似地，決策樹也在把總體盡量分割到不同的組裡去。

更多信息請見：決策樹演算法的簡化(http://www.analyticsvidhya.com/blog/2015/01/decision-tree-simplified/)

Python代碼

#Import Library

#Import other necessary libraries like pandas, numpy...

fromsklearnimporttree

#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset

# Create tree object

model=tree.DecisionTreeClassifier(criterion="gini")# for classification, here you can change the algorithm as gini or entropy (information gain) by default it is gini

# model = tree.DecisionTreeRegressor() for regression

# Train the model using the training sets and check score

model.fit(X,y)

model.score(X,y)

#Predict Output

predicted=model.predict(x_test)

R代碼

library(rpart)

x

# grow tree

fit

summary(fit)

#Predict Output

predicted=predict(fit,x_test)

4、支持向量機

這是一種分類方法。在這個演算法中，我們將每個數據在N維空間中用點標出（N是你所有的特徵總數），每個特徵的值是一個坐標的值。

舉個例子，如果我們只有身高和頭髮長度兩個特徵，我們會在二維空間中標出這兩個變數，每個點有兩個坐標（這些坐標叫做支持向量）。

現在，我們會找到將兩組不同數據分開的一條直線。兩個分組中距離最近的兩個點到這條線的距離同時最優化。

上面示例中的黑線將數據分類優化成兩個小組，兩組中距離最近的點（圖中A、B點）到達黑線的距離滿足最優條件。這條直線就是我們的分割線。接下來，測試數據落到直線的哪一邊，我們就將它分到哪一類去。

更多請見：支持向量機的簡化(http://www.analyticsvidhya.com/blog/2014/10/support-vector-machine-simplified/)

將這個演算法想作是在一個 N 維空間玩 JezzBall。需要對遊戲做一些小變動：

比起之前只能在水平方向或者豎直方向畫直線，現在你可以在任意角度畫線或平面。

遊戲的目的變成把不同顏色的球分割在不同的空間里。

球的位置不會改變。

Python代碼

#Import Library

fromsklearnimportsvm

#Assumed you have, X (predic

tor)andY(target)fortraining data setandx_test(predictor)oftest_dataset

# Create SVM classification object

model=svm.svc()# there is various option associated with it, this is simple for classification. You can refer link, for mo# re detail.

# Train the model using the training sets and check score

model.fit(X,y)

model.score(X,y)

#Predict Output

predicted=model.predict(x_test)

R代碼

library(e1071)

x

# Fitting model

fit

summary(fit)

#Predict Output

predicted=predict(fit,x_test)

5、樸素貝葉斯

在預示變數間相互獨立的前提下，根據貝葉斯定理可以得到樸素貝葉斯這個分類方法。用更簡單的話來說，一個樸素貝葉斯分類器假設一個分類的特性與該分類的其它特性不相關。舉個例子，如果一個水果又圓又紅，並且直徑大約是 3 英寸，那麼這個水果可能會是蘋果。即便這些特性互相依賴，或者依賴於別的特性的存在，樸素貝葉斯分類器還是會假設這些特性分別獨立地暗示這個水果是個蘋果。

樸素貝葉斯模型易於建造，且對於大型數據集非常有用。雖然簡單，但是樸素貝葉斯的表現卻超越了非常複雜的分類方法。

貝葉斯定理提供了一種從P(c)、P(x)和P(x|c) 計算後驗概率 P(c|x) 的方法。請看以下等式：

在這裡，

P(c|x) 是已知預示變數（屬性）的前提下，類（目標）的後驗概率

P(c) 是類的先驗概率

P(x|c) 是可能性，即已知類的前提下，預示變數的概率

P(x) 是預示變數的先驗概率

例子：讓我們用一個例子來理解這個概念。在下面，我有一個天氣的訓練集和對應的目標變數「Play」。現在，我們需要根據天氣情況，將會「玩」和「不玩」的參與者進行分類。讓我們執行以下步驟。

步驟1：把數據集轉換成頻率表。

步驟2：利用類似「當Overcast可能性為0.29時，玩耍的可能性為0.64」這樣的概率，創造 Likelihood 表格。

步驟3：現在，使用樸素貝葉斯等式來計算每一類的後驗概率。後驗概率最大的類就是預測的結果。

問題：如果天氣晴朗，參與者就能玩耍。這個陳述正確嗎？

我們可以使用討論過的方法解決這個問題。於是 P（會玩 | 晴朗）= P（晴朗 | 會玩）* P（會玩）/ P （晴朗）

我們有 P （晴朗 |會玩）= 3/9 = 0.33，P（晴朗） = 5/14 = 0.36, P（會玩）= 9/14 = 0.64

現在，P(會玩 | 晴朗）= 0.33 * 0.64 / 0.36 = 0.60，有更大的概率。

樸素貝葉斯使用了一個相似的方法，通過不同屬性來預測不同類別的概率。這個演算法通常被用於文本分類，以及涉及到多個類的問題。

Python代碼

#Import Library

fromsklearn.naive_bayesimportGaussianNB

#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset

# Create SVM classification object model = GaussianNB() # there is other distribution for multinomial classes like Bernoulli Naive Bayes, Refer link

# Train the model using the training sets and check score

model.fit(X,y)

#Predict Output

predicted=model.predict(x_test)

R代碼

library(e1071)

x

# Fitting model

fit

summary(fit)

#Predict Output

predicted=predict(fit,x_test)

6、KNN（K – 最近鄰演算法）

該演算法可用於分類問題和回歸問題。然而，在業界內，K – 最近鄰演算法更常用於分類問題。K – 最近鄰演算法是一個簡單的演算法。它儲存所有的案例，通過周圍k個案例中的大多數情況劃分新的案例。根據一個距離函數，新案例會被分配到它的 K 個近鄰中最普遍的類別中去。

這些距離函數可以是歐式距離、曼哈頓距離、明式距離或者是漢明距離。前三個距離函數用於連續函數，第四個函數（漢明函數）則被用於分類變數。如果 K=1，新案例就直接被分到離其最近的案例所屬的類別中。有時候，使用 KNN 建模時，選擇 K 的取值是一個挑戰。

更多信息：K – 最近鄰演算法入門（簡化版）

我們可以很容易地在現實生活中應用到 KNN。如果想要了解一個完全陌生的人，你也許想要去找他的好朋友們或者他的圈子來獲得他的信息。

在選擇使用 KNN 之前，你需要考慮的事情：

KNN 的計算成本很高。

變數應該先標準化（normalized），不然會被更高範圍的變數偏倚。

在使用KNN之前，要在野值去除和噪音去除等前期處理多花功夫。

Python代碼

#Import Library

fromsklearn.neighborsimportKNeighborsClassifier

#Assumed you have, X (predictor) and Y (target) for training data set and x_test(predictor) of test_dataset

# Create KNeighbors classifier object model

KNeighborsClassifier(n_neighbors=6)# default value for n_neighbors is 5

# Train the model using the training sets and check score

model.fit(X,y)

#Predict Output

predicted=model.predict(x_test)

R代碼

library(knn)

x

# Fitting model

fit

summary(fit)

#Predict Output

predicted=predict(fit,x_test)

7、K 均值演算法

K – 均值演算法是一種非監督式學習演算法，它能解決聚類問題。使用 K – 均值演算法來將一個數據歸入一定數量的集群（假設有 k 個集群）的過程是簡單的。一個集群內的數據點是均勻齊次的，並且異於別的集群。

還記得從墨水漬里找出形狀的活動嗎？K – 均值演算法在某方面類似於這個活動。觀察形狀，並延伸想像來找出到底有多少種集群或者總體。

K – 均值演算法怎樣形成集群：

K – 均值演算法給每個集群選擇k個點。這些點稱作為質心。

每一個數據點與距離最近的質心形成一個集群，也就是 k 個集群。

根據現有的類別成員，找出每個類別的質心。現在我們有了新質心。

當我們有新質心後，重複步驟 2 和步驟 3。找到距離每個數據點最近的質心，並與新的k集群聯繫起來。重複這個過程，直到數據都收斂了，也就是當質心不再改變。

如何決定 K 值：

K – 均值演算法涉及到集群，每個集群有自己的質心。一個集群內的質心和各數據點之間距離的平方和形成了這個集群的平方值之和。同時，當所有集群的平方值之和加起來的時候，就組成了集群方案的平方值之和。

我們知道，當集群的數量增加時，K值會持續下降。但是，如果你將結果用圖表來表示，你會看到距離的平方總和快速減少。到某個值 k 之後，減少的速度就大大下降了。在此，我們可以找到集群數量的最優值。

Python代碼

#Import Library

fromsklearn.clusterimportKMeans

#Assumed you have, X (attributes) for training data set and x_test(attributes) of test_dataset

# Create KNeighbors classifier object model

k_means=KMeans(n_clusters=3,random_state=)

# Train the model using the training sets and check score

model.fit(X)

#Predict Output

predicted=model.predict(x_test)

R代碼

library(cluster)

fit

8、隨機森林

隨機森林是表示決策樹總體的一個專有名詞。在隨機森林演算法中，我們有一系列的決策樹（因此又名「森林」）。為了根據一個新對象的屬性將其分類，每一個決策樹有一個分類，稱之為這個決策樹「投票」給該分類。這個森林選擇獲得森林裡（在所有樹中）獲得票數最多的分類。

每棵樹是像這樣種植養成的：

如果訓練集的案例數是 N，則從 N 個案例中用重置抽樣法隨機抽取樣本。這個樣本將作為「養育」樹的訓練集。

假如有 M 個輸入變數，則定義一個數字 m

儘可能大地種植每一棵樹，全程不剪枝。

喜歡這篇文章嗎？立刻分享出去讓更多人知道吧！