单变量线性回归

1 模型表示

让我们通过一个例子来开始:这个例子是预测住房价格的,我们要使用一个数据集,数据集包含俄勒冈州波特兰市的住房价格。在这里,我要根据不同房屋尺寸所售出的价格,画出我的数据集。比方说,如果你朋友的房子是1250平方尺大小,你要告诉他们这房子能卖多少钱。那么,你可以做的一件事就是构建一个模型,也许是条直线,从这个数据模型上来看,也许你可以告诉你的朋友,他能以大约220000(美元)左右的价格卖掉这个房子。这就是监督学习算法的一个例子。
《单变量线性回归》

它被称作监督学习是因为对于每个数据来说,我们给出了“正确的答案”,即告诉我们:根据我们的数据来说,房子实际的价格是多少,而且,更具体来说,这是一个回归问题。回归一词指的是,我们根据之前的数据预测出一个准确的输出值,对于这个例子就是价格,同时,还有另一种最常见的监督学习方式,叫做分类问题,当我们想要预测离散的输出值,例如,我们正在寻找癌症肿瘤,并想要确定肿瘤是良性的还是恶性的,这就是0/1离散输出的问题。更进一步来说,在监督学习中我们有一个数据集,这个数据集被称训练集。
以之前的房屋交易问题为例,假使我们回归问题的训练集(Training Set)如下表所示:

《单变量线性回归》

我们将要用来描述这个回归问题的标记如下:
《单变量线性回归》
代表学习算法的解决方案或函数也称为假设(hypothesis)
《单变量线性回归》
这就是一个监督学习算法的工作方式,我们可以看到这里有我们的训练集里房屋价格 我们把它喂给我们的学习算法,学习算法输出一个函数,通常表示为小写 h 表示。h 代表hypothesis(假设),表示一个函数,输入是房屋尺寸大小, 根据输入的x 值来得出y 值,y 值对应房子的价格, 因此h是一个从x 到 y 的函数映射。

选择最初的使用规则代表hypothesis,因而,要解决房价预测问题,我们实际上是要将训练集“喂”给我们的学习算法,进而学习得到一个假设h,然后将我们要预测的房屋的尺寸作为输入变量输入给h,预测出该房屋的交易价格作为输出变量输出为结果。那么,对于我们的房价预测问题,我们该如何表达h ?

一种可能的表达方式为: h θ h_{\theta } hθ(x)= θ 0 \theta _{0} θ0+ θ 1 \theta _{1} θ1x,因为只含有一个特征/输入变量,因此这样的问题叫作单变量线性回归问题。

2 代价函数

《单变量线性回归》
在线性回归中我们有一个像这样的训练集,代表了训练样本的数量,比如 。而我们的假设函数,也就是用来进行预测的函数,是这样的线性函数形式: h θ h_{\theta } hθ(x)= θ 0 \theta _{0} θ0+ θ 1 \theta _{1} θ1x。

接下来要做的便是为我们的模型选择合适的参数(parameters) θ 0 \theta _{0} θ0 θ 1 \theta _{1} θ1,选择的参数决定了我们得到的直线相对于我们的训练集的准确程度,模型所预测的值与训练集中实际值之间的差距(下图中蓝线所指)就是建模误差(modeling error)。
《单变量线性回归》
我们的目标便是选择出可以使得建模误差的平方和能够最小的模型参数。 即使得代价函数最小。
《单变量线性回归》
我们绘制一个等高线图,三个座标分别为 θ 0 \theta _{0} θ0 θ 1 \theta _{1} θ1 和J( θ 0 \theta _{0} θ0+ θ 1 \theta _{1} θ1):

《单变量线性回归》
则可以看出在三维空间中存在一个使得J( θ 0 \theta _{0} θ0+ θ 1 \theta _{1} θ1)最小的点。

代价函数也被称作平方误差函数,有时也被称为平方误差代价函数。我们之所以要求出误差的平方和,是因为误差平方代价函数,对于大多数问题,特别是回归问题,都是一个合理的选择。还有其他的代价函数也能很好地发挥作用,但是平方误差代价函数可能是解决回归问题最常用的手段了。

也许这个函数J( θ 0 \theta _{0} θ0+ θ 1 \theta _{1} θ1)有点抽象,可能你仍然不知道它的内涵,接下来我们要更进一步解释代价函数J的工作原理,并尝试更直观地解释它在计算什么,以及我们使用它的目的。

3 代价函数的直观理解I

参考视频: 2 – 3 – Cost Function – Intuition I (11 min).mkv

让我们通过一些例子来获取一些直观的感受,看看代价函数到底是在干什么。
《单变量线性回归》
《单变量线性回归》

4 代价函数的直观理解II

