吴恩达-机器学习笔记(第三周)

六、逻辑回归(Logistic Regression)

6.1 分类问题

  本节介绍了分类问题最流行的一种算法——逻辑回归 (Logistic Regression)算法。逻辑回归的本质:输出在(0,1)内的连续值。
  “回归”体现在输出连续值,但因为逻辑回归解决的是分类问题,将输出值限定在(0,1)之间。对于二分类问题来说,问题有两种结果:1代表正向类,0代表负向类。逻辑回归算法设定一个阈值(如0.5),则输出[ 0.5, 1)代表正向类,(0, 0.5)代表负向类。
  注,逻辑回归输出的值表示概率(和输出范围在(0,1)结合理解),e.g.输出为0.7代表结果为1的概率为70%,结果为0的概率为30%。

6.2 假说表示

  本节解决如何让输出为在(0,1)上的连续值。
  引入逻辑函数(logistic function),将input放入逻辑函数中,使得其输出限定在(0,1)。常用的逻辑函数有:Sigmoid函数——g(z)={1\over(1+e^{-z})}

  因此逻辑回归模型h_θ (x)=g(θ^T X)

python代码实现:

import numpy as np
def sigmoid(z):
   return 1 / (1 + np.exp(-z))

h_θ (x)表示对于input,其output =1的概率,即h_θ (x)=P(y=1|x;θ)

6.3 判定边界

  本节介绍了决策边界(decision boundary)的概念。


  看上图,我们可以看到,当z >= 0时,hθ (x) >= 0.5, 预测 y=1。z < 0时,hθ (x) < 0.5, 预测 y=0。其中,z=θT x= θ01 x12 x2+... ,
  假设现有参数θ 是向量[-3,1,1]。则z=θT x=-3+x1+x2

  则z >= 0 表示x1+x2 >= 3的区域(预测 y=1),z < 0表示x1+x2<3的区域(预测 y=0)。

  我们同样可以用非常复杂的模型来适应非常复杂形状的判定边界。比如用二次模型hθ (x)=g(θ01 x12 x23 x12 4 x22 ),θ=[-1 0 0 1 1]表示下图的分布情况。

6.4 代价函数

  和之前一样,要训练参数θ,需要求代价函数。若用线性回归模型的代价函数(所有模型误差的平方和),得到的代价函数将是一个非凸函数(non-convexfunction),这意味会受许多局部最小值的干扰。


  为此我们必须重新定义代价函数:
J(θ)={1\over m }∑_{i=1}^mCost(h_θ (x^{(i) }),y^{(i)})

其中

hθ (x)与 Cost(hθ (x),y)之间的关系如下图所示:

  如果真实值y=1,则hθ (x)离1越近 cost越小,且当hθ (x)=1时,-log1=0,代表无误差,若y=0,则hθ (x)离0越近 cost越小,且当hθ (x)=0时,-log(1-0)=0,代表无误差。

将Cost函数合并:

  Cost(h_θ (x),y)=-y×log(h_θ (x))-(1-y)×log(1-h_θ (x))

则代价函数J为:

J(θ)={1\over m} ∑_{i=1}^m[-y^{(i)} log(h_θ (x^{(i)}))-(1-y^{(i)} )log(1-h_θ (x^{(i)} ))]
=-{1\over m} ∑_{i=1}^m[y^{(i)} log(h_θ (x^{(i)}))+(1-y^{(i)} )log(1-h_θ (x^{(i)} ))]
  如此得到的代价函数J(θ)会是一个凸函数,并且没有局部最优值。

Python代码实现:
import numpy as np
def cost(theta, X, y):
  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)))
  return np.sum(first - second) / (len(X))

代价函数求导过程

  上图看出,得到的梯度下降算法形式上与线性回归的一样,但由于这里给input套了一层Sigmoid函数,所以实际上是不一样的。另,依然要重视在梯度下降前进行特征缩放

6.5 简化的成本函数和梯度下降

逻辑回归的代价函数:


即:
J(θ)=-{1\over m }∑_{i=1}^m[y^{(i)} log(h_θ (x^{(i)} ))+(1-y^{(i)} )log(1-h_θ (x^{(i)} ))]

h_θ (x)=g(θ^T X)

  再次复习,在逻辑回归里 比以前多学了两个概念,(1)逻辑函数(Sigmoid)是在hypothesis里用到的;(2)改变过的Cost是用来算误差的,在训练θ时用到。

目的-->min J(θ)
梯度下降法,就是利用下面的模板,反复更新θ


  形式上,与求线性回归时一样,每一步 θ都减去 预测误差乘以
x_j^{(i)}
的平均值(α是设置的学习率)。但要时刻记得逻辑回归的时候在input外套了一层Sigmoid函数,所以输出的不再是一个的连续的可直接代表实际事物的数值(如线性回归里 预测出的房价),而是分为哪一类的概率。

6.6 高级优化

  本节介绍一些高级优化算法(e.g.共轭梯度法 BFGS (变尺度法) 、L-BFGS (限制变尺度法) ),来提高逻辑回归的效率。
  由于这些算法比较复杂,暂时能学会用就行,不要求搞明白原理。(打代码要用到时在查用法吧,这里有个印象就好了)。

