作者 | Prince Grover

译者 | 张建军

整理 | 阿司匹林、Jane

出品 | AI科技大本营

【导读】没有一个损失函数可以适用于所有类型的数据。损失函数的选择取决于许多因素，包括是否有离群点，机器学习算法的选择，运行梯度下降的时间效率，是否易于找到函数的导数，以及预测结果的置信度。这篇文章的目的就是希望帮助大家了解不同的损失函数。

机器学习中的所有算法都依赖于最小化或最大化某一个函数，我们称之为“目标函数”。最小化的这组函数被称为“损失函数”。损失函数是衡量预测模型预测期望结果表现的指标。寻找函数最小值的最常用方法是“梯度下降”。把损失函数想象成起伏的山脉，梯度下降就像从山顶滑下，目的是到达山脉的最低点。

损失函数可以大致分为两类：分类损失（Classification Loss）和回归损失（Regression Loss）。在这篇文章中，将重点放在回归损失上。此外，文中出现的所有代码和图可以在下面这个 notebook 中找到。

参阅地址：

https://nbviewer.jupyter.org/github/groverpr/Machine-Learning/blob/master/notebooks/05_Loss_Functions.ipynb

回归函数预测实数值，分类函数预测标签

1.均方误差，二次损失，L2 损失

（Mean Square Error, Quadratic Loss, L2 Loss）

均方误差（MSE）是最常用的回归损失函数。MSE 是目标变量与预测值之间距离平方之和。

下面是一个MSE函数的图，其中真实目标值为 100，预测值在 -10,000 至 10,000之间。预测值（X轴）= 100 时，MSE 损失（Y轴）达到其最小值。损失范围为 0 至 ∞。

MSE损失（Y轴）与预测值（X轴）关系图

2.平均绝对误差，L1损失

（Mean Absolute Error, L1 Loss）

平均绝对误差（MAE）是另一种用于回归模型的损失函数。MAE是目标变量和预测变量之间差异绝对值之和。因此，它在一组预测中衡量误差的平均大小，而不考虑误差的方向。（如果我们也考虑方向，那将被称为平均偏差（Mean Bias Error, MBE），它是残差或误差之和）。损失范围也是 0 到 ∞。

MAE损失（Y轴）与预测值（X轴）关系图

▌MSE vs MAE （L2损失 vs L1损失）

使用平方误差更容易求解，但使用绝对误差对离群点更加鲁棒。但是，知其然更要知其所以然！

每当我们训练机器学习模型时，我们的目标就是找到最小化损失函数的点。当然，当预测值正好等于真实值时，这两个损失函数都达到最小值。

下面让我们快速过一遍两个损失函数的 Python 代码。我们可以编写自己的函数或使用 sklearn 的内置度量函数：

让我们来看看两个例子的 MAE 值和 RMSE 值（RMSE，Root Mean Square Error，均方根误差，它只是 MSE 的平方根，使其与 MAE 的数值范围相同）。在第一个例子中，预测值接近真实值，观测值之间误差的方差较小。第二个例子中，有一个异常观测值，误差很高。

左：误差彼此接近  右：有一个误差和其他误差相差很远

▌我们从中观察到什么？我们该如何选择使用哪种损失函数？

由于 MSE 对误差（e）进行平方操作（y - y_predicted = e），如果 e> 1，误差的值会增加很多。如果我们的数据中有一个离群点，e 的值将会很高，将会远远大于 |e|。这将使得和以 MAE 为损失的模型相比，以 MSE 为损失的模型会赋予更高的权重给离群点。在上面的第二个例子中，以 RMSE 为损失的模型将被调整以最小化这个离群数据点，但是却是以牺牲其他正常数据点的预测效果为代价，这最终会降低模型的整体性能。

MAE  损失适用于训练数据被离群点损坏的时候（即，在训练数据而非测试数据中，我们错误地获得了不切实际的过大正值或负值）。

直观来说，我们可以像这样考虑：对所有的观测数据，如果我们只给一个预测结果来最小化 MSE，那么该预测值应该是所有目标值的均值。但是如果我们试图最小化 MAE，那么这个预测就是所有目标值的中位数。我们知道中位数对于离群点比平均值更鲁棒，这使得 MAE 比 MSE 更加鲁棒。

使用 MAE 损失（特别是对于神经网络）的一个大问题是它的梯度始终是相同的，这意味着即使对于小的损失值，其梯度也是大的。这对模型的学习可不好。为了解决这个问题，我们可以使用随着接近最小值而减小的动态学习率。MSE 在这种情况下的表现很好，即使采用固定的学习率也会收敛。MSE 损失的梯度在损失值较高时会比较大，随着损失接近 0 时而下降，从而使其在训练结束时更加精确（参见下图）。

