|
对于所有这些(u,v)∈rn×rn-和任何解b∈θ,我们有supb∈θp(b,u,v)=nxi=1n-ui(bTi)-vi(bTi)o=-bnxi=1{uiti+viti}=0。因此,对于所有(u,v)∈rn×rn-使得参数bp(b,u,v)=0,p(b,u,v)的最优值为0。结合似然共轭表达式(E.1)和规则条件(E.2)给出了所有(u,v)使得参数bp(b,u,v)=0的拉格朗日对偶函数g(u,v)。对偶GT回归问题(5.2)的形式如下。对于(5.2)isL(u,v,b)=-n(1-κ)+nxi=1ui-log(-vi)-bnxi=1{tiui+tivi},具有有限阶条件uil(u,v,b)=ui-bti=0,vil(u,v,b)=-vi-bti=0,i∈{1,..,n},等价地,在求解uiand vi,ui=bTi,vi=-bTi,i∈{1,..,n}时,(e.3)代入(5.2)的约束条件后,我们得到了(5.2)的矩量法表示。在定理5(i)的证明中证明了解Bb∈θ的存在性。该Hessian矩阵为-∑ni=1{titi+titi/(bti)},其概率为负,引理6为1。因此,GT回归问题(4.1)存在唯一的解,概率接近1。程序(5.2)的解(bu,bv)的存在性从解Bb∈θ的存在性到ML问题(4.1)的有限阶条件和对偶问题(5.2)的动量表示法,当对i∈{1,...,n}设置Bui=bbti,Bvi=-1/(bbti)时。我们现在证明,对于所有的b∈θ,映射(u,v)7→L(u,v,b)在Rn×Rn-中存在最大值。对偶问题(5.2)的序条件为ui,uiL(u,v,b)=1,ui,viL(u,v,b)=0 vi,uiL(u,v,b)=0,vi,viL(u,v,b)=vi。因此,(u,v)7→L(u,v,b)的Hessian矩阵对所有的b∈θ是正的,因此,映射(u,v)7→L(u,v,b)是严格凸的,有唯一解(bu,bv)。当使用(E.3)且BEI=BBTI和Bηi=BBTI,i∈1,...,n}时,程序(5.2)isL(bu,bv,bb)=-n(1-κ)+NXI=1BE+log(Bηi)-NXI=1BE-1=nκ-NXI=1BE-log(Bηi)的值,即ML问题(4.1)在解E.3下的值。定理9的证明。设kbk1,BW=PJKL=1BWLBL。类似于定理8(i)的证明,自适应Lasso GT回归问题的一个等价公式Mismismax(b,e,η)∈θ×rn×rn+nκ-nxi=1ei-log(ηi)-λnkbk1,bw服从ei=bTi,ηi=bTi,i∈{1,..,n},且设(u,v)∈Rn×Rn-表示拉格朗日乘子向量,对应的拉格朗日对偶函数可写成:asg(u,v)=sup(e,η)∈Rn×Rn+nxi=1uiei+viηi-nκ+ei-log(ηi)+supb∈θ(-nxi=1ui(b>ti)-nxi=1vi(b>ti)-nxi=1vi(b>ti)-λnkbk1,bw),其中的项是pni=1{-nκ+ei/2-log(ηi)}的凸共轭项,在定理8(i)的证明中导出。对于第二项(bTi)+NXI=1VI(bTi)+λNKBK1,bw,b∈RJK,它在b中是凸的,但不光滑。为了计算F(b)的次梯度,我们计算了kbk1,BW的次梯度。若权值满足bwl>0forbbl6=0和bwl=0,则加权范数kbk1,bw=max{sb:sl∈{-bwl,bwl}}可以写成2n个线性函数的最大值。由主动函数atb的梯度的所有凸组合给出的kbk1,bwis的子函数(Boyd和Vandenberghe,2008)。我们通过定义s=(s,..,sJK),sl∈{-bwl,bwl}来识别一个主动函数sb,使得sb=kbk1,bw。对于每个L选择sl=bwlif bl>0,并且sl=-bwlif bl<0。如果BL=0,则选择SL=-BWLOR或SL=BWL。因此,takezl=bwlif bl>0-bwlif bl<0-bwlor bwlif bl=0,l∈{1,...,JK}。kbk1,bwis:kbk1,bw=nz:zl≤bwl,l∈{1,...,JK},zb=kbk1,bwo。因此,F(b)的子梯度为:F(b)=(Nxi=1uiti.l+nxi=1viti.l+λnzl,l∈{1,...,JK}),其中zl≤bwl,l∈{1,...,JK}。
|