【机器学习】模型过拟合的处理方法整理

港控/mmm° 2022-06-02 01:42 297阅读 0赞

**过拟合标准定义**：给定一个假设空间H，一个假设  h  属于H，如果存在其他的假设  h’  属于H,使得在训练样例上  h  的错误率比  h’  小，但在整个实例分布上  h’  比  h  的错误率小，则假设  h  过度拟合训练数据。 —-《Machine Learning》Tom M.Mitchell

因为在统计学习中，假设数据满足独立同分布（i.i.d，independently and identically distributed），即当前已产生的数据可以对未来的数据进行推测与模拟，因此都是使用历史数据建立模型，即使用已经产生的数据去训练，然后使用该模型去拟合未来的数据。

但是一般独立同分布的假设往往不成立，即数据的分布可能会发生变化（distribution drift），并且可能当前的数据量过少，不足以对整个数据集进行分布估计，因此往往需要防止模型过拟合，提高模型泛化能力。而为了达到该目的的最常见方法便是：正则化，即在对模型的目标函数（objective function）或代价函数（cost function）加上正则项。   
  在对模型进行训练时，有可能遇到训练数据不够，即训练数据无法对整个数据的分布进行估计的时候，或者在对模型进行过度训练（overtraining）时，常常会导致模型的过拟合（overfitting）。如下图所示：   
![这里写图片描述][20151026205032330]   
  通过上图可以看出，随着模型训练的进行，模型的复杂度会增加，此时模型在训练数据集上的训练误差会逐渐减小，但是在模型的复杂度达到一定程度时，模型在验证集上的误差反而随着模型的复杂度增加而增大。此时便发生了过拟合，即模型的复杂度升高，但是该模型在除训练集之外的数据集上却不work。

## **出现过拟合的一些原因** ##

（1）**建模样本抽取错误**，包括（但不限于）样本数量太少，抽样方法错误，抽样时没有足够正确考虑业务场景或业务特点，等等导致抽出的样本数据不能有效足够代表业务逻辑或业务场景；   
（2）**样本里的噪音数据干扰过大**，大到模型过分记住了噪音特征，反而忽略了真实的输入输出间的关系；   
（3）**建模时的“逻辑假设”到了模型应用时已经不能成立了**。任何预测模型都是在假设的基础上才可以搭建和应用的，常用的假设包括：假设历史数据可以推测未来，假设业务环节没有发生显著变化，假设建模数据与后来的应用数据是相似的，等等。如果上述假设违反了业务场景的话，根据这些假设搭建的模型当然是无法有效应用的。   
（4）**参数太多、模型复杂度高**   
（5）**决策树模型**。如果我们对于决策树的生长没有合理的限制和修剪的话，决策树的自由生长有可能每片叶子里只包含单纯的事件数据(event)或非事件数据（no event），可以想象，这种决策树当然可以完美匹配（拟合）训练数据，但是一旦应用到新的业务真实数据时，效果是一塌糊涂。   
（6）**神经网络模型**。   
a.由于对样本数据,可能存在隐单元的表示不唯一,即产生的分类的决策面不唯一.随着学习的进行, BP算法使权值可能收敛过于复杂的决策面,并至极致.   
b.权值学习迭代次数足够多(Overtraining),拟合了训练数据中的噪声和训练样例中没有代表性的特征.

## **解决方案及原理** ##

**（1）权值衰减. 主要应用在神经网络模型中**   
它在每次迭代过程中以某个小因子降低每个权值,这等效于修改E的定义,加入一个与网络权值的总量相应的   
惩罚项,此方法的动机是保持权值较小,避免weight decay,从而使学习过程向着复杂决策面的反方向偏。   
**（2）适当的stopping criterion**   
![70][]   
**在二次误差函数的情况下，关于早停止和权值衰减类似结果的原因说明**。椭圆给出了常数误差函数的轮廓线，Wml表示误差函数的最小值。如果权向量的起始点为原点，按照局部负梯度的方向移动，那么它会沿着曲线给出的路径移动。通过对训练过程早停止，我们找到了一个权值向量w。定性地说，它类似于使用检点的权值衰减正则化项，然后最小化正则化误差函数的方法得到的权值。   
**（3）验证数据**   
一个最成功的方法是在训练数据外再为算法提供一套验证数据,应该使用在验证集合上产生最小误差   
的迭代次数,不是总能明显地确定验证集合何时达到最小误差.   
Typically 30% of training patterns;Validation set error is checked each epoch;   
Stop training if validation error goes up   
**（4）交叉验证**   
交叉验证方法在可获得额外的数据提供验证集合时工作得很好,但是小训练集合的过度拟合问题更为严重.   
原理图   
![20160425090721762][]   
**（5）添加正则项**。L1正则更加容易产生稀疏解、L2正则倾向于让参数w趋向于0.   
**（6）针对树模型**   
**a.在树过于大之前便停止生长**   
每个叶中至少需要多少个数据（threshold）   
如何判断这个阈值（threshold）是重点【可以考虑用假设检验/P-值】   
**b.等树生长到足够大之后进行修剪**   
修剪枝叶，直到任何改动都会降低正确率