▌决定使用哪种损失函数？

如果离群点是会影响业务、而且是应该被检测到的异常值，那么我们应该使用 MSE。另一方面，如果我们认为离群点仅仅代表数据损坏，那么我们应该选择 MAE 作为损失。

推荐阅读：下面这篇文章中有一项很好的研究，比较了在存在和不存在离群点的情况下使用 L1损失和 L2损失的回归模型的性能。请记住，L1 和 L2 损失分别是 MAE 和 MSE 的另一个名称而已。

参阅地址：

http://rishy.github.io/ml/2015/07/28/l1-vs-l2-loss/

L1损失对异常值更加稳健，但其导数并不连续，因此求解效率很低。L2损失对异常值敏感，但给出了更稳定的闭式解（closed form solution）（通过将其导数设置为0）

两种损失函数的问题：可能会出现这样的情况，即任何一种损失函数都不能给出理想的预测。例如，如果我们数据中 90％ 的观测数据的真实目标值是 150，其余 10％ 的真实目标值在 0-30 之间。那么，一个以 MAE 为损失的模型可能对所有观测数据都预测为 150，而忽略 10％ 的离群情况，因为它会尝试去接近中值。同样地，以 MSE 为损失的模型会给出许多范围在 0 到 30 的预测，因为它被离群点弄糊涂了。这两种结果在许多业务中都是不可取的。

在这种情况下怎么做？一个简单的解决办法是转换目标变量。另一种方法是尝试不同的损失函数。这是我们的第三个损失函数——Huber Loss——被提出的动机。

3.Huber Loss，平滑的平均绝对误差

Huber Loss 对数据离群点的敏感度低于平方误差损失。它在0处也可导。基本上它是绝对误差，当误差很小时，误差是二次形式的。误差何时需要变成二次形式取决于一个超参数，(delta)，该超参数可以进行微调。当  ? ~ 0时， Huber Loss 接近 MAE，当  ? ~ ∞（很大的数）时，Huber Loss 接近 MSE。

Huber Loss（Y轴）与预测值（X轴）关系图。真值= 0

delta 的选择非常重要，因为它决定了你认为什么数据是离群点。大于 delta 的残差用 L1 最小化（对较大的离群点较不敏感），而小于 delta 的残差则可以“很合适地”用 L2 最小化。

▌为什么使用Huber Loss？

使用 MAE 训练神经网络的一个大问题是经常会遇到很大的梯度，使用梯度下降时可能导致训练结束时错过最小值。对于 MSE，梯度会随着损失接近最小值而降低，从而使其更加精确。

在这种情况下，Huber Loss 可能会非常有用，因为它会使最小值附近弯曲，从而降低梯度。另外它比 MSE 对异常值更鲁棒。因此，它结合了 MSE 和 MAE 的优良特性。但是，Huber Loss 的问题是我们可能需要迭代地训练超参数delta。

4. Log-Cosh Loss

Log-cosh 是用于回归任务的另一种损失函数，它比 L2 更加平滑。Log-cosh 是预测误差的双曲余弦的对数。

Log-cosh Loss（Y轴）与预测值（X轴）关系图。真值= 0

优点： log(cosh(x)) 对于小的 x 来说，其大约等于 (x ** 2) / 2，而对于大的 x 来说，其大约等于 abs(x) - log(2)。这意味着 'logcosh' 的作用大部分与均方误差一样，但不会受到偶尔出现的极端不正确预测的强烈影响。它具有Huber Loss 的所有优点，和 Huber Loss 不同之处在于，其处处二次可导。

为什么我们需要二阶导数？许多机器学习模型的实现（如XGBoost）使用牛顿方法来寻找最优解，这就是为什么需要二阶导数（Hessian）的原因。对于像 XGBoost 这样的机器学习框架，二阶可导函数更有利。

XGBoost中使用的目标函数。注意其对一阶和二阶导数的依赖性。

但 Log-chsh Loss 并不完美。它仍然存在梯度和 Hessian 问题，对于误差很大的预测，其梯度和 hessian 是恒定的。因此会导致 XGBoost 中没有分裂。

Huber 和 Log-cosh 损失函数的 Python 代码：

5.Quantile Loss（分位数损失）

在大多数真实预测问题中，我们通常想了解我们预测的不确定性。了解预测值的范围而不仅仅是单一的预测点可以显着改善许多业务问题的决策过程。

