利用 Scikit Learn的Python数据预处理实战指南

简介:



简而言之,预处理是指在你将数据“喂给”算法之前进行的一系列转换操作。在Python中,scikit-learn库在sklearn.preprocessing下有预装的功能。有更多的选择来进行预处理,这将是我们要探索的。


读完本文,你将具备数据预处理的基本技能并对其有更深入的理解。为了方便起见,我附上了进一步学习机器学习算法的一些资源,并且为更好地掌握这些概念,设计了几个小练习。


可用数据集

本文中,我使用了部分的贷款预测数据,缺失观测值的数据已被移除(需要数据的读者朋友,请在评论区留下电邮地址,我们会把数据发给你——译者注)。


备注:贷款预测问题中,测试集数据是训练集的子集。


现在,让我们从导入重要的包和数据集开始。


导入pandas

>> import pandas as pd

导入训练用的数据集

>> X_train=pd.read_csv('X_train.csv')

>> Y_train=pd.read_csv('Y_train.csv')

导入测试数据集

>> X_test=pd.read_csv('X_test.csv')

>> Y_test=pd.read_csv('Y_test.csv')

对我们的数据集进行仔细观察。

>> print (X_train.head())

     Loan_ID Gender Married Dependents Education Self_Employed

15   LP001032   Male      No          0  Graduate            No   

248  LP001824   Male     Yes          1  Graduate            No   

590  LP002928   Male     Yes          0  Graduate            No   

246  LP001814   Male     Yes          2  Graduate            No   

388  LP002244   Male     Yes          0  Graduate            No   

 
 

     ApplicantIncome  CoapplicantIncome  LoanAmount  Loan_Amount_Term

15              4950                0.0       125.0             360.0   

248             2882             1843.0       123.0             480.0   

590             3000             3416.0        56.0             180.0   

246             9703                0.0       112.0             360.0   

388             2333             2417.0       136.0             360.0   

 
 

     Credit_History Property_Area  

15              1.0         Urban  

248             1.0     Semiurban  

590             1.0     Semiurban  

246             1.0         Urban  

388             1.0         Urban


特征缩放

特征缩放是用来限制变量范围的方法,以让它们能在相同的尺度上进行比较。这是在连续变量上操作的。让我们输出数据集中所有连续变量的分布。


>> import matplotlib.pyplot as plt

>> X_train[X_train.dtypes[(X_train.dtypes=="float64")|(X_train.dtypes=="int64")]

                        .index.values].hist(figsize=[11,11])




理解以上图示后,我们推断ApplicantIncome(申请人收入) 和CoapplicantIncome(共同申请人收入) 有相似的尺度范围(0-50000$),LoanAmount(贷款额度) 以千为单位,范围在0 到 600$之间,而Loan_Amount_Term(贷款周期)与其它变量完全不同,因为它的单位是月份,而其它变量单位是美元。


如果我们尝试应用基于距离的算法,如KNN,在这些特征上,范围最大的特征会决定最终的输出结果,那么我们将得到较低的预测精度。我们可通过特征缩放解决这个问题。让我们实践一下。


