机器学习技法 之 聚合模型(Aggregation Model)

本文深入解析聚合模型原理,探讨模型融合提升性能之道。从简单平均到条件加权,揭示算法如Bagging、AdaBoost如何利用数据多样性增强预测力。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

聚合模型实际上就是将许多模型聚合在一起,从而使其分类性能更佳。

aggregation models: mix or combine hypotheses (for better performance)

下面举个例子:
你有 TTT 朋友,他们对于股票涨停的预测表现为 g1,⋯ ,gTg_1,\cdots ,g_Tg1,,gT。 常见的聚合(aggregation)方法有:

  1. select the most trust-worthy friend from their usual performance
    根据他们的平常表现,选出最值得信任的朋友
    G(x)=gt∗(x) with t∗=argmin⁡t∈{1,2,…,T}Eval (gt−) G(\mathbf{x})=g_{t_{*}}(\mathbf{x}) \text { with } t_{*}=\operatorname{argmin}_{t \in\{1,2, \ldots, T\}} E_{\text {val }}\left(g_{t}^{-}\right) G(x)=gt(x) with t=argmint{1,2,,T}Eval (gt)
  2. mix the predictions from all your friends uniformly
    将所有朋友的预测取平均值
    G(x)=sign⁡(∑t=1T1⋅gt(x)) G(\mathbf{x})=\operatorname{sign}\left(\sum_{t=1}^{T} 1 \cdot g_{t}(\mathbf{x})\right) G(x)=sign(t=1T1gt(x))
  3. mix the predictions from all your friends non-uniformly
    将所有朋友的预测值取加权平均值
    G(x)=sign⁡(∑t=1Tαt⋅gt(x)) with αt≥0 G(\mathbf{x})=\operatorname{sign}\left(\sum_{t=1}^{T} \alpha_{t} \cdot g_{t}(\mathbf{x})\right) \text { with } \alpha_{t} \geq 0 G(x)=sign(t=1Tαtgt(x)) with αt0
  4. combine the predictions conditionally
    根据当前状态 x\mathbf{x}x 确定权重后结合。
    G(x)=sign⁡(∑t=1Tqt(x)⋅gt(x)) with qt(x)≥0G(\mathbf{x})=\operatorname{sign}\left(\sum_{t=1}^{T} q_{t}(\mathbf{x}) \cdot g_{t}(\mathbf{x})\right) \text { with } q_{t}(\mathbf{x}) \geq 0G(x)=sign(t=1Tqt(x)gt(x)) with qt(x)0

学到这里,可能有一种感觉,与模型选择比较相近,并根据直观印象,取平均获得是分类器一定比最好的差,比最差的好。所以会感觉 aggregation 用处不大,那现在看一下, aggregation 的真正的用处是什么?

以下图为例:
在这里插入图片描述
左侧第一个图中,实际上是使用三条竖线或横线实现了二分类,虽然竖线或横线是很弱的一种分类器,但是如此结合便获得了一个较强的分类器,其分类效果好于任何一个分类器独自分类的结果。

右侧第一个图中,是许多直线的取平均值获得的,这种状态存在于数据样本较少时,可以获取一种与SVM类似的效果,虽然这么多直线对于训练样本(采样数据)的分类效果一样,但是对于测试样本(全局数据)可能有更好的分类效果。

所以说真正的 aggregation 并不只是单纯的取平均而已,其可能是为了弥补当前分类器的不足(分类器分类性能较弱,分类器的泛化能力较弱)。即合理的聚合(aggregation)代表了更好的性能(performance)。

Blending

均值融合(uniform blending)

用于分类:

数学表达如下:

G(x)=sign⁡(∑t=1Tgt(x)) G(\mathbf{x})=\operatorname{sign} \left( \sum_{t=1}^{T} g_{t}(\mathbf{x}) \right) G(x)=sign(t=1Tgt(x))

TTT 个人,每人一票。当 gtg_{t}gt 预测值相近,那么性能不变。当 gtg_{t}gt 多样民主时,少数服从多数(majority can correct minority)

在多分类中的数学表达为:

