三种集成学习算法原理及核心公式推导

简介: 本文主要介绍3种集成学习算法的原理及重要公式推导部分,包括随机森林(Random Forest)、自适应提升(AdaBoost)、梯度提升(Gradient Boosting)。仅对重点理论和公式推导环节做以简要介绍。

640.png集成学习3大流派


01 集成学习流派


在经典机器学习场景下,当单个学习模型性能不足以有效满足算法精度时,人们开始向集成学习模型发力——其思想和出发点很直观,就是三个臭皮匠赛过诸葛亮。


进一步地,根据这三个臭皮匠在致力于赛过诸葛亮期间的协作模式不同,集成学习又细分为bagging和boosting两大学派,其中前者是并行模式,意味着三个臭皮匠各搞各的然后将最后结果进行融合以期带来提升,这里bagging是一个合成词汇,本意为booststrap aggregating;后者是串行模型,先让一个臭皮匠去探探底,根据反馈结果再派第二个、第三个,颇有些车轮战的味道。进一步的,车轮战还可以根据具体战术落脚点的不同,而分为Adaboost和GB,其中AdaBoost算法中各个臭皮匠的迭代重点在于不断弥补自己的过失/错误的地方,而GB算法中各个臭皮匠的迭代重点在于不断弥补自己与诸葛亮/理想型的差距,或者说残差。


当然,除了bagging和boosting两大流派之外,其实集成学习还有第三大流派——stacking,从其名字可以看出端倪是有些堆栈的意味,其实质就是将前一轮学习之后的输出/标签组织变换后作为下一轮学习的输入/特征,然后再次训练。其实就是深度学习那一套……


02 随机森林


随机森林是典型的bagging流派集成学习算法,个人非常喜欢这个算法的名字:森林二字可以窥探出这个算法是源于决策树,当以多棵树作为弱学习器,那么组成的集成学习算法自然叫做森林最合适;而随机二字则形象的体现其在构建每棵决策树时,实际上是采取了随机采样的方式来确保各弱学习器结果的多样性。至于为什么要保证多样性,稍后的公式推导可以说明这一点。


这里,既然说随机森林是bagging流派的典型代表,那么言外之意就是还存在其他bagging算法,实际也正是如此。广义来讲,bagging流派集成学习算法可分为4类,主要差别主要源于采样方式的不同:


  • 仅对样本维度(体现为行采样)进行采样,且采样是有放回的,意味着每个弱学习器的K个采样样本中可能存在重复的样本,此时对应bagging算法,这里bagging=bootstrap aggregating。发现,这个具体算法的名字与bagging流派的名字重合,这并不意外,因为这是bagging中一种经典的采样方式,因而以其作为流派的名字。当然,bagging既是一种算法也是流派名,那就要看是狭义还是广义的bagging来加以区分了


  • 仅对样本维度采样,但采样是无放回的,意味着对于每个弱学习器的K个采样样本是完全不同的,由于相当于是每执行一次采样,则该样本就被舍弃掉(pass),所以此时算法叫做pasting


  • 前两者的随机性均来源于样本维度的随机采样,也就是行方向上,那么对于列方向上进行随机采样是否也可以呢?比如说每个弱学习器都选用所有的样本参与训练,但各学习器选取参与训练的特征是不一致的,训练出的算法自然也是具有随机性的,可以满足集成的要求。此时,对应算法叫做subspaces,中文译作子空间,具体说是特征维度的子空间,名字还是比较传神的


  • 发现,既有样本维度的随机,也有特征维度的随机,那么自然想到有没有兼顾二者的随机呢,也就是说每个弱学习器既执行行采样、也有列采样,得到的弱学习器其算法随机性应该更强。当然,这种算法被称作是patches,比如前文已经提到的随机森林就属于这种。实际上,随机森林才是最为广泛使用的bagging流派集成学习算法


