try ai
科普
编辑
分享
反馈
  • 最小二乘法

最小二乘法

SciencePedia玻尔百科
核心要点
  • 最小二乘法通过最小化每个数据点到直线的垂直距离(残差)的平方总和来确定数据集的最佳拟合线。
  • 最小二乘线的一个基本特性是它必然穿过数据的质心(x 和 y 的平均值点)。
  • 决定系数 R2R^2R2 通过衡量模型能够解释的因变量方差的比例来量化拟合优度。
  • 该方法的一个主要弱点是对离群点的高度敏感性,因为误差的平方会使较大的偏差对最终结果产生不成比例的影响。

引言

在几乎所有的科学和工程领域,数据很少能讲述一个简单明了的故事。相反,它通常以图表上散乱的点云形式出现,暗示着某种趋势,但被噪声和随机变异所掩盖。根本的挑战在于穿透这片混乱,找到那个单一的、潜在的关系。最小二乘法为这个问题提供了一个强大且在数学上优雅的解决方案,为确定一组数据的“最佳拟合”线提供了一种明确的方法。它是统计分析的基石,使我们能够建立关系模型、做出预测,并将杂乱的观测数据转化为清晰、可操作的知识。

本文深入探讨了这项基础技术,探索其理论基础和广泛的实际效用。在接下来的章节中,您将发现使该方法如此有效的核心思想。关于​​原理与机制​​的章节将解析最小化平方误差的数学原理,揭示所得直线的优雅几何特性,并解释如何量化拟合质量。随后,关于​​应用与跨学科联系​​的章节将展示该方法如何应用于从化学、生物学到材料科学和工程等不同领域,以解决现实世界的问题,展示其非凡的灵活性和强大功能。

原理与机制

想象一下,你是一位站在河边的环境科学家,正在看一张散点图。图上的每个点代表一个地点,将污染物浓度与某种鱼类的种群数量配对。你看到了一个趋势——一个向下倾斜的点云——但它很杂乱。你如何穿过这片点云画出那条唯一的“最佳”直线,以捕捉这种关系的本质?这是最小二乘法诞生之初就要回答的核心问题。

“最佳”的衡量标准:最小化平方误差

“最佳”到底是什么意思?我们的直觉可能会认为是一条穿过这些点“中间”的线。最小二乘法使这个想法变得精确。对于你可能画出的任何一条线,你都可以测量每个数据点的“误差”。这并非你测量中的失误,而是你的观测值与该线预测值之间的差异。出于非常充分的理由,惯例是将这个误差测量为从观测数据点 (xi,yi)(x_i, y_i)(xi​,yi​) 到线上具有相同 xxx 坐标的点 (xi,y^i)(x_i, \hat{y}_i)(xi​,y^​i​) 的​​垂直距离​​。这个垂直差距 yi−y^iy_i - \hat{y}_iyi​−y^​i​ 被称为​​残差​​。

为什么是垂直距离?因为在许多实验中,比如研究污染物影响的实验,我们认为 xxx 变量(污染物浓度)是可以高精度控制或观测的,而 yyy 变量(鱼类种群)是带有某种随机性或“噪声”的结果。我们试图在给定 xxx 的情况下预测 yyy,所以我们关心的是 yyy 方向上的误差。

现在,我们有了一组这样的垂直误差线,有些是正的(点在线上方),有些是负的(点在线下方)。我们可以直接将它们相加,但正负误差会相互抵消,得出一个具有误导性的小总和。我们也可以将它们的绝对值相加,但这会在数学上带来一些麻烦。

独立发展出此方法的 Carl Friedrich Gauss 和 Adrien-Marie Legendre 的天才之处在于,在将每个误差相加之前先将其平方。这个简单的举动带来了深远的影响。我们试图最小化的量是​​平方误差和​​(Sum of Squared Errors, SSE),有时也称为残差平方和:

S=∑i=1n(yi−y^i)2=∑i=1n(yi−(mxi+b))2S = \sum_{i=1}^{n} (y_i - \hat{y}_i)^2 = \sum_{i=1}^{n} (y_i - (mx_i + b))^2S=∑i=1n​(yi​−y^​i​)2=∑i=1n​(yi​−(mxi​+b))2