## **Tips** ##

（1）增加样本的全面性和数量；   
（2）控制模型的复杂度；   
（3）不要过度训练   
（4）模型融合本质上也是一种提高泛化能力的方法

### 防止过拟合用到一些方法 early stopping、数据集扩增（Data augmentation）、正则化（Regularization）、Dropout 详解 ###

### Early stopping ###

对模型进行训练的过程即是对模型的参数进行学习更新的过程，这个参数学习的过程往往会用到一些迭代方法，如梯度下降（Gradient descent）学习算法。Early stopping便是一种迭代次数截断的方法来防止过拟合的方法，即在模型对训练数据集迭代收敛之前停止迭代来防止过拟合。   
  Early stopping方法的具体做法是，在每一个Epoch结束时（一个Epoch集为对所有的训练数据的一轮遍历）计算validation data的accuracy，当accuracy不再提高时，就停止训练。这种做法很符合直观感受，因为accurary都不再提高了，在继续训练也是无益的，只会提高训练的时间。那么该做法的一个重点便是怎样才认为validation accurary不再提高了呢？并不是说validation accuracy一降下来便认为不再提高了，因为可能经过这个Epoch后，accuracy降低了，但是随后的Epoch又让accuracy又上去了，所以不能根据一两次的连续降低就判断不再提高。一般的做法是，在训练的过程中，记录到目前为止最好的validation accuracy，当连续10次Epoch（或者更多次）没达到最佳accuracy时，则可以认为accuracy不再提高了。此时便可以停止迭代了（Early Stopping）。这种策略也称为“No-improvement-in-n”，n即Epoch的次数，可以根据实际情况取，如10、20、30……

### 数据集扩增 ###

在数据挖掘领域流行着这样的一句话，“有时候往往拥有更多的数据胜过一个好的模型”。因为我们在使用训练数据训练模型，通过这个模型对将来的数据进行拟合，而在这之间又一个假设便是，训练数据与将来的数据是独立同分布的。即使用当前的训练数据来对将来的数据进行估计与模拟，而更多的数据往往估计与模拟地更准确。因此，更多的数据有时候更优秀。但是往往条件有限，如人力物力财力的不足，而不能收集到更多的数据，如在进行分类的任务中，需要对数据进行打标，并且很多情况下都是人工得进行打标，因此一旦需要打标的数据量过多，就会导致效率低下以及可能出错的情况。所以，往往在这时候，需要采取一些计算的方式与策略在已有的数据集上进行手脚，以得到更多的数据。   
  通俗得讲，数据机扩增即需要得到更多的符合要求的数据，即和已有的数据是独立同分布的，或者近似独立同分布的。一般有以下方法：

*  从数据源头采集更多数据
 *  复制原有数据并加上随机噪声
 *  重采样
 *  根据当前数据集估计数据分布参数，使用该分布产生更多数据等

### 正则化方法 ###

正则化方法是指在进行目标函数或代价函数优化时，在目标函数或代价函数后面加上一个正则项，一般有L1正则与L2正则等。

