误差预算先找结果依赖哪些输入
不确定度描述结果范围,预算能帮助选择最有效的改进步骤。
读图时保留这项条件:标准不确定度和扩展不确定度不同,相关输入不能总按独立处理。
改善一项输入,观察总预算减少多少
分别拖动四项相对标准不确定度,看平方后贡献与合成结果。
7.3.1 区分离散、偏差与测量不确定度
一个定量结果需要同时给出估计值及其可信程度。重复采集的数值接近,说明该重复层面的离散较小;若每次都使用不充分的弛豫条件或错误标准含量,数值仍可能共同偏离参照。精密度(precision)描述重复结果的接近程度,偏差(bias)涉及估计值与合适参照之间的系统差异,两者需要分别检查。
展开说明与推导
测量不确定度(measurement uncertainty)依据所用信息表征赋予被测量值的分散程度。它由测量模型与输入资料建立,既可能包括重复变化,也可能包括标准赋值、称量、响应和处理的不充分知识。误差描述某个结果相对参照或真值的差,不确定度则不能简单当成已经知道的误差绝对值。
本节补充分析本章结果所需的统计基础。JCGM 的 GUM 文档提供测量模型、标准不确定度和传播框架;Thompson、Ellison 与 Wood(2002)讨论单实验室方法验证。定量 NMR 的具体控制可参考 Malz 与 Jancke(2005)的验证研究。文献中某种方法达到的表现,不会自动成为当前仪器、基质和操作的结果。
测量模型把输出量与称量、标准含量、积分比及必要校正联系起来。不确定度评估应覆盖这个模型中的实际输入。只对最终结果做标准差计算,会遗漏所有在重复期间保持相同的输入。
展开说明与推导
发现显著偏差后,优先定位原因并采用有依据的校正,再评估校正的不确定度。不能任意扩大误差条来使错误方法看起来与参照相容,也不能把未知的所有偏差直接设为零。对尚未可靠控制的重叠、反应或损失,应说明方法范围和未解决的限制,再决定结果是否足以回答当前问题。
7.3.2 重复层面决定标准差的含义
对 \(n\) 个相同定义的观测 \(x_i\),算术平均为 \(\bar x=\sum_i x_i/n\),样本标准差为 \(s\)。在独立、相同分布且适用的条件下,平均值的标准不确定度可以用 \(s/\sqrt n\) 估计。单次结果的离散与平均值的离散不同,报告平均值时要说明平均了多少次以及怎样重复。
展开说明与推导
\(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)\),每个输入有估计值、标准不确定度及可能的相关性。当输入变化较小、模型可在当前范围内适当线性化时,合成标准不确定度由灵敏度系数加权得到。某个输入的不确定度很大,若输出对它不敏感,贡献可能仍小;反之,一个较小输入变化也可能显著影响结果。
展开说明与推导
\(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 个百分点。相对百分比与质量分数的百分点必须分开。
展开说明与推导
面积比本身也含两个积分的相关性。若同一基线变化同时抬高两者,部分影响可能在比值中抵消;若两区域偏移方向不同,则可能放大。教学两输入相对不确定度均为 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)
- 检验方法在指定条件和范围内是否适合任务。
参考文献与延伸阅读
- JCGM. 2008. Evaluation of measurement data — Guide to the expression of uncertainty in measurement. JCGM 100:2008. 标准不确定度与传播框架。
- JCGM. 2008. Evaluation of measurement data — Supplement 1 to the GUM — Propagation of distributions using a Monte Carlo method. JCGM 101:2008. 非线性分布传播。
- Thompson M, Ellison SLR, Wood R. 2002. Harmonized guidelines for single-laboratory validation of methods of analysis. Pure and Applied Chemistry, 74: 835–855. 验证设计。
- Malz F, Jancke H. 2005. Validation of quantitative NMR. Journal of Pharmaceutical and Biomedical Analysis, 38: 813–823. qNMR 原始验证研究。
- 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. 纯度赋值与不确定度。