平方误差完美地实现了两件事:它使所有误差都变为正数,因此不会相互抵消;并且它对远离直线的点给予了更大的惩罚。一个距离远一倍的点对平方和的贡献是四倍。最小二乘线是使这个总平方和尽可能小的唯一一条线。

寻找碗底

所以,我们的目标是:找到使函数 S(m,b)S(m, b)S(m,b) 最小化的斜率 mmm 和截距 bbb。我们该怎么做呢?想象函数 SSS 是一个地貌。由于它是平方和,它在三维空间中形成一个光滑、向上弯曲的碗状曲面,其中两个水平方向代表 mmm 和 bbb 的值,垂直方向是 SSS 的值。找到“最佳拟合线”就等同于找到这个碗状曲面最底部的坐标 (m,b)(m, b)(m,b)。

我们如何找到碗底呢?我们找到那个曲面完全平坦的唯一一点。用微积分的语言来说,这意味着找到 SSS 对 mmm 和 bbb 的偏导数都为零的地方。

∂S∂m=0and∂S∂b=0\frac{\partial S}{\partial m} = 0 \quad \text{and} \quad \frac{\partial S}{\partial b} = 0∂m∂S​=0and∂b∂S​=0

通过微积分推导——这是任何统计学学生的必经之路——可以得到关于 mmm 和 bbb 的一对联立线性方程组,称为​​正规方程​​。对于一组数据点,我们可以计算所有必要的总和(如 ∑xi\sum x_i∑xi​, ∑yi\sum y_i∑yi​, ∑xi2\sum x_i^2∑xi2​, 和 ∑xiyi\sum x_i y_i∑xi​yi​),将它们代入正规方程,然后解出最小化误差的唯一一对 (m,b)(m, b)(m,b)。这正是一位材料科学专业的学生从力-伸长数据中找出一种新型聚合物纤维的刚度(mmm)和初始伸长(bbb)所遵循的程序。

最小化带来的优雅推论

这个最小化过程不仅给了我们一条线,它还赋予了这条线一些卓越而优美的特性。

首先,直接由 ∂S∂b=0\frac{\partial S}{\partial b} = 0∂b∂S​=0 导出的一个正规方程可以简化为一个绝妙的结论:∑(yi−(mxi+b))=0\sum (y_i - (mx_i + b)) = 0∑(yi​−(mxi​+b))=0。换句话说,最小二乘线的全部残差之和​​永远精确地为零​​。正误差(在线上方的点)和负误差(在线下方的点)完美地相互抵消。这意味着,如果一位物理学家在不同伸长量下测量弹簧的力,那么测量力与最佳拟合线预测值之间的差值之和将为零。

其次,同一个方程可以重新整理,证明 yˉ=mxˉ+b\bar{y} = m\bar{x} + byˉ​=mxˉ+b,其中 xˉ\bar{x}xˉ 和 yˉ\bar{y}yˉ​ 分别是 xxx 和 yyy 数据的平均值。这证明了一个绝佳的几何事实:最小二乘回归线​​必然穿过​​数据的​​质心​​,即点 (xˉ,yˉ)(\bar{x}, \bar{y})(xˉ,yˉ​)。在非常真实的意义上,这条线是数据云的枢轴点或“重心”。

分解变异:拟合效果好吗?

找到最佳拟合线是一回事;知道它是否是一个好的拟合是另一回事。我们的线究竟揭示了数据中多少信息?关键在于考察变异。

想象一下,你没有一条拟合线,有人让你预测一个新观测值的 yyy 值。你最好的猜测就是你所见过的所有 yyy 值的平均值 yˉ\bar{y}yˉ​。数据中的总变异可以看作是每个观测值 yiy_iyi​ 与这个平均值之差的平方和,这个量被称为​​总平方和 (SST)​​。

SST=∑i=1n(yi−yˉ)2\text{SST} = \sum_{i=1}^{n} (y_i - \bar{y})^2SST=∑i=1n​(yi​−yˉ​)2

最小二乘法的奇妙之处在于,它将这个总变异分解为两个有意义的部分。第一部分是由我们的回归线“解释”的变异。这是回归线的预测值 y^i\hat{y}_iy^​i​ 与总体均值 yˉ\bar{y}yˉ​ 之差的平方和。这被称为​​回归平方和 (SSR)​​。

SSR=∑i=1n(y^i−yˉ)2\text{SSR} = \sum_{i=1}^{n} (\hat{y}_i - \bar{y})^2SSR=∑i=1n​(y^​i​−yˉ​)2

