支持向量回归模型SVR

支持向量回归模型SVR1.SVM回归模型的损失函数度量    回顾下我们前面SVM分类模型中,我们的目标函数是让12||w||2212||w||22最小,同时让各个训练集中的点尽量远离自己类别一边的的支持向量,即yi(w∙ϕ(xi)+b)≥1yi(w∙ϕ(xi)+b)≥1。如果是加入一个松弛变量ξi≥0ξi≥0,则目标函数是12||w||22+C∑i=1mξi12||w||22+C∑i=1mξi,对应的约束条

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

1. SVM回归模型的损失函数度量

    回顾下我们前面SVM分类模型中,我们的目标函数是让 12||w||22 12||w||22最小,同时让各个训练集中的点尽量远离自己类别一边的的支持向量,即 yi(wϕ(xi)+b)1 yi(w∙ϕ(xi)+b)≥1。如果是加入一个松弛变量 ξi0 ξi≥0,则目标函数是 12||w||22+Ci=1mξi 12||w||22+C∑i=1mξi,对应的约束条件变成: yi(wϕ(xi)+b)1ξi yi(w∙ϕ(xi)+b)≥1−ξi

    但是我们现在是回归模型,优化目标函数可以继续和SVM分类模型保持一致为 12||w||22 12||w||22,但是约束条件呢?不可能是让各个训练集中的点尽量远离自己类别一边的的支持向量,因为我们是回归模型,没有类别。对于回归模型,我们的目标是让训练集中的每个点 (xi,yi) (xi,yi),尽量拟合到一个线性模型 yi =wϕ(xi)+b yi =w∙ϕ(xi)+b。对于一般的回归模型,我们是用均方差作为损失函数,但是SVM不是这样定义损失函数的。

    SVM需要我们定义一个常量 ϵ>0 ϵ>0,对于某一个点 (xi,yi) (xi,yi),如果 |yiwϕ(xi)b|ϵ |yi−w∙ϕ(xi)−b|≤ϵ,则完全没有损失,如果 |yiwϕ(xi)b|>ϵ |yi−w∙ϕ(xi)−b|>ϵ,则对应的损失为 |yiwϕ(xi)b|ϵ |yi−w∙ϕ(xi)−b|−ϵ,这个均方差损失函数不同,如果是均方差,那么只要 yiwϕ(xi)b0 yi−w∙ϕ(xi)−b≠0,那么就会有损失。

    如下图所示,在蓝色条带里面的点都是没有损失的,但是外面的点的是有损失的,损失大小为红色线的长度。

支持向量回归模型SVR

    总结下,我们的SVM回归模型的损失函数度量为:

