LASSO回归与L1正则化 西瓜书「建议收藏」

1.结构风险与经验风险在支持向量机部分,我们接触到松弛变量,正则化因子以及最优化函数,在朴素贝叶斯分类,决策树我们也遇到类似的函数优化问题。其实这就是结构风险和经验风险两种模型选择策略,经验风险负责最小化误差,使得模型尽可能的拟合数据,而结构风险则负责规则化参数,使得参数的形式尽量简洁,从而达到防止过拟合的作用.所以针对常见模型,我们都有下式:                           …

大家好,又见面了,我是你们的朋友全栈君。

1.结构风险与经验风险

在支持向量机部分,我们接触到松弛变量,正则化因子以及最优化函数,在朴素贝叶斯分类,决策树我们也遇到类似的函数优化问题。其实这就是结构风险和经验风险两种模型选择策略,经验风险负责最小化误差,使得模型尽可能的拟合数据,而结构风险则负责规则化参数,使得参数的形式尽量简洁,从而达到防止过拟合的作用.所以针对常见模型,我们都有下式:

                                                        LASSO回归与L1正则化 西瓜书「建议收藏」

第一项经验风险L(yi,f(xi,w))衡量真实值与预测值之间的误差,第二项结构风险Ω(w)规则化项使得模型尽可能简单.而第二项Ω(w)一般是模型复杂度的单调函数,模型越复杂,则规则化项的值越大,这里常引入范数作为规则化项,这也就引入了我们常见的L0范数,L1范数以及L2范数.

2.L0范数,L1范数,L2范数与LASSO回归,岭回归

1)广义定义

                L0范数 LASSO回归与L1正则化 西瓜书「建议收藏」 : 向量中非零元素的个数.

                L1范数 LASSO回归与L1正则化 西瓜书「建议收藏」:向量中各个元素绝对值之和

                L2范数 LASSO回归与L1正则化 西瓜书「建议收藏」:向量各元素平方和然后求平方根

L0,L1范数可以实现稀疏化,而L1系数因为有比L0更好的特性而被广泛应用,L2范数在回归里就是岭回归,也叫均值衰减,常用于解决过拟合,通过对向量各元素平方和再求平方根,使得L2范数最小,从而使得参数W的各个元素都接近于0,与L1范数不同,L2范数规划后w的值会接近于0但不到0,而L1范数规范后则可能令w的一些值为0,所以L1范数规范在特征选择中经常用到,而L2范数在参数规则化时经常用到.在回归模型中,通过添加L1,L2范数引入正则化项,便得到了LASSO回归和岭回归:

2)回归模型

常见线性模型回归:

                                                                    LASSO回归与L1正则化 西瓜书「建议收藏」

LASSOO回归:

                                                                    LASSO回归与L1正则化 西瓜书「建议收藏」

岭回归:

                                                                    LASSO回归与L1正则化 西瓜书「建议收藏」

3.嵌入式选择与LASSO回归

这里主要针对西瓜书第11节的内容,对近端梯度下降PGD法解决L1正则化问题进行讨论.

1)优化目标

令▽表示微分算子,对优化目标:

                                                                     LASSO回归与L1正则化 西瓜书「建议收藏」

若f(x)可导,且▽f 满足L-Lipschitz(利普希茨连续条件),即存在常数L>0使得:

                                                        LASSO回归与L1正则化 西瓜书「建议收藏」

2)泰勒展开

则在Xk处我们可以泰勒展开:

                                                LASSO回归与L1正则化 西瓜书「建议收藏」

上式是严格相等,由L-Lipschitz条件我们可以看到:

                                                                    LASSO回归与L1正则化 西瓜书「建议收藏」

这里给出了一个L的下界,且下界的形式与二阶导函数形式类似,从而泰勒展开式的二阶导便通过L替代,从而严格不等也变成了近似:

                                                   LASSO回归与L1正则化 西瓜书「建议收藏」

3)简化泰勒展开式

接下来我们化简上式:

                                       LASSO回归与L1正则化 西瓜书「建议收藏」

                                LASSO回归与L1正则化 西瓜书「建议收藏」LASSO回归与L1正则化 西瓜书「建议收藏」

                                LASSO回归与L1正则化 西瓜书「建议收藏」

其中φ(xk)是与x无关的const常数.

                                                              LASSO回归与L1正则化 西瓜书「建议收藏」

4)简化优化问题