接下来更深入地学习代价函数的作用,假设你已经认识等高线图,如果你对等高线图不太熟悉的话,这段视频中的某些内容你可能会听不懂,但不要紧,如果你跳过这段视频的话,也没什么关系,不听这节课对后续课程理解影响不大。
《单变量线性回归》

代价函数的样子,等高线图,则可以看出在三维空间中存在一个使得最小的点。
《单变量线性回归》
通过这些图形,希望你能更好地理解这些代价函数所表达的值是什么样的,它们对应的假设是什么样的,以及什么样的假设对应的点,更接近于代价函数的最小值。

当然,我们真正需要的是一种有效的算法,能够自动地找出这些使代价函数取最小值的参数和来。

5 梯度下降

梯度下降是一个用来求函数最小值的算法,我们将使用梯度下降算法来求出代价函数J( θ 0 \theta _{0} θ0+ θ 1 \theta _{1} θ1) 的最小值。

梯度下降背后的思想是:开始时我们随机选择一个参数的组合,计算代价函数,然后我们寻找下一个能让代价函数值下降最多的参数组合。我们持续这么做直到到到一个局部最小值(local minimum),因为我们并没有尝试完所有的参数组合,所以不能确定我们得到的局部最小值是否便是全局最小值(global minimum),选择不同的初始参数组合,可能会找到不同的局部最小值。
《单变量线性回归》
想象一下你正站立在山的这一点上,站立在你想象的公园这座红色山上,在梯度下降算法中,我们要做的就是旋转360度,看看我们的周围,并问自己要在某个方向上,用小碎步尽快下山。这些小碎步需要朝什么方向?如果我们站在山坡上的这一点,你看一下周围,你会发现最佳的下山方向,你再看看周围,然后再一次想想,我应该从什么方向迈着小碎步下山?然后你按照自己的判断又迈出一步,重复上面的步骤,从这个新的点,你环顾四周,并决定从什么方向将会最快下山,然后又迈进了一小步,并依此类推,直到你接近局部最低点的位置。

批量梯度下降(batch gradient descent)算法的公式为:
《单变量线性回归》

其中 α \alpha α 是学习率(learning rate),它决定了我们沿着能让代价函数下降程度最大的方向向下迈出的步子有多大,在批量梯度下降中,我们每一次都同时让所有的参数减去学习速率乘以代价函数的导数。
《单变量线性回归》

在梯度下降算法中,还有一个更微妙的问题,梯度下降中,我们要更新 θ 0 \theta _{0} θ0 θ 1 \theta _{1} θ1) ,当j=0和j=1时,会产生更新,所以你将更新J( θ 0 \theta _{0} θ0)和J( θ 1 \theta _{1} θ1) 。实现梯度下降算法的微妙之处是,在这个表达式中,如果你要更新这个等式,你需要同时更新和,我的意思是在这个等式中,我们要这样更新 θ 0 \theta _{0} θ0:= θ 0 \theta _{0} θ0,并更新 θ 1 \theta _{1} θ1:= θ 1 \theta _{1} θ1

实现方法是:你应该计算公式右边的部分,通过那一部分计算出 θ 0 \theta _{0} θ0 θ 1 \theta _{1} θ1的值,然后同时更新和。

让我进一步阐述这个过程:
《单变量线性回归》

在梯度下降算法中,这是正确实现同时更新的方法。这里不打算解释为什么你需要同时更新,同时更新是梯度下降中的一种常用方法。我们之后会讲到,同步更新是更自然的实现方法。当人们谈到梯度下降时,他们的意思就是同步更新。

2.6 梯度下降的直观理解

上面我们给出了一个数学上关于梯度下降的定义,接下来我们更深入研究一下,更直观地感受一下这个算法是做什么的,以及梯度下降算法的更新过程有什么意义。梯度下降算法如下:
《单变量线性回归》
描述:对 θ \theta θ赋值,使得J( θ \theta θ)按梯度下降最快方向进行,一直迭代下去,最终得到局部最小值。其中 α \alpha α是学习率(learning rate),它决定了我们沿着能让代价函数下降程度最大的方向向下迈出的步子有多大。
《单变量线性回归》

对于这个问题,求导的目的,基本上可以说取这个红点的切线,就是这样一条红色的直线,刚好与函数相切于这一点,让我们看看这条红色直线的斜率,就是这条刚好与函数曲线相切的这条直线,这条直线的斜率正好是这个三角形的高度除以这个水平长度,现在,这条线有一个正斜率,也就是说它有正导数,因此,我得到的新的 θ 1 \theta _{1} θ1,更新后 θ 1 \theta _{1} θ1等于减去一个正数乘以 θ 1 \theta _{1} θ1