G(x)=argmax⁡1≤k≤K∑t=1T[[gt(x)=k]] G(\mathbf{x})=\underset{1 \leq k \leq K}{\operatorname{argmax}} \sum_{t=1}^{T}\left[\kern-0.15em\left[g_{t}(\mathbf{x})=k\right]\kern-0.15em\right] G(x)=1kKargmaxt=1T[[gt(x)=k]]

用于回归:

G(x)=1T∑t=1Tgt(x) G(\mathbf{x})=\frac{1}{T} \sum_{t=1}^{T} g_{t}(\mathbf{x}) G(x)=T1t=1Tgt(x)

gtg_{t}gt 预测值相近,那么性能不变。当 gtg_{t}gt 多样民主时,一些分类结果 gt(x)>f(x)g_{t}(\mathbf{x})>f(\mathbf{x})gt(x)>f(x) ,另一些分类结果 gt(x)<f(x)g_{t}(\mathbf{x})<f(\mathbf{x})gt(x)<f(x),那么理想状态取平均可以获得最佳解。

综合上述两种需求,多样性的 hypotheses 更容易使得融合模型性能更佳。

现在进行理论分析,其性能是否改善,这里以回归模型为例:
这里的取平均是针对全部的 hypothesis 或者说 TTTgtg_tgt 进行的,并针对的是随机的单个样本。
avg⁡((gt(x)−f(x))2)=avg⁡(gt2−2gtf+f2)=avg⁡(gt2)−2Gf+f2=avg⁡(gt2)−G2+(G−f)2=avg⁡(gt2)−2G2+G2+(G−f)2=avg⁡(gt2)−2avg⁡(gt)G+G2+(G−f)2=avg⁡(gt2−2gtG+G2)+(G−f)2=avg⁡((gt−G)2)+(G−f)2 \begin{aligned} \operatorname{avg}\left(\left(g_{t}(\mathrm{x})-f(\mathrm{x})\right)^{2}\right) &=\operatorname{avg}\left(g_{t}^{2}-2 g_{t} f+f^{2}\right) \\ &=\operatorname{avg}\left(g_{t}^{2}\right)-2 G f+f^{2} \\ &=\operatorname{avg}\left(g_{t}^{2}\right)-G^{2}+(G-f)^{2} \\ &=\operatorname{avg}\left(g_{t}^{2}\right)-2 G^{2}+G^{2}+(G-f)^{2} \\ &=\operatorname{avg}\left(g_{t}^{2}\right)-2\operatorname{avg}\left(g_{t}\right)G+G^{2}+(G-f)^{2} \\ &=\operatorname{avg}\left(g_{t}^{2}-2 g_{t} G+G^{2}\right)+(G-f)^{2} \\ &=\operatorname{avg}\left(\left(g_{t}-G\right)^{2}\right)+(G-f)^{2} \end{aligned} avg((gt(x)f(x))2)=avg(gt22gtf+f2)=avg(gt2)2Gf+f2=avg(gt2)G2+(Gf)2=avg(gt2)2G2+G2+(Gf)2=avg(gt2)2avg(gt)G+G2+(Gf)2=avg(gt22gtG+G2)+(Gf)2=avg((gtG)2)+(Gf)2

也就是说,在对全部训练样本 xn\mathbf{x}_nxn 进行分析取全部误差的平均值。这里 用E\mathcal{E}E 表示平均值。举个例子:1N∑n=1N(gt(xn)−f(xn))2=E(gt−f)2\frac{1}{N}\sum_{n = 1}^{N}\left(g_{t}(\mathrm{x}_n)-f(\mathrm{x}_n)\right)^{2} = \mathcal{E}\left(g_{t}-f\right)^{2}N1n=1N(gt(xn)f(xn))2=E(gtf)2

avg⁡(E(gt−f)2)=avg⁡(E(gt−G)2)+E(G−f)2avg⁡(Eout (gt))=avg⁡(E(gt−G)2)+Eout (G)≥+Eout (G) \begin{aligned} \operatorname{avg}\left(\mathcal{E}\left(g_{t}-f\right)^{2}\right) &=\operatorname{avg}\left(\mathcal{E}\left(g_{t}-G\right)^{2}\right) & +\mathcal{E}(G-f)^{2}\\ \operatorname{avg}\left(E_{\text {out }}\left(g_{t}\right)\right) &=\operatorname{avg}\left(\mathcal{E}\left(g_{t}-G\right)^{2}\right) &+E_{\text {out }}(G) \\ & \geq & +E_{\text {out }}(G) \end{aligned} avg(E(gtf)2)avg(Eout (gt))=avg(E(gtG)2)=avg(E(gtG)2)+E(Gf)2+Eout (G)+Eout (G)

