06多因素随机区组试验的统计探究.docVIP

二、多因素随机区组试验资料的统计分析 如果试验中有两个考察因素，因素A有a个水平，因素B有b个水平，共有ab个处理组合，在试验中安排了r个随机区组，那么，此类资料观察值的数学模型为： （i＝1，2，…，a；j＝1，2，…，b；k＝1，2，…，r） 其中为第ij处理组合的效应值，为第k区组的效应值，，为第i个A水平的效应值，为第j个B水平的效应值，为Ai与Bj之间的交互作用。从数学模型看出，对于这类资料的分析可以分两步进行： 先按模型将总变异分解为区组间变异、处理(组合)间变异和误差变异。这是一个两向分类的方差分析，方差分析表如表8.15所示。该表与表8.14有点类似，只是在表8.14中的处理是指一个因素的一个水平()，而这里的处理是指多因素各水平交错配合成的一个组合()。如果这一步中对处理间差异的F测验不显著，就不必进行第二步。分析就此结束。 表8.15 第一步：将总变异分解为区组间变异、处理(组合)间变异和误差变异 变异来源 自由度 平方和 均方 F值 区组间 dfr＝r－1 MSr MSr/ MSe 处理间 dft＝ab－1 MSt MSt/ MSe 误差 dfe＝(r－1) (ab－1) MSe 总变异 dfT＝rab－1 2. 如果第一步对处理间差异的F测验显著，应进一步按模型将处理变异分解为A因素各水平间的变异、B因素各水平间的变异和交互作用(AB)引起的变异。方差分析表如表8.16所示。在第二步中的F测验值如何计算需要视研究目的所确定的模型而定。表8.17列出了总的方差分析表和各种模型中的期望均方。 表8.16 第二步：将处理间变异分解为A间变异、B间变异和交互作用 变异来源 自由度 平方和 均方 F值 A间 dfA＝a－1 MSA 视模 型而 确定 B间 dfB＝b－1 MSB AB互作 dfAB＝(a－1) (b－1) MSAB 处理间 dft ＝ab－1 MSt MSt/ MSe 表8.17 两因素随机区组试验资料的方差分析表 变异 来源 自由度 平方和 均方 F值 期望均方(EMS) 固定模型 随机模型 A固定B随机 区组间 dfr SSr MSr 视 模 型 而 确 定 或 处理间 dft SSt MSt A间 dfA SSA MSA B间 dfB SSB MSB AB互作 dfAB SSAB MSAB 误差 dfe SSe MSe 总变异 dfT SST 例8.5 考察三种生长素Ai (i＝1,2,3)和两种葡萄糖浓度Bj (j＝1,2)所配成的6种培养基对香草兰胚状体发育的影响。将接种好的培养皿放置在四个培养箱(k＝Ⅰ,Ⅱ,Ⅲ,Ⅳ)中，按随机区组排列，每培养箱为一个区组。所得数据如表8.18所示。如果除区组效应不需考察外，其余效应都是固定效应，试对资料进行适当的分析。 表8.18 三种生长素和两种葡萄糖浓度所配成的6种培养基对香草兰胚状体发育的影响 生长素(A) 糖浓度(B) 区组(培养箱k) Ⅰ Ⅱ Ⅲ Ⅳ A1 B1 8 9 6 9 32 262 1024 8 B2 5 6 8 5 24 150 576 6 A2 B1 4 3 3 2 12 38 144 3 B2 4 5 3 8 20 114 400 5 A3 B1 5 9 6 8 28 206 784 7 B2 3 2 5 6 16 74 256 4 29 34 31 38 132 844 3184 5.5 155 236 179 274 844 841 1156 961 1444 4402 4.83 5.70 5.17 6.33 5.50 分析过程分为两个步骤： 第一步：按两向分类资料方差分析的方法，将总变异分解为区组间变异、处理间变异和试验误差。根据表8.18的整理结果可以得到 总自由度　dfT ＝观察值总数－1＝3×2×4－1＝24－1＝23 区组间自由度　dfr ＝区组数－1＝4－1＝3 处理间自由度　dft ＝处理数－1＝3×2－1＝6－1＝5 误差自由度 dfe ＝ dfT －dft－dfr ＝23－3－5＝15 矫正项 C.T. ＝ 观察值总和的平方／观察值总数目＝1322／24＝726 总平方和 SST ＝各观察值平方之和－C.T.＝844－726＝118 区组间平方和 SSr ＝各区组和的平方之和／每区组的观察值数目－C.T. ＝4402／6－726＝7.6667 处理间平方和 SSt ＝各处理和的平方之和／每处理的观察值数目－C.T