煩惱

作為一名資料分析師，你來到這家跨國銀行工作已經半年了。

今天上午，老闆把你叫到辦公室，面色凝重。

你心裡直打鼓，以為自己捅了什麼簍子。幸好老闆的話讓你很快打消了顧慮。

客戶主要分佈在法國、德國和西班牙。

你手裡掌握的資訊，包括他們的年齡、性別、信用、辦卡資訊等。客戶是否已流失的資訊在最後一列（Exited）。

請選擇左側的Python 3.6版本下載安裝。

其次是新建資料夾，起名為demo-customer-churn-ann，並且從這個連結下載資料，放到該資料夾下。

點選介面右上方的New按鈕，新建一個Python 3 Notebook，起名為customer-churn-ann。

準備工作結束，下面我們開始清理資料。

清理

首先，讀入資料清理最常用的pandas和numpy包。

可以看到，資料完整無誤讀入。但是並非所有的列都對我們預測使用者流失有作用。我們一一甄別一下：

• RowNumber：行號，這個肯定沒用，刪除
• CustomerID：使用者編號，這個是順序發放的，刪除
• Surname：使用者姓名，對流失沒有影響，刪除
• CreditScore：信用分數，這個很重要，保留
• Geography：使用者所在國家/地區，這個有影響，保留
• Gender：使用者性別，可能有影響，保留
• Age：年齡，影響很大，年輕人更容易切換銀行，保留
• Tenure：當了本銀行多少年使用者，很重要，保留
• Balance：存貸款情況，很重要，保留
• NumOfProducts：使用產品數量，很重要，保留
• HasCrCard：是否有本行信用卡，很重要，保留
• IsActiveMember：是否活躍使用者，很重要，保留
• EstimatedSalary：估計收入，很重要，保留
• Exited：是否已流失，這將作為我們的標籤資料

在Scikit-learn工具包裡面，專門提供了方便的工具 LabelEncoder ，讓我們可以方便地將類別資訊變成數值。

這樣是不是就完事大吉了呢？

不對，Gender還好說，只有兩種取值方式，要麼是男，要麼是女。我們可以把“是男性”定義為1，那麼女性就取值為0。兩種取值只是描述類別不同，沒有歧義。

而Geography就不同了。因為資料集裡面可能的國家地區取值有3種，所以就轉換成了0（法國）、1（德國）、2（西班牙）。問題是，這三者之間真的有序列（大小）關係嗎？

答案自然是否定的。我們其實還是打算用數值描述分類而已。但是取值有數量的序列差異，就會給機器帶來歧義。它並不清楚不同的取值只是某個國家的程式碼，可能會把這種大小關係帶入模型計算，從而產生錯誤的結果。

沒有。

因為本例中，OneHotEncoder轉換出來的3列數字，實際上是不獨立的。給定其中兩列的資訊，你自己都可以計算出其中的第3列取值。

好比說，某一行的前兩列數字是 (0, 0) ，那麼第三列肯定是1。因為這是轉換規則決定的。3列裡只能有1個是1，其餘都是0。

如果你做過多元線性迴歸，應該知道這種情況下，我們是需要去掉其中一列，才能繼續分析的。不然會落入“虛擬變數陷阱”（dummy variable trap）。

我們刪掉第0列，避免掉進坑裡。

X = np.delete(X, [0], 1)

再次列印第一行：

這樣在後面訓練的時候，他就可以和前面的特徵矩陣一一對應來操作計算了。

既然標籤代表了類別，我們也把它用OneHotEncoder轉換，這樣方便我們後面做分類學習。

onehotencoder = OneHotEncoder()

y = onehotencoder.fit_transform(y).toarray()

此時的標籤變成兩列資料，一列代表顧客存留，一列代表顧客流失。

y

array([[ 0., 1.],

[ 1., 0.],

[ 0., 1.],

...,

[ 0., 1.],

[ 0., 1.],

[ 1., 0.]])

總體的資料已經齊全了。但是我們 不能 把它們 都用來 訓練。

你會發現，許多列的方差比原先小得多。機器學習起來，會更加方便。

資料清理和轉換工作至此完成。

決策樹

如果讀過我的《 貸還是不貸：如何用Python和機器學習幫你決策？ 》一文，你應該有一種感覺——這個問題和貸款審批決策很像啊！既然在該文中，決策樹很好使，我們繼續用決策樹不就好了？

