JMP数据分析论坛

标题: 单因素方差解析及其在药品检测实验室中的应用 [打印本页]
作者: admin    时间: 2018-10-25 09:57
标题: 单因素方差解析及其在药品检测实验室中的应用
单因素方差解析是常见的一种对连续数据的方差解析方法，其中最广泛的是多组数据的差异显著性检验，以找出组间差异。

然而，除上述应用外，日常工作中，还经常遇到下述情况：

在药品检测实验室，对方法精密度的计算和评价；


这些工作都离不开对实验所测数据的方差解析，求出其中所需的表达参数。

在计算上述各种目的所需指标时，需要有周密而完善的实验设计。一个切实可行的实验设计，不仅可以有效排除偶然因素造成的变异性，使实验误差降到最低限度，从而使所获结果更可靠，而且还可以节省人力、时间，降低实验成本，达到高效的实验目的。目前，在进行上述实验或评价中，常见的问题是手工计算繁琐、实验设计不规范、所需计算的方差成分不了解，从而导致许多困惑。

本文拟利用单因素多水平随机效应模型的实验设计，形成一个通用模板，对上述情形所需的参数（指标）通过方差解析的方式进行计算演示，以便当遇到上述情况时，便于大家理解和科学、快速地计算所需指标。

本文作者谭德讲老师等，经作者授权发布。如需转载，请联系原作者授权。

1. 单因素多水平随机效应模型的方差分析及其具体应用

根据研究因素个数的不同，方差解析可以分为单因素方差解析、双因素方差解析等。试验设计中只有一个因素取不同的水平进行试验，而其他因素保持不变，这样的试验称为单因素试验（one factor trial），相应的方差解析就是单因素的方差解析。

根据因素水平的选择，方差解析的模型又可分为固定效应模型（fixed-effects model）、随机效应模型（random-effects model）和混合效应模型（mixed-effects model）。



随机效应模型的实验设计是经典线性模型设计的一种推广，就是把原各固定效应因素视为来自总体的随机效应因素，而在原来各因素的记录数据或观测（如重复次数或多时间点的记录）间建立起一定的相关性模型。


该模型可用来拟合非独立观测的数据，具有计算简便、能选择相关的预测参数（predictors），对因素的压缩程度较小等优点。

本文研究探讨的即是采用单因素多水平随机效应的实验设计模型。详见表1。

表1 单因素多水平随机效应试验模型
[attach]4084[/attach]

假设因素有j个水平，j=1,2,3,…,c ，在每个水平j下进行i次独立重复，i=1,2,3,…,r，得到的试验结果为xj,1,xj,2,...,xj,r, 试验结果服从正态分布N（μj,σ2）; 应特别指出在每个水平下重复次数r=2时，其计算形式较重复次数r≥3有一定的区别。

实验的总样本数为n=cr。

因素在j水平下观测值的样本均值，即组内均值：
[attach]4085[/attach]

样本的总均值为：
[attach]4086[/attach]

总离差平方和SST，描述了所有观察效应之间的差异。
[attach]4087[/attach]

组内离差平方和SSE反映了各组内部因重复性实验而产生的随机抽样误差。
[attach]4088[/attach]

组间离差平方和SSA主要反映了各组样本均值间的差异，既包括了随机抽样误差，也包括不同的实验室所造成的系统误差。
[attach]4089[/attach]
总离差平方和SST的自由度dfT=n-1；组内离差平方和SSE的自由度dfE=n-c；组间离均差平方和SSA的自由度dfA=c-1；
由此得组内均方：
[attach]4090[/attach]

组间均方：
[attach]4091[/attach]

组内方差分量（Sr）为：
[attach]4092[/attach]

组内方差分量（SL）为：
[attach]4093[/attach]

上述的单因素方差解析过程均可由软件计算汇总得“方差分析表”，本文采用JMP软件，如图1。

上述设计和相应的方差计算，对药品检测实验室中常见的分析方法精密度计算、实验室检测能力评价、样本均匀性评估以及标准品协作标定的赋值等具有很好的应用价值。这四方面所需的参数，均可由此单因素的方差解析获得MSE和MSA，然后根据每个试验的试验目的进行进一步的运算。

下面我们来逐一说明。

1.1 应用于确定分析方法的精密度计算

在确定分析方法的精密度试验中，

中间精密度或再现性标准差的计算公式为：
[attach]4094[/attach]

重复性标准差的计算公式：
[attach]4095[/attach]


1.2 应用于样品均匀性的计算

在进行实验室能力评定或标准物质标定过程中，所使用的实验样品均应满足均匀性这一前提条件。

样品间的标准差，亦即组间标准差SL，与能力评定的标准差只要满足下式就表明样本是均匀的。
[attach]4096[/attach]

对样品间标准差的计算一般按照下列规则进行：

由主实验室制备并且包装一轮能力验证的样品或待标标准品的样品，然后独自或指定一个实验室进行均匀性的检验。

实验时，从本批总体中随机选择c份，即在此试验中的因素为随机样本，从每个样品中制备r个独立测试份量，整个测量在独立重复性条件下完成。

样本间标准差的计算公式为：
[attach]4097[/attach]


1.3 应用于实验室的能力评定标准的计算

实验室能力评定实验中所使用的试验样本均匀、稳定，实验室使用已标准化的测量方法，并且该方法的重复性和再现性可获得时，能力评定的标准差也可由上述的单因素方差解析进一步求算获得。

能力评定标准差的计算如下：
[attach]4098[/attach]


1.4 应用于标准物质的协作标定

协作标定（Collaborative Calibration）是指在多个具有同等能力的实验室间，使用法定方法进行测量，各实验室所得的实验数据按统计程序处理后得到标准物质的特性值和不确定度/置信区间。