这就是我梯度下降法的更新规则:《单变量线性回归》
让我们来看看如果太小或太大会出现什么情况:

如果 α \alpha α太小了,即我的学习速率太小,结果就是只能这样像小宝宝一样一点点地挪动,去努力接近最低点,这样就需要很多步才能到达最低点,所以如果 α \alpha α太小的话,可能会很慢,因为它会一点点挪动,它会需要很多步才能到达全局最低点。

如果 α \alpha α太大,那么梯度下降法可能会越过最低点,甚至可能无法收敛,下一次迭代又移动了一大步,越过一次,又越过一次,一次次越过最低点,直到你发现实际上离最低点越来越远,所以,如果 α \alpha α太大,它会导致无法收敛,甚至发散。

现在,我还有一个问题,当我第一次学习这个地方时,我花了很长一段时间才理解这个问题,如果我们预先把放在一个局部的最低点,你认为下一步梯度下降法会怎样工作?

假设你将初始化在局部最低点,在这儿,它已经在一个局部的最优处或局部最低点。结果是局部最优点的导数将等于零,因为它是那条切线的斜率。这意味着你已经在局部最优点,它使得不再改变,也就是新的等于原来的,因此,如果你的参数已经处于局部最低点,那么梯度下降法更新其实什么都没做,它不会改变参数的值。这也解释了为什么即使学习速率保持不变时,梯度下降也可以收敛到局部最低点。

我们来看一个例子,这是代价函数J( θ \theta θ) 。

《单变量线性回归》

我想找到它的最小值,首先初始化我的梯度下降算法,在那个品红色的点初始化,如果我更新一步梯度下降,也许它会带我到这个点,因为这个点的导数是相当陡的。现在,在这个绿色的点,如果我再更新一步,你会发现我的导数,也即斜率,是没那么陡的。随着我接近最低点,我的导数越来越接近零,所以,梯度下降一步后,新的导数会变小一点点。然后我想再梯度下降一步,在这个绿点,我自然会用一个稍微跟刚才在那个品红点时比,再小一点的一步,到了新的红色点,更接近全局最低点了,因此这点的导数会比在绿点时更小。所以,我再进行一步梯度下降时,我的导数项是更小的,更新的幅度就会更小。所以随着梯度下降法的运行,你移动的幅度会自动变得越来越小,直到最终移动幅度非常小,你会发现,已经收敛到局部极小值。

回顾一下,在梯度下降法中,当我们接近局部最低点时,梯度下降法会自动采取更小的幅度,这是因为当我们接近局部最低点时,很显然在局部最低时导数等于零,所以当我们接近局部最低时,导数值会自动变得越来越小,所以梯度下降将自动采取较小的幅度,这就是梯度下降的做法。所以实际上没有必要再另外减小。

这就是梯度下降算法,你可以用它来最小化任何代价函数,不只是线性回归中的代价函数。

7 梯度下降的线性回归

在以前的视频中我们谈到关于梯度下降算法,梯度下降是很常用的算法,它不仅被用在线性回归上和线性回归模型、平方误差代价函数。在这段视频中,我们要将梯度下降和代价函数结合。我们将用到此算法,并将其应用于具体的拟合直线的线性回归算法里。

梯度下降算法和线性回归算法比较如图:
《单变量线性回归》

对我们之前的线性回归问题运用梯度下降法,关键在于求出代价函数的导数,即:
《单变量线性回归》

我们刚刚使用的算法,有时也称为批量梯度下降。指的是在梯度下降的每一步中,我们都用到了所有的训练样本,在梯度下降中,在计算微分求导项时,我们需要进行求和运算,所以,在每一个单独的梯度下降中,我们最终都要计算这样一个东西,这个项需要对所有个训练样本求和。因此,批量梯度下降法这个名字说明了我们需要考虑所有这一”批”训练样本,而事实上,有时也有其他类型的梯度下降法,不是这种”批量”型的,不考虑整个的训练集,而是每次只关注训练集中的一些小的子集。在后面的课程中,我们也将介绍这些方法。

但就目前而言,应用刚刚学到的算法,你应该已经掌握了批量梯度算法,并且能把它应用到线性回归中了,这就是用于线性回归的梯度下降法。

如果你之前学过线性代数,有些同学之前可能已经学过高等线性代数,你应该知道有一种计算代价函数最小值的数值解法,不需要梯度下降这种迭代算法。在后面的课程中,我们也会谈到这个方法,它可以在不需要多步梯度下降的情况下,也能解出代价函数的最小值,这是另一种称为正规方程(normal equations)的方法。实际上在数据量较大的情况下,梯度下降法比正规方程要更适用一些。

点赞