GGG 优于 gtg_tgt 的平均值。

现在假设在分布为 PNP^{N}PN (i.i.d.) 的数据上选取大小为 NNN 的数据集 Dt\mathcal{D}_{t}Dt,并通过 A(Dt)\mathcal{A}\left(\mathcal{D}_{t}\right)A(Dt) 获取 gtg_{t}gt。那么执行无数次可以获取到 gˉ\bar ggˉ,表达式如下:

gˉ=lim⁡T→∞G=lim⁡T→∞1T∑t=1Tgt=EDA(D) \bar{g}=\lim _{T \rightarrow \infty} G=\lim _{T \rightarrow \infty} \frac{1}{T} \sum_{t=1}^{T} g_{t}=\underset{\mathcal{D}}{\mathcal{E}} \mathcal{A}(\mathcal{D}) gˉ=TlimG=TlimT1t=1Tgt=DEA(D)

那么现在用 gˉ\bar{g}gˉ 代替 GGG,之前所求仍然成立,即:

avg⁡(Eout (gt))=avg⁡(E(gt−gˉ)2)+Eout (gˉ) \begin{aligned} \operatorname{avg}\left(E_{\text {out }}\left(g_{t}\right)\right) &=\operatorname{avg}\left(\mathcal{E}\left(g_{t}-\bar{g}\right)^{2}\right) &+E_{\text {out }}(\bar{g}) \\ \end{aligned} avg(Eout (gt))=avg(E(gtgˉ)2)+Eout (gˉ)

其中

  • avg⁡(Eout (gt))\operatorname{avg}\left(E_{\text {out }}\left(g_{t}\right)\right)avg(Eout (gt)) 代表了算法的期望性能(expected performance of A)。
  • Eout (gˉ)E_{\text {out }}(\bar{g})Eout (gˉ) 代表了共识性能(performance of consensus),又叫偏差(bias)
  • avg⁡(E(gt−gˉ)2)\operatorname{avg}\left(\mathcal{E}\left(g_{t}-\bar{g}\right)^{2}\right)avg(E(gtgˉ)2) 代表了共识的期望偏差(expected deviation to consensus),又叫方差(variance)

线性融合(Linear Blending)

用于分类:

数学表达如下:

G(x)=sign⁡(∑t=1Tαt⋅gt(x)) with αt≥0 G(\mathbf{x})=\operatorname{sign}\left(\sum_{t=1}^{T} \alpha_{t} \cdot g_{t}(\mathbf{x})\right) \text { with } \alpha_{t} \geq 0 G(x)=sign(t=1Tαtgt(x)) with αt0

与均值融合相似,有 TTT 个人,但是每人 αt\alpha_tαt 票,而不是都只有一票。

用于回归:
min⁡αt≥01N∑n=1N(yn−∑t=1Tαtgt(xn))2 \min _{\alpha_{t} \geq 0} \frac{1}{N} \sum_{n=1}^{N}\left(y_{n}-\sum_{t=1}^{T} \alpha_{t} g_{t}\left(\mathbf{x}_{n}\right)\right)^{2} αt0minN1n=1N(ynt=1Tαtgt(xn))2

这里重温一下线性回归加非线性转换的结合模型,其数学表达如下:

min⁡wi1N∑n=1N(yn−∑i=1d~wiϕi(xn))2 \min _{w_{i}} \frac{1}{N} \sum_{n=1}^{N}\left(y_{n}-\sum_{i=1}^{\tilde{d}} w_{i} \phi_{i}\left(\mathbf{x}_{n}\right)\right)^{2} wiminN1n=1Nyni=1d~wiϕi(xn)2

可以看出两种非常相似。

所以说线性融合就是线性回归使用假设函数作为非线性转换工具,并且有约束条件。