所测样品的总均值为：
[attach]4099[/attach]

标准不确定度：
[attach]4100[/attach]


即：
[attach]4101[/attach]

置信区间为：
[attach]4102[/attach]
t1-α:c-1 代表t分布，左侧面积百分数为1-α，自由度为c-1。


2. 单因素多水平的随机效应模型的方差实例分析

根据表1的单因素多水平随机效应的实验设计，我们将假定获得的检测数据汇总于表2中。

表2 原始数据汇总表
[attach]4103[/attach]

通过统计软件（基于JMP软件）对表2进行数据的方差解析可得图1的方差分析表。

图1.方差分析表
[attach]4104[/attach]

由表2可知r=3, c=8, n=24。

由图1确定SSA、SSE和SST的自由度分别为：
[attach]4105[/attach]

离均差平方和分别为：
[attach]4106[/attach]

均方分别为：
[attach]4107[/attach]

下面将以此实例作为上述四类试验的试验结果，对于不同实验类型所需参数进行计算演示。


2.1分析方法精密度试验中的重复性标准差和中间精密度或再现性标准差计算

假定欲对某一实验方法的精密度进行评价。

如果表2中的检测结果是在同一实验室中，使用同一样品（均质），采用不同分析人员、或不同仪器等组合为8个条件，每个条件下对样本重复测定3次所获得的，即可求出方法的重复性标准差和中间精密度；

如果表2中的数据是采用8个实验室对同一样品进行检测，每个实验室重复测定3次，即可获得方法的重复性标准差和再现性标准差。计算结果如下：

根据公式10，中间精密度或再现性标准差：
[attach]4108[/attach]

根据公式11，重复性标准差：
[attach]4109[/attach]


2.2 样本均匀性试验中样品间变异的计算

假定表2的原始数据是对在实验室能力评定或标准物质标定中制备的样品进行的均匀性研究结果，即在同一个实验室中，从包装好的样本中随机选择8份，每一份样本制备3个测定份量所测得的原始数据。

根据公式13，样本间标准差为：
[attach]4110[/attach]

根据公式11，样品内变异相当于重复性标准差，为：
[attach]4111[/attach]

将SL与能力评定标准进行比较即可判断该批制备的样品是否具有均匀性。


2.3实验室能力评定标准差的计算

假定表2的原始数据是经过8个经验实验室对某一产品（已通过均匀性、稳定性评估）采用同一方法测得的原始数据，那么，可以使用这些数据求出作为能力评定标准的标准差。


根据公式14，能力评定标准差：
[attach]4112[/attach]

2.4 标准物质协作标定的均值和不确定度（置信区间）计算

假定表2的原始数据是协作标定的结果，即筛选出8个具有同等能力的实验室，每个实验室均采用规定方法，重复测定3次实验数据。

则标准物质的均值根据公式2得：
[attach]4113[/attach]

根据公式15，其标准不确定度为：
[attach]4114[/attach]

置信区间为：
[attach]4115[/attach]



3. 讨论与启发

1） 本文给出的方差解析及各种计算，需要有严格的实验设计和数据探测，并保证其数据的独立性、正态性（离群值检查或对数据进行转换）和方差齐性这些前提。

本实例是一个数据集完整的嵌套式随机效应的单因素方差解析特例。当数据集不完全符合表1的结构，导致数据结构不均衡时，公式更复杂，这时，最好直接使用专业统计软件，计算出组内和组间的方差分量（如图1中的下半部分）；有些需使用其他计算方法，如标准物质的合并计算，需要使用权重法等。

2） 如若所建立的模型中在每个水平的试验重复测定次数仅为r=2时，一般使用极差法，即：
[attach]4116[/attach]

这与当r≥3的计算方式不同，后者使用每个水平的标准差计算，每个水平的标准差记为：
[attach]4117[/attach]
总组内均方MSE见公式17。r=2与r≥3在方差解析的其他过程均一致。
[attach]4118[/attach]

3） 在实际应用中对于四类分析应就具体的试验要求进行实验设计

研究表明，增加实验组数（c）比增加每组的重复次数（r）可更大地增加所得参数的稳健性或效能。

本文的主要目的，是拟通过一个适用于四类试验的模型建立一个通用的计算模式，在实例分析中没有考虑到具体应用中因素水平的个数（c）以及在每个水平下的重复测定次数（r）。

实际实施过程中，由于受现实条件（如成本）限制，需特别权衡成本与结果可靠性问题。一般应结合具体法规和实验的变异大小进行考虑，并尽可能地根据实验设计要求，增加实验组数。如对参加实验室能力认证的独立样本数c≥10，每个样本重复数r≥2次。

本文通过建立一个适用于上述四类分析的通用性实验设计，即单因素多水平随机效应模型的实验设计，利用JMP统计软件对该类型数据进行方差解析，将原本复杂的方差分析过程简化为方差分析表，并进一步结合具体应用中的方差成分需求进行计算演示，使药品检测试验室的分析人员能够用最短的时间快速理解，并系统掌握获得稳健可靠的所需参数，为今后在药品检测实验室开展相应的工作提供了很好的参考依据。

想体验JMP 14的朋友们，欢迎点击这里下载JMP 30天免费试用版。

关注JMP官方微信公众号，了解更多统计分析思路、JMP使用技巧、最新活动资讯等：

作者: Elina3355    时间: 2018-11-26 14:35
得出的这些结果都要和既定标准进行比较才会有意义，但往往就是这些标准是最难确定的
作者: Elina3355    时间: 2018-11-26 14:38
就比如得出的中间精密度或再现性标准差多大才算再现性比较好呢



欢迎光临 JMP数据分析论坛 (http://www.jmpforum.net/) X3