04-02 AdaBoost算法

目录

• AdaBoost算法
• 一、AdaBoost算法学习目标
• 二、AdaBoost算法详解
  • 2.1 Boosting算法回顾
  • 2.2 AdaBoost算法
  • 2.3 AdaBoost算法目标函数优化
• 三、AdaBoost算法流程
  • 3.1 输入
  • 3.2 输出
  • 3.3 强分类器流程
  • 3.4 强回归器流程
• 四、AdaBoost算法优缺点
  • 4.1 优点
  • 4.2 缺点
• 五、小结

更新、更全的《机器学习》的更新网站，更有python、go、数据结构与算法、爬虫、人工智能教学等着你：https://www.cnblogs.com/nickchen121/p/11686958.html

AdaBoost算法

集成学习中弱学习器之间有强依赖关系的，称之为Boosting系列算法，而AdaBoost则是Boosting系列算法中最著名的算法之一。

AdaBoost算法强大之处在于既可以解决分类问题，又可以解决回归问题。

一、AdaBoost算法学习目标

1. AdaBoost算法目标函数优化
2. 强分类器和强回归器流程
3. AdaBoost算法优缺点

二、AdaBoost算法详解

2.1 Boosting算法回顾

Boosting算法的流程是：首先训练处一个弱学习器，根据弱学习器的误差率更新训练样本的权重，然后基于调整权重后的训练集训练第二个弱学习器，直到弱学习器达到事先指定的数目T，停止算法。

对于Boosting算法的流程，可以看到如果我们解决以下4个问题，既可以得到完整的Boosting算法

1. 弱学习器的误差率
2. 训练样本的权重w$w$更新方法
3. 更新样本权重的方法
4. 结合策略

2.2 AdaBoost算法

上面讲到了Boosting算法需要解决的4个问题，因为AdaBoost算法隶属于Boosting算法，那么AdaBoost算法也需要解决这4个问题，其实也可以说成只要是Boosting系列的算法，都需要解决这4个问题。

2.3 AdaBoost算法目标函数优化

AdaBoost算法可以理解成模型是加法模型、目标函数是指数函数、学习算法是前向分步算法时的学习方法。其中加法模型可以理解成强学习器是由之前所有的弱学习器加权平均得到的；前向分步算法则可以理解成弱学习器训练数据的权重通过前一个弱学习器更新。

AdaBoost算法的模型是加法模型，即强学习器的模型为

f(x)=∑k=1KαkGk(x)$$f(x)=\sum _{k=1}^K{\alpha }_k{G}_k(x)$$

其中K$K$是K$K$个弱学习器。

AdaBoost算法的事前向分步算法，即经过k−1$k-1$次迭代后，第k−1$k-1$轮后强学习器为

fk−1(x)=α1G1(x)+α2G2(x)+⋯+αk−1Gk−1(x)=fk−2(x)+αk−1Gk−1(x)(1)(2)$$\begin{array}{}\text{(1)}& f_{k-1}(x)& ={\alpha }_1{G}_1(x)+{\alpha }_2{G}_2(x)+\cdots +{\alpha }_{k-1}{G}_{k-1}(x)\text{(2)}& & =f_{k-2}(x)+{\alpha }_{k-1}{G}_{k-1}(x)\end{array}$$

经过k$k$次迭代后，第k$k$轮后强学习器为

fk(x)=∑i=1kαiGi(x)=fk−1(x)+αkGk(x)$$f_k(x)=\sum _{i=1}^k{\alpha }_i{G}_i(x)=f_{k-1}(x)+{\alpha }_k{G}_k(x)$$

得到第k$k$轮强学习器后，我们知道AdaBoost的目标函数是指数函数，因此我们的目标是使前向分步算法得到的αk${\alpha }_k$和Gk(x)$G_k(x)$使fk(x)$f_k(x)$在训练数据集上的指数损失最小，即AdaBoost的目标函数为