资料:阅读这篇关于KNN的文章获得更好的理解。(https://www.analyticsvidhya.com/blog/2014/10/introduction-k-neighbours-algorithm-clustering/)

让我们在我们的数据集中试试KNN,看看它表现如何。

# 初始化和拟合一个KNN模型

>> from sklearn.neighbors import KNeighborsClassifier

>> knn=KNeighborsClassifier(n_neighbors=5)

>> knn.fit(X_train[['ApplicantIncome', 'CoapplicantIncome','LoanAmount', 

                   'Loan_Amount_Term', 'Credit_History']],Y_train)

# Checking the performance of our model on the testing data set

# 检查我们的模型在测试数据集上的性能

>> from sklearn.metrics import accuracy_score

>> accuracy_score(Y_test,knn.predict(X_test[['ApplicantIncome', 'CoapplicantIncome',

                             'LoanAmount', 'Loan_Amount_Term', 'Credit_History']]))


Out : 0.61458333333333337


我们得到了大约61%的正确预测,这不算糟糕,但在真正实践中,这是否足够?我们能否将该模型部署于实际问题中?为回答该问题,让我们看看在训练集中关于Loan_Status(贷款状态) 的分布。

>> Y_train.Target.value_counts()/Y_train.Target.count()


Out : Y 0.705729

     N 0.294271

Name: Loan_Status, dtype: float64


大约有70%贷款会被批准,因为有较高的贷款批准率,我们就建立一个所有贷款都通过的预测模型,继续操作并检测我们的预测精度。


 >> Y_test.Target.value_counts()/Y_test.Target.count()


Out :  Y 0.635417

      N 0.364583

Name: Loan_Status, dtype: float64


哇!通过猜测,我们获得63%的精度。这意味着,该模型比我们的预测模型得到更高的精度?


这可能是因为某些具有较大范围的无关紧要的变量主导了目标函数。我们可以通过缩小所有特征到同样的范围来消除该问题。Sklearn提供了MinMaxScaler 工具将所有特征的范围缩小到0-1之间,MinMaxScaler 的数学表达式如下所示:


让我们在我们的问题中试试该工具。

# 导入MinMaxScaler并初始化

>> from sklearn.preprocessing import MinMaxScaler

>> min_max=MinMaxScaler()# Scaling down both train and test data set

>> X_train_minmax=min_max.fit_transform(X_train[['ApplicantIncome', 'CoapplicantIncome',

                'LoanAmount', 'Loan_Amount_Term', 'Credit_History']])

>> X_test_minmax=min_max.fit_transform(X_test[['ApplicantIncome', 'CoapplicantIncome',

                'LoanAmount', 'Loan_Amount_Term', 'Credit_History']])


现在,我们已经完成缩放操作,让我们在缩放后的数据上应用KNN并检测其精度。


在我们缩小后的数据集上拟合KNN

>> knn=KNeighborsClassifier(n_neighbors=5)

>> knn.fit(X_train_minmax,Y_train)

# 检查该模型的精度

>> accuracy_score(Y_test,knn.predict(X_test_minmax))


Out : 0.75


太好了!我们的精度从61%提升到了75%。这意味在基于距离的方法中(如:KNN),一些大范围的特征对预测结果有决定性作用。


应当牢记,当使用基于距离的算法时,我们必须尝试将数据缩放,这样较不重要的特征不会因为自身较大的范围而主导目标函数。此外,具有不同度量单位的特征也应该进行缩放,这样给每个特征具有相同的初始权重,最终我们会得到更好的预测模型。


练习1:


尝试利用逻辑回归模型做相同的练习(参数: penalty=’l2′,C=0.01), 并请在评论区留下缩放前后的精度。


特征标准化

在进入这部分内容前,我建议你先完成练习1。


在之前的章节,我们在贷款预测数据集之上操作,并在其上拟合出一个KNN学习模型。通过缩小数据,我们得到了75%的精度,这看起来十分不错。我在逻辑回归模型上尝试了同样的练习, 并得到如下结果:

Before Scaling : 61%

After Scaling : 63%

缩放前:61%

缩放后:63%


缩放后的精度与我们凭猜测得到的预测精度相近,这并不是很了不起的成就。那么,这是怎么回事呢?在精度上,为什么不像用KNN一样有令人满意的提升?


资料:阅读本文(https://www.analyticsvidhya.com/blog/2015/08/comprehensive-guide-regression/)获得对逻辑回归更好的理解。


答案在此:


在逻辑回归中,每个特征都被分配了权重或系数(Wi)。如果某个特征有相对来说比较大的范围,而且其在目标函数中无关紧要,那么逻辑回归模型自己就会分配一个非常小的值给它的系数,从而中和该特定特征的影响优势,而基于距离的方法,如KNN,没有这样的内置策略,因此需要缩放。


我们是否忘了什么?我们的逻辑模型的预测精度和猜测的几乎接近。


现在,我将在此介绍一个新概念,叫作标准化。很多Sklearn中的机器学习算法都需要标准化后的数据,这意味数据应具有零均值和单位方差。


标准化(或Z-score正则化)是对特征进行重新调整,让数据服从基于 μ=0 和 σ=1的标准正态分布,其中μ是均值(平均值)而σ是关于均值的标准偏差。样本的标准分数(也称为z-scores)按如下所示的方法计算:



线性模型中因子如l1,l2正则化和学习器的目标函数中的SVM中的RBF核心假设所有的特征都集中在0周围并且有着相同顺序的偏差。


有更大顺序的方差的特征将在目标函数上起决定作用,因为前面的章节中,有着更大范围的特征产生过此情形。 正如我们在练习1中看到的,没进行任何预处理的数据之上的精度是61%,让我们标准化我们的数据,在其上应用逻辑回归。Sklearn提供了尺度范围用于标准化数据。


# 标准化训练和测试数据

>> from sklearn.preprocessing import scale

>> X_train_scale=scale(X_train[['ApplicantIncome', 'CoapplicantIncome',

                'LoanAmount', 'Loan_Amount_Term', 'Credit_History']])

>> X_test_scale=scale(X_test[['ApplicantIncome', 'CoapplicantIncome',

               'LoanAmount', 'Loan_Amount_Term', 'Credit_History']])

# 在我们的标准化了的数据集上拟合逻辑回归

>> from sklearn.linear_model import LogisticRegression

>> log=LogisticRegression(penalty='l2',C=.01)

>> log.fit(X_train_scale,Y_train)

# 检查该模型的精度

>> accuracy_score(Y_test,log.predict(X_test_scale))


Out : 0.75


我们再次达到缩放后利用KNN所能达到的我们最大的精度。这意味着,当使用l1或l2正则化估计时,标准化数据帮助我们提高预测模型的精度。其它学习模型,如有欧几里得距离测量的KNN、k-均值、SVM、感知器、神经网络、线性判别分析、主成分分析对于标准化数据可能会表现更好。


尽管如此,我还是建议你要理解你的数据和对其将要使用的算法类型。过一段时间后,你会有能力判断出是否要对数据进行标准化操作。


备注:在缩放和标准化中二选一是个令人困惑的选择,你必须对数据和要使用的学习模型有更深入的理解,才能做出决定。对于初学者,你可以两种方法都尝试下并通过交叉验证精度来做出选择。


资料:阅读本文(https://www.analyticsvidhya.com/blog/2015/11/improve-model-performance-cross-validation-in-python-r/)会对交叉验证有更好的理解



练习2


尝试利用SVM模型做相同的练习,并请在评论区留下标准化前后的精度。


资料:阅读本文(https://www.analyticsvidhya.com/blog/2015/10/understaing-support-vector-machine-example-code/)会对SVM有更好的理解。


标签编码

在前面的章节里,我们对连续数字特征做了预处理。但是,我们的数据集还有其它特征,如性别(Gender)、婚否(Married)、供养人(Dependents)、自雇与否(Self-Employed)和教育程度(Education)。所有这些类别特征的值是字符型的。例如,性别(Gender)有两个层次,或者是男性(Male),或者是女性(Female)。让我们把这些特征放进我们的逻辑回归模型中。


#在整个数据集上拟合放逻辑回归模型

>> log=LogisticRegression(penalty='l2',C=.01)

>> log.fit(X_train,Y_train)

#检查模型的精度

>> accuracy_score(Y_test,log.predict(X_test))


Out : ValueError: could not convert string to float: Semiurban


我们得到一个错误信息:不能把字符型转换成浮点型。因此,这里真正在发生的事是像逻辑回归和基于距离的学习模式,如KNN、SVM、基于树的方法等等,在Sklearn中需要数字型数组。拥有字符型值的特征不能由这些学习模式来处理。


Sklearn提供了一个非常有效的工具把类别特征层级编码成数值。LabelEncoder用0到n_classes-1之间的值对标签进行编码。


让我们对所有的类别特征进行编码。


#导入LabelEncoder并初始化

>> from sklearn.preprocessing import LabelEncoder

>> le=LabelEncoder()

#遍历在训练和测试集中所有的公共列

>> for col in X_test.columns.values:

        #只对类别变量编码

       if X_test[col].dtypes=='object':

        #利用全部数据形成一个所有层次的详尽清单

       data=X_train[col].append(X_test[col])

       le.fit(data.values)

       X_train[col]=le.transform(X_train[col])

       X_test[col]=le.transform(X_test[col])


我们所有的类别特征都已编码。用X_train.head()可以查看更新了的数据集。我们将看下性别(Gender)在编码前后的频率分布。

Before : Male 318 

         Female 66 

Name: Gender, dtype: int64


After : 1 318 

        0 66 

Name: Gender, dtype: int64


现在我们已经完成了标签编码,让我们在同时有着类别和连续特征的数据集上运行逻辑回归模型。


#特征标准化

>> X_train_scale=scale(X_train)

>> X_test_scale=scale(X_test)

#拟合逻辑回归模型

>> log=LogisticRegression(penalty='l2',C=.01)

>> log.fit(X_train_scale,Y_train)

#检查模型的精度

>> accuracy_score(Y_test,log.predict(X_test_scale))


Out : 0.75


现在可以用了。但是,精度仍然和我们从数字特征标准化之后用逻辑回归得到的一样。这意味着我们加入的类别特征在我们的目标函数中不是非常显著。


练习3

试试用所有的特征作为非独立变量进行决策树分类,并评论一下你得到的精度。


资料:浏览本文(https://www.analyticsvidhya.com/blog/2016/04/complete-tutorial-tree-based-modeling-scratch-in-python/)中关于决策树的内容以更好地理解。


一位有效编码(One-Hot-Encoding,主要是采用位状态寄存器来对某个状态进行编码,每个状态都有自己独立的寄存器位,并且在任意时候只有一位有效——译者注)。


一位有效编码把每个带有n个可能值的类别特征转换成n个二进制特征,只有一个是有效的。


大多数机器学习算法不是为每个特征设置单个权重就是计算样本之间的距离。如线性模型算法(例如:逻辑回归)属于第一类。


让我们看一看一个来自loan_prediction数据集的例子。特征从属(Feature Dependents)有4个可能的值:0、1、2和3+,这些是编过码的,没有丢掉0、1、2和3的一般性。


在线性分类器中,我们就分配一个权重“W”给这个特征,这将在W*Dependents+K>0或相当于W*Dependents < span="" style="margin: 0px; padding: 0px; max-width: 100%; box-sizing: border-box !important; word-wrap: break-word !important;"> <>

Let f(w)= W*Dependents

让f(w)=W*Dependents

Possible values that can be attained by the equation are 0, W, 2W and 3W. A problem with this equation is that the weight “W” cannot make decision based on four choices. It can reach to a decision in following ways:

由方程获得的可能值是0、W、2W和3W。这个方程的一个问题是权重W不能在4个选择的基础上得到。它可以用下面的方法来决定:

•        所有导致同样的决定(所有的值 < span="" style="margin: 0px; padding: 0px; max-width: 100%; box-sizing: border-box !important; word-wrap: break-word !important;"> <>

•        3:2的层级分配(当f(w)>2W时的决策区间)

•        2:2的层级分配(当f(w)>W时的决策区间)


这里我们可以看到丢失了许多不同的可能决策,比如:“0”和“2W”应该给予相同的标签,“3W”和“W”是额外的。


这个问题可以通过一位有效编码来解决,因为它有效地把特征“从属”的维度从1变成4,这样特征“从属”的每个值将有自己的权重。更新了的决策方程将是f'(w) < K。


这里,f'(w) = W1*D_0 + W2*D_1 + W3*D_2 + W4*D_3


所有4个新变量有布尔型值(0或1)。


同样的事发生在基于距离的方法中,如KNN。没有编码,“0”和“1”从属值之间的距离是1,在“0”和“3+”之间的距离是3,这不是所期望的,因为这两个距离应该类似。在编码后,值将有新特征(列序列是0,1,2,3+):[1,0,0,0]和[0,0,0,1](最初我们找到的在“0”和“3+”之间的距离),现在这个距离将会是√2。


对于基于树的方法,同样的情况(在一个特征中有2个以上的值)可能在一定程度上影响输出,但是如果像随机森林的方法,若有足够深的深度,无需一位有效编码就能够处理类别变量。


现在,让我们看下不同算法中的一位有效编码的实现。


让我们创建一个逻辑回归模型用于分类,而不使用一位有效编码。


#我们使用的是缩放后的变量,因为我们看到在上一节中缩放会影响L1L2的正则化算法

>> X_train_scale=scale(X_train)

>> X_test_scale=scale(X_test)

拟合逻辑回归模型

>> log=LogisticRegression(penalty='l2',C=1)

>> log.fit(X_train_scale,Y_train)

检查模型的精度

>> accuracy_score(Y_test,log.predict(X_test_scale))


Out : 0.73958333333333337


现在,我们对数据进行编码。

>> from sklearn.preprocessing import OneHotEncoder

>> enc=OneHotEncoder(sparse=False)

>> X_train_1=X_train

>> X_test_1=X_test

>> columns=['Gender', 'Married', 'Dependents', 'Education','Self_Employed',

          'Credit_History', 'Property_Area']

>> for col in columns:

        #创建一个包含所有可能分类值的详尽列表

       data=X_train[[col]].append(X_test[[col]])

       enc.fit(data)

        #在训练数据上拟合一位有效编码

       temp = enc.transform(X_train[[col]])

        #用新列名把编了码的特征改为数据帧

       temp=pd.DataFrame(temp,columns=[(col+"_"+str(i)) for i in data[col]

            .value_counts().index])

        #在并排级联中,检索值应该是相同的。

        #设置与X_train数据帧类似的索引值

       temp=temp.set_index(X_train.index.values)

        #把新的一位有效编码了的变量加入训练数据帧

       X_train_1=pd.concat([X_train_1,temp],axis=1)

        #在测试数据上拟合一位有效编码

       temp = enc.transform(X_test[[col]])

        #把它变成数据帧并加上列名

       temp=pd.DataFrame(temp,columns=[(col+"_"+str(i)) for i in data[col]

            .value_counts().index])

        #设置合适的级联索引

       temp=temp.set_index(X_test.index.values)

        #把新的一位有效编码了的变量加入到测试数据帧

       X_test_1=pd.concat([X_test_1,temp],axis=1)


现在,让我们在一位有效编码了的数据上应用逻辑回归模型

#标准化数据集

>> X_train_scale=scale(X_train_1)

>> X_test_scale=scale(X_test_1)

#拟合逻辑回归模型

>> log=LogisticRegression(penalty='l2',C=1)

>> log.fit(X_train_scale,Y_train)

#检查模型的精度

>> accuracy_score(Y_test,log.predict(X_test_scale))

Out : 0.75


到此,我们再次得到最大的精度是0.75,这是我们迄今所能得到的。在这个例子中,逻辑回归正则(C)参数是1,早前我们用的是C=0.01。

原文发布时间为:2017-01-22


本文来自云栖社区合作伙伴“大数据文摘”,了解相关信息可以关注“BigDataDigest”微信公众号

相关文章
|
25天前
|
数据采集 数据可视化 数据挖掘
利用Python自动化处理Excel数据:从基础到进阶####
本文旨在为读者提供一个全面的指南,通过Python编程语言实现Excel数据的自动化处理。无论你是初学者还是有经验的开发者,本文都将帮助你掌握Pandas和openpyxl这两个强大的库,从而提升数据处理的效率和准确性。我们将从环境设置开始,逐步深入到数据读取、清洗、分析和可视化等各个环节,最终实现一个实际的自动化项目案例。 ####
|
22天前
|
存储 数据采集 人工智能
Python编程入门:从零基础到实战应用
本文是一篇面向初学者的Python编程教程,旨在帮助读者从零开始学习Python编程语言。文章首先介绍了Python的基本概念和特点,然后通过一个简单的例子展示了如何编写Python代码。接下来,文章详细介绍了Python的数据类型、变量、运算符、控制结构、函数等基本语法知识。最后,文章通过一个实战项目——制作一个简单的计算器程序,帮助读者巩固所学知识并提高编程技能。
|
2天前
|
数据采集 Web App开发 监控
Python爬虫:爱奇艺榜单数据的实时监控
Python爬虫:爱奇艺榜单数据的实时监控
|
22天前
|
小程序 开发者 Python
探索Python编程:从基础到实战
本文将引导你走进Python编程的世界,从基础语法开始,逐步深入到实战项目。我们将一起探讨如何在编程中发挥创意,解决问题,并分享一些实用的技巧和心得。无论你是编程新手还是有一定经验的开发者,这篇文章都将为你提供有价值的参考。让我们一起开启Python编程的探索之旅吧!
44 10
|
23天前
|
数据采集 分布式计算 大数据
构建高效的数据管道:使用Python进行ETL任务
在数据驱动的世界中,高效地处理和移动数据是至关重要的。本文将引导你通过一个实际的Python ETL(提取、转换、加载)项目,从概念到实现。我们将探索如何设计一个灵活且可扩展的数据管道,确保数据的准确性和完整性。无论你是数据工程师、分析师还是任何对数据处理感兴趣的人,这篇文章都将成为你工具箱中的宝贵资源。
|
1月前
|
算法 Unix 数据库
Python编程入门:从基础到实战
本篇文章将带你进入Python编程的奇妙世界。我们将从最基础的概念开始,逐步深入,最后通过一个实际的项目案例,让你真正体验到Python编程的乐趣和实用性。无论你是编程新手,还是有一定基础的开发者,这篇文章都将为你提供有价值的信息和知识。让我们一起探索Python的世界吧!
|
1月前
|
并行计算 调度 开发者
探索Python中的异步编程:从基础到实战
在Python的世界里,异步编程是一种让程序运行更加高效、响应更快的技术。本文不仅会介绍异步编程的基本概念和原理,还将通过具体代码示例展示如何在Python中实现异步操作。无论你是初学者还是有经验的开发者,都能从中获益,了解如何运用这一技术优化你的项目。
|
1月前
|
传感器 物联网 开发者
使用Python读取串行设备的温度数据
本文介绍了如何使用Python通过串行接口(如UART、RS-232或RS-485)读取温度传感器的数据。详细步骤包括硬件连接、安装`pyserial`库、配置串行端口、发送请求及解析响应等。适合嵌入式系统和物联网应用开发者参考。
54 3
|
1月前
|
数据处理 Python
探索Python中的异步编程:从基础到实战
在Python的世界中,“速度”不仅是赛车手的追求。本文将带你领略Python异步编程的魅力,从原理到实践,我们不单单是看代码,更通过实例感受它的威力。你将学会如何用更少的服务器资源做更多的事,就像是在厨房里同时烹饪多道菜而不让任何一道烧焦。准备好了吗?让我们开始这场技术烹饪之旅。
|
1月前
|
机器学习/深度学习 数据采集 数据可视化
Python数据科学实战:从Pandas到机器学习
Python数据科学实战:从Pandas到机器学习