6.7 多类别分类:一对多

  本节讨论用逻辑回归解决多分类问题。
  我们已经知道如何使用逻辑回归进行二元分类。如,用直线可以将数据集分为正、负两类。


处理多分类问题:

假设有一个训练集,有3个类别。我们可以将一个训练集分成3个二元分类问题:
    类1为正类,类2类3为负类;
    类2为正类,类1类3为负类;
    类3为正类,类1类2为负类.

三分类问题

  得到一系列的模型简记为:
h_θ^{(i)} (x)=p(y=i|x;θ)
其中:
i=(1,2,3....k)

分成3次二分类

  这样我们就有了3个逻辑回归分类器
h_θ^{(i) }(x)
,在预测一个input的时候,会运行所有的分类器,然后找到output最大的那个分类器。比如,当我的input为class1时,要让第1个分类器
h_θ^{(1)}(x)
的输出值(概率)最高,input为class i时,第i个分类器
h_θ^{(i)}(x)
最高。所以第i个分类器的输出最大,就表示是第i类的概率最大。这就是多类别分类问题。

七、正则化(Regularization)

7.1 过拟合的问题

本节讨论了过拟合(over-fitting)的问题。

欠拟合 刚好 过拟合

左图:线性模型,欠拟合,不能很好地适应我们的训练集;
右图:四次方模型,过拟合,能非常好地拟合训练集,但预测新input时,效果可能很差;
中图:最合适。

分类中也如此:


模型x的次数越高,拟合的越好,但越可能过拟合。

如何处理过拟合?

1.丢弃一些不重要的特征。可以是手工选择保留哪些特征,或者使用一些模型选择的算法来帮忙(例如PCA)
2.正则化。 保留所有的特征,但是减少参数的大小(magnitude)。

7.2 代价函数

假设有模型:
h_θ (x)=θ_0+θ_1 x_1+θ_2 x_2^2+θ_3 x_3^3+θ_4 x_4^4
  该模型的过拟合是由高此项产生的,因此,我们要想办法让高次项的系数θ-->0,这样就能很好的拟合了。这就是正则化的基本方法。
  为了减少θ_3和θ_4的大小,我们给θ_3和θ_4设置一点惩罚,操作如下:


  这样选择出的
θ_3和θ_4
就小了许多。
  假如我们有非常多的特征,且不知道该惩罚哪些特征,我们可以惩罚所有特征,并让优化代价函数的程序自动选择惩罚程度,操作如下:
J(θ)={1\over 2m}[∑_{i=1}^m(h_θ (x^{(i)})-y^{(i)})^2+λ∑_{j=1}^nθ_j^2 ]

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

  若λ过大,会使所有θ过小,导致模型
h_θ (x)-→θ_0
,如红色直线所示,造成欠拟合。 因此,正则化要取一个合理的 λ 的值。

7.3 正则化线性回归
正则化线性回归的代价函数为:

    J(θ)={1\over 2m} ∑_{i=1}^m[((h_θ (x^{(i)})-y^{(i)})^2+λ∑_{j=1}^nθ_j^2 )]

梯度下降算法变为:

对上面的算法中j=1,2,...,n时的更新式子进行调整可得:
    θ_j:=θ_j (1-a {λ\over m})-a {1\over m} ∑_{i=1}^m(h_θ (x^{(i)})-y^{(i)})x_j^{(i) }
  可以看出,正则化线性回归的梯度下降算法的变化在于,每次都在原有算法更新规则的基础上令θ值减少了一个额外的值。


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


  图中的矩阵尺寸为 (n+1)*(n+1)。(其实正规方程这块我也没搞太懂,先留个坑)

7.4 正则化的逻辑回归模型
正则化的逻辑回归代价函数:

J(θ)={1\over m }∑_{i=1}^m[-y^{(i)} log(h_θ (x^{(i)} ))-(1-y^{(i)} )log(1-h_θ (x^{(i)} ))]+{λ\over 2m} ∑_{j=1}^nθ_j^2

梯度下降算法
注:

  (1)依旧是,炒冷饭一样的提醒,逻辑回归线性回归的正则化代价函数“形式一样,概念不同”。
  (2)θ_0不参与其中的任何一个正则化。

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

推荐阅读更多精彩内容

  • 以西瓜书为主线,以其他书籍作为参考进行补充,例如《统计学习方法》,《PRML》等 第一章 绪论 1.2 基本术语 ...
    danielAck阅读 4,514评论 0 6
  • 机器学习术语表 本术语表中列出了一般的机器学习术语和 TensorFlow 专用术语的定义。 A A/B 测试 (...
    yalesaleng阅读 1,964评论 0 11
  • https://developers.google.com/machine-learning/crash-cour...
    iOSDevLog阅读 2,656评论 1 11
  • 机器学习系统设计 第九十三课:确定执行的优先级 当你在设计复杂的机器学习系统时所面临的主要问题,并会在如何制定一个...
    带刺的小花_ea97阅读 1,540评论 0 6
  • 前言 我跳过了Octave部分的学习,转而利用这部分时间去研究python如何实现这系列课程的小作业,当作是熟悉一...
    yh_y阅读 1,741评论 0 0