COURSE1 WEEK2
多维特征
在线性回归中,往往特征不止一个,而是具有多维特征
例如,在预测房价的例子中,我们知道更多的信息:
x 1 x_1 x1:房屋的面积
x 2 x_2 x2:卧室的数目
x 3 x_3 x3:楼层数目
x 4 x_4 x4:房屋的年限
因此,我们每一个特征 x ( i ) x^{(i)} x(i) 的表示变成了向量形式, x j ( i ) x^{(i)}_j xj(i) 表示具体的某的 特征( i 行 j 列)
从而,我们的线性模型公式转化为:
f w , b ( x ) = w 1 x 1 + w 2 x 2 + w 3 x 3 + w 4 x 4 + b f_{w,b}(x) = w_1x_1 + w_2x_2 + w_3x_3 + w_4x_4 + b fw,b(x)=w1x1+w2x2+w3x3+w4x4+b
其中, w i w_i wi可以理解为第 i i i 个特征对目标的贡献程度
进而,推广到更一般的形式:
f w , b ( x ) = w 1 x 1 + w 2 x 2 + ⋯ + w n x n + b f_{w,b}(x) = w_1x_1 + w_2x_2 + \dots + w_nx_n + b fw,b(x)=w1x1+w2x2+⋯+wnxn+b
将所有的参数 w i w_i wi组合在一起,形成向量 w ⃗ = [ w 1 , w 2 , … , w n ] \vec {w} = [w_1, w_2, \dots, w_n] w=[w1,w2,…,wn],将所有的特征 x i x_i xi 组合在一起,形成一个向量 x ⃗ = [ x 1 , x 2 , … , x n ] \vec{x} = [\boldsymbol x_1, \boldsymbol x_2, \dots, \boldsymbol x_n] x=[x1,x2,…,xn]
从而,将模型形式可以写为:
f w ⃗ , b ( x ⃗ ) = w ⃗ ⋅ x ⃗ + b f_{\vec w, b}(\vec x) = \vec w \cdot \vec x + b fw,b(x)=w⋅x+b
其中, ⋅ \cdot ⋅ 代表点乘
这种具有多个特征的线性回归模型叫做多元线性回归
向量化
在多元线性回归中,使用向量化的方法,可以使得代码的编写更加简洁,实现更加快速
如果不使用向量化,在代码的编写中,我们需要计算以下结果:
f w ⃗ , b ( x ⃗ ) = ∑ j = 1 n w j x j + b f_{\vec w, b}(\vec x) = \sum _{j=1}^{n}w_jx_j + b fw,b(x)=j=1∑nwjxj+b
当 n n n 较大时,需要较大的计算量,实现较为复杂
f = 0
for j in range(n):f = f + w[j] * x[j]
f = f + b
所谓向量化,就是把数据都看作向量,在每一步的计算中使用向量的计算。例如多元线性回归模型,将 w w w 和 特征 x x x 进行点乘计算
f = np.dot(w, x) + b
向量化的好处:
- 使得代码更加简洁
- 运算速度更快
多元线性回归的梯度下降
与单变量线性回归的梯度下降相似,唯一不同的是,此时要把参数 w w w 当作是一个向量 w ⃗ \vec w w,因此得到参数更新的公式:
w j = w j − α ∂ ∂ w j J ( w ⃗ , b ) b = b − α ∂ ∂ b J ( w ⃗ , b ) w_j = w_j - \alpha \frac{\partial}{\partial w_j}J(\vec w, b) \\ b = b - \alpha \frac{\partial}{\partial b}J(\vec w, b) wj=wj−α∂wj∂J(w,b)b=b−α∂b∂J(w,b)
正规方程
正规方程即最小二乘法。
由于我们要求解损失函数最小的时候对应的参数值,所以不妨将损失函数看作是参数的函数,然后对损失函数求一阶导函数,令一阶导函数等于 0,求解其极小值点,就对应着最优的参数
特点:
- 仅适用于线性回归
- 解决最小化参数问题(同梯度下降算法),但是不需要迭代
- 当特征较多时( > 10000),运行速度较慢
只要特征变量的数目并不大,标准方程是一个很好的计算参数的替代方法。具体地说,只要特征变量数量小于一万,通常使用标准方程法,而不使用梯度下降法。
特征缩放
使用特征缩放的方法,能够使得梯度下降算法的运行速度得到提升
所谓特征缩放,就是在进行模型训练之前,对数据进行归一化操作
例如,以房价预测为例,特征又房屋面积 x 1 x_1 x1 和卧室数量 x 2 x_2 x2,因此:
p r i c e ^ = w 1 x 1 + x 2 x 2 + b \hat {price} = w_1x_1 + x_2x_2 + b price^=w1x1+x2x2+b
其中, x 1 ∈ [ 300 , 2000 ] , x 2 ∈ [ 0 , 5 ] x_1 \in [300,2000],x_2 \in [0, 5] x1∈[300,2000],x2∈[0,5]数据集:
$x_1 = 2000, x_2 = 5, price = $500K$
由此可以看出,我们的参数 w 1 w_1 w1应该较小, w 2 w_2 w2 应该较大
即,对于一个好的模型来说:
- 当特征的可能指较小时,其参数的合理值将相对较大
- 当特征的可能指较大时,其参数的合理值将相对较小
当每个特征的取值范围相差较大时,特征关系与损失函数图像如下:
如右侧的损失函数梯度图,当我们使用梯度下降算法时,如果学习率设置不当,算法会来回左右横跳动,经过很长一段时间才会收敛到最优值
当我们使用特征缩放时,即将 x 1 x_1 x1和 x 2 x_2 x2进行重新标度,归一化到区间 [ 0 , 1 ] [0,1] [0,1]内,保证了两个特征拥有一个可比较的范围,从而使得损失函数梯度图更像一个圆形,便于算法能够快速收敛到最优点
特征缩放方法
特征缩放的目的是将一列数据变化到某个固定区间(范围)中
均值归一化
将数据归一化到区间 [ − 1 , 1 ] [-1,1] [−1,1]内
x = x − μ x m a x − x m i n x = \frac{x - \mu}{x_{max} - x_{min}} x=xmax−xminx−μ
其中, μ \mu μ 是数据 x x x 的均值
Z-score 归一化
即,将数据转化为均值为0,标准差为1的分布
x = x − μ σ x = \frac{x - \mu}{\sigma } x=σx−μ
其中, μ \mu μ 是数据的均值, σ \sigma σ 是标准差
特征缩放的好坏,具体取决于所有特征进行特征缩放后的取值范围是否尽量一致,以保证梯度下降算法的有效进行
梯度下降法则
如何判断梯度下降是否收敛
一般而言,在模型训练阶段,随着迭代次数的进行,损失值如下图所示:
可以看到的是,当迭代次数大于300时,曲线接近平行,下降的趋势非常平缓,此时意味着我们的梯度下降开始收敛了
通常,可以使用 epsilon
法进行自动收敛测试,即设置收敛阈值 ε = 0.001 \varepsilon = 0.001 ε=0.001,当损失值下降幅度小于阈值时,即认为算法开始收敛,但从实际来看,要想确定一个正确的阈值是非常困难的
如何设置学习率
如果学习率设置过大,则最终结果不容易收敛
如果学习率设置太小,则会导致算法运行较长时间
通过绘制损失函数与迭代次数关系的图像,如果损失函数出现时而下降,时而上升,即不是一直下降的趋势,那么则表明学习率的设置可能较大(也可能是代码存在错误
因此,在实际工作中,一般会选择一系列的值不断去尝试,且在尝试的过程中,只对模型的部分数据进行有限的迭代次数,通过对比来选择最优的学习率
多项式回归
特征工程
在实际问题中,使用的模型往往比较复杂,因此有时需要利用特征工程的方法来对模型加入一些重要的特征
例如,在预测房价时,目前存在临街长度 x 1 x_1 x1 和深度 x 2 x_2 x2,因此房价预测模型为:
f w ⃗ , b ( x ⃗ ) = w 1 x 1 + w 2 x 2 + b f_{\vec w, b}(\vec x) = w_1x_1 + w_2x_2 + b fw,b(x)=w1x1+w2x2+b
但是在实际中,根据生活经验,使用房屋面积作为单特征可能会更好的帮助我们进行预测,因此引入第三个变量房屋面积 x 3 x_3 x3,且 x 3 = x 1 x 2 x_3 = x_1x_2 x3=x1x2,从而我们的模型转化为:
f w ⃗ , b ( x ⃗ ) = w 1 x 1 + w 2 x 2 + w 3 x 3 + b f_{\vec w, b}(\vec x) = w_1x_1 + w_2x_2 +w_3x_3 + b fw,b(x)=w1x1+w2x2+w3x3+b
这种方法叫做创建新特征
多元线性回归+特征工程
将多元线性回归与特征工程的思想结合起来,就是多项式回归的算法,这可以使我们获得更好的数据模型
对于通过面积来预测房价的例子,根据数据集的分布情况,可以看出如果使用二次函数来拟合,效果可能会更好,如下图:
但是考虑到二次函数在达到最高点之后会再次下降,而实际情况中房屋面积越大,价格应该是越高,因此对模型进行调整,改为三次函数模型
同时,在加入高次幂时,要记得对使用特征缩放得方法,保证我们的梯度下降算法有效的进行
除此之外,由于观察到随着面积的增长,价格增长的趋势不在那么陡峭,因此也可以考虑使用平方根函数