前言
线性回归作为机器学习中经典方法之一,以其形式简单、易于建模而被广泛应用。虽然说线性回归形式简单但却蕴含着机器学习中一些非常重要的思想。许多更能更为强大的非线性模型可在线性模型的基础上通过引入层级结构或者高纬映射而得到。因此线性模型在机器学习领域是极为重要的模型之一。
而线性回归中多元线性回归较之一元线性回归应用范围更广,应用过程也更加复杂。因此有必要对多元线性回归在机器学习相关领域一些问题的一般过程进行归纳和总结。当然这一系列步骤成立的前提在于:通过一定分析,可以预见该问题可以通过多元线性回归来解决。
多元线性回归解决问题的一般方法
## 问题抽象
当遇到一个具体问题时,首先要对具体的问题抽象成数学语言,并以恰当的数学符号来表示,从而建立起数学模型,使得问题更加直观便于分析。 而一般来说线性回归遇到的问题一般是这样的:给定由 d d d个属性描述的示例 x = ( x 1 ; x 2 ; … ; x d ) \boldsymbol{x}=\left(x_{1} ; x_{2} ; \ldots ; x_{d}\right) x=(x1;x2;…;xd),其中$ x_{i} 是 是 是x 在 第 在第 在第i$个属性上的取值,线性模型(linear model)试图通过属性的线性组合来进行预测的函数,即
f ( x ) = w 1 x 1 + w 2 x 2 + … + w d x d + b f(\boldsymbol{x})=w_{1} x_{1}+w_{2} x_{2}+\ldots+w_{d} x_{d}+b f(x)=w1x1+w2x2+…+wdxd+b
一般向量形式写成
f ( x ) = w T x + b f(x)=w^{T} x+b f(x)=wTx+b
其中 w = ( w 1 ; w 2 ; w 3 . . . ; w d ) w=(w_1;w_2;w_3...;w_d) w=(w1;w2;w3...;wd)。 w w w和 b b b学得之后,模型就得以确定。因此通过该步的数学抽象就可以明确,整个线性回归的整个过程主要目的就是确定合适的 w w w 和 b b b使得训练的模型尽可能地与实际情况相吻合。
数据预处理
其实这一步对所有算法都是适用的。毕竟机器学习算法都是建立在数据的基础之上的,因此对于数据的预处理对所有算法来说都是极其重要的。该过程大体上分为以下几步:
- 数据收集:简单来说就是通过各种手段收集自己所需要的数据,例如爬虫、开源数据集等
- 数据清洗:包括数据格式的转化,数据的清洗(处理噪声数据缺失值数据),以及数据的采样
- 数据等价转换:包括统一数据的度量(这在距离计算时十分重要)、零均值化、属性分解以及合并
属性的分解是,一个属性能够分解成多个属性,只有某些子属性对于输出有显著的影响。那我们就可以只存储这些子属性对于输出具有显著影响,那我们就可以只存储这些子属性,而不用去存储原来的属性。其中数据的合并和属性的分解是对立的。将一些子属性合并成一个新属性后,这个属性对于输出的影响会更加显著。那就将这个属性进行合并。
通过数据预处理阶段后,数据会变得更加规整,便于训练模型时直接使用。
确定假设函数
经过数据预处理之后,一般来说就可以确定出一些与输出相关性强的属性,因而可以开始尝试构建出解决问题的假设模型。比如说如果确定出与输出相关的属性是n维的,那么就可以建立起如下的假设函数:
f ( x ) = w 1 x 1 + w 2 x 2 + … + w n x n + b f(\boldsymbol{x})=w_{1} x_{1}+w_{2} x_{2}+\ldots+w_{n} x_{n}+b f(x)=w1x1+w2x2+…+wnxn+b
其中:
- x = ( x 0 , x 1 , x 2 , … , x n ) T x=\left(x_{0}, x_{1}, x_{2}, \ldots, x_{n}\right)^{T} x=(x0,x1,x2,…,xn)T为输入, n + 1 n+1 n+1维的列向量,其中 x 0 = 1 x_0=1 x0=1
- w = ( b , w 1 , w 2 , … , w n ) T w=\left(b, w_{1}, w_{2}, \ldots, w_{n}\right)^{T} w=(b,w1,w2,…,wn)T为参数, n + 1 n+1 n+1维的列向量, b b b是 b i a s bias bias
构建代价函数LMS
在确定假设函数之后,需要确定模型中的参数,具体如何确定那?显然关键在于衡量 f ( x ) f(x) f(x)与 y y y之间的差别,因此需要建立起一个代价函数来进行衡量。均方误差是回归任务重最常见的性能度量,因此我们以均方误差为例,建立起如下的代价函数:
E ( w , b ) = arg min ( w , b ) ∑ i = 1 n ( f ( x i ) − y i ) 2 E(w,b) =\underset{(w, b)}{\arg \min } \sum_{i=1}^{n}\left(f\left(x_{i}\right)-y_{i}\right)^{2} E(w,b)=(w,b)argmini=1∑n(f(xi)−yi)2
模型训练
对于不同的参数,成本函数有不同的值,而我们所需要的是让代价函数最小化即:
( w ∗ , b ∗ ) = arg min ( w , b ) ∑ i = 1 n ( f ( x i ) − y i ) 2 = arg min ( w , b ) ∑ i = 1 n ( y i − w x i − b ) 2 \begin{aligned} \left(w^{*}, b^{*}\right) &=\underset{(w, b)}{\arg \min } \sum_{i=1}^{n}\left(f\left(x_{i}\right)-y_{i}\right)^{2} \\ &=\underset{(w, b)}{\arg \min } \sum_{i=1}^{n}\left(y_{i}-w x_{i}-b\right)^{2} \end{aligned} (w∗,b∗)=(w,b)argmini=1∑n(f(xi)−yi)2=(w,b)argmini=1∑n(yi−wxi−b)2
求解 w w w和 b b b使得 E ( w , b ) E_(w,b) E(w,b)最小化的过程,称之为线性回归模型的最小二乘参数估计(parameter estimation)。通过最小二乘法我们可以推导出通用的线性回归模型如下:
f ( x ^ i ) = x ^ i T ( X T X ) − 1 X T y f\left(\hat{\boldsymbol{x}}_{i}\right)=\hat{\boldsymbol{x}}_{i}^{\mathrm{T}}\left(\mathbf{X}^{\mathrm{T}} \mathbf{X}\right)^{-1} \mathbf{X}^{\mathrm{T}} \boldsymbol{y} f(x^i)=x^iT(XTX)−1XTy
其中:
- X X X:数据集 D D D所表示的 m × ( d + 1 ) m \times(d+1) m×(d+1)大小的矩阵
- ( X T X ) − 1 \left(\mathbf{X}^{\mathrm{T}} \mathbf{X}\right)^{-1} (XTX)−1为 ( X T X ) \left(\mathbf{X}^{\mathrm{T}} \mathbf{X}\right) (XTX)的逆矩阵
如果训练集的矩阵维数不是太大,并且 X T X \mathbf{X}^{\mathrm{T}} \mathbf{X} XTX刚好为满秩矩阵。我们可以直接套用上边公式得出结果。
当然现实任务中 X T X \mathbf{X}^{\mathrm{T}} \mathbf{X} XTX往往不是满秩矩阵。例如在许多任务中我们会遇到大量的变量,其数目甚至超过样例数导致 X X X的列数多于行数,此时 X T X \mathbf{X}^{\mathrm{T}} \mathbf{X} XTX显然不是满秩的。此时需要解出多个 w ^ \hat{w} w^,它们都能使均方差最小化,具体选择哪一种需要通过学习算法的归纳偏好来决定,常见的做法是引入正则化项。当然如果如果训练集的矩阵维数很大这种方法,显然不太合适,此时我们可以考虑通过梯度下降算法来求出参数。关于梯度下降算法的具体使用,可以参考周志华老师的《机器学习》,此处不再详述。
预测与结果优化
在参数求出来之后,我们便建立起了一个针对该问题的线性回归模型,然后便可以使用该模型进行预测。当然可能此时得出的模型与现实问题特别吻合,可能造成输入结果与实际结果有出入。此时,需要对结果进行优化。
总结
本文尝试从实际问题入手,归纳出了多元线性回归解决实际问题的一般方法的步骤。一般来说,多元线性回归在解决实际问题时,需要经过问题抽象、数据预处理、假设函数的确定、构建代价函数、模型训练和预测与结果优化六个部分。一个实际问题经过这几个步骤之后,一般来说便可以建立起一个比较合适的多元线性模型。
参考文献
-
https://www.cnblogs.com/ordili/p/10262569.html
-
https://zhuanlan.zhihu.com/p/49222906
-
https://blog.csdn.net/tuqinag/article/details/54730360
-
https://baike.baidu.com/item/%E6%95%B0%E6%8D%AE%E9%A2%84%E5%A4%84%E7%90%86/2711288#3](https://baike.baidu.com/item/数据预处理/2711288#3)
-
周志华 《机器学习》[M].北京:清华大学出版社,2016,74-8