这里若通过梯度下降法对f(x)进行最小化,则每一步下降迭代实际上等价于最小化二次函数f(x),从而推广到我们最上面的优化目标,类似的可以得到每一步的迭代公式:

                                                LASSO回归与L1正则化 西瓜书「建议收藏」

令  

                                                                      LASSO回归与L1正则化 西瓜书「建议收藏」

则我们可以先计算z,再求解优化问题:

                                                            LASSO回归与L1正则化 西瓜书「建议收藏」

5)求解

令xi为x的第i个分量,将上式展开可以看到没有xixj(x≠j)的项,即x的各分量互不影响,所以优化目标有闭式解.这里对于上述优化问题求解需要用到Soft Thresholding软阈值函数,其解为:

                                                            LASSO回归与L1正则化 西瓜书「建议收藏」

对于本例,带入求解即得:

                                                            LASSO回归与L1正则化 西瓜书「建议收藏」

因此,PGD能使LASSO和其他基于L1范数最小化的方法得以快速求解.

4.Soft Thresholding软阈值函数证明

1)软阈值函数

上式求解用到了软阈值函数,下面对软阈值函数的解进行求证,从而更好理解上述求解过程.

先看一下软阈值函数:

                                                           LASSO回归与L1正则化 西瓜书「建议收藏」

2)证明

Proof:

对于优化问题:

                                                                    LASSO回归与L1正则化 西瓜书「建议收藏」

此处X,Z均为n维向量.

展开目标函数:

      LASSO回归与L1正则化 西瓜书「建议收藏」

LASSO回归与L1正则化 西瓜书「建议收藏」LASSO回归与L1正则化 西瓜书「建议收藏」

从而优化问题变为求解N个独立的函数:

                                                                        LASSO回归与L1正则化 西瓜书「建议收藏」

这是我们常见的二次函数,对其求导:

                                                                    LASSO回归与L1正则化 西瓜书「建议收藏」

令导数为0:

                                                                            LASSO回归与L1正则化 西瓜书「建议收藏」

看到两边都有x,所以我们要对上述情况进行讨论:

A.z>λ/2时

假设 x<0 ,   所以  sign(x)=-1 ,但 z-λ/2sign(x)>0 ,所以矛盾.

假设 x>0 ,所以  sign(x) = 1,z-λ/2sign(x)>0,所以极小值在x>0 取得1.

此时极小值小于f(0):

                                            LASSO回归与L1正则化 西瓜书「建议收藏」

再看x<0,

                                      LASSO回归与L1正则化 西瓜书「建议收藏」

                                                LASSO回归与L1正则化 西瓜书「建议收藏」

所以f(x)在负无穷到0单调递减,所以最小值在z-λ/2处取得.

B.z<-λ/2时

假设 x<0 ,   所以  sign(x)=-1 ,z-λ/2sign(x)<0 ,所以极值点在 x<0 处取得.

假设 x>0 ,所以  sign(x) = 1,z-λ/2sign(x)<0,所以矛盾.

此时极值小于f(0):

                                              LASSO回归与L1正则化 西瓜书「建议收藏」

再看 x>0 ,

                                        LASSO回归与L1正则化 西瓜书「建议收藏」

                                                  LASSO回归与L1正则化 西瓜书「建议收藏」

所以f(x)在0到正无穷单调递增,所以最小值在z+λ/2处取得.

C.λ/2<z<λ/2时

假设 x<0 ,   所以  sign(x)=-1 ,z-λ/2sign(x)>0 ,所以矛盾.

假设 x>0 ,所以  sign(x) = 1,z-λ/2sign(x)<0,所以矛盾.

所以x>0,x<0均不满足条件.

所以有:

                                                                LASSO回归与L1正则化 西瓜书「建议收藏」

                                                                        LASSO回归与L1正则化 西瓜书「建议收藏」

                                                                  LASSO回归与L1正则化 西瓜书「建议收藏」

                                                                  LASSO回归与L1正则化 西瓜书「建议收藏」

当△x>0时,由条件z<λ/2:

                                                                LASSO回归与L1正则化 西瓜书「建议收藏」

                                                            LASSO回归与L1正则化 西瓜书「建议收藏」

                                                            LASSO回归与L1正则化 西瓜书「建议收藏」

当△x<0时,由条件z<λ/2:

                                                                LASSO回归与L1正则化 西瓜书「建议收藏」

                                                            LASSO回归与L1正则化 西瓜书「建议收藏」

                                                            LASSO回归与L1正则化 西瓜书「建议收藏」

