xgboost系列丨xgboost原理及公式推导

举报
就挺突然 发表于 2021/01/05 11:13:22 2021/01/05
【摘要】 建树过程中如何选择使用哪个特征哪个值来进行分裂?什么时候停止分裂?如何计算叶节点的权值?建完了第一棵树之后如何建第二棵树?为防止过拟合,XGB做了哪些改进树的集成本文主要针对xgboost的论文原文中的公式细节做了详细的推导,对建树过程进行详细分析。对于样本个数为n特征个数为m的数据集 ,其中。树的集成学习方法使用K个增量函数来预测输出:为子模型的预测函数,每个即是一棵树。函数空间即树的搜索...

  • 建树过程中如何选择使用哪个特征哪个值来进行分裂?
  • 什么时候停止分裂?
  • 如何计算叶节点的权值?
  • 建完了第一棵树之后如何建第二棵树?
  • 为防止过拟合,XGB做了哪些改进




树的集成

图片
本文主要针对xgboost的论文原文中的公式细节做了详细的推导,对建树过程进行详细分析。
对于样本个数为n特征个数为m的数据集图片 ,其中图片
树的集成学习方法使用K个增量函数来预测输出:
图片

图片为子模型的预测函数,每个图片即是一棵树。


函数空间 图片即树的搜索空间。其中q为每棵树的结构,q将 图片中每个样本对应到唯一的叶节点上,最终产生T个叶节点, 图片则是该叶节点对应的权重,w即从节点到权重的映射(权重即叶节点的值)。每个 图片对应一个独立的树结构q和该树每个叶节点的权重w。(这里树结构是指每个分裂点和对应的分裂值)。 图片可以看做一个分段函数,q对应的不同的分段,w对应的为该分段的值,图片即分段到值的映射。
对我们的预测函数 图片,目标函数为:
图片
从公式1中可以看出,对于最终的预测函数图片,其参数为一个个的函数图片,因为参数为函数,所以图片无法使用传统的优化方法在欧氏空间中进行优化,而是采用了加法模型来进行训练。
boost的思想是将一系列弱分类器串行的组合起来,在前面的分类器的基础上迭代的优化新的分类器。
图片
首先我们对所有的数据默认预测一个固定值 图片(对应xgboost中参数base_score,注意并不等于base_score,而是经过Sigmoid函数映射后的值),在此基础上根据该预测值与真实y值的损失 ,建立第一棵树图片,之后每次迭代时都是根据其之前所有树做出的预测之和与真实y值的损失来建立新树。也就是每次迭代建树时用新树图片来优化前一个树的损失 。

图片为第t棵树对第i个样本做出的预测。我们每次添加新树的时候,要优化的目标函数为上一个树产生的损失。

因此我们建立第t棵树时有损失函数:
图片
图片为新建的这棵树做出的预测, 图片为之前所有的树预测值之和, 图片即是新建了当前这棵树后模型做出的预测值,求其与真实值 图片之间的损失(注意这里是损失不是残差,这里的 图片可以是log_loss, mse等)。



泰勒展开
图片
gbdt的目标函数与xgboost区别就是带不带正则项,也就是上面式子中的 图片gbdt对损失函数的优化是直接使用了损失函数的负梯度,沿着梯度下降的方向来减小损失,其是也就是一阶泰勒展开。而xgboost在这里使用了二阶泰勒展开,因为包含了损失函数的二阶信息,其优化的速度大大加快。
图片
下面来看一下泰勒展开的推导。首先我们来复习一下泰勒定理:
设n是一个正整数。如果定义在一个包含a的区间上的函数f在a点处n+1次可导,那么对于这个区间上的任意x,则有:图片其中的多项式称为函数在a处的泰勒展开式,剩余的图片是泰勒公式的余项,是图片的高阶无穷小。 
该公式经过变换图片可以得到二阶展开式:
图片
对于式子:
图片
可以这样分析, 图片为预测值 图片和真实值 图片之间的损失, 图片为常量,因此图片是以预测值图片为自变量的函数,当建立新树给出新的预测图片后,相当于在上一次的预测 图片上增加了一个无穷小量
图片
则有
图片
其中真实标签图片是常数, 图片是上次迭代求出的值即这里的 图片图片为无穷小量 图片有了这个对应之后。
图片
因此我们建立第t棵树时有损失函数:
图片
令损失函数的一阶、二阶偏导分别为图片,其中 图片图片
图片
式中图片为常量,优化的是损失函数的最小值,因此常量值可以从损失函数中去掉。上式可简化为:
图片

叶节点权重

图片
式中正则项图片进行展开,得:
图片
其中图片是新建的树的值,对于每个样本来说,就是对应的叶节点的权重 图片。定义 图片为分到叶节点 图片的样本(叶节点总数为T,样本总数为n)

上式是对本次建树时n个样本的损失求和,下面分两步:先对每个叶节点的样本损失求和,再对所有叶节点求和,两者结果一样。
图片
对于叶节点图片上的损失:
图片
对于当前的树结构求 图片使 图片最小,显然这是个一元二次方程求最小值问题。 图片
可以得到叶节点权重图片的最优值:
图片


分裂准则

图片
面是对单个叶节点计算出了最优权重,对于新建的这树(树结构 图片)在此权重下对应的的最小损失为每个叶节点上样本最小损失之和(将上式中的 图片代入):
图片
在树结构图片下产生的最优损失图片可以做为树结构的评价函数,也就是作为树分裂时候的评价指标。
图片为每次分裂时分到左子树上的样本, 图片为每次分裂时分到右子树上的样本,有 图片则在该次分裂后损失的减小量为:
图片
因此将分裂时增益定义为:
图片

我们在建树的过程(也就是求分段函数的过程)包括两步:一是选择分裂依据的特征和特征值(将自变量分段),二是确定叶节点的权重(确定每段对应的函数值)。划分的依据准则是Gain,其实也就是损失函数的解析解,划分后叶节点的权重图片是使函数达到解析解的权重图片
从最优化的角度来看:GBDT采用的是数值优化的思维, 用的最速下降法去求解Loss Function的最优解, 其中用CART决策树去拟合负梯度, 用牛顿法求步长。XGboost用的解析的思维, 对Loss Function展开到二阶近似, 求得解析解, 用解析解作为Gain来建立决策树, 使得Loss Function最优.

图片

除了对目标函数添加正则项外,为了减小过拟合,xgboost还使用了列采样和缩减方法(Shrinkage,即Learning rate)。

损失函数计算

图片
对于二分类问题常使用log损失作为损失函数,下面推导一下log loss的一阶梯度G和海森矩阵H。图片

图片
其中p为预测概率。 图片为预测值,则有:
图片
因此:
图片
即:
图片

声明:本文首发于华为NAIE《网络人工智能园地》微信公众号,如有转载,请注明出处。微信公众号二维码为:

二维码-15CM.jpg

本文首发:https://mp.weixin.qq.com/s?__biz=Mzg2MjExMDIzNw==&mid=2247492915&idx=1&sn=49ad15f532ac4b7e3414d5d739dcd065&chksm=ce0e7241f979fb571fc2dc563b7ea1b024bd074b571021f37481c389b466970b7abb17328be2&token=1640375181&lang=zh_CN#rd


【版权声明】本文为华为云社区用户原创内容,转载时必须标注文章的来源(华为云社区)、文章链接、文章作者等基本信息, 否则作者和本社区有权追究责任。如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。