經檢測，決策樹在咱們的資料集上，表現得還是不錯的。總體的準確率為0.81，召回率為0.80，f1分數為0.81，已經很高了。對10個客戶做流失可能性判斷，它有8次都能判斷正確。

但是，這樣是否足夠？

我們或許可以調整決策樹的引數做優化，嘗試改進預測結果。

或者我們可以採用 深度學習 。

深度

深度學習的使用場景，往往是因為原有的模型經典機器學習模型過於簡單，無法把握複雜資料特性。

我不準備給你講一堆數學公式，咱們動手做個實驗。

請你開啟這個網址。

你會看到如下圖所示的深度學習遊樂場：

右側的圖形，裡面是藍色資料，外圈是黃色資料。你的任務就是要用模型分類兩種不同資料。

你說那還不容易？我一眼就看出來了。

你看出來沒有用。通過你的設定，讓機器也能正確區分，才算數。

圖中你看到許多加減號。咱們就通過操縱它們來玩兒一玩兒模型。

首先，點圖中部上方的"2 HIDDEN LAYERS"左側減號，把中間隱藏層數降低為1。

然後，點選"2 neurons"上面的減號，把神經元數量減少為1。

把頁面上方的Activation函式下拉框開啟，選擇“Sigmoid”。

現在的模型，其實就是經典的邏輯迴歸（Logistic Regression）。

點選左上方的執行按鈕，我們看看執行效果。

由於模型過於簡單，所以機器絞盡腦汁，試圖用一條直線切分二維平面上的兩類節點。

損失(loss)居高不下。訓練集和測試集損失都在0.4左右，顯然不符合我們的分類需求。

下面我們試試增加層數和神經元數量。這次點選加號，把隱藏層數加回到2，兩層神經元數量都取2。

再次點選執行。

經過一段時間，結果穩定了下來，你發現這次電腦用了兩條線，把平面切分成了3部分。

測試集損失下降到了0.25左右，而訓練集損失更是降低到了0.2以下。

模型複雜了，效果似乎更好一些。

再接再厲，我們把第一個隱藏層的神經元數量增加為4看看。

點選執行，不一會兒有趣的事情就發生了。

機器用一條近乎完美的曲線把平面分成了內外兩個部分。測試集和訓練集損失都極速下降，訓練集損失甚至接近於0。

這告訴我們，許多時候模型過於簡單帶來的問題，可以通過加深隱藏層次、增加神經元的方法提升模型複雜度，加以改進。

目前流行的劃分方法，是用隱藏層的數量多少來區分是否“深度”。當神經網路中隱藏層數量達到3層以上時，就被稱為“深度神經網路”，或者“深度學習”。

久聞大名的深度學習，原來就是這麼簡單。

如果有時間的話，建議你自己在這個遊樂場裡多動手玩兒一玩兒。你會很快對神經網路和深度學習有個感性認識。

框架

遊樂場背後使用的引擎，就是Google的深度學習框架Tensorflow。

所謂框架，就是別人幫你構造好的基礎軟體應用。你可以通過呼叫它們，避免自己重複發明輪子，大幅度節省時間，提升效率。

支援Python語言的深度學習的框架有很多，除了Tensorflow外，還有PyTorch, Theano和MXNet等。

我給你的建議是，找到一個你喜歡的軟體包，深入學習使用，不斷實踐來提升自己的技能。 千萬不要 跟別人爭論哪個深度學習框架更好。一來蘿蔔白菜各有所愛，每個人都有自己的偏好；二來深度學習的江湖水很深，言多有失。說錯了話，別的門派可能會不高興喲。

我比較喜歡Tensorflow。但是Tensorflow本身是個底層庫。雖然隨著版本的更迭，介面越來越易用。但是對初學者來說，許多細節依然有些過於瑣碎，不容易掌握。

初學者的耐心有限，挫折過多容易放棄。

幸好，還有幾個高度抽象框架，是建立在Tensorflow之上的。如果你的任務是 應用 現成的深度學習模型，那麼這些框架會給你帶來非常大的便利。

這些框架包括Keras, TensorLayer等。咱們今天將要使用的，叫做TFlearn。

它的特點，就是長得很像Scikit-learn。這樣如果你熟悉經典機器學習模型，學起來會特別輕鬆省力。

實戰

