第7章 · 7.3

7.3 不确定度与方法验证

Measurement Uncertainty and Method Validation
本页范围与学习目标
补充平均、标准差和协方差的统计基础,以内标质量分数建立不确定度预算,区分重复性与偏差,说明方法验证和覆盖区间的适用条件。
先读图解和核心解释;需要细节时再展开。
先看图 · 阅读入口
打开本节交互图

误差预算先找结果依赖哪些输入

不确定度描述结果范围,预算能帮助选择最有效的改进步骤。

误差预算先找结果依赖哪些输入 步骤 1 列出测量输入 面积比、纯度和称量 标明单位与来源 步骤 2 传播不确定度 考虑灵敏度和相关性 不要重复计入同一误差 步骤 3 报告结果范围 给出覆盖因子与条件 说明主要贡献来源 误差预算先找结果依赖哪些输入 步骤 1 列出测量输入 面积比、纯度和称量 标明单位与来源 步骤 2 传播不确定度 考虑灵敏度和相关性 不要重复计入同一误差 步骤 3 报告结果范围 给出覆盖因子与条件 说明主要贡献来源
图 7.3-1 本页概念关系图。箭头表示本页采用的检查顺序。

读图时保留这项条件:标准不确定度和扩展不确定度不同,相关输入不能总按独立处理。

动手看一看 · 教学模型

改善一项输入,观察总预算减少多少

分别拖动四项相对标准不确定度,看平方后贡献与合成结果。

改善一项输入,观察总预算减少多少先平方,再比较方差贡献面积比内标称量内标纯度样品称量
各项相对方差贡献
各项相对方差贡献份额00.250.50.751面积比内标称量内标纯度样品称量输入项
合成相对标准不确定度0.3873%
绝对标准不确定度0.3468 个百分点
扩展不确定度,k=20.6937 个百分点
0.3 %
0.1 %
0.2 %
0.1 %
改善一项输入,观察总预算减少多少先平方,再比较方差贡献面积比内标称量内标纯度样品称量

图 7.3-2 独立、小相对不确定度的乘除模型;质量分数沿用 89.55%,覆盖因子取二,未计入摩尔质量影响。

7.3.1 区分离散、偏差与测量不确定度

一个定量结果需要同时给出估计值及其可信程度。重复采集的数值接近,说明该重复层面的离散较小;若每次都使用不充分的弛豫条件或错误标准含量,数值仍可能共同偏离参照。精密度(precision)描述重复结果的接近程度,偏差(bias)涉及估计值与合适参照之间的系统差异,两者需要分别检查。

展开说明与推导

测量不确定度(measurement uncertainty)依据所用信息表征赋予被测量值的分散程度。它由测量模型与输入资料建立,既可能包括重复变化,也可能包括标准赋值、称量、响应和处理的不充分知识。误差描述某个结果相对参照或真值的差,不确定度则不能简单当成已经知道的误差绝对值。

本节补充分析本章结果所需的统计基础。JCGM 的 GUM 文档提供测量模型、标准不确定度和传播框架;Thompson、Ellison 与 Wood(2002)讨论单实验室方法验证。定量 NMR 的具体控制可参考 Malz 与 Jancke(2005)的验证研究。文献中某种方法达到的表现,不会自动成为当前仪器、基质和操作的结果。

重复结果可以紧密但偏离参照,也可以围绕参照但分散教学示意:竖虚线为可靠参照值重复很接近重复较分散左:离散小,仍有明显偏移右:围绕参照,单次离散较大必须检查偏差、重复层面,以及未包含在重复中的其他输入
图 7.3-3 原创概念示意,点的位置未对应实际实验数据。离散和偏移需要不同资料评估,不能由一种表现代替另一种。
定义 · 先写测量模型

测量模型把输出量与称量、标准含量、积分比及必要校正联系起来。不确定度评估应覆盖这个模型中的实际输入。只对最终结果做标准差计算,会遗漏所有在重复期间保持相同的输入。

展开说明与推导

发现显著偏差后,优先定位原因并采用有依据的校正,再评估校正的不确定度。不能任意扩大误差条来使错误方法看起来与参照相容,也不能把未知的所有偏差直接设为零。对尚未可靠控制的重叠、反应或损失,应说明方法范围和未解决的限制,再决定结果是否足以回答当前问题。