发散一下,其实bagging算法无非就是区分到底是行采样、列采样还是行列采样,那为什么会出来4种呢?原来是行采样在采样执行过程中又细分了是否有放回。那既然行采样可以区分是否有放回,列采样是否也可区分一下衍生出两种具体算法呢?实际上是不可以的,因为特征的有放回意味着特征重复,在机器学习训练过程中重复的特征或者说复制特征参与训练不能带来学习效果的改变;但样本的重复或者说复制则能带来学习效果的差异,因为这直接带来算法训练样本的类别平衡性,自然影响训练结果。


bagging算法原理非常简单易懂,算法出发点也非常朴素,其中涉及到的公式推导也不多,主要需要理解以下三个要点:


  • 有放回采样中,预计有36.8%的样本不会被用于训练,也就意味着这些样本可以天然的作为测试集用于评估集成学习模型效果。36.8%是如何而来呢?设样本数量为N,则对任意样本而言在1次采样中不被抽中的概率为 1 - 1/N,在N次采样中均未被抽中的概率则为


微信截图_20220527161310.png


  • bagging流派的各种算法中,一直在强调每个弱学习器训练效果之间的随机性,或者说差异性,那么为什么要保证这种差异性呢?直观来看,把集成学习比作是一次考试,A同学在参考了周边的BCD3名同学答案的基础上做判断,如果BCD3人的答案都是一样的,正确的都正确,错误的都错误,那么A同学的在综合3人答案之后丝毫不会对最终结果带来任何提升,仍然是BCD三人的正确率。而只有当3人擅长的考点不一样得出的答案也不尽相同时,A综合之后才可能带来提升。具体到其中一道题,用数学语言描述就是:设BCD三人答对的概率分别为Pb、Pc、Pd,如果三人答案的正确与否是独立事件,那么A本着少数服从多数的原则,得到正确答案的概率为PbPc(1-Pd) + (1-Pb)PcPd + Pb(1-Pc)Pd + PbPcPd;反之,如果三人答案的正确与否是相关联的事件,即同时正确或错误,那么A的最终正确率其实是不变的。所以,bagging流派要求集成学习各弱学习器之间要求算法尽可能是随机的,同时也有准确率的最低要求,即不能低于50%,否则会对集成效果带来负作用


  • bagging流派将带来算法的低方差,但不会降低偏差。衡量机器学习模型效果时,往往需要考虑偏差和方差两方面因素,偏差意味着模型训练效果距离理想目标还有很大差距,常常效果很差;方差意味着模型训练效果时好时坏,虽然有时候效果很好,但不够稳定,二者的形象描述示意如下图所示。那么,当我们说bagging集成学习相较于单个弱学习器效果有所提升,具体是指哪方面有所提升呢,换言之是降低了偏差还是方差呢?答案是偏差不变,方差降低。这可以由以下两个公式意会的表达:


微信截图_20220527161321.png


当然,以上两个公式也仅可意会的表达集成学习效果的提升,而不能具体刻画。实际中的bagging算法当然是会带来偏差的降低,只是统计意义下的均值与弱学习器效果相当,具体如何降低偏差那就需要具体问题具体调参了!


640.png

理解机器学习方差和偏差

图片源于《Understanding the Bias-Variance Tradeoff 2012》一文

原文链接:http://scott.fortmann-roe.com/docs/BiasVariance.html


此外,bagging流派由于是并行训练多个模型,而后综合各个弱学习器的效果进行综合决策,那么综合决策的方式其实也是一个值得探讨的问题。简言之,就是hard voting和soft voting两种方式,以二分类问题为例,前者是直接统计所有弱学习器的结果,然后取其中结果最多的那一类作为最终结果;而后者则不是直接统计结果,而是统计各弱学习器的分类概率,并分别计算所有弱学习器中两类的概率之和,以概率之和较大者作为最终结果。具体不再赘述。


03 Adaboost


与bagging流派集成学习思想不同,boosting流派侧重于后人站在前人的肩膀上看问题,Adaboost和GB概莫能外。当然,虽然理论上最后一个弱学习器在综合了所有前人的经验基础上可能会具有最好的性能,但Adaboost也不是简单的以最后的弱学习器结果为准,而是仍然加权考虑所有弱学习器的结果。


