跳到主要内容

统计建模

非线性回归

在参数进入函数的方式已经弯曲时,用非线性最小二乘估计米氏饱和或指数衰减,并用迭代分箱查看求解器走了多远。

更新日期 2026-08-29 回归

这篇能回答什么

非线性回归(工作台目录常写作「拟合工具箱」)问的是:浓度升高以后速率会不会饱和、存量会不会按指数掉下去——而不是「再加一分成绩,总分平均涨多少」。均值函数 f(x,θ) 对参数 θ 不能写成 Xβ,要用高斯-牛顿或莱文伯格-马夸尔特去最小化残差平方和。教学表不打印 Wald z,只把求解器记下的步数打进三个互斥迭代分箱。

适合:已知机理曲线。米氏(Michaelis–Menten)速率 v=V_{max}[S]/(K_m+[S]);指数 y=ae^{bx} 或一阶衰减 y=a(1-e^{-kx})。教学占用合计 265 步:早段 A 65(24.53%),中段 B 72(27.17%),近收敛 C 128(48.3%)。C 大只说明多数记录步已经走在平坦的 RSS 面上,不能写成「曲线已经对、K_m 可信」。

什么时候不该用

斜率在参数里仍是一条直线时,用线性回归,不要为「看起来高级」去套指数。要指定分布族和连接、参数仍进线性预测子,用广义线性模型。对多维列联表的对数期望频数建模、没有指定哪一列是 Y,用对数线性模型。学号、年份编码本身不是浓度。

如何选择方法

先问弯曲发生在哪一层。函数形式对 θ 非线性 → 本页 NLS。对 θ 线性、只是 Y 的分布不是正态 → GLM。没有指定响应、只拆格子关联 → 对数线性。同一人多次随访、要总体平均并改工作相关 → 广义估计方程。把米氏式两边取倒数做成 Lineweaver–Burk 直线,会把误差结构拧歪,不能代替本页。

  • 非线性回归:一条定量 X、一条定量 Y、预指定 f;输出拟合轨迹与迭代分箱。

  • 线性回归(OLS):Y≈Xβ,斜率对参数是平的。

  • 广义线性模型:g(μ)=Xβ,弯的是连接,不是 θ 的进入方式。

  • 对数线性模型:对格子期望取对数,对象是列联表而不是浓度–速率点。

输入、输出与限制

输入

自变量是浓度、时间或剂量,因变量是速率、存量或株高,都必须是定量。函数形式从米氏、指数、logistic 生长里选,或按教材写入。缺失按所用列成对删除。

  • 自变量 X:最少 1 个,最多 1 个定量列。

  • 因变量 Y:最少 1 个,最多 1 个定量列。

  • 函数形式:米氏、指数或用户给出的 f(x,θ)。

输出

结果表项目列写「拟合工具箱」。类别 A/B/C 是早/中/晚迭代分箱的步数与占比,合计 265。没有系数、标准误或 Wald 列,不要把 48.3% 读成某个参数的显著性。图是自变量–因变量散点加一条拟合轨迹。

产品限制

  • 占用合计是求解器记下的步,不必等于工作表行数(年报「annual」只有二十来期)。

  • 初值离盆地太远时,莱文伯格阻尼会把步数堆进晚段分箱,看起来「C 很大」其实在原地踏。

  • 标识列、纯日期序号不能当浓度。

  • 散点图是配对视图;完整分箱人数以表的『取值』列为准。

案例

研究问题。 一条年度产值序列看起来往上弯,研究者想用指数或饱和生长去贴,而不是强行画一条过原点的直线;同时先看求解器的步是堆在早段还是已经挤在近收敛分箱。

教学示例下载时间序列表,练习曲线用「annual」工作表:年份对产值。占用合计 265 步,早段 65、中段 72、晚段 128,与二十期年报行数不必相等。若要练「两个定量列配对」,定量成绩表里的数学成绩对物理成绩通常更接近直线,只适合对照,不适合宣称米氏动力学。

下载教学示例表(时间序列)

在 StatsX 中操作

  1. 步骤 1 · 拖入年报或成绩表

    打开工作台后放入 Excel / CSV / JSON。练指数或饱和生长选「annual」;对照直线可选定量成绩表。体积上限 20 MB。

  2. 步骤 2 · 核对数列而不是学号

    画像里确认产值或成绩为定量、按时间或浓度升序。不要把学号、姓名拖进 X。仍在线性爬坡且机理就是直线时,应离开本页。

  3. 步骤 3 · 找到拟合工具箱

    目录搜索「非线性」或「拟合工具箱」。它估计弯曲的 f(x,θ),不是自动选阶的时序按钮。

  4. 步骤 4 · 指定 X、Y 与曲线

    X 放年份或浓度,Y 放产值或速率。函数选指数或米氏,给一个说得通的初值(V_{max} 靠近最大观测速率),再运行。

  5. 步骤 5 · 先看分箱再看轨迹

    早段是否很快降 RSS,晚段是否只是在阻尼里空转。轨迹与表不一致时引用表。