7.3.2 重复层面决定标准差的含义

对 \(n\) 个相同定义的观测 \(x_i\),算术平均为 \(\bar x=\sum_i x_i/n\),样本标准差为 \(s\)。在独立、相同分布且适用的条件下,平均值的标准不确定度可以用 \(s/\sqrt n\) 估计。单次结果的离散与平均值的离散不同,报告平均值时要说明平均了多少次以及怎样重复。

展开说明与推导
\[s=\sqrt{\frac{\sum_{i=1}^{n}(x_i-\bar x)^2}{n-1}},\qquad u(\bar x)=\frac{s}{\sqrt n}.\]
(7.3-1)

\(x_i,\bar x,s,u\) 使用相同单位,且 \(n\) 至少为二。教学五次独立制样的质量分数数值为 89.2%、89.8%、89.4%、89.6%、89.5%。平均为 89.5%,偏离平均的百分点分别为 −0.3、0.3、−0.1、0.1、0;平方和 0.20,除以四后开方,得到 \(s\approx0.2236\) 个百分点,平均的相应标准不确定度为 0.1000 个百分点。

这五个数的计算只估计当前重复设计涵盖的变化。若五次只是同一谱图重复读取积分,它们不包含称量、溶解、样品间差异或独立采集;若五份样品共用同一瓶标准和同一纯度赋值,标准含量的不确定度也不会因五次平均而除以根号五。重复设计的范围与计算公式一样需要写清。

同一条件、短时间内重复可以检查重复性;跨日期、操作人员或相关仪器条件的设计可检查更多变化。随时间漂移的数据可能彼此相关,连续扫描也可能共享基线和处理偏差。没有独立性时机械使用 \(s/\sqrt n\) 会低估平均值的不确定度。应先检查趋势、批次和相关来源,再决定合适统计模型。

GUM 中的 A 类、B 类区分评估方式。A 类通常由统计分析得到,B 类利用证书、校准、分辨率、历史资料或其他信息。这个分类不等于“随机误差”和“系统误差”的分箱。证书若给扩展不确定度,还要按其说明的覆盖因子等转换为标准不确定度,不能直接与一个标准差平方相加。

习题 7.3-1 · 两种重复设计

甲重复处理同一 FID 五次,乙独立称样、加入标准、溶解和采集五次。两者都得到小标准差,能否认为覆盖相同误差来源?

参考解答

不能。甲主要反映所采用处理的变化;乙包含更多制样和采集变化,但仍可能共享标准赋值、同一天平校准等输入。应按模型登记共同和独立来源,不把同一输入的不确定度重复缩小或重复计入。

7.3.3 传播要考虑灵敏度与相关性

设输出 \(y=f(x_1,\ldots,x_m)\),每个输入有估计值、标准不确定度及可能的相关性。当输入变化较小、模型可在当前范围内适当线性化时,合成标准不确定度由灵敏度系数加权得到。某个输入的不确定度很大,若输出对它不敏感,贡献可能仍小;反之,一个较小输入变化也可能显著影响结果。

展开说明与推导
\[u_c^2(y)=\sum_i c_i^2u^2(x_i)+2\sum_{i<j}c_ic_j\,u(x_i,x_j),\qquad c_i=\frac{\partial f}{\partial x_i}.\]
(7.3-2)

\(c_i\) 是灵敏度系数,单位为输出单位除以相应输入单位;\(u(x_i,x_j)\) 表示两个输入的协方差,单位为它们单位的乘积。输入独立时相应协方差为零。相关性可能来自同一校准、共同基线或共同制样,不能只因数据位于表格不同列就把它们当成独立。

上一节内标模型写成 \(w_A=R\,m_Sw_SM_A/(m_XM_S)\),其中 \(R=(I_A/N_A)/(I_S/N_S)\)。若各连续输入独立且相对不确定度较小,乘除模型的相对方差近似为各输入相对方差之和。报告氢数在指认已确定时为确定计数;若氢数归属尚有歧义,应先解决模型选择,不能用一个很小数值不确定度把错误指认包住。