那么该最优化问题可以写为:
min⁡αt≥01N∑n=1Nerr⁡(yn,∑t=1Tαtgt(xn)) \min _{\alpha_{t} \geq 0} \frac{1}{N} \sum_{n=1}^{N} \operatorname{err}\left(y_{n}, \sum_{t=1}^{T} \alpha_{t} g_{t}\left(\mathbf{x}_{n}\right)\right) αt0minN1n=1Nerr(yn,t=1Tαtgt(xn))

在实际运用中,常常不用约束条件 αt>0\alpha_t > 0αt>0,因为:
 if αt<0⇒αtgt(x)=∣αt∣(−gt(x)) \text { if } \alpha_{t}<0 \Rightarrow \alpha_{t} g_{t}(\mathbf{x})=\left|\alpha_{t}\right|\left(-g_{t}(\mathbf{x})\right)  if αt<0αtgt(x)=αt(gt(x))
也就是说认为 gtg_tgt 的分类错误很高,与预测值常常相反。那么取反便会得到较好性能的分类器。

与模型选择一样,虽然使用训练集获取 gtg_tgt,但是最好使用验证集获取 αt\alpha_tαt

堆叠融合(Stacking or Any Blending)

前面提到的均值融合和线性融合实际上类似于滤波,将预测值乘以一个系数后输出,若将其视为一个模型,那么该模型表达式为 g~(g1,g2,⋯ ,gT)=α1g1+α2g2+⋯+αTgT\tilde g(g_1,g_2,\cdots,g_T) = \alpha_1 g_1 + \alpha_2 g_2 + \cdots + \alpha_T g_Tg~(g1,g2,,gT)=α1g1+α2g2++αTgT。那么 blending 的一般形式便不局限于输入参数的线性组合,可能 g~\tilde gg~ 是也是一个 hypothesis。

Given g1−,g2−,…,gT−g_{1}^{-}, g_{2}^{-}, \ldots, g_{T}^{-}g1,g2,,gT from Dtrain ,\mathcal{D}_{\text {train }},Dtrain , transform (xn,yn)\left(\mathbf{x}_{n}, y_{n}\right)(xn,yn) in Dval \mathcal{D}_{\text {val }}Dval  to (zn=Φ−(xn),yn),\left(\mathbf{z}_{n}=\Phi^{-}\left(\mathbf{x}_{n}\right), y_{n}\right),(zn=Φ(xn),yn), where

学习步骤如下:

  1. 从训练集 Dtrain\mathcal{D}_{\text {train}}Dtrain 中获取 g1−,g2−,…,gT−g_{1}^{-}, g_{2}^{-}, \ldots, g_{T}^{-}g1,g2,,gT,将验证集数据映射到 Z\mathcal ZZ 空间,即 zn=(Φ−(xn),yn)\mathbf{z}_{n}=\left(\Phi^{-}\left(\mathbf{x}_{n}\right), y_{n}\right)zn=(Φ(xn),yn) ,其中映射函数为:Φ−(x)=(g1−(x),…,gT−(x))\Phi^{-}(\mathbf{x})=\left(g_{1}^{-}(\mathbf{x}), \ldots, g_{T}^{-}(\mathbf{x})\right)Φ(x)=(g1(x),,gT(x))
  2. Z\mathcal{Z}Z 空间训练出融合各种模型的模型(函数) g~\tilde{g}g~ === AnyModel ({(zn,yn)})\left(\left\{\left(\mathbf{z}_{n}, y_{n}\right)\right\}\right)({(zn,yn)})
  3. 最终的堆叠融合模型 GANYB(x)=g~(Φ(x))G_{\mathrm{ANYB}}(\mathbf{x})=\tilde{g}(\Phi(\mathbf{x}))GANYB(x)=g~(Φ(x))

优缺点:

  • 很强大(powerful),可以完成有条件的融合(conditional blending)
  • 很容易过拟合(模型复杂度过高)

应用(Blending in Practice)

在这里插入图片描述
在 any blending 的基础上,将原来的 gtg_tgtGGG 结合在一起再做一次融合。

Bagging

blending : 在获取 gtg_tgt 之后,进行聚合;
learning : 在聚合(学习)过程中获取 gtg_tgt

获得多样 gtg_tgt 的方法有:

  • diversity by different models
  • diversity by different parameters: 例如优化方法GD的步长变化多样
  • diversity by algorithmic randomness
  • diversity by data randomness

下面便从数据出发,来满足假设函数的多样性。