err(xi,yi)={
0|yiwϕ(xi)+b|ϵ|yiwϕ(xi)b|ϵ|yiwϕ(xi)b|>ϵ
err(xi,yi)={0|yi−w∙ϕ(xi)−b|≤ϵ|yi−w∙ϕ(xi)+b|−ϵ|yi−w∙ϕ(xi)−b|>ϵ

2. SVM回归模型的目标函数的原始形式

    上一节我们已经得到了我们的损失函数的度量,现在可以可以定义我们的目标函数如下:

min12||w||22s.t|yiwϕ(xi)b|ϵ(i=1,2,...m) min12||w||22s.t|yi−w∙ϕ(xi)−b|≤ϵ(i=1,2,…m)

    和SVM分类模型相似,回归模型也可以对每个样本 (xi,yi) (xi,yi)加入松弛变量 ξi0 ξi≥0, 但是由于我们这里用的是绝对值,实际上是两个不等式,也就是说两边都需要松弛变量,我们定义为 ξi,ξi ξi∨,ξi∧, 则我们SVM回归模型的损失函数度量在加入松弛变量之后变为:

min12||w||22+Ci=1m(ξi+ξi) min12||w||22+C∑i=1m(ξi∨+ξi∧)

s.t.ϵξiyiwϕ(xi)bϵ+ξi s.t.−ϵ−ξi∨≤yi−w∙ϕ(xi)−b≤ϵ+ξi∧

ξi0,ξi0(i=1,2,...,m) ξi∨≥0,ξi∧≥0(i=1,2,…,m)

    依然和SVM分类模型相似,我们可以用拉格朗日函数将目标优化函数变成无约束的形式,也就是拉格朗日函数的原始形式如下:

L(w,b,α,α,ξi,ξi,μ,μ)=12||w||22+Ci=1m(ξi+ξi)+i=1mα(ϵξiyi+wϕ(xi)+b)+i=1mα(yiwϕ(xi)bϵξi)i=1mμξii=1mμξi L(w,b,α∨,α∧,ξi∨,ξi∧,μ∨,μ∧)=12||w||22+C∑i=1m(ξi∨+ξi∧)+∑i=1mα∨(−ϵ−ξi∨−yi+w∙ϕ(xi)+b)+∑i=1mα∧(yi−w∙ϕ(xi)−b−ϵ−ξi∧)−∑i=1mμ∨ξi∨−∑i=1mμ∧ξi∧

    其中  μ0,μ0,αi0,αi0 μ∨≥0,μ∧≥0,αi∨≥0,αi∧≥0,均为拉格朗日系数。

3. SVM回归模型的目标函数的对偶形式

    上一节我们讲到了SVM回归模型的目标函数的原始形式,我们的目标是

minw,b,ξi,ξimaxμ0,μ0,αi0,αi0L(w,b,α,α,ξi,ξi,μ,μ) min⏟w,b,ξi∨,ξi∧max⏟μ∨≥0,μ∧≥0,αi∨≥0,αi∧≥0L(w,b,α∨,α∧,ξi∨,ξi∧,μ∨,μ∧)

    和SVM分类模型一样,这个优化目标也满足KKT条件,也就是说,我们可以通过拉格朗日对偶将我们的优化问题转化为等价的对偶问题来求解如下:

maxμ0,μ0,αi0,αi0minw,b,ξi,ξiL(w,b,α,α,ξi,ξi,μ,μ) max⏟μ∨≥0,μ∧≥0,αi∨≥0,αi∧≥0min⏟w,b,ξi∨,ξi∧L(w,b,α∨,α∧,ξi∨,ξi∧,μ∨,μ∧)

    我们可以先求优化函数对于 w,b,ξi,ξi w,b,ξi∨,ξi∧的极小值, 接着再求拉格朗日乘子 α,α,μ,μ α∨,α∧,μ∨,μ∧的极大值。

    首先我们来求优化函数对于 w,b,ξi,ξi w,b,ξi∨,ξi∧的极小值,这个可以通过求偏导数求得:

Lw=0w=i=1m(αiαi)ϕ(xi) ∂L∂w=0⇒w=∑i=1m(αi∧−αi∨)ϕ(xi)

Lb=0i=1m(αiαi)=0 ∂L∂b=0⇒∑i=1m(αi∧−αi∨)=0

Lξi=0Cαμ=0 ∂L∂ξi∨=0⇒C−α∨−μ∨=0

Lξi=0Cαμ=0 ∂L∂ξi∧=0⇒C−α∧−μ∧=0

 

    好了,我们可以把上面4个式子带入 L(w,b,α,α,ξi,ξi,μ,μ) L(w,b,α∨,α∧,ξi∨,ξi∧,μ∨,μ∧)去消去 w,b,ξi,ξi w,b,ξi∨,ξi∧了。

    看似很复杂,其实消除过程和系列第一篇第二篇文章类似,由于式子实在是冗长,这里我就不写出推导过程了,最终得到的对偶形式为:

maxα,αi=1m(ϵyi)αi+(ϵ+yi)αi)12i=1,j=1m(αiαi)(αjαj)Kij max⏟α∨,α∧∑i=1m(ϵ−yi)αi∧+(ϵ+yi)αi∨)−12∑i=1,j=1m(αi∧−αi∨)(αj∧−αj∨)Kij

s.t.i=1m(αiαi)=0 s.t.∑i=1m(αi∧−αi∨)=0

0<αi<C(i=1,2,...m) 0<αi∨<C(i=1,2,…m)

0<αi<C(i=1,2,...m) 0<αi∧<C(i=1,2,…m)


     对目标函数取负号,求最小值可以得到和SVM分类模型类似的求极小值的目标函数如下:

minα,α12i=1,j=1m(αiαi)(αjαj)Kiji=1m(ϵyi)αi+(ϵ+yi)αi min⏟α∨,α∧12∑i=1,j=1m(αi∧−αi∨)(αj∧−αj∨)Kij−∑i=1m(ϵ−yi)αi∧+(ϵ+yi)αi∨

s.t.i=1m(αiαi)=0 s.t.∑i=1m(αi∧−αi∨)=0

0<αi<C(i=1,2,...m) 0<αi∨<C(i=1,2,…m)

0<αi<C(i=1,2,...m) 0<αi∧<C(i=1,2,…m)

    对于这个目标函数,我们依然可以用第四篇讲到的SMO算法来求出对应的 α,α α∨,α∧,进而求出我们的回归模型系数 w,b w,b

4. SVM回归模型系数的稀疏性

    在SVM分类模型中,我们的KKT条件的对偶互补条件为:  αi(yi(wϕ(xi)+b)1)=0 αi∗(yi(w∙ϕ(xi)+b)−1)=0,而在回归模型中,我们的对偶互补条件类似如下:

αi(ϵ+ξi+yiwϕ(xi)b)=0 αi∨(ϵ+ξi∨+yi−w∙ϕ(xi)−b)=0