第二部分是线无法“解释”的变异——即剩余的误差。这正是我们一开始就最小化的平方误差和 (SSE)。

SSE=∑i=1n(yi−y^i)2\text{SSE} = \sum_{i=1}^{n} (y_i - \hat{y}_i)^2SSE=∑i=1n​(yi​−y^​i​)2

事实证明,这三个量通过一个极其简洁的恒等式联系在一起,该恒等式构成了方差分析 (ANOVA) 的基石:

SST=SSR+SSE\text{SST} = \text{SSR} + \text{SSE}SST=SSR+SSE

总变异 = 已解释变异 + 未解释变异。这个强大的方程使我们能够量化拟合优度。比率 SSR/SST\text{SSR}/\text{SST}SSR/SST,即我们所知的 R2R^2R2,告诉我们 yyy 的总方差中可以由 xxx 预测的比例。R2R^2R2 为 0.9 意味着鱼类种群数量变化的 90% 可以由污染物浓度来解释。

挑战极限:完美拟合与不可能的拟合

当我们将该方法推向极限时,会发生什么?

  • ​​完美拟合:​​ 假设我们只有两个不同的数据点。我们的直觉告诉我们,最佳拟合线就是穿过这两个点的直线。最小二乘机制完美地证实了这一点。当你用两个点 (x1,y1)(x_1, y_1)(x1​,y1​) 和 (x2,y2)(x_2, y_2)(x2​,y2​) 进行正规方程的计算时,结果恰好是我们熟悉的斜率公式 m=y2−y1x2−x1m = \frac{y_2 - y_1}{x_2 - x_1}m=x2​−x1​y2​−y1​​ 和相应的截距。在这种情况下,SSE 为零,因为这条线精确地穿过了两个点。

  • ​​过拟合:​​ 这个想法可以进一步延伸。如果你有 nnn 个数据点(具有不同的 xxx 值),数学上可以证明,你总能找到一个唯一的 n−1n-1n−1 次多项式,精确地穿过每一个点。如果你用最小二乘法来拟合这样一个多项式,算法会找到它,并且平方误差和将恰好为零。这听起来很棒,但这是一个陷阱!这个“完美”的模型只是记住了数据,包括其随机噪声。在预测新数据点时,它的表现可能会非常糟糕。这是机器学习中一个被称为​​过拟合​​的关键概念。

  • ​​不可能的拟合:​​ 如果分析师犯了一个错误,将所有化学标准品都配制在完全相同的浓度下,比如说 x=3.0x=3.0x=3.0?数据点将在图上形成一条垂直线。现在的“最佳拟合”线是什么?单一斜率的概念本身就变得毫无意义了。最小二乘法反映了这种模糊性。正规方程变得线性相关,不再提供唯一的 (m,b)(m, b)(m,b) 解,而是提供了满足简单关系 3m+b=yˉ3m+b = \bar{y}3m+b=yˉ​ 的无限多个解。所有这些线都围绕点 (3.0,yˉ)(3.0, \bar{y})(3.0,yˉ​) 旋转,并且每一条线都产生完全相同的最小平方误差和。数学并没有崩溃;它正确地告诉我们这个问题是病态的(ill-posed)。

阿喀琉斯之踵:对离群点的敏感性

最小二乘法最大的优点——源自误差平方的简洁数学特性——同时也是其最大的弱点。考虑一个工程师正在确定一个组件的热阻。大多数测量都是准确的,但由于瞬时传感器故障,有一个读数严重偏离。因为最小二乘法极度“厌恶”大误差(记住,它会对其进行平方),它会急剧改变直线的斜率,使其远离真实关系,只为了减小那一个巨大的平方误差。单个离群点就像一个引力恶霸,施加不成比例的影响,严重扭曲最终结果。在任何真实世界的应用中,数据很少是完美的,记住这种敏感性至关重要。虽然存在更稳健的方法,但它们牺牲了最小二乘法的数学优雅性。

从最小化垂直线段平方这一简单前提出发,最小二乘法发展成为一个丰富而强大的理论体系,它拥有优雅的几何特性、一套自我评估的框架以及清晰的局限性。它之所以成为科学和工程的基石,不仅因为它行之有效,更因为其原理揭示了我们将杂乱数据转化为清晰知识的过程中所蕴含的深刻而优美的结构。

