单因子方差分析

单因子方差分析可检验组均值的差异

单因子方差分析是一种统计方法，可用于检验三个或更多总体均值相等的原假设 (H₀) 与至少有一个均值不同的备择假设 (H_a)。对于 k 个均值，我们可用正式的统计假设表示法写为：

$ H_0:\mu_1=\mu_2=\cdots=\mu_k $

$ H_a:\mathrm{并非所有的均值都相等} $

其中 $\mu_i$ 是因子第 i 个水平的均值。

您可能在想，什么情况下我需要确定多个总体的均值是相同还是不同呢？一种常见的情况是，您怀疑某个特定的过程变量是该过程的重要结果的驱动因素。例如，您可能怀疑不同的生产批次、操作员或原材料批次会影响生产过程的输出（即质量测量值）。

为了验证您的想法，您可以使用此自变量（即因子）的三个或更多变异（即水平）来运行过程，然后从每次运行的结果中抽取观测值样本。如果使用方差分析比较各组观测值的均值时发现了差异，那么（假设您已正确完成所有操作！），您就有证据来证明您的猜想是正确的 — 您所检验的因子可能对结果有影响！

单因子方差分析示例

让我们来详细地研究一个单因子方差分析示例。假设您在一家制造粘合剂凝胶的公司工作，这些凝胶装在小罐中出售。凝胶的粘度很重要：太稠会难以施涂；太稀会粘合性不足。您最近收到了一些客户对产品不满意的反馈，他们抱怨粘合剂的粘度不像以前那样稳定。老板要求您调查此事。

您认为第一步最好是检查五个最近生产批次的平均粘度。如果发现批次之间存在差异，似乎就可以确认问题的确存在。另外，它还有助于您对可能导致批次之间不一致的因素建立假设。

使用一种仪器来测量粘度，将其主轴浸入粘合剂罐中并进行旋转。该检验会产生一个测量值，称为扭矩阻抗。然后分别从最近 5 个批次中随机选择 5 个罐子进行检验，取得每个罐子的扭矩阻抗测量值，并绘制数据。

从数据图中可以观测到，第 3 批罐子的扭矩测量值看起来低于其他批次样本的扭矩测量值。计算所有测量值的均值后，发现第 3 批的平均扭矩是 26.77，远低于其他 4 个批次，其他批次的均值都是 30 左右。

表 1：检验 5 个批次粘合剂得出的平均扭矩测量值

方差分析表

方差分析结果通常显示在方差分析表中。方差分析表包含：

• 来源：变异的来源，如所检验的因子（本例为批次）、误差和总计。
• DF：每个变异来源的自由度。
• 平方和：每个变异来源的平方和 (SS) 以及所有来源的总计。
• 均方：平方和除以与其相关联的自由度。
• F 比：因子（批次）的均方除以误差的均方。
• 概率 > F：p 值。

表 2：方差分析表（结果取自我们的扭矩测量值）

我们将在下文中解释这张表的组成部分是如何得出的。在此表中，目前应关注的一个关键元素是 p 值。p 值用于可评估“所有均值都相同”这个原假设的有效性。在这个例子中，p 值（概率 > F）为 0.0012，可以证明并非所有均值都相同。我们的样本提供了证明：在 5 个批次中，有一个或多个批次之间的平均扭矩阻抗值存在差异。

基于上述结果，您决定对第 3 批做进一步的检验。这时，您可以在报告中这样写：分别从 5 个最近的生产批次中抽取了 5 罐产品并测量了其扭矩。方差分析发现，观测值支持批次之间的平均扭矩存在差异 (p = 0.0012)。数据图显示第 3 批的平均扭矩 (26.77) 比其他 4 个批次低。因此，我们将对第 3 批做进一步的评估。

单因子方差分析计算

现在，我们来深入地思考一下扭矩测量值的示例。请回忆一下，我们有 5 批材料，从每个批次中随机选择了 5 罐进行检验。这称为单因子设计。批次这 1 个因子有 5 个水平。每个水平又重复（检验） 5 次。检验结果如下所示。

表 3：按批次列出的扭矩测量值

为了探索上面的方差分析表（表 2）中的计算结果，让我们首先建立以下定义：

$n_i$ = 待处理的观测值数量 $i$（在我们的例子中是批次 $i$）

$N$ = 观测值总数

$Y_{ij}$ = 第 i 次处理的第 j 个观测值

$\overline{Y}_i$ = 第 i 次处理的样本均值

$\overline{\overline{Y}}$ = 所有观测值的均值（总均值）

平方和

记住了这些定义之后，让我们来处理方差分析表的“平方和”列。平方和为我们提供了一种方法，可通过关注每个数据点与该数据集中所有数据点均值之间的差异来量化数据集的变异性。下面的公式将总体变异性分成两部分：由模型或因子水平引起的变异性，以及由随机误差引起的变异性。  

$$ \sum_{i=1}^{a}\sum_{j=1}^{n_i}(Y_{ij}-\overline{\overline{Y}})^2\;=\;\sum_{i=1}^{a}n_i(\overline{Y}_i-\overline{\overline{Y}})^2+\sum_{i=1}^{a}\sum_{j=1}^{n_i}(Y_{ij}-\overline{Y}_i)^2 $$

$$ SS（总计）\;     =     \;SS（因子）\;     +     \;SS（误差）$$

虽然该方程式看起来很复杂，但只要分别关注每个元素就会更容易掌握这个方程式。下面的表 4 列出了公式的每个组成部分，然后将其构建到组成平方和的平方项中。第 1 列的 ($Y_{ij}$) 数据包含了我们在上面表 3 中收集到的扭矩测量值。

表 4：平方和计算

自由度 (DF)

自由度 (DF) 是一个与平方和相关联的值。自由度表示用来计算每个平方和的独立信息的数量。对于因子水平为 k（本例为 5 个批次）且总观测次数为N（每批次 5 个罐子，总共 25 个）的单因子设计，其自由度如下：

表 5：确定自由度

均方 (MS) 和 F 比

我们将每个平方和除以相应的自由度即可得到均方值。当原假设为真（即均值相等）时，MS（因子）和 MS（误差）均为误差方差的估计值且大小大致相同。它们的比率（即 F 比）将接近 1。当原假设不为真时，MS（因子）将大于MS（误差）且其比率将大于 1。在我们的粘合剂检验例子中，计算出的 F 比是 6.90，这充分证明了原假设的均值不相等。

表 6：计算均方和 F 比