教学预算取 \(u_r(R)=0.30\%\)、\(u_r(m_S)=0.10\%\)、\(u_r(w_S)=0.20\%\)、\(u_r(m_X)=0.10\%\),摩尔质量影响在本算例中暂略。独立近似给 \(u_r(w_A)=\sqrt{0.30^2+0.10^2+0.20^2+0.10^2}\%=0.3873\%\)。沿用 \(w_A=0.8955\),绝对标准不确定度约为 0.003468,即 0.3468 个百分点。相对百分比与质量分数的百分点必须分开。

面积比标准称量标准纯度样品称量分别贡献60%、6.67%、26.67%和6.67%的相对方差教学独立输入预算:比较方差贡献面积比 R标准质量 mS标准含量 wS样品质量 mX60.00%6.67%26.67%6.67%总相对方差:0.30² + 0.10² + 0.20² + 0.10²条形按贡献比例绘制;输入需使用标准不确定度
图 7.3-4 原创预算图。四项贡献按平方后的相对不确定度比较,不能直接用 0.30:0.10:0.20:0.10 当作方差比例。
展开说明与推导

面积比本身也含两个积分的相关性。若同一基线变化同时抬高两者,部分影响可能在比值中抵消;若两区域偏移方向不同,则可能放大。教学两输入相对不确定度均为 1%,相关系数为 0.8,商的相对不确定度约为 \(\sqrt{1^2+1^2-2\times0.8}\%=0.6325\%\),独立假设则给 1.4142%。相关系数应有数据或模型依据,不能为了减小报告数字任意设定。

基线影响可以从面积的定义检查:近似常数偏移乘以积分区间宽度,形成额外面积。两个峰共用同一基线偏移,如果区间宽度和真实面积不同,相对影响仍然不同;弱峰或宽区间常更敏感。可以合理改变基线处理与区间边界,观察输出变化并联系背景模型。重叠拟合给出的参数标准误差还依赖拟合模型,若峰数或背景选错,那个小标准误差不会包括模型遗漏。

线性传播在接近零、强非线性、非对称分布或参数边界附近可能不足。JCGM 101 给出通过分布传播的 Monte Carlo 框架;此类计算还需要合理输入分布、相关性与模型。重复运行程序很多次只能减少模拟误差,不能改善输入资料的真实性。实际采用时保存模型、分布假设和计算检查,不把模拟次数当作测量准确度。

习题 7.3-2 · 找主要贡献

本算例中仅把标准称量相对不确定度从 0.10% 降到 0.05%,合成值是多少?为何改善有限?

参考解答

新值为 \(\sqrt{0.30^2+0.05^2+0.20^2+0.10^2}\%\approx0.3775\%\),只从 0.3873% 小幅下降。面积比与标准含量的方差贡献更大;应优先检查主导项,同时避免把同一个来源同时放入面积比预算和重复制样预算。

7.3.4 方法验证针对样品和量程展开

方法验证检查方法能否在指定条件下完成指定任务。选择性、工作范围、精密度、偏差或回收、稳定性和对条件变化的响应各有用途。Malz 与 Jancke 的研究展示了对采集和处理的系统控制;它的结论属于相应验证设计。当前方法若改变溶剂、抑制序列或基质,就需要检查这些变化是否影响原有量关系。

展开说明与推导

选择性可通过空白、单组分参照、加标及其他可靠信号检查。一个峰在纯待测物中干净,在实际样品中仍可能被赋形剂或副产物污染。范围检验应覆盖拟报告的含量范围,检查响应和残差,不能只用一个相关系数很高的拟合替代逐点偏差检查。强峰端的接收非线性与弱峰端的基线误差可能同时存在。

若样品本身不均匀,取自同一溶液的重复测量只能说明这份溶液的表现,不能直接代表整批材料。制样前的取样范围、均匀化方式与份额进入被测量定义。需要报告批次平均组成时,可能还要设计不同位置或不同容器的取样。把这些变化与仪器噪声区分,可以避免大量重复扫描掩盖真正主导的取样差异。

加标回收可用加入前后目标量的增量与已知加入量比较。若加标在提取后进行,它主要检查提取后测量与基质响应,不能独立证明原样中的目标已被完全提取。若基质中的原有目标与新加标准行为不同,回收也有解释边界。报告加标阶段、量和计算方式,避免只给一个“回收良好”的形容。

与参照方法比较时,参照值也有不确定度。两结果的差应结合双方的不确定度和可能共享的校准来源解释。若共用相同标准赋值,二者接近可能部分来自共同输入;若独立方法在合理区间内相容,则提供更有区分能力的检查。不能把参照显示的末位当作无误差真值,再将任何末位差异都判成当前方法失败。