那应该怎么做呢,在前面提到有共识便是一个模型的期望表现:

 consensus gˉ= expected gt from Dt∼PN \text { consensus } \bar { g } = \text { expected } g _ { t } \text { from } \mathcal { D } _ { t } \sim P ^ { N }  consensus gˉ= expected gt from DtPN

其优于单个的 gtg_tgt

其由两个部分组成,一个是无穷多个 gtg_tgt,另一个则是丰富的样本数据。对于第一个问题这里提供有限个但相当多个 gtg_tgt,第二个问题只能从手中的数据入手,来创造多样的样本数据集 Dt\mathcal{D}_tDt

拔靴法(Bootstrap Aggregation)

Bagging 实际上就是指 Bootstrap Aggregation,拔靴法实际上是从手中的数据重采样来获得仿真的 Dt\mathcal{D}_tDt。其实现方法是:

  1. 在原有的大小为 NNN 的数据集 D\mathcal{D}D 上,有放回的采样 N′N^\primeN 次获得仿真数据集 D~t→\tilde \mathcal{D} _ t \rightarrowD~t 这一步便是 Bootstrap 操作
  2. 通过 A(D~t)\mathcal A (\tilde \mathcal{D} _ t)A(D~t) 获取 gtg_tgt,再使用均值融合获得:G=Uniform⁡({gt})G = \operatorname {Uniform}(\{g_t\})G=Uniform({gt})

拔靴法(bootstrap aggregation)是一种简单的基于基算法(base algorithm A\mathcal AA)的融合算法(meta algorithm)。方法合理前提是:数据集的多样性和基算法 A\mathcal AA 对于随机数据敏感。

Adaptive Boosting

Adaptive Boosting (AdaBoost )实际上是从 Bagging 的核心 bootstrap 出发实现的一种融合算法。具体实现如下:

加权基算法(Weighted Base Algorithm)

数据集的构造相当于对于不同样本的权重不同,也就是说重采样(Re-sample)过程相当于重赋予权重(Re-weighting)过程:

假设重采样如下:

D={(x1,y1),(x2,y2),(x3,y3),(x4,y4)}⟹ bootstrap D~t={(x1,y1),(x1,y1),(x2,y2),(x4,y4)} \begin{aligned} \mathcal { D } = \left\{ \left( \mathbf { x } _ { 1 } , y _ { 1 } \right) , \left( \mathbf { x } _ { 2 } , y _ { 2 } \right) , \left( \mathbf { x } _ { 3 } , y _ { 3 } \right) , \left( \mathbf { x } _ { 4 } , y _ { 4 } \right) \right\} \\ \stackrel { \text { bootstrap } } { \Longrightarrow } \tilde { \mathcal { D } } _ { t } = \left\{ \left( \mathbf { x } _ { 1 } , y _ { 1 } \right) , \left( \mathbf { x } _ { 1 } , y _ { 1 } \right) , \left( \mathbf { x } _ { 2 } , y _ { 2 } \right) , \left( \mathbf { x } _ { 4 } , y _ { 4 } \right) \right\} \end{aligned} D={(x1,y1),(x2,y2),(x3,y3),(x4,y4)} bootstrap D~t={(x1,y1),(x1,y1),(x2,y2),(x4,y4)}

原来的误差计算如下:
Ein0/1(h)=14∑(x,y)∈D~t[[y≠h(x)]] E _ { \mathrm { in } } ^ { 0 / 1 } ( h ) = \frac { 1 } { 4 } \sum _ { ( \mathbf { x } , y ) \in \tilde { D } _ { t } } \left[\kern-0.15em\left[ y \neq h ( \mathbf { x } ) \right]\kern-0.15em\right] Ein0/1(h)=41(x,y)D~t[[y=h(x)]]

现在则是:

Einu(t)(h)=14∑n=14un(t)⋅[[yn≠h(xn)]] E _ { \mathrm { in } } ^ { \mathrm { u }^{(t)} } ( h ) = \frac { 1 } { 4 } \sum _ { n = 1 } ^ { 4 } u _ { n } ^ { ( t ) } \cdot \left[\kern-0.15em\left[ y _ { n } \neq h \left( \mathbf { x } _ { n } \right) \right]\kern-0.15em\right] Einu(t)(h)=41n=14un(t)[[yn=h(xn)]]