工作台上传时间序列或成绩表
把教学表示例拖进上传区。
列画像核产值或成绩为定量
画像里确认目标列为定量,标识列不要进入曲线。
目录中点选非线性回归或拟合工具箱
在算法目录选定非线性回归(拟合工具箱)。
为非线性回归指定自变量与因变量
X 放浓度或年份,Y 放速率或产值,并指定曲线后再运行。

如何阅读输出

下表是迭代分箱,不是系数表。合计 265 步。A 早段 65(24.53%)对应 RSS 仍在大幅下降的步;B 中段 72(27.17%)对应步长已经缩短;C 晚段 128(48.3%)对应靠近驻点或被阻尼钉住的步。占比分母是记下的步数,不是年数,也不是学生人数。

拟合工具箱结果表
项目类别统计量取值
拟合工具箱类别 A频数65
拟合工具箱类别 B频数72
拟合工具箱类别 C频数128
拟合工具箱合计有效样本 N265
拟合工具箱类别 A占比(%)24.53
拟合工具箱类别 B占比(%)27.17
拟合工具箱类别 C占比(%)48.3

左右滑动可查看完整表格

教学示例表在 StatsX 工作台中的计算结果口径。

非线性回归迭代分箱结果表
工作台里的分箱人数应与上表同一口径。
自变量与因变量散点及拟合轨迹
散点轴标注自变量、因变量;橙线是工作台给出的拟合轨迹。教学图可能看起来接近直线,仍不要把它读成 OLS 斜率表,分箱含义以表为准。

C 独占约一半步(48.3%)很常见:求解器在盆地附近会连记许多小步。这不能翻译成「米氏已经拟合成功」。早段若几乎为空,多半是初值已经贴在终值上,或阻尼一开始就把更新压没了。反过来,早段很大而晚段几乎走不到,要怀疑函数选错或初值飞出定义域(K_m 被推成负数)。

阅读顺序:合计步数是否对得上这次运行 → 三箱占比是否畸形堆在一箱 → 散点有没有饱和或指数弯曲 → 只有在轨迹和机理一致时,才去讨论 V_{max}、K_m 或衰减常数。表头没有 Wald,就不要补一列「显著」。换初值或换指数/米氏,分箱占比可以整份改写。

理论与公式

非线性最小二乘寻找 \hatθ=arg min_θ S(θ),其中

米氏与指数是本页最常用的 f。米氏在高浓度饱和;指数在半对数纸上接近直线,但误差仍在原始尺度上相加。

高斯-牛顿用雅可比 J=∂f/∂θ 做局部线性化;莱文伯格-马夸尔特在法矩阵对角加上阻尼 λ,λ 大时步子像梯度下降。教学分箱按「这一步 RSS 相对下降」切开:下降猛进早段,下降微弱进晚段。渐近协方差要在收敛以后用法矩阵的逆再乘残差方差另算,本页占用表不提供那一列。Lineweaver–Burk 把米氏改写成 1/v 对 1/s 的直线,低速率点会被放大,因此本页坚持在原始尺度上做 NLS。

注意事项

  • 晚段分箱大 ≠ 收敛成功;看 RSS 是否还在降,以及参数有没有撞到边界。

  • 初值决定你掉进哪一口井;V_{max} 请从最大观测速率附近起,不要从 0 起。

  • 年份当 X 时,最好先平移到从 0 起,否则指数里的 b 会病态。

  • 本页不管聚类;同一株作物反复量株高,应改 GEE 或混合模型。

  • 图上的直线感不能反证「所以该用 OLS」——先看机理,再看分箱。

常见问题

参考文献

  • Bates, D. M., Watts, D. G. Nonlinear Regression Analysis and Its Applications. Wiley.

  • Seber, G. A. F., Wild, C. J. Nonlinear Regression. Wiley.

  • Michaelis, L., Menten, M. L. Die Kinetik der Invertinwirkung. Biochemische Zeitschrift, 1913.

  • 何晓群, 刘文卿. 应用回归分析[M]. 中国人民大学出版社.