第8章:遗憾界
编辑:赵志民,詹好
本章前言
本章的内容围绕学习理论中的遗憾(Regret)概念展开(有的教材里也翻译为“悔”)。通常,我们使用超额风险(Excess Risk)来评估批量学习的分类器性能,而用遗憾来评估在线学习的分类器性能。二者的不同在于,前者衡量的是整个学习过程结束后所得到的分类器性能,可以理解为学习算法最终输出的模型与假设空间内最优模型的风险之差;而后者衡量的是算法运行过程中,所产生的模型与假设空间内最优模型的损失之差的和。
8.1 【概念解释】超额风险与遗憾的区别
8.1介绍了遗憾这一评估指标的基本概念,我们在此基础上梳理一下其与超额风险这一评估指标的区别。
超额风险这一评估指标被定义为:
ER=E(x,y)∼D[l(wT+1,(x,y))]−w∈WminE(x,y)∼D[l(w,(x,y))]
其中,ER 指的是excess risk,等式右边的前半部分 E(x,y)∼D[l(wT+1,(x,y))] 指的是模型 wT+1 的风险,等式右边的后半部分 minw∈WE(x,y)∼D[l(w,(x,y))] 指的是假设空间内的最优模型的风险。值得注意的是,这里的评估是在整个数据集上进行的,也正是因为如此,我们必须要引入期望的操作。
而遗憾这一评估指标,被定义为:
regret=t=1∑Tft(wt)−w∈Wmint=1∑Tft(w)
其中,ft(wt) 指的是:
t=1∑Tl(wt,(xt,yt))−w∈Wmint=1∑Tl(w,(xt,yt))
由于wt的计算过程与样本(xt,yt) 无关,而是与(x1,y1),...,(xt−1,yt−1) 有关,因此可以直接使用 l(w,(xt,yt)) 来衡量性能。
由此,我们可以总结出二者之间的两个主要区别:首先,超额风险引入了期望,而遗憾没有;其次,超额风险是在所有数据上进行的一次性计算,而遗憾是对多次损失的一个求和。同时,由于在线学习不依赖于任何分布假设,因此适用于非独立同分布样本或固定分布的情形。
8.2 【案例分享】Maler 算法
在8.2.3节的170页末尾,作者提到了Maler算法(Multiple Sub-algorithms & Learning Rates)(详细证明参考:Adaptivity and Optimality: A Universal Algorithm for Online Convex Optimization),这是一个能够自适应选择最优专家的在线学习算法,并在不同类型的损失函数上实现最优的遗憾界限:
- 一般凸函数:R(T)≤OT)
- 指数凹函数:R(T)≤O(dlogT)
- 强凸函数:R(T)≤O(logT)
这里T表示时间总步数,d表示特征空间的维度。
下面,我们简要补充Maler算法的原理和实现。
假设和定义
-
假设 1(梯度有界性):所有损失函数 ft(x) 的梯度被 G 所有界:
∀t>0,x∈Dmax∥∇ft(x)∥≤G
-
假设 2(行动集的直径有界性):行动集 D 的直径被 D 所有界:
x1,x2∈Dmax∥x1−x2∥≤D
-
定义 1(凸函数):函数 f:D→R 是凸的,如果:
f(x1)≥f(x2)+∇f(x2)⊤(x1−x2),∀x1,x2∈D
-
定义 2(强凸函数):函数 f:D→R 是 λ-强凸的,如果:
f(x1)≥f(x2)+∇f(x2)⊤(x1−x2)+2λ∥x1−x2∥2,∀x1,x2∈D
-
定义 3(指数凹函数):函数 f:D→R 是 α-指数凹的(简称 α-exp-concave),如果:
exp(−αf(x))是凹的
元算法(Maler)
输入:学习率 ηc,η1,η2,…,专家的先验权重 π1c,π1η1,s,π1η2,s…,以及 π1η1,l,π1η2,l,…。
- 对于每个回合 t=1,…,T:
-
从凸专家算法(专家 1)获取预测 xtc,从指数凹专家算法(专家 2)和强凸专家算法(专家 3)分别获取 xtη,l 和 xtη,s。
-
执行:
xt=πtcηc+∑η(πtη,sη+πtη,lη)πtcηcxtc+∑η(πtη,sηxtη,s+πtη,lηxtη,l)
-
观察梯度 gt 并发送给所有专家算法。
-
对所有的 η 更新权重:
πt+1c=Φtπtce−ct(xtc),πt+1η,s=Φtπtη,se−stη(xtη,s),πt+1η,l=Φtπtη,le−ltη(xtη,l)
其中:
Φt=η∑(πtη,se−stη(xtη,s)+πtη,le−ltη(xtη,l))+πtce−ct(xtc)
凸专家算法(专家 1)
- x1c=0
- 对于每个回合 t=1,…,T:
- 将 xtc 发送给元算法
- 从元算法接收梯度 gt
- 更新:
xt+1c=ΠDId(xtc−ηcGtD∇ct(xtc))
其中 ∇ct(xtc)=ηcgt
指数凹专家算法(专家 2)
- 输入:学习率 η
- x1η,l=0,β=21min{4GlD1,1},Gl=25D7,Σ1=β2D21Id
- 对于每个回合 t=1,…,T:
- 将 xtη,l 发送给元算法
- 从元算法接收梯度 gt
- 更新:
Σt+1xt+1η,l=Σt+∇ltη(xtη,l)∇ltη(xtη,l)⊤=ΠDΣt+1(xtη,l−β1Σt+1−1∇ltη(xtη,l))
其中 ∇ltη(xtη,l)=ηgt+2η2gtgt⊤(xtη,l−xt)
强凸专家算法(专家 3)
- 输入:学习率 η
- x1η,s=0
- 对于每个回合 t=1,…,T:
- 将 xtη,s 发送给元算法
- 从元算法接收梯度 gt
- 更新:
xt+1η,s=ΠDId(xtη,s−2η2G2t1∇stη(xtη,s))
其中 ∇stη(xtη,s)=ηgt+2η2G2(xtη,s−xt)
8.3 【证明补充】随机多臂赌博机的遗憾界
172页中定理8.3给出了随机多臂赌博机的遗憾界,我们在此基础上对公式(8.42)至(8.47)证明过程进行补充。
首先,(8.42)给出当μ∗(p)+p2lnt≤μi(q)+q2lnt成立时,必然有三种可能情况中的一种成立。但这三种情况并不是互斥的,因此显得不直观,这里将第二种情况做了细微调整,即:
μ∗(p)+p2lnt≤μ∗,μ∗≤μi(q)+q2lnt,μi(q)+q2lnt≤μi(p)
此时,构造(8.44)和(8.45)的逻辑更加顺畅。我们令l=⌈(2lnT)/Δi2⌉,则(8.45)转化为:
P(μ∗≤μi+q2lnt)=0,q≥l
代入(8.44),可得:
E[niT]≤⌈Δi22lnT⌉+2t=1∑T−1p=1∑t−1q=l∑t−1t−4≤Δi22lnT+1+2t=1∑T−1p=1∑tq=1∑tt−4≤Δi22lnT+1+2T→+∞limt=1∑T−1t−2
根据p-级数判别法,当p=2>1时,级数收敛,因此limT→+∞∑t=1T−1t−2是有界的。至于该级数的具体值,对定理的结论没有影响,因此我们可以将其视为一个常数,然后带入后续推导中。为了证明的完整性,我们对此进行简要说明。
limT→+∞∑t=1T−1t−2的取值在数学界被称为Basel问题,推导过程涉及诸多前置定理,感兴趣的读者可以查看这个讲义:The Basel Problem - Numerous Proofs。此处提供另一种在微积分变换中常见的缩放方法:
t=1∑T−1t−2≤1+∫1T−1x21dx=1+(−x1)∣1T−1=2−T1
对不等式两边同时取极限,可得:
T→+∞limt=1∑T−1t−2≤2
代入(8.46),同样可得类似(8.47)的结论。
这里继续沿用书中给出的limT→+∞∑t=1Tt−2=6π2,代入(8.46)得到遗憾界(8.47):
E[regret]≤i=1∑KΔi22lnT+O(1)
此时(8.46)变为:
E[niT]≤i=∗∑KΔi2lnT+(1+3π2)Δi=O(KlogT)
观察(8.47)可知,求和公式中的每一项符合对钩函数的构造,即:
f(x)=Ax+xB,x>0,A>0,B>0
这里x=Δi,A=1+3π2,B=2lnT,因此无论Δi过大或过小时,都会导致遗憾界的上界变大。另外,遗憾界跟摇臂的个数K呈线性关系,当K越大时,遗憾界也越大。
8.4 【概念解释】线性赌博机
176页的8.3.2节介绍了线性赌博机的概念,我们在此基础上对参数估计部分进行补充。
为了估计线性赌博机的参数,我们将原问题转化为岭回归问题,即(8.52):
f(w)=(Y−wTX)T(Y−wTX)+λwTw
为了求得最优解w∗,我们令f′(w)=0,可推导出(8.53):
∂w∂f(w)=−2XT(Y−wTX)+2λw→XTY→w∗=0=(XTX+λI)w=(XTX+λI)−1XTY
相比于每次传入新数据(xt,yt)时从头计算wt,这里巧妙地利用了 Sherman-Morrison-Woodbury 公式将任何形如(A+uvT)−1的矩阵逆转化为可逆矩阵A和列向量u,v之间的运算,在O(d2)的时间复杂度内完成参数的更新。
8.5 【证明补充】Sherman-Morrison-Woodbury (或 Woodbury) 公式
177页的 Sherman-Morrison-Woodbury 公式变种是矩阵求逆中的一个重要工具,它可以通过已知矩阵的逆来快速计算被低秩修正的矩阵的逆。
该公式如下所示:
(A+UCV)−1=A−1−A−1U(C−1+VA−1U)−1VA−1
其中,A 是一个 n×n 的矩阵,C 是 k×k 的矩阵,U 和 V 是 n×k 的矩阵,(8.54)中C为单位矩阵。
证明
该公式可以通过验证 A+UCV 与其假设的逆(公式右侧)的乘积是否为单位矩阵来证明。我们对以下乘积进行计算:
(A+UCV)[A−1−A−1U(C−1+VA−1U)−1VA−1]
逐步推导如下:
====={I+UCVA−1}−{U(C−1+VA−1U)−1VA−1+UCVA−1U(C−1+VA−1U)−1VA−1}I+UCVA−1−(U+UCVA−1U)(C−1+VA−1U)−1VA−1I+UCVA−1−UC(C−1+VA−1U)(C−1+VA−1U)−1VA−1I+UCVA−1−UCVA−1I
8.6 【证明补充】单样本的近似梯度
第181页的引理8.2给出了单样本条件下的梯度近似公式,本节将提供该引理的完整证明过程。
Eu∈S[f(x+δu)u]=dδ∇Ev∈B[f(x+δv)]
其中:
- d 为空间的维数;
- δ 为任意正数;
- B 为单位球的空间,即 B={v∈Rd∣∥v∥≤1};
- S 为单位球的表面,即 S={u∈Rd∣∥u∥=1}。
证明
为了证明上述等式,我们将分三个步骤进行推导。
1. 表达左边的期望
首先,考虑左边的期望:
Eu∈S[f(x+δu)u]=Vold−1(S)1∫Sf(x+δu)udS(u)
其中,Vold−1(S) 表示 (d−1) 维单位球面的体积,dS(u) 为球面上的微分面积元素。
进行变量替换,令 w=δu。此时:
- 当 u∈S 时,w∈δS;
- 球面上的微分面积元素变化为 dS(u)=δd−1dS(w),因为每个维度按 δ 缩放,(d−1) 维体积按 δd−1 缩放。
将变量替换代入期望的表达式:
Eu∈S[f(x+δu)u]=Vold−1(S)1∫Sf(x+δu)udS(u)=Vold−1(S)⋅δd−11∫δSf(x+w)δwdS(w)
简化后得到:
Eu∈S[f(x+δu)u]=Vold−1(δS)1∫δSf(x+w)∥w∥wdS(w)
2. 表达右边的期望及其梯度
接下来,考虑右边的期望:
Ev∈B[f(x+δv)]=Vold(B)1∫Bf(x+δv)dv
其中,Vold(B) 表示 d 维单位球的体积,dv 为体积上的微分元素。
同样进行变量替换,令 w=δv。则:
- 当 v∈B 时,w∈δB;
- 微分体积元素变化为 dv=δddw,因为每个维度按 δ 缩放,体积按 δd 缩放。
代入后得到:
Ev∈B[f(x+δv)]=Vold(B)⋅δd1∫δBf(x+w)dw=Vold(δB)1∫δBf(x+w)dw
为了计算 ∇Ev∈B[f(x+δv)],令:
F(x)=Ev∈B[f(x+δv)]=Vold(δB)1∫δBf(x+w)dw
梯度作用在积分上,由于 x 和 w 是独立变量,可以将梯度算子移入积分内部:
∇F(x)=Vold(δB)1∫δB∇xf(x+w)dw
注意到:
∇xf(x+w)=∇wf(x+w)
这是因为 x 和 w 的关系是通过相加连接的,故梯度对 x 的作用等同于对 w 的作用。
根据散度定理,有:
∫δB∇wf(x+w)dw=∫δSf(x+w)n(w)dS(w)
其中,δS 是半径为 δ 的球面,n(w) 为点 w 处的单位外法向量。因此:
∇F(x)=Vold(δB)1∫δSf(x+w)∥w∥wdS(w)
3. 关联两边的表达式
将步骤 1 和步骤 2 的结果进行对比,可以得到:
Eu∈S[f(x+δu)u]=Vold−1(δS)Vold(δB)∇Ev∈B[f(x+δv)]
为了确定系数,我们需要利用 d 维球的体积与表面积之间的关系。
d 维球的体积与半径 δ 的关系为:
Vold(δB)=δd⋅Vold(B)
而球面的表面积与半径 δ 的关系为:
Vold−1(δS)=δd−1⋅Vold−1(S)
结合这两个关系,可以得到:
Vold(δB)=∫0δVold−1(rS)dr=∫0δVold−1(S)rd−1dr=dVold−1(S)⋅δd=dδ⋅Vold−1(δS)
带入上述等式中,得证:
Eu∈S[f(x+δu)u]=dδ∇Ev∈B[f(x+δv)]
8.7 【证明补充】凸赌博机的在线梯度下降
182页中引理8.3给出了凸赌博机的随机版本在线梯度下降,我们在此给出完整的证明过程。
设 f1,f2,…,fT:W→R 为一列凸且可微的函数,ω1,ω2,…,ωT∈W 的定义满足 ω1 为任意选取的点,且 ωt+1=ΠW(ωt−ηgt),其中 η>0,且 g1,…,gT 是满足 E[gt∣ωt]=∇ft(ωt) 的随机向量变量,且 ∥gt∥≤l,其中 l>0。则当 η=lTΛ 时,有:
t=1∑TE[ft(ωt)]−ω∈Wmint=1∑Tft(ω)≤lΛT
证明:
设 ω⋆ 为在 W 中使 ∑t=1Tft(ω) 最小化的点。由于 ft 是凸且可微的,我们可以使用梯度界定 ft(ωt) 和 ft(ω⋆) 之间的差异:
ft(ω⋆)−ft(ωt)≥∇ft(ωt)⊤(ω⋆−ωt)=E[gt∣ωt]⊤(ω⋆−ωt)
对该不等式取期望,得到:
E[ft(ωt)−ft(ω⋆)]≤E[gt⊤(ωt−ω⋆)]
我们使用 ∥ωt−ω⋆∥2 作为潜在函数。注意到 ∥ΠW(ω)−ω⋆∥≤∥ω−ω⋆∥,因此:
∥ωt+1−ω⋆∥2=∥ΠW(ωt−ηgt)−ω⋆∥2≤∥ωt−ηgt−ω⋆∥2=∥ωt−ω⋆∥2+η2∥gt∥2−2η(ωt−ω⋆)⊤gt≤∥ωt−ω⋆∥2+η2l2−2η(ωt−ω⋆)⊤gt
整理后得到:
gt⊤(ωt−ω⋆)≤2η∥ωt−ω⋆∥2−∥ωt+1−ω⋆∥2+η2l2
因此,我们有:
t=1∑TE[ft(ωt)]−t=1∑Tft(ω⋆)=t=1∑TE[ft(ωt)−ft(ω⋆)]≤t=1∑TE[gt⊤(ωt−ω⋆)]≤t=1∑TE[2η∥ωt−ω⋆∥2−∥ωt+1−ω⋆∥2+η2l2]=2ηE[∥ω1−ω⋆∥2]−E[∥ωT+1−ω⋆∥2]+2Tηl2≤2ηE[∥ω1−ω⋆∥2]+2Tηl2≤2ηΛ2+2Tηl2
代入 η=lTΛ 可得最终结果。
8.8 【证明补充】凸赌博机的缩减投影误差
182页中引理8.4给出了凸赌博机的缩减投影误差,我们在此给出完整的证明过程。
设 f1,f2,…,fT:W→R 为一列凸且可微的函数且 ∀ω∈W,i∈[T] 满足 ∣fi(ω)∣≤c,有:
ω∈(1−α)Wmint=1∑Tft(ω)−ω∈Wmint=1∑Tft(ω)≤2αcT
证明
显然,(1−α)W⊆W。因此,有:
ω∈(1−α)Wmint=1∑Tft(ω)=ω∈Wmint=1∑Tft((1−α)ω)
由于每个ft是凸函数,且0∈W,则我们有:
ω∈Wmint=1∑Tft((1−α)ω)≤ω∈Wmint=1∑Tαft(0)+(1−α)ft(ω)=ω∈Wmint=1∑Tα(ft(0)−ft(ω))+ft(ω)
最后,由于对于任意ω∈W和t∈{1,…,T},我们有∣ft(ω)∣≤c,因此可以得出:
t=1∑Tω∈Wminα(ft(0)−ft(ω))+ft(ω)≤ω∈Wmint=1∑T2αc+ft(ω)=2αcT+ω∈Wmint=1∑Tft(ω)
进行适当移项即可得原不等式。
8.9 【证明补充】凸赌博机的遗憾界
182页中定理8.5给出了凸赌博机的遗憾界,在证明开始时,作者对η,α,δ的取值进行了限定。我们可以发现这些取值不是很直观,证明给出的解释也较为分散,部分取值与证明略有出入,因此我们在此进行补充。
对于步长η,在缩放(8.87)中 E[∑t=1Tf^t(zt)]−minw∈(1−α)W∑t=1Tf^t(w) 时,为使用引理8.3创造条件,因此采用步长η=l′TΛ。根据(8.89)的推导,我们可令Λ=Λ2且l′=δdc,此时,将η=(dc/δ)TΛ2带入到更新公式(8.76)中即可得到(8.88)。
对于缩减系数α与扰动系数δ,可以一同考虑这两个系数的取值。观察(8.91)第一个不等式的形式,我们发现这是一个关于δ的对钩函数:
f(δ)=Aδ+δB+C
假设α的取值与δ无关,那么:
A=3lT,B=dcΛ2T,C=2αcT
令f′(δ)=0,可得:
δ∗=T−1/43ldcΛ2
此时,f(δ)的最小值为:
f(δ∗)=O(T3/4)
如果我们想加速收敛,则可将α的取值与δ相关联。根据上面的结论,当迭代次数T足够大时,必然有δ→0。因此,不妨取α=Λ1δ,代入(8.91)中并利用对钩函数f(δ)的性质,得到:
δ∗=T−1/43(lΛ1+c)dcΛ1Λ2f(δ∗)=O(T3/4)
进一步地,可以发现,δ∗的取值并不唯一,这是因为(8.91)的第二个不等式缩放并非必需。如果取δ∗=T−1/43lΛ1+2cdcΛ1Λ2,同样可以得到更紧致的遗憾界,并保证定理的结论不变。