其中 u1=2,u2=1,u3=0,u4=1u_1 = 2,u_2 = 1,u_3 = 0,u_4 = 1u1=2,u2=1,u3=0,u4=1

那么袋中的每一个 gtg_tgt 都是通过最小化加权误差(bootstrap-weighted error)获得的。

所以加权基算法(Weighted Base Algorithm)的数学表达为:

Einu(h)=1N∑n=1Nun⋅err⁡(yn,h(xn)) E _ { \mathrm { in } } ^ { \mathrm { u } } ( h ) = \frac { 1 } { N } \sum _ { n = 1 } ^ { N } u _ { n } \cdot \operatorname { err } \left( y _ { n } , h \left( \mathbf { x } _ { n } \right) \right) Einu(h)=N1n=1Nunerr(yn,h(xn))

那么通过重新赋值获取多样的 gtg_tgt 是另一种可行的方法:

假如两个 gtg_tgt 的获取方法如下:
gt←argmin⁡h∈H(∑n=1Nun(t)[[yn≠h(xn)]])gt+1←argmin⁡h∈H(∑n=1Nun(t+1)[[yn≠h(xn)]]) \begin{aligned} g _ { t } & \leftarrow \underset { h \in \mathcal { H } } { \operatorname { argmin } } \left( \sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t ) } \left[ \kern-0.15em \left[ y _ { n } \neq h \left( \mathbf { x } _ { n } \right) \right]\kern-0.15em\right] \right) \\ g _ { t + 1 } & \leftarrow \underset { h \in \mathcal { H } } { \operatorname { argmin } } \left( \sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } \left[ \kern-0.15em \left[ y _ { n } \neq h \left( \mathbf { x } _ { n } \right) \right]\kern-0.15em\right] \right) \end{aligned} gtgt+1hHargmin(n=1Nun(t)[[yn=h(xn)]])hHargmin(n=1Nun(t+1)[[yn=h(xn)]])

什么时候两个人分类器会很不一样呢?就是当 gtg_tgt 对于权重 un(t)u _ { n } ^ { ( t ) }un(t) 的性能很好,但是 gtg_tgt 对于权重 un(t+1)u _ { n } ^ { ( t + 1) }un(t+1) 的性能很差,最差的 ggg 则是随机值也就是说有 50% 的概率会预测准确。即:

∑n=1Nun(t+1)[[yn≠gt(xn)]]∑n=1Nun(t+1)=12 \frac {\sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } \left[ \kern-0.15em \left[y _ { n } \neq g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] } { \sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } } = \frac { 1 } { 2 } n=1Nun(t+1)n=1Nun(t+1)[[yn=gt(xn)]]=21

所以现在希望的效果是:

∑n=1Nun(t+1)[[yn≠gt(xn)]]∑n=1Nun(t+1)=□t+1□t+1+◯t+1=12, where □t+1=∑n=1Nun(t+1)[[yn≠gt(xn)]]◯t+1=∑n=1Nun(t+1)[[yn=gt(xn)]] \frac {\sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } \left[ \kern-0.15em \left[y _ { n } \neq g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] } { \sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } } = \frac { \square_ { t + 1 } } { \square_ { t + 1 } + \bigcirc_{ t + 1 } } = \frac { 1 } { 2 } , \text { where } \\ \square_ { t + 1 } = \sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } \left[ \kern-0.15em \left[y _ { n } \neq g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right]\\ \bigcirc_{ t + 1 } = \sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } \left[ \kern-0.15em \left[y _ { n } = g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] n=1Nun(t+1)n=1Nun(t+1)[[yn=gt(xn)]]=t+1+t+1t+1=21, where t+1=n=1Nun(t+1)[[yn=gt(xn)]]t+1=n=1Nun(t+1)[[yn=gt(xn)]]

那么通过重新放缩权重(re-scaling (multiplying) weights)便可以实现,即:

对于 gtg_tgt 分类错误的样本:

un(t+1)=◯t⋅un(t) u _ { n } ^ { ( t + 1 ) } = \bigcirc_{ t } \cdot u _ { n } ^ { ( t ) } un(t+1)=tun(t)

对于 gtg_tgt 分类正确的样本:

un(t+1)=□t⋅un(t) u _ { n } ^ { ( t + 1 ) } = \square_{ t } \cdot u _ { n } ^ { ( t ) } un(t+1)=tun(t)

那么在实际中如何实现呢?这里提出放缩系数。

放缩系数(Scaling Factor)

错误率 ϵt\epsilon _ { t }ϵt 定义如下:
ϵt=∑n=1Nun(t)[[yn≠gt(xn)]]∑n=1Nun(t) \epsilon _ { t } = \frac {\sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t ) } \left[ \kern-0.15em \left[y _ { n } \neq g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] } { \sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t) } } ϵt=n=1Nun(t)n=1Nun(t)[[yn=gt(xn)]]

