源 | 程序猿的进击之路
序言
本文尽可能的不涉及到繁杂的数学公式,把面试中常问的模型核心点,用比较通俗易懂但又不是专业性的语言进行描述。希望可以帮助大家在找工作时提纲挈领的复习最核心的内容,或是在准备的过程中抓住每个模型的重点。
篇幅较长,可下载PDF慢慢研读
后台回复“20180704”
或者
文末点击阅读原文
实战环境说明:
Python 2.7
Sklearn 0.19.0
graphviz 0.8.1 决策树可视化
1.1 原理
顾名思义,决策树就是用一棵树来表示我们的整个决策过程。这棵树可以是二叉树(比如CART只能是二叉树),也可以是多叉树(比如ID3、C4.5可以是多叉树或二叉树)。
根节点包含整个样本集,每个叶节点都对应一个决策结果(注意,不同的叶节点可能对应同一个决策结果),每一个内部节点都对应一次决策过程或者说是一次属性测试。从根节点到每个叶节点的路径对应一个判定测试序列。
举个例子:
就像上面这个例子,训练集由三个特征:outlook(天气),humidity(湿度),windy(是否有风)。那么我们该如何选择特征对训练集进行划分那?连续型特征(比如湿度)划分的阈值又是如何确定的那?
决策树的生成就是不断的选择最优的特征对训练集进行划分,是一个递归的过程。递归返回的条件有三种:
(1)当前节点包含的样本属于同一类别,无需划分
(2)当前属性集为空,或所有样本在属性集上取值相同,无法划分
(3)当前节点包含样本集合为空,无法划分
1.2 ID3、C4.5、CART
这三个是非常著名的决策树算法。简单粗暴来说,ID3使用信息增益作为选择特征的准则;C4.5使用信息增益比作为选择特征的准则;CART使用Gini指数作为选择特征的准则。
ID3
熵表示的是数据中包含的信息量大小。熵越小,数据的纯度越高,也就是说数据越趋于一致,这是我们希望的划分之后每个子节点的样子。
信息增益 = 划分前熵 - 划分后熵。信息增益越大,则意味着使用属性a来进行划分所获得的“纯度提升”越大。也就是说,用属性a来划分训练集,得到的结果中纯度比较高。
ID3仅仅适用于二分类问题。ID3仅仅能够处理离散属性。
C4.5
C4.5克服了ID3仅仅能够处理离散属性的问题,以及信息增益偏向选择取值较多特征的问题,使用信息增益比来选择特征。信息增益比 = 信息增益 /划分前熵 选择信息增益比最大的作为最优特征。
C4.5处理连续特征是先将特征取值排序,以连续两个值中间值作为划分标准。尝试每一种划分,并计算修正后的信息增益,选择信息增益最大的分裂点作为该属性的分裂点。
CART
CART与ID3,C4.5不同之处在于CART生成的树必须是二叉树。也就是说,无论是回归还是分类问题,无论特征是离散的还是连续的,无论属性取值有多个还是两个,内部节点只能根据属性值进行二分。
CART的全称是分类与回归树。从这个名字中就应该知道,CART既可以用于分类问题,也可以用于回归问题。
回归树中,使用平方误差最小化准则来选择特征并进行划分。每一个叶子节点给出的预测值,是划分到该叶子节点的所有样本目标值的均值,这样只是在给定划分的情况下最小化了平方误差。要确定最优化分,还需要遍历所有属性,以及其所有的取值来分别尝试划分并计算在此种划分情况下的最小平方误差,选取最小的作为此次划分的依据。由于回归树生成使用平方误差最小化准则,所以又叫做最小二乘回归树。
分类树种,使用Gini指数最小化准则来选择特征并进行划分;
Gini指数表示集合的不确定性,或者是不纯度。基尼指数越大,集合不确定性越高,不纯度也越大。这一点和熵类似。另一种理解基尼指数的思路是,基尼指数是为了最小化误分类的概率。
1.3 信息增益 vs 信息增益比
之所以引入了信息增益比,是由于信息增益的一个缺点。那就是:信息增益总是偏向于选择取值较多的属性。信息增益比在此基础上增加了一个罚项,解决了这个问题。
1.4 Gini指数 vs 熵
既然这两个都可以表示数据的不确定性,不纯度。那么这两个有什么区别那?
Gini指数的计算不需要对数运算,更加高效
Gini指数更偏向于连续属性,熵更偏向于离散属性
1.5 剪枝
决策树算法很容易过拟合(overfitting),剪枝算法就是用来防止决策树过拟合,提高泛华性能的方法。
剪枝分为预剪枝与后剪枝。
预剪枝是指在决策树的生成过程中,对每个节点在划分前先进行评估,若当前的划分不能带来泛化性能的提升,则停止划分,并将当前节点标记为叶节点。
后剪枝是指先从训练集生成一颗完整的决策树,然后自底向上对非叶节点进行考察,若将该节点对应的子树替换为叶节点,能带来泛化性能的提升,则将该子树替换为叶节点。
那么怎么来判断是否带来泛化性能的提升那?最简单的就是留出法,即预留一部分数据作为验证集来进行性能评估。
1.6 总结
决策树算法主要包括三个部分:特征选择、树的生成、树的剪枝。常用算法有ID3、C4.5、CART。
特征选择。特征选择的目的是选取能够对训练集分类的特征。特征选择的关键是准则:信息增益、信息增益比、Gini指数。
决策树的生成。通常是利用信息增益最大、信息增益比最大、Gini指数最小作为特征选择的准则。从根节点开始,递归的生成决策树。相当于是不断选取局部最优特征,或将训练集分割为基本能够正确分类的子集。
决策树的剪枝。决策树的剪枝是为了防止树的过拟合,增强其泛化能力。包括预剪枝和后剪枝。
要说随机森林就要先说Bagging,要说Bagging就要先说集成学习。
2.1 集成学习方法
集成学习(ensemble learning)通过构建并组合多个学习器来完成学习任务。集成学习通过将多个学习器进行结合,常获得比单一学习器显著优越的泛化性能。
根据个体学习器是否是同类型的学习器(由同一个算法生成,比如C4.5,BP等),分为同质和异质。同质的个体学习器又叫做基学习器,而异质的个体学习器则直接成为个体学习器。
原则:要获得比单一学习器更好的性能,个体学习器应该好而不同。即个体学习器应该具有一定的准确性,不能差于弱学习器,并且具有多样性,即学习器之间有差异。
根据个体学习器的生成方式,目前集成学习分为两大类:
个体学习器之间存在强依赖关系、必须串行生成的序列化方法。代表是Boosting。
个体学习器之间不存在强依赖关系、可同时生成的并行化方法。代表是Bagging和随机森林(Random Forest)
2.2 Bagging
前面提到,想要集成算法获得性能的提升,个体学习器应该具有独立性。虽然“独立”在现实生活中往往无法做到,但是可以设法让基学习器尽可能的有较大的差异。Bagging给出的做法就是对训练集进行采样,产生出若干个不同的子集,再从每个训练子集中训练一个基学习器。由于训练数据不同,我们的基学习器可望具有较大的差异。
Bagging是并行式集成学习方法的代表,采样方法是自助采样法,用的是有放回的采样。初始训练集中大约有63.2%的数据出现在采样集中。
Bagging在预测输出进行结合时,对于分类问题,采用简单投票法;对于回归问题,采用简单平均法。
Bagging优点:
高效。Bagging集成与直接训练基学习器的复杂度同阶。
Bagging能不经修改的适用于多分类、回归任务。
包外估计。使用剩下的样本作为验证集进行包外估计(out-of-bag estimate)
Bagging主要关注降低方差。(low variance)
2.3 随机森林(Random Forest)
2.3.1 原理
随机森林(Random Forest)是Bagging的一个变体。Ramdon Forest在以决策树为基学习器构建Bagging集成的基础上,进一步在决策树的训练过程中引入随机属性选择。
原来决策树从所有属性中,选择最优属性。Ramdom Forest的每一颗决策树中的每一个节点,先从该节点的属性集中随机选择K个属性的子集,然后从这个属性子集中选择最优属性进行划分。K控制了随机性的引入程度,是一个重要的超参数。
预测:
分类:简单投票法
回归:简单平均法
2.3.2 优缺点
优点:
由于每次不再考虑全部的属性,而是一个属性子集,所以相比于Bagging计算开销更小,训练效率更高。
由于增加了属性的扰动,随机森林中基学习器的性能降低,使得在随机森林在起始时候性能较差,但是随着基学习器的增多,随机森林通常会收敛于更低的泛化误差,相比于Bagging。
两个随机性的引入,使得随机森林不容易陷入过拟合,具有很好的抗噪声能力
对数据的适应能力强,可以处理离散和连续的,无需要规范化
可以得到变量的重要性, 基于oob误分类率和基于Gini系数的变化
缺点:
在噪声较大的时候容易过拟合
AdaBoost是Boosting的代表,前面我们提到Boosting是集成学习中非常重要的一类串行化学习方法。
3.1 Boosting
Boosting是指个体学习器之间存在强依赖关系,必须串行序列化生成的集成学习方法。他的思想来源是三个臭皮匠顶个诸葛亮。Boosting意为提升,意思是希望将每个弱学习器提升为强学习器。
工作机制如下:
先从初始训练集中学习一个基学习器
根据基学习器的表现对训练样本分布进行调整,使得先前基学习器做错的训练样本在后续收到更多关注
基于调整后的样本分布来训练下一个基学习器
如此反复,直到基学习器数目达到T,最终将这T个基学习器进行加权结合
对训练样本分布调整,主要是通过增加误分类样本的权重,降低正确分类样本的权重。
Boosting关注的主要是降低偏差。(low bias)
3.2 AdaBoost原理
面临两个问题:
(1)在每一轮,如何改变训练数据的概率分布或者权值分布。(也可以重采样,但是AdaBoost没这么做)
(2)如何将弱分类器组合成强分类器
AdaBoost的做法:
(1)提高那些被前一轮弱分类器错误分类样本的权值,降低那些被正确分类的样本的权值。这样一来,那些没有得到正确分类的数据,由于其权值的加大而受到后一轮弱分类器的关注。
(2)采用加权多数表决。具体的,加大分类错误率低的分类器的权值,使其在表决中起较大作用,减少分类误差率大的弱分类器的权值,使其在表决中起较小作用。
弱分类器被线性组合成为一个强分类器。
训练目标:
最小化指数损失函数。
三部分组成:
(1)分类器权重更新公式
(2)样本分布(也就是样本权重)更新公式
(3)加性模型。 最小化指数损失函数
3.3 AdaBoost优缺点
优点:
不改变所给的训练数据,而不断改变训练数据的权值分布,使得训练数据在基本分类器的学习中起不同的作用。这是AdaBoost的一个特点。
利用基本分类器的加权线性组合构建最终分类器,是AdaBoost的另一个特点。
AdaBoost被实践证明是一种很好的防止过拟合的方法,但至今为什么至今没从理论上证明。
缺点:
AdaBoost只适用于二分类问题。
GBDT(Gradient Boosting Decision Tree)又叫MART(Multiple Additive Regression Tree)。是一种迭代的决策树算法。
本文从以下几个方面进行阐述:
Regression Decision Tree(DT)
Gradient Boosting(GB)
Shrinkage(算法的一个重要演进分支,目前大部分源码都是基于该版本实现)
GBDT适用范围
与随机森林的对比
4.1 DT:回归树 Regression Decision Tree
GDBT中的树全部都是回归树,核心就是累加所有树的结果作为最终结果。只有回归树的结果累加起来才是有意义的,分类的结果加是没有意义的。
GDBT调整之后可以用于分类问题,但是内部还是回归树。
这部分和决策树中的是一样的,无非就是特征选择。回归树用的是最小化均方误差,分类树是用的是最小化基尼指数(CART)
4.2 GB:梯度迭代 Gradient Boosting
首先Boosting是一种集成方法。通过对弱分类器的组合得到强分类器,他是串行的,几个弱分类器之间是依次训练的。GBDT的核心就在于,每一颗树学习的是之前所有树结论和的残差。
Gradient体现在:无论前面一颗树的cost function是什么,是均方差还是均差,只要它以误差作为衡量标准,那么残差向量都是它的全局最优方向,这就是Gradient。
4.3 Shrinkage
Shrinkage(缩减)是GBDT算法的一个重要演进分支,目前大部分的源码都是基于这个版本的。
核心思想在于:Shrinkage认为每次走一小步来逼近结果的效果,要比每次迈一大步很快逼近结果的方式更容易防止过拟合。也就是说,它不信任每次学习到的残差,它认为每棵树只学习到了真理的一小部分,累加的时候只累加一小部分,通过多学习几棵树来弥补不足。
具体的做法就是:仍然以残差作为学习目标,但是对于残差学习出来的结果,只累加一小部分(step*残差)逐步逼近目标,step一般都比较小0.01-0.001,导致各个树的残差是渐变而不是陡变的。
本质上,Shrinkage为每一颗树设置了一个weight,累加时要乘以这个weight,但和Gradient没有关系。这个weight就是step。跟AdaBoost一样,Shrinkage能减少过拟合也是经验证明的,目前还没有理论证明。
4.4 GBDT适用范围
GBDT可以适用于回归问题(线性和非线性),相对于logistic regression仅能用于线性回归,GBDT适用面更广。
GBDT也可用于二分类问题(设定阈值,大于为正,否则为负)和多分类问题。
4.5 GBDT和随机森林
GBDT和随机森林的相同点:
都是由多棵树组成
最终的结果都由多棵树共同决定
GBDT和随机森林的不同点:
组成随机森林的可以是分类树、回归树;组成GBDT只能是回归树
组成随机森林的树可以并行生成(Bagging);GBDT只能串行生成(Boosting)
对于最终的输出结果而言,随机森林使用多数投票或者简单平均;而GBDT则是将所有结果累加起来,或者加权累加起来
随机森林对异常值不敏感,GBDT对异常值非常敏感
随机森林对训练集一视同仁权值一样,GBDT是基于权值的弱分类器的集成
随机森林通过减小模型的方差提高性能,GBDT通过减少模型偏差提高性能
TIP
1. GBDT相比于决策树有什么优点
泛化性能更好!GBDT的最大好处在于,每一步的残差计算其实变相的增大了分错样本的权重,而已经分对的样本则都趋向于0。这样后面就更加专注于那些分错的样本。
2. Gradient体现在哪里?
可以理解为残差是全局最优的绝对方向,类似于求梯度。
3. re-sample
GBDT也可以在使用残差的同时引入Bootstrap re-sampling,GBDT多数实现版本中引入了这个选项,但是是否一定使用有不同的看法。原因在于re-sample导致的随机性,使得模型不可复现,对于评估提出一定的挑战,比如很难确定性能的提升是由于feature的原因还是sample的随机因素。
LR原理
参数估计
LR的正则化
为什么LR能比线性回归好?
LR与MaxEnt的关系
5.1 LR模型原理
首先必须给出Logistic分布:
u是位置参数,r是形状参数。对称点是(u,1/2),r越小,函数在u附近越陡峭。
然后,二分类LR模型,是参数化的logistic分布,使用条件概率来表示:
然后,一个事件的几率(odds):指该事件发生与不发生的概率比值:
对数几率:
那么对于逻辑回归而言,Y=1的对数几率就是:
最终,输出Y=1的对数几率是输入x的线性函数表示的模型,这就是逻辑回归模型。
5.2 参数估计
在统计学中,常常使用极大似然估计法来估计参数。即找到一组参数,使得在这组参数下,我们数据的似然度(概率)最大。
似然函数:
对数似然函数:
对应的损失函数:
5.3 最优化方法
逻辑回归模型的参数估计中,最后就是对J(W)求最小值。这种不带约束条件的最优化问题,常用梯度下降或牛顿法来解决
使用梯度下降法求解逻辑回归参数估计
求J(W)梯度:g(w):
更新Wk:
$$ W{k+1} = Wk - lambda*g(W_k) $$
5.4 正则化
正则化为了解决过拟合问题。分为L1和L2正则化。主要通过修正损失函数,加入模型复杂性评估。
正则化是符合奥卡姆剃刀原理:在所有可能的模型中,能够很好的解释已知数据并且十分简单的才是最好的模型。
p表示范数,p=1和p=2分别应用L1和L2正则。
L1正则化。向量中各元素绝对值之和。又叫做稀疏规则算子(Lasso regularization)。关键在于能够实现特征的自动选择,参数稀疏可以避免非必要的特征引入的噪声。
L2正则化。使得每个元素都尽可能的小,但是都不为零。在回归里面,有人把他的回归叫做岭回归(Ridge Regression),也有人叫他“权值衰减”(weight decay)
一句话总结就是:L1会趋向于产生少量的特征,而其他的特征都是0,而L2会选择更多的特征,这些特征都会接近于0.
5.5 逻辑回归 vs 线性回归
首先,逻辑回归比线性回归要好。
两者都属于广义线性模型。
线性回归优化目标函数用的最小二乘法,而逻辑回归用的是最大似然估计。
逻辑回归只是在线性回归的基础上,将加权和通过sigmoid函数,映射到0-1范围内空间。
线性回归在整个实数范围内进行预测,敏感度一致,而分类范围,需要在[0,1]。而逻辑回归就是一种减小预测范围,将预测值限定为[0,1]间的一种回归模型。逻辑曲线在z=0时,十分敏感,在z>>0或z<<0处,都不敏感,将预测值限定为(0,1)。逻辑回归的鲁棒性比线性回归要好。
5.6 逻辑回归模型 vs 最大熵模型 MaxEnt
简单粗暴的说:逻辑回归跟最大熵模型没有本质区别。逻辑回归是最大熵对应为二类时的特殊情况,也就是说,当逻辑回归扩展为多类别的时候,就是最大熵模型。
最大熵原理:学习概率模型的时候,在所有可能的概率模型(分布)中,熵最大的模型是最好的模型。
虽然咱们的目标是尽可能的不涉及到公式,但是提到SVM就没有办法不涉及到公式推导,因为面试中只要问到SVM,最基本也是最难的问题就是:SVM的对偶问题数学公式推导。所以想学好机器学习,还是要塌下心来,不仅仅要把原理的核心思想掌握了,公式推导也要好好学习才行。:)
6.1 SVM原理
简单粗暴的说:SVM的思路就是找到一个超平面将数据集进行正确的分类。对于在现有维度不可分的数据,利用核函数映射到高纬空间使其线性可分。
支持向量机SVM是一种二分类模型。它的基本模型是定义在特征空间上的间隔最大的线性分类器,间隔最大使它有别于感知机。SVM的学习策略是间隔最大化,可形式化为求解凸二次规划问题。
SVM分为:
线性可分支持向量机。当训练数据线性可分时,通过硬间隔最大化,学习到的一个线性分类器。
线性支持向量机。当训练数据近似线性可分时,通过软间隔最大化,学习到的一个线性分类器。
非线性支持向量机。当训练数据线性不可分,通过使用核技巧及软间隔最大化,学习非线性支持向量机。
上图中,X表示负例,O表示正例。此时的训练数据可分,线性可分支持向量机对应着将两类数据正确划分并且间隔最大的直线。
6.1.1 支持向量与间隔
支持向量:在线性可分的情况下,训练数据样本集中的样本点中与分离超平面距离最近的样本点的实例称为支持向量(support vector)。
函数间隔定义如下:
yi表示目标值,取值为+1、-1. 函数间隔虽然可以表示分类预测的准确性以及确信度。但是有个不好的性质:只要成倍的改变W和B,虽然此时的超平面并没有改变,但是函数间隔会变大。所以我们想到了对超平面的法向量W施加一些约束,比如规范化,使得间隔确定,这就引出了几何间隔:
支持向量学习的基本思想就是求解能够正确划分训练数据集并且几何间隔最大的分类超平面。
6.1.2 对偶问题
为了求解线性可分支持向量机的最优化问题:
将它作为原始最优化问题,应用拉格朗日对偶性,通过求解对偶问题得到原始问题的最优解,这就是线性可分支持向量机的对偶算法。
本来的算法也可以求解SVM,但是之所以要用对偶问题来求解,优点是:
一是对偶问题往往更容易求解
二是自然引入核函数,进而推广到非线性分类问题
说点题外话,这也是面试中会被问到的一个问题:原始问题既然可以求解,为什么还要转换为对偶问题。答案就是上述这两点。由于篇幅的愿意,此处就不在展开数学公式的推导了,大家可以参考《机器学习》与《统计学习方法》
6.1.3 核函数
对于在原始空间中不可分的问题,在高维空间中是线性可分的。
对于线性不可分的问题,使用核函数可以从原始空间映射到高纬空间,使得问题变得线性可分。核函数还可以使得在高维空间计算的内积在低维空间中通过一个函数来完成。
常用的核函数有:高斯核、线性核、多项式核。
6.1.4 软间隔
线性可分问题的支持向量机学习方法,对现行不可分训练数据是不适用的。所以讲间隔函数修改为软间隔,对于函数间隔,在其上加上一个松弛变量,使其满足大于等于1。约束条件变为:
6.2 优缺点
缺点:
时空开销比较大,训练时间长
核函数的选取比较难,主要靠经验
优点:
在小训练集上往往得到比较好的结果
使用核函数避开了高纬空间的复杂性
泛化能力强
七、利用sklearn进行实战
使用sklearn用决策树来进行莺尾花数据集的划分问题。代码上没有固定随机种子,所以每次运行的结果会稍有不同。
数据集三维可视化:
在Sepal length和Sepal width二维上的可视化:
原文链接:https://www.toutiao.com/i6571035846022005261/?tt_from=weixin&utm_campaign=client_share&from=singlemessage×tamp=1530418475&app=news_article_lite&utm_source=weixin&iid=36909045601&utm_medium=toutiao_android&wxshare_count=2&pbid=6570804297075213832