应用与跨学科联系

既然我们已经熟悉了最小二乘法的数学核心,我们就可以开始真正的冒险:看它如何付诸实践。学习原理就像学习一门新语言的语法;真正的乐趣来自于阅读它所能讲述的诗篇和故事。而这些故事何其精彩!最小二乘法是一条贯穿人类几乎所有量化研究领域的线索,从平凡到宏伟。它是一种与自然“对话”并探寻其秘密的通用工具。

最佳猜测的艺术:寻找趋势

在最基本的层面上,最小二乘法是我们穿过一堆散乱数据点绘制直线的最佳数学工具。这听起来可能很简单,但却是经验科学的基石。想象一下,你是一位分析化学家,试图测量一座老房子油漆碎片中的铅浓度。你使用一种叫做原子吸收光谱法的技术,其中样品吸收的光量与其中元素的浓度成正比。问题是,你的仪器并不完美;每次测量都带有少量不可避免的“噪声”。

为了确定你的未知浓度,你首先准备几个已知浓度的样品并测量它们的吸光度。你将结果绘制成图,看到一个总体趋势——铅越多,吸光度越高——但这些点并不在一条完美的直线上。真实的关系是什么?最小二乘法给了我们最民主的答案:它画出一条唯一的直线,最小化到所有数据点的总垂直距离平方和,为我们提供了一把坚固的校准“尺子”,用来测量我们的未知样品。同样的想法也让分析师能够找到每日温度和冰淇淋销量之间的最佳拟合关系,为根据现实世界数据预测商业结果提供了一个简单而强大的模型。

超越直线:变换的力量

当然,自然界并不总是那么直截了当。许多关系并非线性。这是否意味着我们闪亮的新工具毫无用处?完全不是!这时就需要一点科学创造力了。通常,一个非线性关系可以通过从不同角度审视它来“变换”成线性关系。

考虑物理学中的波义耳定律,该定律指出,对于恒温下的气体,压力 PPP 与体积 VVV 成反比,即 P=k/VP = k/VP=k/V。如果你绘制 PPP 对 VVV 的图,你会得到一条曲线。然而,如果你巧妙地决定绘制 PPP 对体积的倒数 x=1/Vx = 1/Vx=1/V 的图,关系就变成了 P=kxP = kxP=kx,一条穿过原点的完美直线!然后我们可以使用最小二乘法从实验数据中找到常数 kkk 的最佳拟合值。

这种强大的线性化思想无处不在。在生物学中,种群通常呈指数增长,遵循像 y=αexp⁡(βx)y = \alpha \exp(\beta x)y=αexp(βx) 这样的模型。这是另一条曲线。但是,如果我们对两边取自然对数,我们得到 ln⁡(y)=ln⁡(α)+βx\ln(y) = \ln(\alpha) + \beta xln(y)=ln(α)+βx。突然之间,这成了一个关联 ln⁡(y)\ln(y)ln(y) 和 xxx 的线性方程。我们可以再次引入最小二乘机制来估计参数,将一个困难的非线性问题转变为我们熟悉的直线拟合练习。该方法本身具有极好的灵活性;它既能为 x=my+bx = my + bx=my+b 找到最佳拟合线,也能为更常见的 y=mx+by = mx + by=mx+b 找到最佳拟合线,这完全取决于我们希望预测哪个变量。

不完美的智慧:从噪声中提取信号

一个关键问题可能一直困扰着你:为什么不直接找到一条精确穿过每个数据点的曲线呢?完美的拟合难道不比近似的拟合更好吗?答案是响亮的不,理解这一点是理解最小二乘法深层智慧的关键。

想象一下,你有一组关于平滑变化现象的带噪声的测量数据。你可以使用一个高次多项式来“连接这些点”,迫使你的曲线完美地穿过每个点。但这样做,你不仅拟合了潜在的信号,还拟合了随机、无意义的噪声。这样的曲线通常会在数据点之间表现出剧烈的振荡,这种行为被称为龙格现象。它对数据的特质学习得太好,结果成了一个糟糕的真实潜在过程预测器。