放缩系数的定义如下:
⋆t=1−ϵtϵt \mathbf { \star } _ { t } = \sqrt { \frac { 1 - \epsilon _ { t } } { \epsilon _ { t } } } t=ϵt1ϵt

那么:

[[yn≠gt(xn)]]un(t+1)←un(t)⋅⋆t[[yn=gt(xn)]]un(t+1)←un(t)/⋆t \begin{aligned} \left[ \kern-0.15em \left[y _ { n } \neq g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] \quad u ^ { ( t + 1 ) } _ n &\leftarrow u ^ { ( t ) } _ n \cdot \mathbf { \star } _ { t } \\ \left[ \kern-0.15em \left[y _ { n } = g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] \quad u ^ { ( t + 1 ) } _ n &\leftarrow u ^ { ( t ) } _n / \mathbf { \star } _ { t } \end{aligned} [[yn=gt(xn)]]un(t+1)[[yn=gt(xn)]]un(t+1)un(t)tun(t)/t

Linear Aggregation on the Fly

有了上述的前提,便可以设计一个由数据多样化创造的融合算法,而AdaBoost 除了上述一些前提外,还有一步,那就是 Linear Aggregation on the Fly,在学习中获得线性融合的参数 αt\alpha_tαt,即:

αt=ln⁡(⋆t) \alpha_t = \ln(\mathbf { \star } _ { t }) αt=ln(t)

  1. ϵt→0\epsilon _ { t } \rightarrow 0ϵt0 时,⋆t→inf⁡,ln⁡(⋆t)→inf⁡\mathbf { \star } _ { t } \rightarrow \inf, \ln(\mathbf { \star } _ { t }) \rightarrow \inftinf,ln(t)inf,也就是说当无错误时,给予无穷权重,即当前 gtg_tgt 完全可以完成任务。
  2. ϵt=12\epsilon _ { t } = \frac{1}{2}ϵt=21 时,⋆t=1,ln⁡(⋆t)=0\mathbf { \star } _ { t } = 1, \ln(\mathbf { \star } _ { t }) = 0t=1,ln(t)=0也就是说当错误率为1/2时,不给予权重,即 gtg_tgt 与随机数的性能一样无用。
  3. ϵt→1\epsilon _ { t } \rightarrow 1ϵt1 时,⋆t=0,ln⁡(⋆t)→−inf⁡\mathbf { \star } _ { t } = 0, \ln(\mathbf { \star } _ { t }) \rightarrow -\inft=0,ln(t)inf 也就是说当全错误时,给予负无穷权重,即只需要将分类结果取反便可以获得非常高准确率的 gtg_tgt

AdaBoost 实现步骤

u(1)=[1N,⋯ ,1N]u^{(1)} = \left[\frac{1}{N},\cdots,\frac{1}{N}\right]u(1)=[N1,,N1]
for t=1,⋯ ,Tt = 1,\cdots,Tt=1,,T

  1. A(D,u(t))\mathcal { A } \left( \mathcal { D } , \mathbf { u } ^ { ( t ) } \right)A(D,u(t)) 获取gtg _ { t }gt , 其中 A\mathcal { A }A 用于优化权重为 u(t)\mathbf { u } ^ { ( t ) }u(t) 的加权误差。

  2. u(t)\mathbf { u } ^ { ( t ) }u(t) 更新 u(t+1)\mathbf { u } ^ { ( t+1 ) }u(t+1)
    [[yn≠gt(xn)]]un(t+1)←un(t)⋅⋆t[[yn=gt(xn)]]un(t+1)←un(t)/⋆t\begin{aligned} \left[ \kern-0.15em \left[y _ { n } \neq g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] \quad u ^ { ( t + 1 ) } _ n &\leftarrow u ^ { ( t ) } _ n \cdot \mathbf { \star } _ { t } \\ \left[ \kern-0.15em \left[y _ { n } = g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] \quad u ^ { ( t + 1 ) } _ n &\leftarrow u ^ { ( t ) } _n / \mathbf { \star } _ { t } \end{aligned} [[yn=gt(xn)]]un(t+1)[[yn=gt(xn)]]un(t+1)un(t)tun(t)/t

    其中:ϵt=∑n=1Nun(t+1)[[yn≠gt(xn)]]∑n=1Nun(t+1)\epsilon _ { t } = \frac {\sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } \left[ \kern-0.15em \left[y _ { n } \neq g _ { t } \left( \mathbf { x } _ { n } \right) \right] \kern-0.15em \right] } { \sum _ { n = 1 } ^ { N } u _ { n } ^ { ( t + 1 ) } }ϵt=n=1Nun(t+1)n=1Nun(t+1)[[yn=gt(xn)]]⋆t=1−ϵtϵt\mathbf { \star } _ { t } = \sqrt { \frac { 1 - \epsilon _ { t } } { \epsilon _ { t } } }t=ϵt1ϵt

  3. 计算线性融合系数 αt=ln⁡(⋆t)\alpha_t = \ln(\mathbf { \star } _ { t })αt=ln(t)

  4. 获得最终hypothesis: G(x)=sign⁡(∑t=1Tαtgt(x))G ( \mathbf { x } ) = \operatorname { sign } \left( \sum _ { t = 1 } ^ { T } \alpha _ { t } g _ { t } ( \mathbf { x } ) \right)G(x)=sign(t=1Tαtgt(x))

理论证明(Theoretical Guarantee)

AdaBoost 的 VC bound 如下:
Eout(G)≤Ein(G)+O(O(dvc(H)⋅Tlog⁡T)⏟dvc of all possible G⋅log⁡NN) E _ { \mathrm { out } } ( G ) \leq E _ { \mathrm { in } } ( G ) + O ( \sqrt { \underbrace { O \left( d _ { \mathrm { vc } } ( \mathcal { H } ) \cdot T \log T \right) }_{d_{\mathbf{vc}} \text{ of all possible } G} \cdot \frac { \log N } { N } } ) Eout(G)Ein(G)+O(dvc of all possible GO(dvc(H)TlogT)NlogN)

原作者有证明最多经过 T=log⁡(N)T= \log(N)T=log(N) 次迭代,便可以实现 Ein(G)=0E_{\text{in}}(G) = 0Ein(G)=0,只要基模型比随机数性能优越即可 ϵt≤ϵ<12\epsilon _ { t } \leq \epsilon < \frac { 1 } { 2 }ϵtϵ<21

也就是说,如果基模型 ggg 很弱(weak),但是总比随机数优秀,那么由AdaBoost + A\mathcal AA 获取的 GGG 也会很强(strong)。

决策树桩(Decision Stump)

数学表达如下:
hs,i,θ(x)=s⋅sign⁡(xi−θ) h _ { s , i , \theta } ( \mathbf { x } ) = s \cdot \operatorname { sign } \left( x _ { i } - \theta \right) hs,i,θ(x)=ssign(xiθ)

一共有三个参数 特征(feature) iii,阈值(threshold)θ\thetaθ ,方向(direction)sss。其实现的功能便是一个分界点,特征(feature) iii 表达的是在第 iii 维的分解点,阈值(threshold)θ\thetaθ 代表在本维度的分界点位于 θ\thetaθ,方向(direction)sss 代表了分界点两边的样本类型。

这是一个弱模型,但是将其作为 AdaBoost 的基模型便可以实现高精度预测了,并且效率很高,时间复杂度为:O(d⋅Nlog⁡N)O(d \cdot N \log N)O(dNlogN)

若使用 Decision Stump 作为 AdaBoost 的基模型,假设一个简单的数据集如下分布:
在这里插入图片描述
那么其学习过程的一种状态可能如下表达:

在这里插入图片描述

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FlameAlpha

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值