*  L1正则   
      L1正则是基于L1范数，即在目标函数后面加上参数的L1范数和项，即参数绝对值和与参数的积项，即：   
      
    其中代表原始的代价函数，是样本的个数，就是正则项系数，权衡正则项与项的比重。后面那一项即为L1正则项。   
      在计算梯度时，的梯度变为：   
      
    其中，是符号函数，那么便使用下式对参数进行更新：   
      
    对于有些模型，如线性回归中（L1正则线性回归即为Lasso回归），常数项的更新方程不包括正则项，即：   
      
    其中，梯度下降算法中，，而在梯度上升算法中则相反。   
      从上式可以看出，当为正时，更新后会变小；当为负时，更新后会变大；因此L1正则项是为了使得那些原先处于零（即）附近的参数往零移动，使得部分参数为零，从而降低模型的复杂度（模型的复杂度由参数决定），从而防止过拟合，提高模型的泛化能力。   
      其中，L1正则中有个问题，便是L1范数在0处不可导，即在0处不可导，因此在为0时，使用原来的未经正则化的更新方程来对进行更新，即令，这样即： 
 *  L2正则   
      L2正则是基于L2范数，即在目标函数后面加上参数的L2范数和项，即参数的平方和与参数的积项，即：   
      
    其中代表原始的代价函数，是样本的个数，与L1正则化项前面的参数不同的是，L2项的参数乘了，是为了便于计算以及公式的美感性，因为平方项求导有个2，就是正则项系数，权衡正则项与项的比重。后面那一项即为L2正则项。   
      L2正则化中则使用下式对模型参数进行更新：   
      
    对于有些模型，如线性回归中（L2正则线性回归即为Ridge回归，岭回归），常数项的更新方程不包括正则项，即：   
      
    其中，梯度下降算法中，，而在梯度上升算法中则相反。   
      从上式可以看出，L2正则项起到使得参数变小加剧的效果，但是为什么可以防止过拟合呢？一个通俗的理解便是：更小的参数值意味着模型的复杂度更低，对训练数据的拟合刚刚好（奥卡姆剃刀），不会过分拟合训练数据，从而使得不会过拟合，以提高模型的泛化能力。   
      在这里需要提到的是，在对模型参数进行更新学习的时候，有两种更新方式，mini-batch （部分增量更新）与 full-batch（全增量更新），即在每一次更新学习的过程中（一次迭代，即一次epoch），在mini-batch中进行分批处理，先使用一部分样本进行更新，然后再使用一部分样本进行更新。直到所有样本都使用了，这次epoch的损失函数值则为所有mini batch的平均损失值。设每次mini batch中样本个数为，那么参数的更新方程中的正则项要改成：   
      
      
    而full-batch即每一次epoch中，使用全部的训练样本进行更新，那么每次的损失函数值即为全部样本的误差之和。更新方程不变。
 *  总结   
      正则项是为了降低模型的复杂度，从而避免模型区过分拟合训练数据，包括噪声与异常点（outliers）。从另一个角度上来讲，正则化即是假设模型参数服从先验概率，即为模型参数添加先验，只是不同的正则化方式的先验分布是不一样的。这样就规定了参数的分布，使得模型的复杂度降低（试想一下，限定条件多了，是不是模型的复杂度降低了呢），这样模型对于噪声与异常点的抗干扰性的能力增强，从而提高模型的泛化能力。还有个解释便是，从贝叶斯学派来看：加了先验，在数据少的时候，先验知识可以防止过拟合；从频率学派来看：正则项限定了参数的取值，从而提高了模型的稳定性，而稳定性强的模型不会过拟合，即控制模型空间。   
      另外一个角度，过拟合从直观上理解便是，在对训练数据进行拟合时，需要照顾到每个点，从而使得拟合函数波动性非常大，即方差大。在某些小区间里，函数值的变化性很剧烈，意味着函数在某些小区间里的导数值的绝对值非常大，由于自变量的值在给定的训练数据集中的一定的，因此只有系数足够大，才能保证导数的绝对值足够大。如下图（引用知乎）：   
    ![这里写图片描述][20151026205322532]   
      另外一个解释，规则化项的引入，在训练（最小化cost）的过程中，当某一维的特征所对应的权重过大时，而此时模型的预测和真实数据之间距离很小，通过规则化项就可以使整体的cost取较大的值，从而，在训练的过程中避免了去选择那些某一维（或几维）特征的权重过大的情况，即过分依赖某一维（或几维）的特征（引用知乎）。   
      L2与L1的区别在于，L1正则是拉普拉斯先验，而L2正则则是高斯先验。它们都是服从均值为0，协方差为。当时，即没有先验）没有正则项，则相当于先验分布具有无穷大的协方差，那么这个先验约束则会非常弱，模型为了拟合所有的训练集数据， 参数可以变得任意大从而使得模型不稳定，即方差大而偏差小。越大，标明先验分布协方差越小，偏差越大，模型越稳定。即，加入正则项是在偏差bias与方差variance之间做平衡tradeoff（来自知乎）。下图即为L2与L1正则的区别：   
    ![这里写图片描述][20151026204740369]   
    上图中的模型是线性回归，有两个特征，要优化的参数分别是w1和w2，左图的正则化是L2，右图是L1。蓝色线就是优化过程中遇到的等高线，一圈代表一个目标函数值，圆心就是样本观测值（假设一个样本），半径就是误差值，受限条件就是红色边界（就是正则化那部分），二者相交处，才是最优参数。可见右边的最优参数只可能在坐标轴上，所以就会出现0权重参数，使得模型稀疏。   
      其实拉普拉斯分布与高斯分布是数学家从实验中误差服从什么分布研究中得来的。一般直观上的认识是服从应该服从均值为0的对称分布，并且误差大的频率低，误差小的频率高，因此拉普拉斯使用拉普拉斯分布对误差的分布进行拟合，如下图：   
    ![这里写图片描述][20151026204837282]   
    而拉普拉斯在最高点，即自变量为0处不可导，因为不便于计算，于是高斯在这基础上使用高斯分布对其进行拟合，如下图：   
    ![这里写图片描述][20151026204901473]   
    具体参见：[正态分布的前世今生][Link 1]