(αk,Gk(x))=argmin����������α,G∑i=1me−yifk(x)=argmin����������α,G∑i=1me[−yi(fk−1(xi)+αG(xi))]=argmin����������α,G∑i=1me[−yi(fk−1(xi))]e[−yi(αG(xi))](3)(4)(5)$$\begin{array}{}\text{(3)}& ({\alpha }_k,G_k(x))& =\underset{\alpha ,G}{\underbrace{\arg min}}\sum _{i=1}^m{e}^{-y_i{f}_k(x)}\text{(4)}& & =\underset{\alpha ,G}{\underbrace{\arg min}}\sum _{i=1}^m{e}^{[-y_i(f_{k-1}(x_i)+\alpha G(x_i))]}\text{(5)}& & =\underset{\alpha ,G}{\underbrace{\arg min}}\sum _{i=1}^m{e}^{[-y_i(f_{k-1}(x_i))]}{e}^{[-y_i(\alpha G(x_i))]}\end{array}$$

由于e[−yi(fk−1(xi))]$e^{[-y_i(f_{k-1}(x_i))]}$的值不依赖α,G$\alpha ,G$，因此他与最小化无关，它仅仅依赖于随着每一轮迭代而变化的fk−1(x)$f_{k-1}(x)$，因此可以把e[−yi(fk−1(xi))]$e^{[-y_i(f_{k-1}(x_i))]}$看做w¯¯¯¯ki${\overline{w}}_{ki}$，即目标函数变为

(αk,Gk(x))=argmin����������α,G∑i=1mw¯¯¯¯kie[−yi(αG(xi))]$$({\alpha }_k,G_k(x))=\underset{\alpha ,G}{\underbrace{\arg min}}\sum _{i=1}^m{\overline{w}}_{ki}{e}^{[-y_i(\alpha G(x_i))]}$$

现在的目标就是最优化AdaBoost的目标函数得到能使目标函数最小化的α∗k${\alpha }_k^{\ast }$和G∗k(x)$G_k^{\ast }(x)$。

首先，对于任意的α>0$\alpha >0$，G∗k(x)$G_k^{\ast }(x)$表示第k$k$轮能够使得加训练数据分类误差率最小的基本分类器，分类误差率为

ek=∑mi=1w¯¯¯¯kiI(yi≠Gk(xi))∑mi=1w¯¯¯¯ki=∑i=1mw¯¯¯¯kiI(yi≠Gk(xi))=∑yi≠Gk(xi)w¯¯¯¯ki$$e_k=\frac{\sum _{i=1}^m{\overline{w}}_{ki}I(y_i\ne G_k(x_i))}{\sum _{i=1}^m{\overline{w}}_{ki}}=\sum _{i=1}^m{\overline{w}}_{ki}I(y_i\ne G_k(x_i))=\sum _{y_i\ne G_k(x_i)}{\overline{w}}_{ki}$$

G∗k(x)$G_k^{\ast }(x)$为

G∗k(x)=argmin����������G∑i=1mw¯¯¯¯kiI(yi≠G(xi))$$G_k^{\ast }(x)=\underset{G}{\underbrace{argmin}}\sum _{i=1}^m{\overline{w}}_{ki}I(y_i\ne G(x_i))$$

G∗k(x)$G_k^{\ast }(x)$即为学习器的Gk(x)$G_k(x)$，把Gk(x)$G_k(x)$代入目标函数对α$\alpha$求导并使导数为0，可以把上述的目标函数优化成

(αk,Gk(x))=argmin����������α,G∑i=1mw¯¯¯¯kie[−yi(αG(xi))]=argmin����������α,G∑yi=Gk(xi)w¯¯¯¯kie−α+∑yi≠Gk(xi)w¯¯¯¯kieα=argmin����������α,G(1−ek)e−α+ekeα(6)(7)(8)$$\begin{array}{}\text{(6)}& ({\alpha }_k,G_k(x))& =\underset{\alpha ,G}{\underbrace{\arg min}}\sum _{i=1}^m{\overline{w}}_{ki}{e}^{[-y_i(\alpha G(x_i))]}\text{(7)}& & =\underset{\alpha ,G}{\underbrace{\arg min}}\sum _{y_i=G_k(x_i)}{\overline{w}}_{ki}{e}^{-\alpha }+\sum _{y_i\ne G_k(x_i)}{\overline{w}}_{ki}{e}^{\alpha }\text{(8)}& & =\underset{\alpha ,G}{\underbrace{\arg min}}(1-e_k)e^{-\alpha }+e_k{e}^{\alpha }\end{array}$$

既得最小的α$\alpha$为

α∗k=12log1−ekek$${\alpha }_k^{\ast }=\frac{1}{2}\log \frac{1-e_k}{e_k}$$

