SAS讲义 第二十五课方差分析

第二十五课 方差分析当影响观察结果的影响因素（原因变量或分组变量）的水平数大于2或原因变量的个数大于1个，一元时常用F 检验（也称一元方差分析），多元时用多元方差分析（最常用Wilks ’∧检验）。

一、 方差分析概述方差分析（analysis of variance ）又称变异数分析，可简记为ANOV A ，主要用于检验计量资料中的两个或两个以上均值间差别显著性的方法。

当欲比较几组均值时，理论上抽得的几个样本，都假定来自正态总体，且有一个相同的方差，仅仅均值可以不相同。

还需假定每一个观察值都由若干部分累加而成，也即总的效果可分成若干部分，而每一部分都有一个特定的含义，称之谓效应的可加性。

所谓的方差是离均差平方和除以自由度，在方差分析中常简称为均方MS （mean square ）。

1. 方差分析的基本思想根据效应的可加性，将总的离均差平方和分解成若干部分，每一部分都与某一种效应相对应，总自由度也被分成相应的各个部分，各部分的离均差平方除以相应部分的自由度得出各部分的均方，然后列出方差分析表算出F 值，作出统计推断。

方差分析的关键是总离均差平方和的分解，分解越细致，各部分的含义就越明确，对各种效应的作用就越了解，统计推断就越准确。

方差分析表的一般形式见表25.1所示：表25.1 方差分析表形式变异来源 source 离差平方和SS 自由度 df 均方 MS F 统计量FP 概率值P 效应S 1 SS 1 df 1 MS 1= SS 1/df 1 F 1(df 1, df e )= MS 1/ MS e P 1 效应S 2 SS 2 df 2 MS 2= SS 2/df 2F 2(df 2, df e )= MS 2/ MS eP 2 …… …… …… …… ……效应S m SS m df m MS m = SS m /df m F m (df m , df e )= MS m / MS eP m 误差S e SS e df e MS e = SS e /df e总变异S TSS T = SS 1+ SS 2+…+ SS m + SS edf T =df 1+ df 2+…+ df m + df eMS T = SS T /df TF T (df T , df e )= MS T / MS eP T表中变异来源一栏，可分为总变异（total ），误差（residual ）,各个效应（effect ）相对应的项。

效应项与试验设计或统计分析的目的有关，一般有：主效应（包括各种因素），交互影响项（因素间的多级交互影响），协变量（来自回归的变异项），等等。

当分析和确定了各个效应项S 后，根据原始观察资料可计算出各个离均差平方和SS ，再根据相应的自由度df ，由公式MS=SS/df ，求出均方MS ，最后由相应的均方，求出各个变异项的F 值，F 值实际上是两个均方之比值，通常情况下，分母的均方是误差项的均方。

根据F 值的分子、分母均方的自由度f 1和f 2，在确定显著性水平为α情况下，由),(21f f F 临界值表查得单侧αF 界限值。

当αF F <时，则α>P ，不拒绝原假设0H ，说明不拒绝这个效应项的效应为0的原假设，也即这个效应项是可能对总变异没有实质影响的；如果αF F ≥,则α≤P ，拒绝原假设0H ，说明拒绝这个效应项的效应为0的原假设，也即这个效应项是很可能对总变异有实质影响的。

2. 方差分析的试验设计为了确定方差分析表中各个有关效应项，需要在试验设计阶段就作出安排，再根据设计要求进行试验，得出原始观察值，按原来设计方案算出方差分析表中的各项。

在试验设计阶段常需要作主要四个方面的考虑： 1) 研究的主要变量方差分析的主要变量，也称响应变量或因变量（dependent variable ），它是我们试验所要观察的主要指标。

一次试验时可以有多个观察指标，方差分析时也可以同时对多个因变量进行分析。

2) 因素和水平试验的因素（factor ）可以是品种、人员、方法、时间、地区等等，因素所处的状态叫水平（level ）。

在每一个因素下面可以分成若干水平。

例如，某工厂的原料来自四个不同地区，那么用不同地区的原料生产的产品质量是否一致呢？所要比较的地区就是因素，四个地区便是地区这一因素的四个水平。

当某个主要因素的各个水平间的主要因变量的均值呈现统计显著性时，必要时可作两两水平间的比较，称为均值间的两两比较。

3) 因素间的交互影响多因素的试验设计，有时需要分析因素间的交互影响（interaction ），2个因素间的交互影响称为一级交互影响，例如因素A 与因素B 的一级交互影响可记为A ×B ，3个因素间的交互影响称为二级交互影响，例如因素A 与因素B 与因素C 的二级交互影响可记为A ×B ×C 。

当交互影响项呈现统计不显著时，表明各个因素独立，当呈现统计显著时，就需要列出这个交互影响项的效应，以助于作出正确的统计推断。

二、 单因素方差分析单因素方差分析（one factor ANOV A 或one-way ANOV A ）或称为完全随机设计的方差分析（completely random design ANOV A ）。