前文提到,Adaboost算法在吸取前一个弱学习器训练效果的基础上,重点关注错误的样本从而针对性的执行后续训练过程,具体说就是加大前一轮训练错的样本权重。


以二分类为例,Adaboost算法核心算法原理(公式推导):

  • 损失函数:指数函数,当训练结果与真实标签一致时(乘积为正)损失较小,不一致时(乘积为负)损失较大


微信截图_20220527161401.png


  • m轮训练后的集成学习模型:前m个弱学习器的加权求和


微信截图_20220527161412.png


其中Gm(X)为第m轮的弱学习器,fm(X)和fm-1(X)分别为第m轮和第m-1轮训练后的集成学习模型,α为第m轮弱学习器的加权系数


  • m轮训练后的模型损失


微信截图_20220527161421.png


首先分析各弱学习器对集成学习模型的贡献权重。直观来看每个弱学习器的权重当然与其学习效果有关:学习效果好的话语权更重,意味着权重更大。那么具体如何计算呢?以第m轮训练的弱学习器权重α为参数,对集成学习模型效果损失函数求导得到如下公式:


微信截图_20220527161441.png


其中em表示第m轮弱学习器的模型训练错误率,根据错误样本的权重之和与总样本的权重之和比值得出。当然在样本权重和进行归一化的基础上,分母为1。

然后分析每轮训练弱学习器时的样本权重迭代方式。再看损失函数,每个样本的损失可看做是两部分的乘积形式,其中前一部分与第m轮训练的弱学习器Gm(X)无关,而后一轮与其直接相关,因而可将前一部分看做是样本加权系数——这也正是不断调整每一轮训练模型的样本权重的根源所在。记第m轮训练弱学习器的第i个样本的权重系数为:


微信截图_20220527161452.png


则可进一步发现权重的更新策略为:


微信截图_20220527161505.png


也就是说,下一轮的样本权重由前一轮样本权重和前一轮的弱学习器训练结果有关,进一步地,当前一轮训练结果与真实标签一致时样本权重前的系数为小于1的值,体现为样本权重降低;否则权重增大。当然,在每轮更新后还需对样本权重加和进行归一化处理。


弱学习器对集成学习最终结果的权重、每轮训练的样本权重,这也是Adaboost算法的灵魂所在,也正因如此,所以在sklearn库中,Adaboost算法的包名叫做_weight_boosting.py,形象的表达了Adaboost是基于权重的提升算法。

另外,与随机森林指定决策树作为弱学习器不同,Adaboost算法不限制具体的弱学习器,所以可以任何经典学习算法为基础进行训练,但默认的弱学习器其实也是决策树:对于分类任务,弱学习器为max_depth=1的决策树,俗称决策树桩;对于回归任务,弱学习器为max_depth=3 的决策树。


04 GB


与Adaboost算法类似,GB(梯度提升)集成学习算法也是基于多个弱学习器的训练效果的加权进行最终判决,且每轮训练也基于前一轮训练效果进行针对性的更新迭代。但与Adaboost聚焦于前一轮训练错误的样本机制不同,GB聚焦于前一轮训练后的残差,相当于是通过集成学习来了个算法接力,以使得最终学习效果不断逼近真实水平。


首先,仍然给出GB的损失函数,其是一个MSE形式:


微信截图_20220527161515.png


其次,给出第m轮后的集成学习模型,仍然是各弱学习器的加权求和:


微信截图_20220527161525.png


我们知道,在机器学习算法优化中,梯度下降法可用于求解最优参数,通过将参数不断朝着自身梯度的反方向进行迭代,最终可不断逼近最优结果。按照这一思想,每轮不断优化迭代的集成学习器fm(X)实际上也可看做是不断迭代优化的参数;另一方面,对于GB的MSE形式损失函数,以fm(X)为参数进行求导得到其梯度为:


微信截图_20220527161532.png


对比fm(X)的梯度和其迭代公式,可知第m轮训练的弱学习器实际上即为拟合当前残差的弱学习器,这也正是GB的核心思想。实际上,GB也被称作是函数空间的梯度下降法。


这是GB集成学习算法的出发点,也是理解不断拟合残差的核心。后续相关公式推导不再给出。


