单因素方差分析（One Way ANOVA）

来源：互联网发布：shake it 编辑：程序博客网时间：2024/04/28 08:57

单因素方差分析是指对单因素试验结果进行分析，检验因素对试验结果有无显著性影响的方法。

单因素方差分析是两个样本平均数比较的引伸，它是用来检验多个平均数之间的差异，从而确定因素对试验结果有无显著性影响的一种统计方法。

因素：影响研究对象的某一指标、变量。
水平：因素变化的各种状态或因素变化所分的等级或组别。
单因素试验：考虑的因素只有一个的试验叫单因素试验。

例如，将抗生素注入人体会产生抗生素与血浆蛋白质结合的现象，以致减少了药效。下表列出了5种常用的抗生素注入到牛的体内时，抗生素与血浆蛋白质结合的百分比。现需要在显著性水平 $α = 0.05$ 下检验这些百分比的均值有无显著的差异。设各总体服从正态分布，且方差相同。

青霉素四环素链霉素红霉素氯霉素29.627.35.821.629.224.332.66.217.432.828.530.811.018.325.032.034.88.319.024.2

　　在这里，试验的指标是抗生素与血浆蛋白质结合的百分比，抗生素为因素，不同的5种抗生素就是这个因素的五个不同的水平。假定除抗生素这一因素外，其余的一切条件都相同。这就是单因素试验。试验的目的是要考察这些抗生素与血浆蛋白质结合的百分比的均值有无显著的差异。即考察抗生素这一因素对这些百分比有无显著影响。这就是一个典型的单因素试验的方差分析问题

与通常的统计推断问题一样，方差分析的任务也是先根据实际情况提出原假设 $H 0$ 与备择假设 $H 1$ ，然后寻找适当的检验统计量进行假设检验。本节将借用上面的实例来讨论单因素试验的方差分析问题。

　　在上例中，因素A（即抗生素）有s（=5）个水平 $A_1,A_2,\cdots,A_5$ ，在每一个水平 $A_j(j=1,2,\cdots,s)$ 下进行了 $n j = 4$ 次独立试验，得到如上表所示的结果。这些结果是一个随机变量。表中的数据可以看成来自s个不同总体（每个水平对应一个总体）的样本值，将各个总体的均值依次记为 $\mu_1,\mu_2,\cdots,\mu_s$ ,则按题意需检验假设

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

　　 $H_1:\mu_1,\mu_2,\cdots,\mu_s$ 不全相等

为了便于讨论，现在引入总平均 $μ$

　　 $\mu=\frac{1}{n}\sum_{j=1}^s n_j \mu_j$ 其中： $n=\sum_{j=1}^s n_j$

再引入水平 $A j$ 的效应 $δ j$

$\delta_j=\mu_j-\mu(j=1,2\ldots,s)$

显然有 $n_1\delta_1+n_2\delta_2+\cdots+n_s\delta_s=0$ ， $δ j$ 表示水平 $A j$ 下的总体平均值与总平均的差异。

利用这些记号，本例的假设就等价于假设

　　 $H_0:\delta_1=\delta_2=\cdots=\delta_s=0$

　　 $H_1:\delta_1,\delta_2,\cdots,\delta_s$ 不全为零

因此，单因素方差分析的任务就是检验s个总体的均值 $μ j$ 是否相等，也就等价于检验各水平 $A j$ 的效应 $δ j$ 是否都等于零。

　　2. 检验所需的统计量

　　假设各总体服从正态分布，且方差相同，即假定各个水平 $A_j(j=1,2,\cdots,s)$ 下的样本 $x_{1j},x_{2j},\cdots,x_{n_jj}$ 来自正态总体 $N (μ j,σ 2)$ ， $μ j$ 与 $σ 2$ 未知，且设不同水平 $A j$ 下的样本之间相互独立，则单因素方差分析所需的检验统计量可以从总平方和的分解导出来。下面先引入：

　　水平 $A j$ 下的样本平均值:

　　 ${\overline x}_{\bullet j}=\frac{1}{n}\sum_{i=1}^{n_j}x_{ij}$

　　数据的总平均:

　　 $\overline{x}=\frac{1}{n}\sum_{j=1}^s\sum_{i=1}^{n_j}x_{ij}=\frac{1}{n}\sum_{j=1}^sn_j{\overline x}_{\bullet j}$

　　总平方和:

　　 $S_T=\sum_{j=1}^s \sum_{i=1}^{n_j}{(x_{ij}-\overline x)}^2$

总平方和 $S T$ 反映了全部试验数据之间的差异，因此 $S T$ 又称为总变差。将其分解为

　　 $S T = S E + S A$

其中:

　　 $S_E=\sum_{j=1}^s \sum_{i=1}^{n_j}{(x_{ij}\overline x}_{\bullet j})}^2$

　　 $S_A=\sum_{j=1}^s \sum_{i=1}^{n_j}{({\overline x}_{\bullet j\overline x)}^2=\sum_{j=1}^s n_j({\overline x}_{\bullet j}-\overline x)^2)$

上述 $S E$ 的各项 $(x_{ij}\overline x}_{\bullet j})^2$ 表示了在水平 $A j$ 下，样本观察值与样本均值的差异，这是由随机误差所引起的，因此 $S E$ 叫做误差平方和。 $S A$ 的各项 $n_j({\overline x}_{\bullet j\overline x)^2$ 表示了在水平 $A j$ 下的样本平均值与数据总平均的差异，这是由水平 $A j$ 以及随机误差所引起的，因此 $S A$ 叫做因素A的效应平方和。