最后看样本的权重更新，利用fk(x)=fk−1(x)+αkGk(x)$f_k(x)=f_{k-1}(x)+{\alpha }_k{G}_k(x)$和w¯¯¯¯ki=e[−yifk−1(xi)]${\overline{w}}_{ki}=e^{[-y_i{f}_{k-1}(x_i)]}$可得

w¯¯¯¯k+1,i=e[−yifk(xi)]=e[−yi(fk−1(xi))−yi(αkGk(xi))]=w¯¯¯¯kie[−yiαkGk(x)](9)(10)(11)$$\begin{array}{}\text{(9)}& {\overline{w}}_{k+1,i}& =e^{[-y_i{f}_k(x_i)]}\text{(10)}& & =e^{[-y_i(f_{k-1}(x_i))-y_i({\alpha }_k{G}_k(x_i))]}\text{(11)}& & ={\overline{w}}_{ki}{e}^{[-y_i{\alpha }_k{G}_k(x)]}\end{array}$$

w¯¯¯¯k+1,i${\overline{w}}_{k+1,i}$即接下来要讲到的AdaBoost算法的训练数据权重的更新公式。

三、AdaBoost算法流程

AdaBoost算法既可以解决分类问题，又可以解决回归问题。对于分类问题，此处我们讲述的AdaBoost算法流程主要是针对二分类问题，二分类问题和多分类问题的区别主要在于弱分类器的系数上，本文会介绍AdaBoost SAMME算法如何计算弱分类器的系数；对于回归问题，由于AdaBoost用来解决回归问题的变种有很多，本文只对AdaBoost R2算法做一个介绍。

3.1 输入

m$m$个样本n$n$个特征的训练数据集T={(x1,y1),(x2,y2),⋯,(xm,ym)}$T=\{(x_1,y_1),(x_2,y_2),\cdots ,(x_m,y_m)\}$。

针对二分类问题，yi∈Y={1,−1}$y_i\in Y=\{1,-1\}$。

3.2 输出

最终强学习器G(x)$G(x)$。

3.3 强分类器流程

1. 初始化训练数据的权重

D1=(w11,⋯,w1i,⋯,w1m),w1i=1m,i=1,2,⋯,m$$D_1=(w_{11},\cdots ,w_{1i},\cdots ,w_{1m}),w_{1i}=\frac{1}{m},i=1,2,\cdots ,m$$

2. 生成弱分类器

Gk(x),k=1,2,⋯,K$$G_k(x),k=1,2,\cdots ,K$$

3. 计算弱分类器Gk(x)$G_k(x)$在训练集上的分类误差率为

ek=∑i=1mP(Gk(xi)≠yi)=∑i=1mwkiI(Gk(xi)≠yi)(12)(13)$$\begin{array}{}\text{(12)}& e_k& =\sum _{i=1}^{m}P(G_k(x_i)\ne y_i)\text{(13)}& & =\sum _{i=1}^m{w}_{k_i}I(G_k(x_i)\ne y_i)\end{array}$$

4. 计算Gk(x)$G_k(x)$的权重系数

αk=12log1−ekek二分类问题αk=12log1−ekek+log(R−1)多分类问题(14)(15)$$\begin{array}{}\text{(14)}& & {\alpha }_k=\frac{1}{2}\log \frac{1-e_k}{e_k}\text{二分类问题}\text{(15)}& & {\alpha }_k=\frac{1}{2}\log \frac{1-e_k}{e_k}+\log (R-1)\text{多分类问题}\end{array}$$

二分类问题的权重系数中，可以看出如果分类误差率ek$e_k$越大，则对应的弱分类器的权重系数αk${\alpha }_k$越小，即误差率小的弱分类器权重系数越大。

多分类问题使用的是AdaBoost SAMME算法，其中R$R$为类别数，如果R=2$R=2$，则该多元分类的权重系数将变成二元分类的权重系数。
5. 更新训练数据的权重

Dk+1=(wk+1,1,⋯,wk+1,i,⋯,wk+1,m)wk+1,i=wkiZke−αkyiGk(xi),i=1,2,⋯,m(16)(17)$$\begin{array}{}\text{(16)}& & D_{k+1}=(w_{k+1,1},\cdots ,w_{k+1,i},\cdots ,w_{k+1,m})\text{(17)}& & w_{k+1,i}=\frac{w_{ki}}{Z_k}{e}^{-{\alpha }_k{y}_i{G}_k(x_i)},i=1,2,\cdots ,m\end{array}$$