640.png





目录
相关文章
|
2天前
|
机器学习/深度学习 算法 PyTorch
【从零开始学习深度学习】43. 算法优化之Adam算法【RMSProp算法与动量法的结合】介绍及其Pytorch实现
【从零开始学习深度学习】43. 算法优化之Adam算法【RMSProp算法与动量法的结合】介绍及其Pytorch实现
|
1天前
|
机器学习/深度学习 存储 算法
【机器学习】深入探索机器学习:线性回归算法的原理与应用
【机器学习】深入探索机器学习:线性回归算法的原理与应用
|
2天前
|
存储 缓存 算法
LRU(Least Recently Used)算法原理
LRU(Least Recently Used)算法原理
3 0
|
2天前
|
机器学习/深度学习 传感器 算法
基于Mediapipe深度学习算法的手势识别系统【含python源码+PyqtUI界面+原理详解】-python手势识别 深度学习实战项目
基于Mediapipe深度学习算法的手势识别系统【含python源码+PyqtUI界面+原理详解】-python手势识别 深度学习实战项目
|
2天前
|
机器学习/深度学习 算法 PyTorch
【从零开始学习深度学习】42. 算法优化之AdaDelta算法【基于AdaGrad算法的改进】介绍及其Pytorch实现
【从零开始学习深度学习】42. 算法优化之AdaDelta算法【基于AdaGrad算法的改进】介绍及其Pytorch实现
|
2天前
|
机器学习/深度学习 算法 PyTorch
【从零开始学习深度学习】41. 算法优化之RMSProp算法【基于AdaGrad算法的改进】介绍及其Pytorch实现
【从零开始学习深度学习】41. 算法优化之RMSProp算法【基于AdaGrad算法的改进】介绍及其Pytorch实现
|
1天前
|
算法 JavaScript 决策智能
基于禁忌搜索算法的TSP路径规划matlab仿真
**摘要:** 使用禁忌搜索算法解决旅行商问题(TSP),在MATLAB2022a中实现路径规划,显示优化曲线与路线图。TSP寻找最短城市访问路径,算法通过避免局部最优,利用禁忌列表不断调整顺序。关键步骤包括初始路径选择、邻域搜索、解评估、选择及禁忌列表更新。过程示意图展示搜索效果。
|
1天前
|
机器学习/深度学习 算法
基于BP神经网络和小波变换特征提取的烟草香型分类算法matlab仿真,分为浓香型,清香型和中间香型
```markdown 探索烟草香型分类:使用Matlab2022a中的BP神经网络结合小波变换。小波分析揭示香气成分的局部特征,降低维度,PCA等用于特征选择。BP网络随后处理这些特征,以区分浓香、清香和中间香型。 ```
|
2天前
|
机器学习/深度学习 算法
m基于PSO-GRU粒子群优化长门控循环单元网络的电力负荷数据预测算法matlab仿真
摘要: 在MATLAB 2022a中,对比了电力负荷预测算法优化前后的效果。优化前为"Ttttttt111222",优化后为"Tttttttt333444",明显改进体现为"Tttttttttt5555"。该算法结合了粒子群优化(PSO)和长门控循环单元(GRU)网络,利用PSO优化GRU的超参数,提升预测准确性和稳定性。PSO模仿鸟群行为寻找最优解,而GRU通过更新门和重置门处理长期依赖问题。核心MATLAB程序展示了训练和预测过程,包括使用'adam'优化器和超参数调整,最终评估并保存预测结果。
6 0
|
3天前
|
算法 安全
基于龙格库塔算法的SIR病毒扩散预测matlab仿真
该程序使用龙格库塔算法实现SIR模型预测病毒扩散,输出易感、感染和康复人群曲线。在MATLAB2022a中运行显示预测结果。核心代码设置时间区间、参数,并定义微分方程组,通过Runge-Kutta方法求解。SIR模型描述三类人群动态变化,常微分方程组刻画相互转化。模型用于预测疫情趋势,支持公共卫生决策,但也存在局限性,如忽略空间结构和人口异质性。