所以在0处取极小值,同时也是最小值.

综合以上三种情况:

                                                                  LASSO回归与L1正则化 西瓜书「建议收藏」

3)对应西瓜书的L1正则化与LASSO回归

这里的解对应的优化问题是:

                                                                    LASSO回归与L1正则化 西瓜书「建议收藏」

而我们PGD优化的问题是:

                                                                LASSO回归与L1正则化 西瓜书「建议收藏」

对上式同乘2/L不影响极值点位置的选取,所以我们的PGD优化问题变成:

                                                                LASSO回归与L1正则化 西瓜书「建议收藏」

带入综合三种情况的到的最终解:

                                                                LASSO回归与L1正则化 西瓜书「建议收藏」

西瓜书上11.14也就得证了~

总结:

终于看完了西瓜书11章特征选择与稀疏学习,发现从头至尾都在提到用LASSO解决问题,所以就结合第六章的正则化和之前的模型评价,对正则化范数以及LASSO重新认识了一下,书中解决LASSO的大致方法就是通过利普希茨连续条件得到L,带入到优化函数中对函数简化变形,简易优化函数,然后通过软阈值函数得到最后的解.LASSO大致就是这些了,有问题欢迎大家交流~

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/128098.html原文链接:https://javaforall.net

赞 (0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 学1个月爬虫就月赚6000?别被骗了,老师傅告诉你爬虫的真实情况!

    学1个月爬虫就月赚6000?别被骗了,老师傅告诉你爬虫的真实情况!爬虫的四个水平,你在哪一个层次?巅峰爬虫是什么样子?

    2022年5月29日
  • 【技术分享】Python3升级pip3[通俗易懂]

    【技术分享】Python3升级pip3[通俗易懂]python3升级pip相关指令错误命令1:错误命令2:正确命令:因为python默认升级pip的命令是:pipinstall–upgradepip,那python3应该如何升级呢(尤其是当系统中同时含有python2.7和python3的时候)错误命令1:pip3install–upgradepip错误命令2:pip3install–upgradepip3正确命令:python3-mpipinstall–upgradepip因提示permiss

    2022年10月22日
  • JDK8新特性_JDK8面试题

    JDK8新特性_JDK8面试题JDK8新特性lambda表达式要求:接口必须是函数式接口,一个接口如果只有一个抽象方法,那他就是函数式接口,我们可以在一个接口上使用_Functionallnterface_注解package com.atguigu.java;import java.util.Comparator;public class LambaTest { public static void main(String[] args) { Runnable r1 = new Runnable()

    2022年8月8日
  • cs架构的软件中服务器作用,cs架构(cs架构基本原理)

    cs架构的软件中服务器作用,cs架构(cs架构基本原理)用最简单的话,让我明白区别就给分。不要复制的!CS架构,就是你的电脑,需要装个软件,才能连接服务器。而BS架构,就是你的电脑,只需要用浏览器,就可以连接服务器了。1.CS(Client/Server):客户端—-服务器结构。C/S结构在技术上很成熟,它的主要特点是交互性强、具有安全的存取模式、网络通信量低、响应速度快、利于处理大量数据.软件中的CS架构指什么?C/S=Client/Server…

    2022年6月16日
  • 传感器尺寸、像素、DPI分辨率、英寸、毫米的关系

    传感器尺寸、像素、DPI分辨率、英寸、毫米的关系虽然网上有很多这种资料,但是太过于复杂,每个人的说法都不一样,看的让人云里雾里的,我总结了一下,不知道对不对!1.1英寸=25.4mm2.传感器尺寸:传感器的尺寸是指传感器的大小,一般描述大小有两种形式,以IMX386感光元件为例,其传感器尺寸1/2.9英寸,是指传感器对角线为1/2.9英寸;还可以描述成传感器尺寸4.97mm×6.2mm,是指水平(竖直)长(宽)为4.97(6.2)m…

    2022年6月13日
  • Android MD5加密

    Android MD5加密概述在网络中传输明文是一件非常危险的事情,所以通常将密码加密后传至服务器,由服务器保存密文在登录判定时只需比较密文是否相同即可。MD5加密是一种常用的加密算法,全称为“Message-DigestAlgorithm5”,即消息摘要算法,由MD2、MD3、MD4演变过来的,是一种单向加密算法,是一种不可逆的加密方式MD5优点:压缩性:任意长度的数据,算出的MD5值长度都是固定的。容易计算:从原数…

    2022年7月11日

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号