其中Zk$Z_k$是规范因子

Zk=∑i=1mwkie−αkyiGk(xi)$$Z_k=\sum _{i=1}^m{w}_{ki}{e}^{-{\alpha }_k{y}_i{G}_k(x_i)}$$

从wk+1,i$w_{k+1,i}$的计算公式中可以看出，如果第i$i$个样本分类错误，则yiGk(xi)<0$y_i{G}_k(x_i)<0$，导致样本的权重在第k+1$k+1$个弱分类器中变大，反之，则样本权重在第k+1$k+1$个弱分类器中变小。
6. 结合策略

f(x)=∑k=1KαkGk(X)线性模型G(x)=sign(∑k=1KαkGk(X))最终强分类器G(x)(18)(19)$$\begin{array}{}\text{(18)}& & f(x)=\sum _{k=1}^K{\alpha }_k{G}_k(X)\text{线性模型}\text{(19)}& & G(x)=sign(\sum _{k=1}^K{\alpha }_k{G}_k(X))\text{最终强分类器}G(x)\end{array}$$

3.4 强回归器流程

1. 初始化训练数据的权重

D1=(w11,⋯,w1i,⋯,w1m),w1i=1m,i=1,2,⋯,m$$D_1=(w_{11},\cdots ,w_{1i},\cdots ,w_{1m}),w_{1i}=\frac{1}{m},i=1,2,\cdots ,m$$

2. 生成弱分类器

Gk(x),k=1,2,⋯,K$$G_k(x),k=1,2,\cdots ,K$$

3. 计算弱回归器Gk(x)$G_k(x)$在训练集上的最大误差

Ek=max|yi−Gk(xi)|,i=1,2,⋯,m$$E_k=max|y_i-G_k(x_i)|,i=1,2,\cdots ,m$$

4. 计算每个样本之间的相对误差

eki=|yi−Gk(xi)|Ek$$e_{ki}=\frac{|y_i-G_k(x_i)|}{E_k}$$

此处也可以使用均方误差，即eki=(yi−Gk(xi))2E2k$e_{ki}=\frac{(y_i-G_k(x_i){)}^2}{E_k^2}$
5. 计算第k$k$弱回归器的误差率和权重系数

ek=∑i=1mwkieki误差率αk=ek1−ek权重系数(20)(21)$$\begin{array}{}\text{(20)}& & e_k=\sum _{i=1}^m{w}_{ki}{e}_{ki}\text{误差率}\text{(21)}& & {\alpha }_k=\frac{e_k}{1-e_k}\text{权重系数}\end{array}$$

6. 更新训练数据的权重

wk+1,i=wkiZkα1−ekik$$w_{k+1,i}=\frac{w_{ki}}{Z_k}{\alpha }_k^{1-e_{ki}}$$

其中Zk$Z_k$是规范因子

Zk=∑i=1mwkiα1−ekik$$Z_k=\sum _{i=1}^m{w}_{ki}{\alpha }_k^{1-e_{ki}}$$

7. 结合策略

G(x)=Gk∗(x)$$G(x)=G_{k^{\ast }}(x)$$

其中Gk∗(x)$G_{k^{\ast }}(x)$是所有ln1αk,k=1,2,⋯,K$\ln \frac{1}{{\alpha }_k},k=1,2,\cdots ,K$的中位数对应序号k∗$k^{\ast }$对应的弱回归器

四、AdaBoost算法优缺点

4.1 优点

1. 不容易过拟合
2. 分类精准度高
3. 由于弱分类器既可以是分类器又可以是回归器，使用灵活

4.2 缺点

1. 由于是对训练数据加权，有可能会赋予训练数据中的异常值较高的权重，影响模型的准确度

五、小结

AdaBoost算法并没有使用较深的数学知识，而是推导过程涉及较为复杂的逻辑。如果看完一遍还不是很理解，需要自己多多揣摩。

AdaBoost算法目前是一个比较流行的Boosting算法，他的弱学习器既可以是回归器，又可以是分类器，这也是AdaBoost较为强大的一点。虽然理论上任何学习器都可以作为AdaBoost的弱学习器，但是AdaBoost算法用的较多的弱学习器一般还是决策树和神经网络。

相信有了第一个集成算法AdaBoost的基础，对于接下来的第二个用的较为广泛的Boosting系列算法你也能很快熟悉他，即梯度提升树(gradient boosting decision tree，GBDT)