吴恩达机器学习——第七章 正则化(Regularization)

7.1 过拟合问题

到现在为止,我们已经学习了几种不同的学习算法,包括线性回归和逻辑回归,它们能够有效地解决许多问题,但是当将它们应用到某些特定的机器学习应用时,会遇到过拟合(over-fitting)的问题,可能会导致它们效果很差。

接下来,我将会解释什么是过度拟合问题,并且会谈论一种称为正则化(regularization)的技术,它可以改善或者减少过度拟合问题。

如果我们有非常多的特征,我们通过学习得到的假设可能能够非常好地适应训练集(代价函数可能几乎为0),但是可能会不能推广到新的数据。

下图是一个回归问题的例子:

第一个模型是一个线性模型,欠拟合,不能很好地适应我们的训练集;第三个模型是一个四次方的模型,过于强调拟合原始数据,而丢失了算法的本质:预测新数据。我们可以看出,若给出一个新的值使之预测,它将表现的很差,是过拟合,虽然能非常好地适应我们的训练集但在新输入变量进行预测时可能会效果不好;而中间的模型似乎最合适。

分类问题中也存在这样的问题:

就以多项式理解,x的次数越高,拟合的越好,但相应的预测的能力就可能变差。

问题是,如果我们发现了过拟合问题,应该如何处理?

1.丢弃一些不能帮助我们正确预测的特征。可以是手工选择保留哪些特征,或者使用一些模型选择的算法来帮忙(例如PCA

2.正则化。保留所有的特征,但是减少参数的大小(magnitude)。



7.2 代价函数

上面的回归问题中如果我们的模型是:h_\theta (x)=\theta _0+\theta _1x_1+\theta _2{x_2}^2+\theta _3{x_3}^3 +\theta _4{x_4}^4

我们可以从之前的事例中看出,正是那些高次项导致了过拟合的产生,所以如果我们能让这些高次项的系数接近于0的话,我们就能很好的拟合了。

所以我们要做的就是在一定程度上减小这些参数\theta 的值,这就是正则化的基本方法。我们决定要减小\theta _3\theta _4的大小,我们要做的便是修改代价函数,在其中\theta _3\theta _4设置一点惩罚。这样做的话,我们在尝试最小化代价时也需要将这个惩罚纳入考虑中,并最终导致选择较小一些的\theta _3\theta _4

修改后的代价函数如下:min\frac{1}{2m}[\sum\nolimits_{i=1}^m(h_\theta (x^{(i)})-y^{(i)} ] )^2+{1000\theta _3}^2+{10000\theta _4}^2]

通过这样的代价函数选择出的\theta _3\theta _4对预测结果的影响就比之前要小许多。假如我们有非常多的特征,我们并不知道其中哪些特征我们要惩罚,我们将对所有的特征进行惩罚,并且让代价函数最优化的软件来选择这些惩罚的程度。这样的结果是得到了一个较为简单的能防止过拟合问题的假设:J(\theta )=\frac{1}{2m}[\sum\nolimits_{i=1}^m(h_\theta (x^{(i)})-y^{(i)} ] )^2+\lambda \sum\nolimits_{j=1}^n\theta _j^2]

其中\lambda 又称为正则化参数(Regularization Parameter)。注:根据惯例,我们不对\theta _0进行惩罚。经过正则化处理的模型与原模型的可能对比如下图所示:

如果选择的正则化参数\lambda 过大,则会把所有的参数都最小化了,导致模型变成h_\theta (x)=\theta _0,也就是上图中红色直线所示的情况,造成欠拟合。

那为什么增加的一项\lambda \sum\nolimits_{j=1}^n\theta _j^2可以使\theta 的值减小呢?

因为如果我们令\lambda 的值很大的话,为了使Cost Function尽可能的小,所有的\theta 的值(不包括\theta _0)都会在一定程度上减小。

但若\lambda 的值太大了,那么\theta (不包括\theta _0)都会趋近于0,这样我们所得到的只能是一条平行于x轴的直线。

所以对于正则化,我们要取一个合理的\lambda 的值,这样才能更好的应用正则化。

回顾一下代价函数,为了使用正则化,让我们把这些概念应用到到线性回归和逻辑回归中去,那么我们就可以让他们避免过度拟合了。


7.3 正则化线性回归

对于线性回归的求解,我们之前推导了两种学习算法:一种基于梯度下降,一种基于正规方程。