### Dropout ###

正则是通过在代价函数后面加上正则项来防止模型过拟合的。而在神经网络中，有一种方法是通过修改神经网络本身结构来实现的，其名为Dropout。该方法是在对网络进行训练时用一种技巧（trick），对于如下所示的三层人工神经网络：   
![这里写图片描述][20151026204937425]   
对于上图所示的网络，在训练开始时，随机得删除一些（可以设定为一半，也可以为1/3，1/4等）隐藏层神经元，即认为这些神经元不存在，同时保持输入层与输出层神经元的个数不变，这样便得到如下的ANN：   
![这里写图片描述][20151026204952319]   
然后按照BP学习算法对ANN中的参数进行学习更新（虚线连接的单元不更新，因为认为这些神经元被临时删除了）。这样一次迭代更新便完成了。下一次迭代中，同样随机删除一些神经元，与上次不一样，做随机选择。这样一直进行瑕疵，直至训练结束。

Dropout方法是通过修改ANN中隐藏层的神经元个数来防止ANN的过拟合。具体可参见[这里][Link 2]。

============================================================

## 整理：过拟合的定义       ##

指在模型参数拟合过程中的问题，由于训练数据包含抽样误差，训练时，复杂的模型将抽样误差也考虑在内，将抽样误差也进行了很好的拟合。具体表现就是训练集上效果好，在测试集上效果差。模型泛化能力弱。

## 解决过拟合的方法 ##

发生过拟合主要有两方面的原因，下面就从数据和模型复杂度来考虑。

### 数据 ###

1.  获取更多的数据：从数据源头获取；根据当前数据集估计数据分布参数，使用分布产生；数据增强，通过一定规则获取数据，例如图像平移，翻转，缩放，切割等       
2.  考虑数据降维，例如PCA，SVD等，选取重要特征，
3.  重新清洗数据，增大数据的训练集

### 模型 ###

1.  使用合适的模型：过拟合的主要原因为，数据太少+模型太复杂。对于神经网络来讲，可以减少网络层数及神经元个数；
2.  训练时间的把控，Early stopping；
3.  Dropout，这是一个很高效的方法，在训练时，每次随机（如50%）忽略隐层的某些节点，这样，我们相当于随机从2^H个模型中采样选择模型，类似于bagging，此外，不同模型之间权值共享，相当于一种权值正则方法，实际效果比L2更好。
4.  限制权值，即正则化（regularization），L1或L2；
5.  增加噪声，噪声随着网络传播，按照权值的平方放大，并传播到输出层，对误差cost产生影响。在输入中加高斯噪声，会达到与L2类似的效果；也可在权值和网络响应上加噪声。
6.  考虑集成的方法，结合多种模型：训练多个模型，以每个模型的平均输出作为结果。如Bagging随机森林，训练一堆互不关联的决策树。Boosting如GBDT，XGBoost。
7.  Batch normalization
8.  决策树的剪枝
9.  通过贝叶斯的相关方法

更多内容：[https://www.zhihu.com/question/59201590/answer/167392763][https_www.zhihu.com_question_59201590_answer_167392763]

[20151026205032330]: /images/20220602/3e75295302ae4a779bf8f2c6b3e004ad.png
[70]: /images/20220602/07019664e1f449ae991a7378402a01a7.png
[20160425090721762]: /images/20220602/f2d0c2b8fa4746808758246e39e98e74.png
[20151026205322532]: /images/20220602/00b6d3538811446eab9d3e0270072fd8.png
[20151026204740369]: /images/20220602/08775f6fb4f94d1f8211ccfd901e2950.png
[20151026204837282]: /images/20220602/6f672b7593d4496980cff88f5bb1d0c3.png
[20151026204901473]: /images/20220602/db59220837ea443c9dd1eef4b1dab46f.png
[Link 1]: http://emma.memect.com/t/05400bfdb821bac2abba403afc4ad2e9d2e7d8ea845aa7227f2433ffbe7a9684/intro-normal-distribution.pdf
[20151026204937425]: /images/20220602/bf2d21e1fa9741a48843f3a3321d93c7.png
[20151026204952319]: /images/20220602/f26a939dcd2a43579e698615a4e2861b.png
[Link 2]: http://papers.nips.cc/paper/4824-imagenet-classification-w
[https_www.zhihu.com_question_59201590_answer_167392763]: https://www.zhihu.com/question/59201590/answer/167392763