大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全系列IDE稳定放心使用
1. 背景
决策树是一种基本的分类与回归方法。决策树模型具有分类速度快,模型容易可视化的解释,但是同时是也有容易发生过拟合,虽然有剪枝,但也是差强人意。
提升方法(boosting)在分类问题中,它通过改变训练样本的权重(增加分错样本的权重,减小分队样本的的权重),学习多个分类器,并将这些分类器线性组合,提高分类器性能。boosting数学表示为:
f ( x ) = w 0 + ∑ m = 1 M w m ϕ m ( x ) f(x) = w_0 + \sum\limits_{m=1}^M w_m \phi_m(x) f(x)=w0+m=1∑Mwmϕm(x)
其中w是权重, ϕ \phi ϕ是弱分类器的集合,可以看出最终就是基函数的线性组合。
于是决策树与boosting结合产生许多算法,主要有提升树、GBDT等。本文主要是GBDT学习笔记。
1.1 Gradient Boosting
Gradient Boosting是一种Boosting的方法,它主要的思想是,每一次建立模型是在之前建立模型损失函数的梯度下降方向。损失函数是评价模型性能(一般为拟合程度+正则项),认为损失函数越小,性能越好。而让损失函数持续下降,就能使得模型不断改性提升性能,其最好的方法就是使损失函数沿着梯度方向下降(讲道理梯度方向上下降最快)。
Gradient Boost是一个框架,里面可以套入很多不同的算法。
1.2 提升树-boosting tree
以决策树为基函数的提升方法称为提升树,其决策树可以是分类树OR回归树。提升树模型可以表示为决策树的加法模型。
f M ( x ) = ∑ m = 1 M T ( x ; Θ m ) f_M(x) = \sum\limits_{m=1}^M T(x;\Theta_m) fM(x)=m=1∑MT(x;Θm)
其中, T ( x ; Θ m ) 表 示 决 策 树 , T(x;\Theta_m)表示决策树, T(x;Θm)表示决策树, Θ m \Theta_m Θm表示树的参数,M为树的个数。
回归问题提升树算法
输入:训练数据集$T={(x_1,y_1),(x_2,y_2),···,(x_N,y_N)}, x_i \in \chi = R^n, y_i \in \gamma, \ i=1,2,···,N ; ; ;\gamma$为输出空间。
输出:提升树 f M ( x ) f_M(x) fM(x)
-
初始化 f 0 ( x ) = 0 f_0(x)=0 f0(x)=0
-
对于 m = 1 , 2 , . . . M m=1,2,…M m=1,2,...M:
-
计算残差(后一棵树拟合前一颗树残差):
r m i = y i − f m − 1 ( x i ) r_{mi} = y_i – f_{m-1}(x_i) rmi=yi−fm−1(xi)
-
拟合残差学习一个回归树,得到 T ( x ; Θ m ) T(x;\Theta_m) T(x;Θm)
-
更新 f m ( x ) = f m − 1 ( x ) + T ( x ; Θ m ) f_m(x) = f_{m-1}(x) + T(x;\Theta_m) fm(x)=fm−1(x)+T(x;Θm)
-
-
M次迭代之后得到提升树:
f M ( x ) = ∑ m = 1 M T ( x ; Θ m ) f_M(x) = \sum\limits_{m=1}^M T(x;\Theta_m) fM(x)=m=1∑MT(x;Θm)
2 Gradient Boosting Decision Tree
提升树的学习优化过程中,损失函数平方损失和指数损失时候,每一步优化相对简单,但对于一般损失函数优化的问题,Freidman提出了Gradient Boosting算法,其利用了损失函数的负梯度在当前模型的值
− [ ∂ L ( y , f ( x i ) ) ∂ f ( x i ) ] f ( x ) = f m − 1 ( x ) -[\frac{\partial L(y,f(x_i))}{\partial f(x_i)}]_{f(x)=f_{m-1}(x)} −[∂f(xi)∂L(y,f(xi))]f(x)=fm−1(x)
作为回归问题提升树算法的残差近似值,去拟合一个回归树。
2.1 函数空间的数值优化
优化目标是使得损失函数最小,(N是样本集合大小):
F ∗ ( x ) = arg min ρ ∑ i = 1 N L ( y i , ρ ) F^*(x)=\underset{\rho}{\arg\min}\sum^N_{i=1}\mathcal{L}(y_i, \rho) F∗(x)=ρargmini=1∑NL(yi,ρ)
GBDT是一个加法模型: f m ( x ) f_m(x) fm(x)是每一次迭代学习的到树模型。
F ^ ( x ) = F M ( x ) = ∑ m = 1 M f m ( x ) \hat{F}(x) = F_M(x) = \sum\limits_{m=1}^M f_m(x) F^(x)=FM(x)=m=1∑Mfm(x)
对于其每一步迭代:
f m ( x ) = − ρ m g m ( x ) f_m(x) = -\rho_m g_m(x) fm(x)=−ρmgm(x)
其中
g m ( x ) = [ ∂ ϕ ( F ( x ) ) ∂ F ( x ) ] F ( x ) = F m − 1 ( x ) ϕ ( F ( x ) ) = E y [ L ( y , F ( x ) ) ∣ x ] , F m − 1 ( x ) = ∑ i = 0 m − 1 f i ( x ) g_m(x) = [\frac{\partial \phi(F(x))}{\partial F(x)}]_{F(x) = F_{m-1}(x)} \\ \phi(F(x)) = E_y[L(y,F(x))|x], F_{m-1}(x) = \sum_{i=0}^{m-1} f_i(x) gm(x)=[∂F(x)∂ϕ(F(x))]F(x)=Fm−1(x)ϕ(F(x))=Ey[L(y,F(x))∣x],Fm−1(x)=i=0∑m−1fi(x)
其实 L ( y , F ( x ) ) L(y, F(x)) L(y,F(x))就是损失函数, ϕ ( F ( x ) ) \phi(F(x)) ϕ(F(x))是当前x下的损失期望, g m ( x ) g_m(x) gm(x)是当前x下的函数梯度。最终 f m ( x ) f_m(x) fm(x)学习的是损失函数在函数空间上的负梯度。
对于权重 ρ m \rho_m ρm通过线性搜索求解(这也是后面算法改进的点):
ρ m = arg min ρ E y , x L ( y , F m − 1 ( x ) − ρ ∗ g m ( x ) ) \rho_m = \arg \min_{\rho} E_{y,x} L(y, F_{m-1}(x) – \rho *g_m(x)) ρm=argρminEy,xL(y,Fm−1(x)−ρ∗gm(x))
理解:每一次迭代可以看做是采用梯度下降法对最优分类器 F ∗ ( x ) F^*(x) F∗(x)的逐渐比较,每一次学习的模型 f m ( x ) f_m(x) fm(x)是梯度,进过M步迭代之后,最后加出来的模型就是最优分类器的一个逼近模型,所以 f m ( x i ) f_m(x_i) fm(xi)使用单步修正方向 − g m ( x i ) -g_m(x_i) −gm(xi):
− g m ( x i ) = g m ( x ) = [ ∂ L ( F ( x ) ) ∂ F ( x ) ] F ( x ) = F m − 1 ( x ) -g_m(x_i) = g_m(x) = [\frac{\partial L(F(x))}{\partial F(x)}]_{F(x) = F_{m-1}(x)} −gm(xi)=gm(x)=[∂F(x)∂L(F(x))]F(x)=Fm−1(x)
这里的梯度变量是函数,是在函数空间上求解(这也是后面XGBoost改进的点),注意以往算法梯度下降是在N维的参数空间的负梯度方向,变量是参数。这里的变量是函数,更新函数通过当前函数的负梯度方向来修正模型,是它更优,最后累加的模型近似最优函数。
2.2 算法
输入:训练数据集 T = { ( x 1 , y 1 ) , ( x 2 , y 2 ) , ⋅ ⋅ ⋅ , ( x N , y N ) } T=\{(x_1,y_1),(x_2,y_2),···,(x_N,y_N)\} T={
(x1,y1),(x2,y2),⋅⋅⋅,(xN,yN)},$ x_i \in \chi = R^n , , , y_i \in \gamma={-1,+1}, \ i=1,2,···,N $;
输出:回归树 f M ( x ) f_M(x) fM(x)
-
初始化
f 0 ( x ) = a r g min c ∑ i = 1 N L ( y i , c ) f_0(x) = arg \min\limits_c \sum\limits_{i=1}^N L(y_i,c) f0(x)=argcmini=1∑NL(yi,c)
-
对m=1,2,…M
-
对i=1,2,…,N,计算
r m i = − [ ∂ L ( y , f ( x i ) ) ∂ f ( x i ) ] f ( x ) = f m − 1 ( x ) r_{mi}= -[\frac{\partial L(y,f(x_i))}{\partial f(x_i)}]_{f(x)=f_{m-1}(x)} rmi=−[∂f(xi)∂L(y,f(xi))]f(x)=fm−1(x)
-
对 r m i r_{mi} rmi拟合一颗回归树,得到第m棵树的叶结点区域 R m j , j = 1 , 2 , . . . J R_{mj}, \ j=1,2,…J Rmj, j=1,2,...J,即一棵由J个叶子节点组成的树。
-
对 j = 1 , 2 , . . . J j=1,2,…J j=1,2,...J,计算
c m j = a r g min c ∑ x i ∈ R m j L ( y i , f m − 1 ( x i ) + c ) c_{mj}=arg \min\limits_c \sum\limits_{x_i \in R_{mj} } L(y_i, f_{m-1}(x_i) + c) cmj=argcminxi∈Rmj∑L(yi,fm−1(xi)+c)
2.2,2.3这一步相当于回归树递归在遍历所有切分变量j和切分点s找到最优j,s,然后在每个节点区域求最优的c。参考回归树生成算法
-
更新 f m ( x ) = f m − 1 ( x ) + ∑ j = 1 J c m j I ( x ∈ R m j ) f_m(x)=f_{m-1}(x) + \sum\limits_{j=1}^J c_{mj} I(x \in R_{mj}) fm(x)=fm−1(x)+j=1∑JcmjI(x∈Rmj)
-
-
得到回归树
f ^ ( x ) = f M ( x ) = ∑ m = 1 M f m ( x ) = ∑ m = 1 M ∑ j = 1 J c m j I ( x ∈ R m j ) \hat{f}(x) = f_M(x) = \sum\limits_{m=1}^M f_m(x) = \sum\limits_{m=1}^M \sum\limits_{j=1}^J c_{mj}I(x \in R_{mj}) f^(x)=fM(x)=m=1∑Mfm(x)=m=1∑Mj=1∑JcmjI(x∈Rmj)
算法1步获得使得损失函数最小的常数估计值,是一个只有根节点的树。在2.1步计算损失函数的负梯度在当前模型的值,将它作为残差估计。在2.2步估计回归树的叶结点区域,来拟合残差的近似值。在2.3步利用线性搜索估计回归树叶结点区域的值,使损失函数最小化。2.4更新回归树。第3步获得输出的最终模型。
Shrinkage
Shrinkage的思想认为,每次走一小步逐渐逼近结果的效果,要比每次迈一大步很快逼近结果的方式更容易避免过拟合。即它不完全信任每一个棵残差树,它认为每棵树只学到了真理的一小部分,累加的时候只累加一小部分,通过多学几棵树弥补不足。
数学方程对比:
- 之前: f m ( x ) = f m − 1 ( x ) + ∑ j = 1 J c m j I ( x ∈ R m j ) f_m(x)=f_{m-1}(x) + \sum\limits_{j=1}^J c_{mj} I(x \in R_{mj}) fm(x)=fm−1(x)+j=1∑JcmjI(x∈Rmj)
- Shrinkage: f m ( x ) = f m − 1 ( x ) + s t e p ∗ ∑ j = 1 J c m j I ( x ∈ R m j ) f_m(x)=f_{m-1}(x) +step* \sum\limits_{j=1}^J c_{mj} I(x \in R_{mj}) fm(x)=fm−1(x)+step∗j=1∑JcmjI(x∈Rmj)
Shrinkage仍然以残差作为学习目标,但对于残差学习的结果,只累加一小部分,step一般取值0.001-0.01(非gradient的step),使得各个树的残差是渐变而不是陡变的,即将大步切成了小步。Shrinkage能减少过拟合发生也是经验证明的,目前还没有看到从理论的证明。
总结
原始的boosting算法开始时,为每一个样本赋上一个权重值。在每一步训练中得到的模型,会使得数据点的估计有对有错,在每一步结束后,增加分错的点的权重,减少分对的点的权重,这样使得某些点如果老是被分错,那么就会被“严重关注”,也就被赋上一个很高的权重。然后等进行了N次迭代(由用户指定),将会得到N个简单的分类器(basic learner),然后我们将它们组合起来(比如说可以对它们进行加权、或者让它们进行投票等),得到一个最终的模型。
那么GBDT算法中并未有权重的改变,哪里有boosting思想 ?
Gradient Boosting与Boosting区别在于,每一计算的是为了减少上一次的残差,下一个模型主要在残差减少的梯度方上建立模型,使得残差往梯度方向上减少。
虽然不同,但是GBDT算法会更关注那些梯度比较大的样本,和Boosting思想类似。
附录
CSDN原文:http://blog.csdn.net/shine19930820/article/details/65633436
公众号:百川NLP
参考资料
- 《统计学习方法》
- 《The Elements of Statistical Learning 》
- 《Machine Learning A Probabilistic Perspective》
- http://www.lai18.com/content/1406280.html
相似算法:
- LightGBM原理: A Highly Efficient Gradient Boosting Decision Tree
- XGBoost原理- A Scalable Tree Boosting System
- 数据挖掘十大算法 』笔记二:SVM-支持向量机
- 《Top 10 algorithms in data mining》
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/185669.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...