閒話就說這麼多，下面咱們繼續寫程式碼吧。

寫程式碼之前，請回到終端下，執行以下命令，安裝幾個軟體包：

pip install tensorflow

pip install tflearn

執行完畢後，回到Notebook裡。

我們呼叫tflearn框架。

import tflearn

然後，我們開始搭積木一樣，搭神經網路層。

首先是輸入層。

net = tflearn.input_data(shape=[None, 11])

注意這裡的寫法，因為我們輸入的資料，是特徵矩陣。而經過我們處理後，特徵矩陣現在有11列，因此shape的第二項寫11。

shape的第一項，None，指的是我們要輸入的特徵矩陣行數。因為我們現在是搭建模型，後面特徵矩陣有可能一次輸入，有可能分成組塊輸入，長度可大可小，無法事先確定。所以這裡填None。tflearn會在我們實際執行訓練的時候，自己讀入特徵矩陣的尺寸，來處理這個數值。

下面我們搭建隱藏層。這裡我們要使用深度學習，搭建3層。

net = tflearn.fully_connected(net, 6, activation='relu')

net = tflearn.fully_connected(net, 6, activation='relu')

net = tflearn.fully_connected(net, 6, activation='relu')

activation剛才在深度學習遊樂場裡面我們遇到過，代表啟用函式。如果沒有它，所有的輸入輸出都是線性關係。

Relu函式是啟用函式的一種。它大概長這個樣子。

如果你想了解啟用函式的更多知識，請參考後文的學習資源部分。

隱藏層裡，每一層我們都設定了6個神經元。其實至今為之，也不存在最優神經元數量的計算公式。工程界的一種做法，是把輸入層的神經元數量，加上輸出層神經元數量，除以2取整。咱們這裡就是用的這種方法，得出6個。

搭好了3箇中間隱藏層，下面我們來搭建輸出層。

net = tflearn.fully_connected(net, 2, activation='softmax')

net = tflearn.regression(net)

這裡我們用兩個神經元做輸出，並且說明使用迴歸方法。輸出層選用的啟用函式為softmax。處理分類任務的時候，softmax比較合適。它會告訴我們每一類的可能性，其中數值最高的，可以作為我們的分類結果。

積木搭完了，下面我們告訴TFlearn，以剛剛搭建的結構，生成模型。

model = tflearn.DNN(net)

有了模型，我們就可以使用擬合功能了。你看是不是跟Scikit-learn的使用方法很相似呢？

model.fit(X_train, y_train, n_epoch=30, batch_size=32, show_metric=True)

注意這裡多了幾個引數，我們來解釋一下。

n_epoch

batch_size

show_metric

以下就是電腦輸出的最終訓練結果。其實中間執行過程看著更激動人心，你自己試一下就知道了。

Training Step: 7499 | total loss: [1m[32m0.39757[0m[0m | time: 0.656s

| Adam | epoch: 030 | loss: 0.39757 - acc: 0.8493 -- iter: 7968/8000

Training Step: 7500 | total loss: [1m[32m0.40385[0m[0m | time: 0.659s

| Adam | epoch: 030 | loss: 0.40385 - acc: 0.8487 -- iter: 8000/8000

--

我們看到訓練集的損失(loss)大概為0.4左右。

開啟終端，我們輸入

tensorboard --logdir=/tmp/tflearn_logs/

然後在瀏覽器裡輸入 http://localhost:6006/

可以看到如下介面：

這是模型訓練過程的視覺化圖形，可以看到準確度的攀升和損失降低的曲線。

開啟GRAPHS標籤頁，我們可以檢視神經網路的結構圖形。

我們搭積木的過程，在此處一目瞭然。

評估

訓練好了模型，我們來嘗試做個預測吧。

看看測試集的特徵矩陣第一行。

X_test[0]

array([ 1.75486502, -0.57369368, -0.55204276, -1.09168714, -0.36890377,

1.04473698, 0.8793029 , -0.92159124, 0.64259497, 0.9687384 ,

1.61085707])

我們就用它來預測一下分類結果。

y_pred = model.predict(X_test)

打印出來看看：

y_pred[0]

array([ 0.70956731, 0.29043278], dtype=float32)

模型判斷該客戶不流失的可能性為0.70956731。

我們看看實際標籤資料：

y_test[0]

array([ 1., 0.])

客戶果然沒有流失。這個預測是對的。

相關推薦