试验设计时按受试对象的抽取或分组的随机程度不同可细分为以下两类：● 完全随机设计——从符合条件的总体中完全随机地抽取所需数目的受试对象，再将全部受试对象完全随机地分配到k 组中去。

此时，受试对象与试验因素间无直接联系。

● 组内完全随机设计——按试验因素的k 个水平将全部受试对象划分成k 个子总体，再分别从k 个子总体中完全随机地抽取所需数目的受试对象。

此时，试验因素的各水平决定了受试对象各自应该归属的组别。

设因素A 有k 个水平k A A A ,,,21 ，在每一个水平下考察的指标可以看成一个总体，现有k 个水平，故有k 个总体，并假定：① 每一总体均服从正态分布； ② 每一总体的方差相同；③ 从每一总体中抽取的样本相互独立。

我们要比较各个总体的均值是否一致，就是要检验各总体的均值是否相同，设第i 个总体的均值为i μ，那么就是要检验如下原假设：k H μμμ=== 210:其备选假设为：k H μμμ,,,:211 不全相同。

设从第i 个总体获得容量为i n 的样本观察值为i in i i y y y ,,,21 ，k i ,,2,1 =，各样本间还是相互独立的。

样本观察值ij y 可看成是来自均值为i μ的总体，这样ij y 就是其均值i μ与随机误差ij ε迭加而产生的。

上面我们已经假定在i A 水平下的ij y 服从),(2σμi N 分布，则有),0(~2σεN ij 。

因此，我们有单因素方差分析的统计模型：⎪⎩⎪⎨⎧==+=),0(,,,2,1,,,2,1,2σεεμN n j k i y ij i ij i ij 且都服从相互独立各 (25.1)为了能更仔细地描述数据，常在方差分析模型中引人一般平均与效应的概念。

称各个i μ的加权平均i ki i n n μμ∑==11(25.2)为总平均，其中∑==ki inn 1。

称k i a i i ,,2,1,=-=μμ(25.3)为因素A 在第i 水平的主效应，也简称为i A 的效应，同时也表明第i 个总体的均值是一般平均与其效应的迭加。

容易看出效应间有如下关系式：01=∑=iki i an(25.4)此时，单因素方差分析的统计模型可改写成包含效应的形式：⎪⎪⎩⎪⎪⎨⎧===++=∑=),0(,0,,2,1,,,2,1,21σεεμN a n n j k i a y ij k i i i iij i i ij 且都服从相互独立各(25.5)所要检验的原假设也可改写成：0:210====k a a a H现在，我们知道造成各ij y 间差异的原因可能有两个：一个可能是假设0H 不真，即各水平下总体均值i μ（或水平效应i a ）不同，因此从各总体中获得的样本观察值也就有差异了；另一可能是0H 为真，差异是由于随机误差引起的。

为了进一步定量分析这些差异，我们需要把这些差异表达出来。

由（25.1）可推导出：∙∙+=i i i y εμ(25.6)其中i n j iji n yy i/1∑=∙=，i n j ij i n i/1∑=∙=εε。

即组内样本观察值的平均值等于组内总体均值加上组内随机误差的平均值。

还可由（25.5）推导出：εμ+=y (25.7)其中n yy k i n j iji/11∑∑===，∑∑===k i n j ij i11εε。

即所有样本观察值的平均值等于总平均（各组均值的加权平均）加上所有随机误差的平均值。

这样，每一个观察值ij y 与总平均y 的偏差可以分解成两部分：)()(y y y y y y i i ij ij -+-=-∙∙(25.8)其中∙-i ij y y 称为组内偏差，由（25.1）和（25.6）代入得到：∙∙∙-=+-+=-i ij i i ij i i ij y y εεεμεμ)()((3.2.9)说明组内偏差仅仅反映了随机误差。

而y y i -∙称为组间偏差，由（25.6）、（25.7）和（25.3）代入得到：εεεμεμ-+=+-+=-∙∙∙i i i i i a y y )()((25.9)说明第i 组间偏差除了反映随机误差外还反映了第i 个水平的效应i a 。

各ij y 间总的差异大小可用总偏差平方和T S 表示：∑∑==-=k i n j ij T iy y S 112)((25.10)由（25.9）随机误差引起的数据间的差异可以用组内偏差平方和表示，也称误差偏差平方和e S ：∑∑==∙-=k i n j i ij e iy y S 112)((25.11)由于组间偏差除了随机误差外，还反映了效应的差异，故由于效应不同引起的数据差异可以用组间偏差平方和表示，也称因素A 的偏差平方和A S ：∑=∙-=ki i i A y y n S 12)((25.12)将表示总差异的平方和进行分解：Ae i ki i ki n j i ij i k i n j i ij ki n j i ki n j i ij k i n j i i ij k i n j ij T S S y y n y y y y y y y y y y y y y y y y S iii iii +=-+-=--+-+-=-+-=-=∙===∙∙==∙==∙==∙==∙∙==∑∑∑∑∑∑∑∑∑∑∑∑∑2111211112112112112)()()()(2)()()()( (25.13)其中0)(1=-∑=∙in j i ijy y。