检出限关注能否可靠区分信号与空白,定量限关注达到规定定量性能的最低范围,两者不同。低含量峰能在图上看到,不保证面积足以准确量化。限值应通过适合的空白、低水平样品及完整流程评估;从一张高浓度谱的信噪比简单外推,可能遗漏低水平时重叠、基线和污染的主导作用。

方法 · 小幅改变关键条件

在合理范围内改变重复时间、积分边界、温度或相关制样条件,检查输出的稳定性。对变化敏感的因素进入方法控制或不确定度评估。若变化使样品化学状态发生改变,应先处理模型适用范围,而非只将它当作随机波动。

7.3.5 扩展区间须注明覆盖依据

扩展不确定度(expanded uncertainty)通常写为 \(U=k\,u_c\),\(k\) 为覆盖因子。在适当近似正态分布和足够有效自由度等条件下,取二可给接近 95% 的覆盖;样本很少、分布明显非对称或模型受边界限制时,不能无条件这样解释。区间的概率或覆盖含义应与所用统计框架一致。

展开说明与推导

沿用教学预算,若采用 \(k=2\),得到 \(U\approx0.006936\),即 0.6936 个百分点,可写质量分数 89.55%,扩展不确定度 0.69 个百分点,并说明覆盖因子二及条件假设。也可将无量纲结果写为 \(0.8955\pm0.0069\)。应保留被测量定义和预算信息,避免把 ±0.69% 含混地解释成相对百分比。

警示 · 平均不能缩小共同校准项

增加独立样本可以降低部分重复变化,不能使所有标准、共同响应和模型不确定度都除以根号样本数。若要改善最终结果,先查预算的主要贡献,以及它是否已经在别处计入。

展开说明与推导

质量分数的物理范围有限,测量估计值在接近上限时可能受到噪声或偏差影响。若计算超过一,应检查模型、样品量定义与参照,不能只把结果截成一百并删除不确定度。复杂边界问题可参考 Duewer 等的纯度赋值讨论。最终报告说明验证范围、主要限制和区间来源,使他人知道结果能支持哪一种比较或决策。

习题 7.3-3 · 窄区间是否足够

计算得到很小标准差,但预算没有标准纯度、积分重叠和外标响应传递项。取 \(k=2\) 后能否称为完整约95%测量区间?

参考解答

不能。覆盖因子只放大已经评估的合成标准不确定度,不能自动补入遗漏来源。应先完善测量模型、检查影响并建立合理输入,避免重复计入,再按分布与有效自由度确定区间。否则只描述有限重复层面的离散。

关键术语

精密度 (precision)
指定重复条件下结果的接近程度。
偏差 (bias)
估计结果相对合适参照的系统差异。
标准不确定度 (standard uncertainty)
以标准差形式表达的不确定度。
合成标准不确定度 (combined standard uncertainty)
输入经模型传播后得到的标准不确定度。
协方差 (covariance)
描述两个输入共同变化的量。
灵敏度系数 (sensitivity coefficient)
输出对一个输入局部变化的导数。
覆盖因子 (coverage factor)
将合成标准不确定度转为扩展量的乘数。
方法验证 (method validation)
检验方法在指定条件和范围内是否适合任务。

参考文献与延伸阅读

  1. JCGM. 2008. Evaluation of measurement data — Guide to the expression of uncertainty in measurement. JCGM 100:2008. 标准不确定度与传播框架。
  2. JCGM. 2008. Evaluation of measurement data — Supplement 1 to the GUM — Propagation of distributions using a Monte Carlo method. JCGM 101:2008. 非线性分布传播。
  3. Thompson M, Ellison SLR, Wood R. 2002. Harmonized guidelines for single-laboratory validation of methods of analysis. Pure and Applied Chemistry, 74: 835–855. 验证设计。
  4. Malz F, Jancke H. 2005. Validation of quantitative NMR. Journal of Pharmaceutical and Biomedical Analysis, 38: 813–823. qNMR 原始验证研究。
  5. Duewer DL, Lippa KA, et al. 2023. Methods for the SI Value Assignment of the Purity of Organic Compounds (IUPAC Technical Report). Pure and Applied Chemistry. 纯度赋值与不确定度。