αi(ϵ+ξiyi+wϕ(xi)+b)=0 αi∧(ϵ+ξi∧−yi+w∙ϕ(xi)+b)=0

    根据松弛变量定义条件,如果 |yiwϕ(xi)b|<ϵ |yi−w∙ϕ(xi)−b|<ϵ,我们有 ξi=0,ξi=0 ξi∨=0,ξi∧=0,此时 ϵ+ξi+yiwϕ(xi)b0,ϵ+ξiyi+wϕ(xi)+b0 ϵ+ξi∨+yi−w∙ϕ(xi)−b≠0,ϵ+ξi∧−yi+w∙ϕ(xi)+b≠0这样要满足对偶互补条件,只有 αi=0,αi=0 αi∨=0,αi∧=0

    我们定义样本系数系数

βi=αiαi βi=αi∧−αi∨

    根据上面 w w的计算式 w=i=1m(αiαi)ϕ(xi) w=∑i=1m(αi∧−αi∨)ϕ(xi),我们发现此时 βi=0 βi=0,也就是说 w w不受这些在误差范围内的点的影响。对于在边界上或者在边界外的点, αi0,αi0 αi∨≠0,αi∧≠0,此时 βi0 βi≠0

5. SVM 算法小结

    这个系列终于写完了,这里按惯例SVM 算法做一个总结。SVM算法是一个很优秀的算法,在集成学习和神经网络之类的算法没有表现出优越性能前,SVM基本占据了分类模型的统治地位。目前则是在大数据时代的大样本背景下,SVM由于其在大样本时超级大的计算量,热度有所下降,但是仍然是一个常用的机器学习算法。

    SVM算法的主要优点有:

    1) 解决高维特征的分类问题和回归问题很有效,在特征维度大于样本数时依然有很好的效果。

    2) 仅仅使用一部分支持向量来做超平面的决策,无需依赖全部数据。

    3) 有大量的核函数可以使用,从而可以很灵活的来解决各种非线性的分类回归问题。

    4)样本量不是海量数据的时候,分类准确率高,泛化能力强。

    SVM算法的主要缺点有:

    1) 如果特征维度远远大于样本数,则SVM表现一般。

    2) SVM在样本量非常大,核函数映射维度非常高时,计算量过大,不太适合使用。

    3)非线性问题的核函数的选择没有通用标准,难以选择一个合适的核函数。

    4)SVM对缺失数据敏感。

   

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

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/146058.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • goland2021永久激活码破解方法

    goland2021永久激活码破解方法,https://javaforall.cn/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

  • centos7.4安装docker_pythondocker

    centos7.4安装docker_pythondocker前言当我们在一台电脑上搭建了python3.6的环境,下次换台电脑,又得重新搭建一次,设置环境变量等操作。好不容易安装好,一会提示pip不是内部或外部命令,一会又提示pip:commandno

  • Netbeans安装教程

    Netbeans安装教程准备工作:电脑安装完JDK能翻墙的浏览器开始:假如没有安装JDK,可以看我写的如何安装JDK1.8.0_101在Oracle官网上下载NetBeans选择相应的版本下载,按照箭头指示方向,即点击AcceotLicenseAgreement然后选择适合自己电脑的版本Windowsx64右面的下载地址建议将安装包放在D盘,新建Netbeans…

    2022年10月22日
  • postgresql数据库报“connections on Unix domain socket “/tmp/.s.PGSQL.5432”?”

    postgresql数据库报“connections on Unix domain socket “/tmp/.s.PGSQL.5432”?”使用postgresql数据库的时候经常遇到的问题:[postgres@test~]$psqlpsql:couldnotconnecttoserver:Nosuchfileordirectory Istheserverrunninglocallyandaccepting connectionsonUnixdomainsocket"/tmp/.s…

  • 京东活动+自动运行脚本+签到2021.05.26更新

    京东活动+自动运行脚本+签到2021.05.26更新先下载Node.js运行环境https://nodejs.org/en下载14.15.4版本。WIN7下载12.20的版本https://nodejs.org/dist/latest-v12.x/node-v12.20.1-x64.msi。只需要一路安装就可以。检测是否安装成功:点击开始-运行-cmd(win+R),打开dos,输入“node–version”检查Node.js版本:只要有显示就说明没问题了。每天大概300-400个京豆,有时候现金红包会有6元左右。下载脚本包:https://..

  • 文末彩蛋 | 这个 Request URL 长得好不一样

    文末彩蛋 | 这个 Request URL 长得好不一样有朋友拿到一个网站请求的链接问这要怎么解密?很明显这不是加密的数据,这是一张图片base64后的结果,第一次写爬虫朋友遇到这样的请求,可能需要琢磨一下这是什么东西。如…

    2022年10月19日

发表回复

您的电子邮箱地址不会被公开。

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