当我们有兴趣预测一个区间而不仅仅是预测一个点时，Quantile Loss 函数就很有用。最小二乘回归的预测区间是基于这样一个假设：残差（y - y_hat）在独立变量的值之间具有不变的方差。我们不能相信线性回归模型，因为它违反了这一假设。当然，我们也不能仅仅认为这种情况一般使用非线性函数或基于树的模型就可以更好地建模，而简单地抛弃拟合线性回归模型作为基线的想法。这时，Quantile Loss 就派上用场了。因为基于 Quantile Loss 的回归模型可以提供合理的预测区间，即使是对于具有非常数方差或非正态分布的残差亦是如此。

让我们看一个有效的例子，以更好地理解为什么基于 Quantile Loss 的回归模型对异方差数据表现良好。

▌Quantile 回归 vs 普通最小二乘（Ordinary Least Square, OLS）回归

左：线性关系b/w X1和Y，残差的方差恒定。右：线性关系b/w X2和Y，但Y的方差随着X2增加而变大（异方差）。

橙线表示两种情况下的OLS估计

Quantile回归：虚线表示基于0.05和0.95 分位数损失函数的回归估计

如上所示的 Quantile 回归代码在下面这个 notebook 中。

参阅地址：

https://github.com/groverpr/Machine-Learning/blob/master/notebooks/09_Quantile_Regression.ipynb）

▌了解 Quantile Loss 函数

基于 Quantile 回归的目的是，在给定预测变量的某些值时，估计因变量的条件“分位数”。Quantile Loss 实际上只是 MAE 的扩展形式（当分位数是第 50 个百分位时，Quantile Loss 退化为 MAE）。

Quantile Loss 的思想是根据我们是打算给正误差还是负误差更多的值来选择分位数数值。损失函数根据所选 quantile (γ) 的值对高估和低估的预测值给予不同的惩罚值。举个例子，γ= 0.25 的 Quantile Loss 函数给高估的预测值更多的惩罚，并试图使预测值略低于中位数。

γ 是给定的分位数，其值介于0和1之间

Quantile Loss（Y轴）与预测值（X轴）关系图。真值为Y= 0

我们也可以使用这个损失函数来计算神经网络或基于树的模型的预测区间。下图是sklearn实现的梯度提升树回归。

使用quantile loss的预测区间（梯度提升回归）

http://scikit-learn.org/stable/auto_examples/ensemble/plot_gradient_boosting_quantile.html

上图显示的是 sklearn 库的 GradientBoostingRegression 中的 quantile loss 函数计算的 90％ 预测区间。上限的计算使用了 γ = 0.95，下限则是使用了 γ = 0.05。

比较研究

“Gradient boosting machines, a tutorial” 中提供了一个很好的比较研究。为了演示上述所有的损失函数的性质，研究人员创造了一个人工数据集，数据集从 sinc(x) 函数中采样，其中加入了两种人造模拟噪声：高斯噪声分量和脉冲噪声分量。脉冲噪声项是用来展示结果的鲁棒效果的。以下是使用不同损失函数来拟合 GBM（Gradient Boosting Machine, 梯度提升回归）的结果。

连续损失函数：（A）MSE损失函数; （B）MAE损失函数; （C）Huber损失函数; （D）Quantile损失函数。用有噪声的 sinc(x)数据来拟合平滑GBM的示例：（E）原始 sinc(x)函数; （F）以MSE和MAE为损失拟合的平滑GBM; （G）以Huber Loss拟合的平滑GBM， = {4,2,1}; （H）以Quantile Loss拟合的平滑GBM，

模拟实验中的一些观察结果

• 以 MAE 为损失的模型预测较少受到脉冲噪声的影响，而以 MSE 为损失的模型的预测由于脉冲噪声造成的数据偏离而略有偏差。

• 以 Huber Loss 为损失函数的模型，其预测对所选的超参数不太敏感。

• Quantile Loss 对相应的置信水平给出了很好的估计。

一张图中画出所有损失函数

原文链接：

https://heartbeat.fritz.ai/5-regression-loss-functions-all-machine-learners-should-know-4fb140e9d4b0

--【完】--

寻找2018最佳AI应用案例

AI科技大本营希望找到在汽车、金融、教育、医疗、安防、零售、家居、文娱、工业等 9 大行业的最佳 AI 应用案例，记录 AI 时代影响人类发展的变革性产品/解决方案。

如果您有优秀的 AI 产品/技术解决方案，欢迎【扫码提交】，参与评选。

点击「阅读原文」，查看案例分享者特别奖励