相比之下,一个简单的最小二乘拟合——比如用一个低次多项式——并不试图做到完美。它接受每个数据点都略有偏差。它找到一条平滑的曲线,穿行于这些点之间,扮演一个明智的仲裁者,平衡数据的矛盾证词,以找到最可信的趋势。这就是偏差-方差权衡的本质:最小二乘拟合引入了少量的“偏差”(它不完美匹配数据),以换取“方差”的大幅减少(它更稳定,对任何单个点的噪声都不那么敏感)。它是在噪声中分离信号的大师。

规则的变通:方法的推广

标准的最小二乘法建立在几个关键假设之上,其中最重要的一条是数据点的误差是相互独立的。当这个假设不成立时会发生什么?这正是该方法展示其真正力量和适应性的地方。

在进化生物学中,研究人员经常比较不同物种的性状。例如,他们可能会问哺乳动物的体重是否与奔跑速度相关。一种简单的方法是收集 80 个物种的数据并进行标准的最小二乘回归。问题在于,这些数据点并非独立的。狮子和老虎共享一个近期的共同祖先,因此它们彼此之间的相似性比它们与树懒的相似性要高。标准最小二乘法对这种共享历史视而不见,很容易被愚弄,从而找到一些仅仅是进化谱系树伪影的虚假相关性。

解决方案是一种优美的扩展方法,称为​​系统发育广义最小二乘法(PGLS)​​。该方法通过将系统发育树——物种的“家谱”——直接整合到数学计算中来修改标准程序。它实质上是告诉算法:“当亲缘关系近的物种相似时,不要太惊讶。”非独立性不再是一个需要忽略的问题,而是需要利用的关键信息,从而得出更可靠的科学结论。

同样,在控制理论和系统辨识中,一个关键假设是测量噪声与预测变量不相关。在为一个动态系统建模时,如果一个时间步的输出依赖于前一个时间步的输出,而噪声本身是序列相关的(“有色噪声”),那么这个假设就可能失效。这种相关性会系统性地偏置参数估计,有时会带来灾难性的后果——例如,导致工程师根据模型拟合结果,将一个物理上稳定的过程判断为数学上不稳定。这凸显了一个重要的教训:最小二乘法是一个强大的工具,但像任何工具一样,必须在深刻理解其操作假设的情况下使用。

现代前沿:逆问题的引擎

在其最先进的形式中,最小二乘框架成为解决所谓​​逆问题​​的强大引擎——这相当于科学界的侦探工作。在这种情况下,我们观察到结果,必须反向推断原因。

一个惊人的例子来自现代材料科学。当材料开裂时,裂纹尖端周围的应力和应变由一套源自断裂力学的复杂方程组描述。实验人员可以使用一种称为数字图像相关法(DIC)的技术来捕获裂纹扩展周围位移场的高分辨率图——即材料变形的精确方式。逆问题是:给定这张包含数千个位移数据点的图,导致这种变形的潜在应力强度因子(KIK_IKI​ 和 KIIK_{II}KII​)必须是什么?最小二乘法提供了解决这个问题的机制。它找到 KIK_IKI​ 和 KIIK_{II}KII​ 的值,当这些值代入理论方程时,生成的位移场与实验中测量的位移场最佳匹配。这是理论与数据的卓越结合,使我们能够测量那些原本看不见的材料基本属性。

更进一步,考虑增材制造(3D 打印)的挑战。随着打印部件的冷却,内部应力——或称“固有应变”——会产生,常常导致最终部件翘曲。如果我们能测量最终的翘曲形状,能否推断出导致这种变形的固有应变模式?这个逆问题通常是“病态的”,意味着许多不同的应变模式可能导致相似的形状,并且少量的测量噪声可能导致截然不同、不符合物理规律的解。

在这里,最小二乘原理被扩展为所谓的​​吉洪诺夫正则化​​。我们设计一个目标函数来同时最小化两件事:熟悉的预测形状与测量形状之间的平方误差和,以及一个量化解的物理不合理性的第二个“惩罚”项。例如,我们可能会惩罚那些在层与层之间剧烈振荡的应变模式。最终的解是一个平衡:它是在与观测数据一致的前提下,物理上最合理的解释。这不再仅仅是曲线拟合;这是一个将物理直觉直接编码到我们数据分析中的框架。

从在图上画线到重建制造过程中隐藏的物理原理,最小二乘法的发展历程证明了其深远的实用性。它是一项统一的原则,为在整个科学和工程领域内处理不确定数据提供了一种通用语言,并不断提醒我们,有时,“最佳猜测”就是我们能找到的最强大的答案。