深入理解GBDT二分类算法
- 我的个人微信公众号: Microstrong
- 微信公众号ID: MicrostrongAI
- 微信公众号介绍: Microstrong(小强)同学主要研究机器学习、深度学习、计算机视觉、智能对话系统相关内容,分享在学习过程中的读书笔记!期待您的关注,欢迎一起学习交流进步!
- 我的知乎主页: https://www.zhihu.com/people/MicrostrongAI/activities
- Github: https://github.com/Microstrong0305
- 个人博客: https://blog.csdn.net/program_developer
- 本文首发在我的微信公众号里,地址:https://mp.weixin.qq.com/s/4T6OD2Hr5MzNI0j-mT3e-g,如有公式和图片不清楚,可以在我的微信公众号里阅读。
目录:
- GBDT分类算法简介
- GBDT二分类算法
2.1 逻辑回归的对数损失函数
2.2 GBDT二分类原理 - GBDT二分类算法实例
- 手撕GBDT二分类算法
4.1 用Python3实现GBDT二分类算法
4.2 用sklearn实现GBDT二分类算法 - GBDT分类任务常见的损失函数
- 总结
- Reference
本文的主要内容概览:
1. GBDT分类算法简介
GBDT无论用于分类还是回归,一直使用的是CART回归树。GBDT不会因为我们所选择的任务是分类任务就选用分类树,这里的核心原因是GBDT每轮的训练是在上一轮训练模型的负梯度值基础之上训练的。这就要求每轮迭代的时候,真实标签减去弱分类器的输出结果是有意义的,即残差是有意义的。如果选用的弱分类器是分类树,类别相减是没有意义的。对于这样的问题,可以采用两种方法来解决:
- 采用指数损失函数,这样GBDT就退化成了Adaboost,能够解决分类的问题;
- 使用类似于逻辑回归的对数似然损失函数,如此可以通过结果的概率值与真实概率值的差距当做残差来拟合;
下面我们就通过二分类问题,去看看GBDT究竟是如何做分类的。
2. GBDT二分类算法
2.1 逻辑回归的对数损失函数
逻辑回归的预测函数为:
hθ(x)=11+e−θTxh_{\theta}(x) = \frac{1}{1+ e^{-\theta^{T} x}} hθ(x)=1+e−θTx1
函数hθ(x)h_{\theta}(x)hθ(x) 的值有特殊的含义,它表示结果取 111 的概率,因此对于输入 xxx 分类结果为类别 111 和类别 000 的概率分别为:
P(Y=1∣x;θ)=hθ(x)P(Y=1|x;\theta) = h_{\theta}(x) P(Y=1∣x;θ)=hθ(x)
P(Y=0∣x;θ)=1−hθ(x)P(Y=0|x;\theta) =1 - h_{\theta}(x) P(Y=0∣x;θ)=1−hθ(x)
下面我们根据上式,推导出逻辑回归的对数损失函数 L(θ)L(\theta)L(θ) 。上式综合起来可以写成:
P(Y=y∣x;θ)=(hθ(x))y(1−hθ(x))1−yP(Y=y|x;\theta) = (h_{\theta}(x))^{y}(1-h_{\theta}(x))^{1-y} P(Y=y∣x;θ)=(hθ(x))y(1−hθ(x))1−y
然后取似然函数为:
l(θ)=∏i=1NP(Y=yi∣xi;θ)=∏i=1N(hθ(xi))yi(1−hθ(xi))1−yil(\theta) = \prod_{i=1}^{N}P(Y=y_{i}|x_{i};\theta)=\prod_{i=1}^{N}(h_{\theta}(x_{i}))^{y_{i}}(1-h_{\theta}(x_{i}))^{1- y_{i}} l(θ)=i=1∏NP(Y=yi∣xi;θ)=i=1∏N(hθ(xi))yi(1−hθ(xi))1−yi
因为 l(θ)l(\theta)l(θ) 和 logl(θ)logl(\theta)logl(θ) 在同一 θ\thetaθ 处取得极值,因此我们接着取对数似然函数为:
L(θ)=∑i=1N[yiloghθ(xi)+(1−yi)log(1−hθ(xi))]L(\theta) = \sum_{i=1}^{N}\left[ {y_{i}logh_{\theta}(x_{i})} + (1-y_{i})log(1-h_{\theta}(x_{i})) \right] L(θ)=i=1∑N[yiloghθ(xi)+(1−yi)log(1−hθ(xi))]
最大似然估计就是求使L(θ)L(\theta)L(θ) 取最大值时的 θ\thetaθ 。这里对 L(θ)L(\theta)L(θ) 取相反数,可以使用梯度下降法求解,求得的 θ\thetaθ 就是要求的最佳参数:
J(θ)=−1NL(θ)=−1N∑i=1N[yiloghθ(xi)+(1−yi)log(1−hθ(xi))]J(\theta) = - \frac{1}{N}L(\theta) =- \frac{1}{N} \sum_{i=1}^{N}\left[ {y_{i}logh_{\theta}(x_{i})} + (1-y_{i})log(1-h_{\theta}(x_{i})) \right] J(θ)=−N1L(θ)=−N1i=1∑N[yiloghθ(xi)+(1−yi)log(1−hθ(xi))]
2.2 GBDT二分类原理
逻辑回归单个样本 (xi,yi)(x_{i},y_{i})(xi,yi) 的损失函数可以表达为:
L(θ)=−yilogyi^−(1−yi)log(1−yi^)L(\theta) = -y_i \log \hat{y_i} - (1-y_i) \log(1-\hat{y_i}) L(θ)=−yilogyi^−(1−yi)log(1−yi^)
其中, yi^=hθ(x)\hat{y_{i}} = h_{\theta}(x)yi^=hθ(x) 是逻辑回归预测的结果。假设GBDT第 MMM 步迭代之后当前学习器为 F(x)=∑m=0Mhm(x)F(x) = \sum_{m=0}^M h_m(x)F(x)=∑m=0Mhm(x) ,将 yi^\hat{y_i}yi^ 替换为 F(x)F(x)F(x)带入上式之后,可将损失函数写为:
L(yi,F(xi))=yilog(1+e−F(xi))+(1−yi)[F(xi)+log(1+e−F(xi))]L(y_{i},F(x_{i})) = y_i \log \left( 1+e^{-F(x_i)} \right) + (1-y_i) \left[ F(x_i) + \log \left( 1+e^{-F(x_i)} \right) \right]\\ L(yi,F(xi))=yilog(1+e−F(xi))+(1−yi)[F(xi)+log(1+e−F(xi))]
其中,第 mmm 棵树对应的响应值为(损失函数的负梯度,即伪残差):
rm,i=−∣∂L(yi,F(xi))∂F(xi)∣F(x)=Fm−1(x)=yi−11+e−F(xi)=yi−yi^r_{m,i} = - \left|\frac{\partial L(y_{i},F(x_{i}))}{\partial F(x_{i})} \right|_{F(x)=F_{m-1}(x)} = y_i - \frac{1}{1+e^{-F(x_i)}}= y_i - \hat{y_i} rm,i=−∣∣∣∣∂F(xi)∂L(yi,F(xi))∣∣∣∣F(x)=Fm−1(x)=yi−1+e−F(xi)1=yi−yi^
对于生成的决策树,计算各个叶子节点的最佳残差拟合值为:
cm,j=argminc∑xi∈Rm,jL(yi,Fm−1(xi)+c)c_{m,j}=\underset{c}{ \arg \min } \sum_{x_{i} \in R_{m,j}} L(y_{i},F_{m-1}(x_{i})+c) cm,j=cargminxi∈Rm,j∑L(yi,Fm−1(xi)+c)
由于上式没有闭式解(closed form solution),我们一般使用近似值代替:
cm,j=∑xi∈Rm,jrm,i∑xi∈Rm,j(yi−rm,i)(1−yi+rm,i)c_{m,j} = \frac{\sum_{x_{i} \in R_{m,j}}^{}{r_{m,i}}}{\sum_{x_{i} \in R_{m,j}}^{}{(y_{i}-r_{m,i})(1-y_{i}+ r_{m,i}) }} cm,j=∑xi∈Rm,j(yi−rm,i)(1−yi+rm,i)∑xi∈Rm,jrm,i
补充近似值代替过程:
假设仅有一个样本:
L(yi,F(x))=−(yiln11+e−F(x)+(1−yi)ln(1−11+e−F(x)))L(y_{i},F(x)) = -(y_{i}ln\frac{1}{1+e^{-F(x)}}+(1-y_{i})ln(1-\frac{1}{1+e^{-F(x)}})) L(yi,F(x))=−(yiln1+e−F(x)1+(1−yi)ln(1−1+e−F(x)1))
令 Pi=11+e−F(x)P_{i} = \frac{1}{1+e^{-F(x)}}Pi=1+e−F(x)1 ,则 ∂Pi∂F(x)=Pi(1−Pi)\frac{\partial P_{i}}{\partial F(x)} = P_{i}(1-P_{i})∂F(x)∂Pi=Pi(1−Pi)
求一阶导:
∂L(yi,F(x))∂F(x)=∂L(yi,F(x))∂Pi⋅∂Pi∂F(x)=−(yiPi−1−yi1−Pi)⋅(Pi⋅(1−Pi))=Pi−yi\begin{aligned} \frac{\partial L(y_{i},F(x))}{\partial F(x)} &= \frac{\partial L(y_{i},F(x))}{\partial P_{i}} \cdot \frac{\partial P_{i}}{\partial F(x)} \\ &= -(\frac{y_{i}}{P_{i}} - \frac{1-y_{i}}{1-P_{i}}) \cdot (P_{i}\cdot(1-P_{i})) \\ & = P_{i} - y_{i} \end{aligned} ∂F(x)∂L(yi,F(x))=∂Pi∂L(yi,F(x))⋅∂F(x)∂Pi=−(Piyi−1−Pi1−yi)⋅(Pi⋅(1−Pi))=Pi−yi
求二阶导:
∂2L(yi,F(x))∂F(x)2=(Pi−yi)′=Pi(1−Pi)\begin{aligned} \frac{\partial ^{2}L(y_{i},F(x))}{\partial F(x)^{2}} & = (P_{i} - y_{i})^{'}\\ & = P_{i}(1-P_{i}) \end{aligned} ∂F(x)2∂2L(yi,F(x))=(Pi−yi)′=Pi(1−Pi)
对于 L(yi,F(x)+c)L(y_{i},F(x)+c)L(yi,F(x)+c) 的泰勒二阶展开式:
L(yi,F(x)+c)=L(yi,F(x))+∂L(yi,F(x))∂F(x)⋅c+12∂2L(yi,F(x))∂F(x)2c2L(y_{i},F(x)+c) = L(y_{i},F(x)) + \frac{\partial L(y_{i},F(x))}{\partial F(x)}\cdot c + \frac{1}{2} \frac{\partial ^{2}L(y_{i},F(x))}{\partial F(x)^{2}} c^{2} L(yi,F(x)+c)=L(yi,F(x))+∂F(x)∂L(yi,F(x))⋅c+21∂F(x)2∂2L(yi,F(x))c2
L(yi,F(x)+c)L(y_{i},F(x)+c)L(yi,F(x)+c) 取极值时,上述二阶表达式中的ccc为:
c=−b2a=−∂L(yi,F(x))∂F(x)2(12∂2L(yi,F(x))∂F(x)2)=−∂L(yi,F(x))∂F(x)∂2L(yi,F(x))∂F(x)2)⇒一阶、二阶导代入yi−PiPi(1−Pi)⇒ri=yi−Piri(yi−ri)(1−yi+ri)\begin{aligned} c &= - \frac{b}{2a} = - \frac{\frac{\partial L(y_{i},F(x))}{\partial F(x)}}{2( \frac{1}{2} \frac{\partial ^{2}L(y_{i},F(x))}{\partial F(x)^{2}})} \\ &= - \frac{\frac{\partial L(y_{i},F(x))}{\partial F(x)}}{ \frac{\partial ^{2}L(y_{i},F(x))}{\partial F(x)^{2}})} \overset{一阶、二阶导代入}{\Rightarrow} \frac{y_{i} - P_{i}}{P_{i}(1-P_{i})} \\ & \overset{r_{i} = y_{i} - P_{i}}{\Rightarrow} \frac{r_{i}}{(y_{i}-r_{i})(1-y_{i}+ r_{i})} \end{aligned} c=−2ab=−2(21∂F(x)2∂2L(yi,F(x)))∂F(x)∂L(yi,F(x))=−∂F(x)2∂2L(yi,F(x)))∂F(x)∂L(yi,F(x))⇒一阶、二阶导代入Pi(1−Pi)yi−Pi⇒ri=yi−Pi(yi−ri)(1−yi+ri)ri
GBDT二分类算法完整的算法过程如下:
(1)初始化第一个弱学习器 F0(x)F_{0}(x)F0(x) :
F0(x)=logP(Y=1∣x)1−P(Y=1∣x)F_{0}(x)=log\frac{P(Y=1|x)}{1-P(Y=1|x)} F0(x)=log1−P(Y=1∣x)P(Y=1∣x)
其中, P(Y=1∣x)P(Y=1|x)P(Y=1∣x) 是训练样本中 y=1y=1y=1 的比例,利用先验信息来初始化学习器。
(2)对于建立 MMM 棵分类回归树 m=1,2,...,Mm=1,2,...,Mm=1,2,...,M :
a) 对 i=1,2,...,Ni = 1,2,...,Ni=1,2,...,N ,计算第 mmm 棵树对应的响应值(损失函数的负梯度,即伪残差):
rm,i=−[∂L(yi,F(xi))∂F(x)]F(x)=Fm−1(x)=yi−11+e−F(xi)r_{m,i} = - \left[ \frac{\partial L(y_{i},F(x_{i}))}{\partial F(x)} \right]_{F(x)=F_{m-1}(x) } = y_i - \frac{1}{1+e^{-F(x_i)}} rm,i=−[∂F(x)∂L(yi,F(xi))]F(x)=Fm−1(x)=yi−1+e−F(xi)1
b) 对于i=1,2,…,Ni = 1,2 , \dots , Ni=1,2,…,N ,利用CART回归树拟合数据 (xi,rm,i)\left(x_{i}, r_{m, i}\right)(xi,rm,i) ,得到第 mmm 棵回归树,其对应的叶子节点区域为 Rm,jR_{m, j}Rm,j ,其中 j=1,2,…,Jmj=1,2, \dots, J_{m}j=1,2,…,Jm ,且 JmJ_{m}Jm 为第mmm棵回归树叶子节点的个数。
c) 对于JmJ_{m}Jm 个叶子节点区域 j=1,2,…,Jmj=1,2, \dots, J_{m}j=1,2,…,Jm,计算出最佳拟合值:
cm,j=∑xi∈Rm,jrm,i∑xi∈Rm,j(yi−rm,i)(1−yi+rm,i)c_{m,j} = \frac{\sum_{x_{i} \in R_{m,j}}^{}{r_{m,i}}}{\sum_{x_{i} \in R_{m,j}}^{}{(y_{i}-r_{m,i})(1-y_{i}+ r_{m,i}) }} cm,j=∑xi∈Rm,j(yi−rm,i)(1−yi+rm,i)∑xi∈Rm,jrm,i
d) 更新强学习器 Fm(x)F_{m}(x)Fm(x) :
Fm(x)=Fm−1(x)+∑j=1Jmcm,jI(x∈Rm,j)F_{m}(x)=F_{m-1}(x)+\sum_{j=1}^{J_{m}} c_{m, j} I\left(x \in R_{m,j}\right) Fm(x)=Fm−1(x)+j=1∑Jmcm,jI(x∈Rm,j)
(3)得到最终的强学习器 FM(x)F_{M}(x)FM(x) 的表达式:
FM(x)=F0(x)+∑m=1M∑j=1Jmcm,jI(x∈Rm,j)F_{M}(x)=F_{0}(x)+\sum_{m=1}^{M} \sum_{j=1}^{J_{m}} c_{m,j} I\left(x \in R_{m, j}\right) FM(x)=F0(x)+m=1∑Mj=1∑Jmcm,jI(x∈Rm,j)
从以上过程中可知,除了由损失函数引起的负梯度计算和叶子节点的最佳残差拟合值的计算不同,二元GBDT分类和GBDT回归算法过程基本相似。那么二元GBDT是如何做分类呢?
将逻辑回归的公式进行整理,我们可以得到logp1−p=θTx\log \frac{p}{1-p} = \theta^{T}xlog1−pp=θTx ,其中p=P(Y=1∣x)p=P(Y=1|x)p=P(Y=1∣x),也就是将给定输入 xxx 预测为正样本的概率。逻辑回归用一个线性模型去拟合 Y=1∣xY=1|xY=1∣x 这个事件的对数几率(odds)logp1−p\log \frac{p}{1-p}log1−pp。二元GBDT分类算法和逻辑回归思想一样,用一系列的梯度提升树去拟合这个对数几率,其分类模型可以表达为:
P(Y=1∣x)=11+e−FM(x)P(Y=1|x) = \frac{1}{1+e^{-F_{M}(x)}} P(Y=1∣x)=1+e−FM(x)1
3. GBDT二分类算法实例
(1)数据集介绍
训练集如下表所示,一组数据的特征有年龄和体重,把身高大于1.5米作为分类边界,身高大于1.5米的令标签为1,身高小于等于1.5米的令标签为0,共有4组数据。
测试数据如下表所示,只有一组数据,年龄为25、体重为65,我们用在训练集上训练好的GBDT模型预测该组数据的身高是否大于1.5米?
(2)模型训练阶段
参数设置:
- 学习率:learning_rate = 0.1
- 迭代次数:n_trees = 5
- 树的深度:max_depth = 3
1)初始化弱学习器:
F0(x)=logP(Y=1∣x)1−P(Y=1∣x)=log22=0F_{0}(x)=log\frac{P(Y=1|x)}{1-P(Y=1|x)}=log \frac{2}{2} = 0 F0(x)=log1−P(Y=1∣x)P(Y=1∣x)=log22=0
2)对于建立MMM棵分类回归树m=1,2,…,Mm = 1,2 , \dots , Mm=1,2,…,M:
由于我们设置了迭代次数:n_trees=5
,这就是设置了M=5M=5M=5。
首先计算负梯度,根据上文损失函数为对数损失时,负梯度(即伪残差、近似残差)为:
rm,i=−[∂L(yi,F(xi))∂F(x)]F(x)=Fm−1(x)=yi−11+e−F(xi)r_{m,i} = - \left[ \frac{\partial L(y_{i},F(x_{i}))}{\partial F(x)} \right]_{F(x)=F_{m-1}(x) } = y_i - \frac{1}{1+e^{-F(x_i)}} rm,i=−[∂F(x)∂L(yi,F(xi))]F(x)=Fm−1(x)=yi−1+e−F(xi)1
我们知道梯度提升类算法,其关键是利用损失函数的负梯度的值作为回归问题提升树算法中的残差的近似值,拟合一个回归树。这里,为了称呼方便,我们把负梯度叫做残差。
现将残差的计算结果列表如下:
此时将残差作为样本的标签来训练弱学习器 F1(x)F_{1}(x)F1(x) ,即下表数据:
接着寻找回归树的最佳划分节点,遍历每个特征的每个可能取值。从年龄特征值为555开始,到体重特征为707070结束,分别计算分裂后两组数据的平方损失(Square Error),SElSE_{l}SEl 为左节点的平方损失, SErSE_{r}SEr 为右节点的平方损失,找到使平方损失和 SEsum=SEl+SErSE_{sum} = SE_{l}+SE_{r}SEsum=SEl+SEr 最小的那个划分节点,即为最佳划分节点。
例如:以年龄777为划分节点,将小于777的样本划分为到左节点,大于等于777的样本划分为右节点。左节点包括 x0x_{0}x0,右节点包括样本 x1,x2,x3x_{1},x_{2},x_{3}x1,x2,x3, SEl=0SE_{l}=0SEl=0, SEr=0.667SE_{r}=0.667SEr=0.667,SEsum=0.667SE_{sum}=0.667SEsum=0.667 ,所有可能的划分情况如下表所示:
以上划分点的总平方损失最小为0.0000.0000.000,有两个划分点:年龄212121和体重606060,所以随机选一个作为划分点,这里我们选年龄212121。现在我们的第一棵树长这个样子:
我们设置的参数中树的深度max_depth=3
,现在树的深度只有222,需要再进行一次划分,这次划分要对左右两个节点分别进行划分,但是我们在生成树的时候,设置了三个树继续生长的条件:
- 深度没有到达最大。树的深度设置为3,意思是需要生长成3层;
- 点样本数 >= min_samples_split;
- 此节点上的样本的标签值不一样。如果值一样说明已经划分得很好了,不需要再分;(本程序满足这个条件,因此树只有2层)
最终我们的第一棵回归树长下面这个样子:
此时我们的树满足了设置,还需要做一件事情,给这棵树的每个叶子节点分别赋一个参数ccc,来拟合残差。
c1,j=∑xi∈R1,jr1,i∑xi∈R1,j(yi−r1,i)(1−yi+r1,i)c_{1,j} = \frac{\sum_{x_{i} \in R_{1,j}}^{}{r_{1,i}}}{\sum_{x_{i} \in R_{1,j}}^{}{(y_{i}-r_{1,i})(1-y_{i}+ r_{1,i}) }} c1,j=∑xi∈R1,j(yi−r1,i)(1−yi+r1,i)∑xi∈R1,jr1,i
根据上述划分结果,为了方便表示,规定从左到右为第1,21,21,2个叶子结点,其计算值过程如下:
(x0,x1∈R1,1),c1,1=−2.0(x_{0},x_{1}\in R_{1,1}), \qquad c_{1,1}=-2.0 (x0,x1∈R1,1),c1,1=−2.0
(x2,x3∈R1,2),c1,2=2.0(x_{2}, x_{3}\in R_{1,2}), \qquad c_{1,2}=2.0 (x2,x3∈R1,2),c1,2=2.0
此时的第一棵树长下面这个样子:
接着更新强学习器,需要用到学习率:learning_rate=0.1
,用lr
表示。更新公式为:
F1(x)=F0(x)+lr∗∑j=12c1,jI(x∈R1,j)F_{1}(x)=F_{0}(x)+lr * \sum_{j=1}^{2} c_{1, j} I\left(x \in R_{1,j}\right) F1(x)=F0(x)+lr∗j=1∑2c1,jI(x∈R1,j)
为什么要用学习率呢?这是Shrinkage的思想,如果每次都全部加上拟合值 ccc ,即学习率为111,很容易一步学到位导致GBDT过拟合。
重复此步骤,直到 m>5m>5m>5 结束,最后生成555棵树。
下面将展示每棵树最终的结构,这些图都是我GitHub上的代码生成的,感兴趣的同学可以去运行一下代码。https://github.com/Microstrong0305/WeChat-zhihu-csdnblog-code/tree/master/Ensemble%20Learning/GBDT_GradientBoostingBinaryClassifier
第一棵树:
第二棵树:
第三棵树:
第四棵树:
第五棵树:
3)得到最后的强学习器:
F5(x)=F0(x)+∑m=15∑j=12cm,jI(x∈Rm,j)F_{5}(x)=F_{0}(x)+\sum_{m=1}^{5} \sum_{j=1}^{2} c_{m,j} I\left(x \in R_{m, j}\right) F5(x)=F0(x)+m=1∑5j=1∑2cm,jI(x∈Rm,j)
(3)模型预测阶段
- F0(x)=0F_{0}(x) = 0F0(x)=0
- 在 F1(x)F_{1}(x)F1(x) 中,测试样本的年龄为252525,大于划分节点212121岁,所以被预测为2.00002.00002.0000。
- 在 F2(x)F_{2}(x)F2(x) 中,测试样本的年龄为252525,大于划分节点212121岁,所以被预测为1.81871.81871.8187。
- 在 F3(x)F_{3}(x)F3(x) 中,测试样本的年龄为252525,大于划分节点212121岁,所以被预测为1.68261.68261.6826。
- 在 F4(x)F_{4}(x)F4(x) 中,测试样本的年龄为252525,大于划分节点212121岁,所以被预测为1.57691.57691.5769。
- 在 F5(x)F_{5}(x)F5(x) 中,测试样本的年龄为252525,大于划分节点212121岁,所以被预测为1.49271.49271.4927。
最终预测结果为:
F(x)=0.0000+0.1∗(2.0000+1.8187+1.6826+1.5769+1.4927)=0.8571F(x) = 0.0000 + 0.1 * (2.0000+1.8187+1.6826+1.5769+1.4927)=0.8571 F(x)=0.0000+0.1∗(2.0000+1.8187+1.6826+1.5769+1.4927)=0.8571
P(Y=1∣x)=11+e−F(x)=11+e−0.8571=0.7021P(Y=1|x) = \frac{1}{1+e^{-F(x)}}= \frac{1}{1+e^{-0.8571}} = 0.7021 P(Y=1∣x)=1+e−F(x)1=1+e−0.85711=0.7021
4. 手撕GBDT二分类算法
本篇文章所有数据集和代码均在我的GitHub中,地址:https://github.com/Microstrong0305/WeChat-zhihu-csdnblog-code/tree/master/Ensemble%20Learning
4.1 用Python3实现GBDT二分类算法
需要的Python库:
pandas、PIL、pydotplus、matplotlib
其中pydotplus
库会自动调用Graphviz
,所以需要去Graphviz
官网下载graphviz-2.38.msi
安装,再将安装目录下的bin
添加到系统环境变量,最后重启计算机。
由于用Python3实现GBDT二分类算法代码量比较多,我这里就不列出详细代码了,感兴趣的同学可以去我的GitHub中看一下,地址:https://github.com/Microstrong0305/WeChat-zhihu-csdnblog-code/tree/master/Ensemble%20Learning/GBDT_GradientBoostingBinaryClassifier
4.2 用sklearn实现GBDT二分类算法
import numpy as np
from sklearn.ensemble import GradientBoostingClassifier'''
调参:
loss:损失函数。有deviance和exponential两种。deviance是采用对数似然,exponential是指数损失,后者相当于AdaBoost。
n_estimators:最大弱学习器个数,默认是100,调参时要注意过拟合或欠拟合,一般和learning_rate一起考虑。
learning_rate:步长,即每个弱学习器的权重缩减系数,默认为0.1,取值范围0-1,当取值为1时,相当于权重不缩减。较小的learning_rate相当于更多的迭代次数。
subsample:子采样,默认为1,取值范围(0,1],当取值为1时,相当于没有采样。小于1时,即进行采样,按比例采样得到的样本去构建弱学习器。这样做可以防止过拟合,但是值不能太低,会造成高方差。
init:初始化弱学习器。不使用的话就是第一轮迭代构建的弱学习器.如果没有先验的话就可以不用管由于GBDT使用CART回归决策树。以下参数用于调优弱学习器,主要都是为了防止过拟合
max_feature:树分裂时考虑的最大特征数,默认为None,也就是考虑所有特征。可以取值有:log2,auto,sqrt
max_depth:CART最大深度,默认为None
min_sample_split:划分节点时需要保留的样本数。当某节点的样本数小于某个值时,就当做叶子节点,不允许再分裂。默认是2
min_sample_leaf:叶子节点最少样本数。如果某个叶子节点数量少于某个值,会同它的兄弟节点一起被剪枝。默认是1
min_weight_fraction_leaf:叶子节点最小的样本权重和。如果小于某个值,会同它的兄弟节点一起被剪枝。一般用于权重变化的样本。默认是0
min_leaf_nodes:最大叶子节点数
'''gbdt = GradientBoostingClassifier(loss='deviance', learning_rate=0.1, n_estimators=5, subsample=1, min_samples_split=2, min_samples_leaf=1, max_depth=3, init=None, random_state=None, max_features=None, verbose=0, max_leaf_nodes=None, warm_start=False)train_feat = np.array([[1, 5, 20],[2, 7, 30],[3, 21, 70],[4, 30, 60],])
train_label = np.array([[0], [0], [1], [1]]).ravel()test_feat = np.array([[5, 25, 65]])
test_label = np.array([[1]])
print(train_feat.shape, train_label.shape, test_feat.shape, test_label.shape)gbdt.fit(train_feat, train_label)
pred = gbdt.predict(test_feat)total_err = 0
for i in range(pred.shape[0]):print(pred[i], test_label[i])err = (pred[i] - test_label[i]) / test_label[i]total_err += err * err
print(total_err / pred.shape[0])
用sklearn中的GBDT库实现GBDT二分类算法的难点在于如何更好的调节下列参数:
用sklearn实现GBDT二分类算法的GitHub地址:https://github.com/Microstrong0305/WeChat-zhihu-csdnblog-code/tree/master/Ensemble%20Learning/GBDT_Classification_sklearn
5. GBDT分类任务常见的损失函数
对于GBDT分类算法,其损失函数一般有对数损失函数和指数损失函数两种:
(1)如果是指数损失函数,则损失函数表达式为:
L(y,F(x))=exp(−yF(x))L(y,F(x)) = exp(-yF(x)) L(y,F(x))=exp(−yF(x))
其负梯度计算和叶子节点的最佳负梯度拟合可以参看Adaboost算法过程。
(2)如果是对数损失函数,分为二元分类和多元分类两种,本文主要介绍了GBDT二元分类的损失函数。
6. 总结
在本文中,我们首先简单介绍了如何把GBDT回归算法变成分类算法的思路;然后从逻辑回归的对数损失函数推导出GBDT的二分类算法原理;其次不仅用Python3实现GBDT二分类算法,还用sklearn实现GBDT二分类算法;最后介绍了GBDT分类任务中常见的损失函数。GBDT可以完美的解决二分类任务,那么它对多分类任务是否有效呢?如果有效,GBDT是如何做多分类呢?这些问题都需要我们不停的探索和挖掘GBDT的深层原理。让我们期待一下GBDT在多分类任务中的表现吧!
7. Reference
【1】Friedman J H. Greedy function approximation: a gradient boosting machine[J]. Annals of statistics, 2001: 1189-1232.
【2】GBDT详细讲解&常考面试题要点,地址:https://mp.weixin.qq.com/s/M2PwsrAnI1S9SxSB1guHdg
【3】机器学习算法GBDT的面试要点总结-上篇,地址:https://www.cnblogs.com/ModifyRong/p/7744987.html
【4】Gradient Boosting Decision Tree,地址:http://gitlinux.net/2019-06-11-gbdt-gradient-boosting-decision-tree/
【5】GBDT算法用于分类问题 - hunter7z的文章 - 知乎,地址:https://zhuanlan.zhihu.com/p/46445201
【6】当我们在谈论GBDT:Gradient Boosting 用于分类与回归 - 余文毅的文章 - 知乎
,地址: https://zhuanlan.zhihu.com/p/25257856
【7】 GBDT原理与Sklearn源码分析-分类篇,地址:https://blog.csdn.net/qq_22238533/article/details/79192579
【8】 《推荐系统算法实践》,黄美灵著。
【9】《百面机器学习》,诸葛越主编、葫芦娃著。
【10】GBDT模型,地址:https://www.jianshu.com/p/0bc32c8e4ca8
【11】代码实战之GBDT,地址:https://louisscorpio.github.io/2018/01/19/%E4%BB%A3%E7%A0%81%E5%AE%9E%E6%88%98%E4%B9%8BGBDT/
【推荐阅读】
【1】Regression Tree 回归树
【2】深入理解提升树(Boosting Tree)算法
【3】梯度提升(Gradient Boosting)算法
【4】深入理解GBDT回归算法
深入理解GBDT二分类算法相关推荐
- python 二分类的实例_深入理解GBDT二分类算法
我的个人微信公众号:Microstrong 微信公众号ID:MicrostrongAI 微信公众号介绍:Microstrong(小强)同学主要研究机器学习.深度学习.计算机视觉.智能对话系统相关内容, ...
- 深入解析GBDT二分类算法(附代码实现)
目录: GBDT分类算法简介 GBDT二分类算法 2.1 逻辑回归的对数损失函数 2.2 GBDT二分类原理 GBDT二分类算法实例 手撕GBDT二分类算法 4.1 用Python3实现GBDT二分类 ...
- 深入理解GBDT多分类算法
目录: GBDT多分类算法 1.1 Softmax回归的对数损失函数 1.2 GBDT多分类原理 GBDT多分类算法实例 手撕GBDT多分类算法 3.1 用Python3实现GBDT多分类算法 3.2 ...
- python 多分类算法_深入理解GBDT多分类算法
我的个人微信公众号:Microstrong 微信公众号ID:MicrostrongAI 微信公众号介绍:Microstrong(小强)同学主要研究机器学习.深度学习.计算机视觉.智能对话系统相关内容, ...
- 「二分类算法」提供银行精准营销解决方案(各个模型汇总分析)baseline
完整代码见 Github:「二分类算法」提供银行精准营销解决方案 赛事详情 1.比赛页面:「二分类算法」提供银行精准营销解决方案 2.数据与测评算法 训练集下载链接 测试集下载链接 字段说明 测评算法 ...
- python复现感知机的二分类算法
本博客复现感知机的二分类算法 import pandas as pd import numpy as npclass Perceptron:def __init__(self, w=None, b=0 ...
- 感知机(Perceptro)二分类算法原理学习小结记录
感知机是非常古老的一个二分类算法,原理非常简单,直观地来说,对于平面上的两类数据点,找到一条直线来将其划分开来,对于空间中的点,用一个平面来将两部分完全划分开,这样的数据都是线性可分的,也就是在平面或 ...
- 「二分类算法」提供银行精准营销解决方案 代码存档
「二分类算法」提供银行精准营销解决方案 代码存档 参考文章: (1)「二分类算法」提供银行精准营销解决方案 代码存档 (2)https://www.cnblogs.com/starcrm/p/1180 ...
- 基于deap脑电数据集的脑电情绪识别二分类算法(附代码)
想尝试一下脑电情绪识别的各个二分类算法. 代码主要分为三部分:快速傅里叶变换处理(fft).数据预处理.以及各个模型处理. 采用的模型包括:决策树.SVM.KNN三个模型(模型采用的比较简单,可以直接 ...
最新文章
- UVa11855 Buzzwords(hash)
- Django通过URL传递参数的4种方法
- 计算机技能大赛文字录入试题,《文字录入与排版》试题.doc
- python实现地牢迷宫生成
- 阿里云发布第七代云服务器ECS,整机算力提升160%
- JavaWeb学习笔记——Tomcat配置
- juniper CLI 基本操作
- vue引用自定义.css文件 - 语法篇
- oracle中execute函数,Oracle 动态SQL execute immediate写法 zt
- L2-016 愿天下有情人都是失散多年的兄妹(DFS)
- linux acpidtd 进程,clover引导卡在该图上了,求资助
- 【反编译】apk反编译之修改贴吧发贴回帖时显示的地理位置
- 如何1分钟制作自己的软件
- 常见文件扩展名和它们的说明(转)
- python数据库管理系统ess_【简说Python WEB】数据库
- 什么情况下钽电容可以代替铝电解电容?
- 搞懂这 9 个步骤,DNS 访问原理就明明白白了
- URLDecoder.decode与URLEncoder.encode
- 官宣!微软跨平台 UI 框架 .NET MAUI 6 正式发布!
- java 计算两个日期相差月数_Java简单计算两个日期月数差的方法