正则化线性回归的代价函数为:J(\theta )=\frac{1}{2m}[\sum\nolimits_{i=1}^m(h_\theta (x^{(i)})-y^{(i)} ] )^2+\lambda \sum\nolimits_{j=1}^n\theta _j^2]

如果我们要使用梯度下降法令这个代价函数最小化,因为我们未对进行正则化,所以梯度下降算法将分两种情形:

\theta _0:=\theta _0-a\frac{1}{m} \sum\nolimits_{i=1}^m (h_\theta (x^{(i)})-y^{(i)})x_0^{(i)})

\theta _j:=\theta _j-a[\frac{1}{m} \sum\nolimits_{i=1}^m (h_\theta (x^{(i)})-y^{(i)})x_j^{(i)}+\frac{\lambda }{m} \theta _j]

对上面的算法中j=1,2,...n时的更新式子进行调整可得:

\theta _j:=\theta _j(1-a\frac{\lambda }{m}) -a\frac{1}{m} \sum\nolimits_{i=1}^m (h_\theta (x^{(i)})-y^{(i)})x_j^{(i)}

可以看出,正则化线性回归的梯度下降算法的变化在于,每次都在原有算法更新规则的基础上令\theta 值减少了一个额外的值(j=0)。

我们同样也可以利用正规方程来求解正则化线性回归模型,方法如下所示:

图中的矩阵尺寸为(n+1)*(n+1)


7.4  正则化的逻辑回归模型

针对逻辑回归问题,我们在之前的课程已经学习过两种优化算法:我们首先学习了使用梯度下降法来优化代价函数J(\theta ),接下来学习了更高级的优化算法,这些高级优化算法需要你自己设计代价函数J(\theta )

自己计算导数同样对于逻辑回归,我们也给代价函数增加一个正则化的表达式,得到代价函数:J(\theta )=\frac{1}{m} \sum_{i=1}^m[-y^{(i)}log(h_\theta (x^{(i)}))-(1-y^{(i)})log(1-h^\theta (x^{(i)}))]+\frac{\lambda }{2m}  \sum_{j=1}^n{\theta _j}^2

python代码:

import numpy as np

def costReg(theta, X, y, learningRate):

               theta = np.matrix(theta)

               X = np.matrix(X)

               y = np.matrix(y)

               first = np.multiply(-y, np.log(sigmoid(X*theta.T)))

               second = np.multiply((1 - y), np.log(1 - sigmoid(X*theta.T)))

               reg = (learningRate / (2 * len(X))* np.sum(

                          np.power(theta[:,1:theta.shape[1]],2))

               return np.sum(first - second) / (len(X)) + reg

要最小化该代价函数,通过求导,得出梯度下降算法为:

\theta _0:=\theta _0-a\frac{1}{m} \sum\nolimits_{i=1}^m (h_\theta (x^{(i)})-y^{(i)})x_0^{(i)})

\theta _j:=\theta _j-a[\frac{1}{m} \sum\nolimits_{i=1}^m (h_\theta (x^{(i)})-y^{(i)})x_j^{(i)}+\frac{\lambda }{m} \theta _j]

for j=1,2,...n

注:看上去同线性回归一样,但是知道h_\theta (x)=g(\theta ^TX),所以与线性回归不同。


参考:https://github.com/fengdu78/Coursera-ML-AndrewNg-Notes

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,445评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,889评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,047评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,760评论 1 276
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,745评论 5 367
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,638评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,011评论 3 398
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,669评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,923评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,655评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,740评论 1 330
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,406评论 4 320
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,995评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,961评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,197评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,023评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,483评论 2 342

推荐阅读更多精彩内容

  • 其中第五章是一些软件教程,这里就跳过了,继续跟着吴恩达老师学习接下来的内容。 6.1 分类问题 在分类问题中,你要...
    Colleen_oh阅读 1,555评论 1 5
  • 多变量线性回归(Linear Regression with Multiple Variables) 多维特征 n...
    钢的炼金术士阅读 979评论 0 1
  • 以西瓜书为主线,以其他书籍作为参考进行补充,例如《统计学习方法》,《PRML》等 第一章 绪论 1.2 基本术语 ...
    danielAck阅读 4,484评论 0 6
  • 最近在整理 Coursera 上 Andrew Ng 的 Machine Learning 课程的笔记,总共近 8...
    kamidox阅读 1,282评论 2 2
  • 今天是我正式操作公众号项目,对于一个网络小白难度可想而知,好在有心理准备,单单是一个菜单设置就弄了多半个小时,又是...
    林强_a818阅读 71评论 0 0