KKT与套索回归的无约束公式


20

L1惩罚回归(又名套索)以两种形式表示。设两个目标函数为

Q1=12||YXβ||22Q2=12||YXβ||22+λ||β||1.
然后将两种不同的制剂是
argminβQ1
||β||1t,
和等同
argminβQ2.
使用Karush-Kuhn-Tucker(KKT)条件,很容易看出第一种配方的平稳性条件等同于采用第二种配方的梯度并将其设置为0。我找不到,也找不到,是第一种配方的补充松弛条件λ(||β||1t)=0如何通过第二种配方的溶液保证得到满足。

Answers:


16

就第一配方中的每个t值而言,第二配方中存在一个λ,使得两个配方具有相同的最小化值这两个配方是等效的β

理由如下:

考虑套索公式: 设最小化子为β,设b=| |。| β*| | 1。我的主张是,如果在第一个公式中设置t=b,那么第一个公式的解也将为β。这是证明:

f(β)=12||YXβ||22+λ||β||1
βb=||β||1t=bβ

min12||YXβ||22 s.t.||β||1b
β^||β^||1<||β||1=bf(β^)<f(β)ββ

由于,因此在解点处满足互补松弛条件。t=bβ

因此,给定带有的套索公式,您可以使用等于套索解决方案的范数的构造约束公式。相反,给定带的约束公式,您会找到一个,从而套索的解等于约束公式的解。λtl1tλ

(如果您了解次梯度,则可以通过求解方程来找到,其中中的X ŤÝ - X β *= λ Ž * Ž *&Element; &PartialD;&| | β * | | 1λXT(yXβ)=λzz||β||1)


1
优秀。一旦您看到了解决方案,您总是会因为没有自己到达那里而感到愚蠢。我假设你的意思是,在发现矛盾,假设我们找到一个β使得| | β | | 1 < | | β * | | 1 = bβ^||β^||1<||β||1=b
goodepic'4

认为弗拉金答案是正确的
bdeonovic 2014年

2
你能否解释为什么f(β^)<f(β)
goofd

这证明了第一个公式的解也必须具有b的l1-范数。如何证明这两个解决方案确实相同?
broncoAbierto

1
此外,套索并不总是有一个独特的解决方案,所以我们不能引用最小化。arxiv.org/pdf/1206.0313.pdf。我们可以,但是,请参阅一套极小的,并表明一些ββ *必须属于这个组。β^β
broncoAbierto

3

我认为elexhobby提出这一证明的想法是一个好主意,但我认为这不是完全正确的。

在示出了用于第一制剂溶液的存在β,使得β< β *通向矛盾,我们只能假设的必要性β= β *,不即β = β *β^β^<ββ^=ββ^=β

我建议改为进行如下操作:

为了方便起见,让我们分别用P 2表示第一和第二个公式。让我们假设P 2有一个独特的解决方案,β *,与β *= b。让P 1有一个解决方案ββ *。然后,我们有ββ *(不能因为约束更大),因此˚F βP1P2P2ββ=bP1β^ββ^β。如果 ˚F β< ˚F β *然后 β *不是解决 P 2,这与我们的假设。如果 ˚F β= ˚F β * β = β *,因为我们认为该解决方案是唯一的。f(β^)f(β)f(β^)<f(β)βP2f(β^)=f(β)β^=β

但是,套索可能有多种解决方案。由引理1 arxiv.org/pdf/1206.0313.pdf我们知道,所有这些解决方案都具有相同的范数(和相同的最小值,当然)。我们将该范数设置为P 1的约束,然后继续。1P1

让我们表示由该组解决方案,以P 2,与β = b &ForAll; β &Element; 小号。设P 1有一个解决方案β&NotElement 小号。然后,我们有ββ &ForAll; β &Element; 小号,因此˚F β˚F β &ForAll; β &Element; 小号。如果fSP2β=b βSP1β^Sβ^ββSf(β^)f(β)βS为一些 β &Element; 小号(并且因此对于所有的)然后 β&Element; 小号,这与我们的假设。如果 ˚F β< ˚F β 为一些 β &Element; 小号然后小号不在设定解决方案,以 P 2。因此,对 P 1的每个解都在 S中,即对 P 1的任何解f(β^)=f(β)βSβ^Sf(β^)<f(β)βSSP2P1SP1 is also a solution to P2. It would remain to prove that the complementary holds too.

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.