田间试验的布置与管理.DOCVIP

PAGE PAGE 19 第六章 方差分析 PAGE 20 第六章 方差分析 第一节 方差分析的基本原理 上章介绍了1个或两个样本平均数的假设测验方法。本章将介绍k(k≥3)个样本平均数的假设测验方法，即方差分析(analysis of variance)。方差分析就是将总变异剖分为各个变异来源的相应部分，从而发现各变异原因在总变异中相对重要程度的一种统计分析方法。其中，扣除了各种试验原因所引起的变异后的剩余变异提供了试验误差的无偏估计，作为假设测验的依据。因而，方差分析象上章的t测验一样也是通过将试验处理的表面效应与其误差的比较来进行统计推断的，只不过这里采用均方来度量试验处理产生的变异和误差引起的变异而已。方差分析是科学的试验设计和分析中的一个十分重要的工具。本章将在介绍方差分析基本原理和方法的基础上进一步介绍数学模型和基本假定。 一、自由度和平方和的分解 方差是平方和除以自由度的商。要将一个试验资料的总变异分解为各个变异来源的相应变异，首先必须将总自由度和总平方和分解为各个变异来源的相应部分。因此，自由度和平方和的分解是方差分析的第一步。下面先从简单的类型说起。设有k组数据，每组皆具n个观察值，则该资料共有nk个观察值，其数据分组如表6.1。 表6.1 每组具n个观察值的k组数据的符号表 组别 观察值(，i=1，2，…，k；j=1，2，…，n) 总和 平均 均方 1 … … 2 … … … … i … … … … k … … 在表6.1中，总变异是nk个观察值的变异，故其自由度，而其平方和则为： （6·1） (6·1)中的C称为矫正数： (6·2) 这里，可通过总变异的恒等变换来阐明总变异的构成。对于第i组的变异，有 EMBED Equation.3 总变异为第1，2，…，k组的变异相加，利用上式总变异(6·1)可以剖分为: （6·3） 即 总平方和=组内(误差)平方和+处理平方和 组间变异由k个的变异引起，故其自由度,组间平方和为： （6·4） 组内变异为各组内观察值与组平均数的变异，故每组具有自由度和平方和；而资料共有组，故组内自由度,组内平方和为： (6·5) 因此，得到表6.1类型资料的自由度分解式为： (6·6) 总自由度DFT=组间自由度DFt+组内自由度DFe 求得各变异来源的自由度和平方和后，进而可得: (6·7) 若假定组间平均数差异不显著(或处理无效)时，(6·7)中与是的两个独立估值，均方用表示，也用表示，两者可以互换。其中组内均方也称误差均方，它是由多个总体或处理所提供的组内变异(或误差)的平均值。 [例6.1] 以A、B、C、D 4种药剂处理水稻种子，其中A为对照，每处理各得4个苗高观察值(cm)，其结果如表6.2，试分解其自由度和平方和。 表6.2 水稻不同药剂处理的苗高(cm) 药 剂 苗高观察值 总和 平均 A 18 21 20 13 72 18 B 20 24 26 22 92 23 C 10 15 17 14 56 14 D 28 27 29 32 116 29 T=336 21 根据(6·6)进行总自由度的剖分： 总变异自由度DFT=(nk-1)=(4?4)-1=15 药剂间自由度DFt=(k-1)=4-1=3 药剂内自由度DFe=k(n-1)=4?(4-1)=12 根据(6·3)进行总平方和的剖分： 或 或 药剂A内： 药剂B内： 药剂C内： 药剂D内： 所以 误差平方和也可直接计算。 进而可得均方： 以上药剂内均方系4种药剂内变异的合并均方值，它是表6.2资料的试验误差估