实验六 拉丁方试验设计.docVIP

实验六 拉丁方实验设计 实验目的 了解拉丁方实验设计的基本方法及数据的分析方法。 实验工具 Spss中的Analyze (General linear Model (Univariate。 知识准备 一、拉丁方设计的概念 将k个不同符号排成k列，使得每一个符号在每一行、每一列都只出现一次的方阵，叫做k×k拉丁方。利用拉丁方阵进行实验设计的方法叫做拉丁方设计。最初设计方案时用拉丁字母组成的方阵来表示后来尽管方阵中的元素改用了阿拉伯数字或其它的符号人们仍称这种方案为拉丁方。如图为×4拉丁方，它的每一行和每一列都是一个区组或一次重复，而每一个处理在每一行或每一列都只出现一次，因此，它的处理数、重复数、行数、列数都等于。 重复数＝处理数＝数＝横行数每个处理在横行的区组内或的区组内都能出现一次，从两个方向都可看成重复，排列呈方形两个方向的排列都是随机的，从两个方向进行局部控制，试验精确度较高。缺点：处理数＝重复数，若处理过多，重复随之增多，使验工作量过大。一般不宜超过个处理。若处理数过少，方差分析时的自由度过小，影响分析结果的精确性。，缺乏灵活性。（1），。在进行拉丁方设计时，首先要根据实验处理数k从标准方表中选定一个k×k的标准方。例如处理数为5，需要选一个5×5的标准方，如图。随后我们要对选定的标准方的行、列和处理进行随机化排列。本例处理数5，因此根据随机数字表任选一页中的一行，除去0、6和重复数字，满5个为一组，要得到这样的3组5位数。假设得到的3随机数字为14325，53124，41235。 （2）列随机第一组5个数字14325调整列顺序，即把第4列调至第2列，第2列调至第4列，其余列不动。如图（3）行随机第二组5个数字53124调整行顺序，即把第5行调至第1行，第3行调至第2行，第1行调至第3行，第2行调至第4行，第4行调至第5行。如图（4）处理随机将处理的编号按第三组5个数字41235的顺序进行随机排列。即4号=A，1号=B，2号=C，3号=D，5号=E。因此经过随机重排的拉丁方中A处理用4，B处理用1，C处理用2，D处理用3，E处理用5。如图 其中，μ为总均值； 为横行效应；为纵列效应；为处理效应；为实验误差，满足正态分布N(0,σ2)。如果是固定效应模型，则满足；满足。 实验的结果中，处理数=横行区组数=纵列区组数=重复次数。这样，实验有个处理，便有个观测值。方差分析时，从总变异方差中除分解出处理间方差和误差项方差外，还可分解出纵横两个区组的方差，这就使误差项方差进一步减小。所以拉丁方实验的精确度比随机区组实验更高。 1.平方和与自由度的分解 总平方和为：，， 。 横行平方和为： ，，为横行数据的和。 纵列平方和为：，，为纵列数据的和。 处理平方和为：，，为各处理的数据和。 误差平方和为：，。 总平方和 = 横行平方和 + 纵行平方和 + 处理平方和 + 误差平方和 总自由度 = 横行自由度 + 纵行自由度 + 处理自由度 + 误差自由度 2．F检验 ，， 对于给定的显著性水平，根据计算得到的F统计量的值和临界值做对比可以得到结论。 如果方差分析的结果有显著的差异，进一步可以进行多重比较。多重比较的方法见前面的内容。 实验背景 对玉米的五个品种进行拉丁方实验，产量（kg）如下表所示，试作方差分析。 表6.12 纵列横行 Ⅱ Ⅲ Ⅳ Ⅴ Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ A 14 D 19 B 23 C 21 E 23 E 22 B 21 A 15 D 18 C 16 D 20 A 16 C 20 E 24 B 23 C 18 E 23 D 18 B 21 A 17 B 25 C 18 E 23 A 17 D 20 实验过程 （1）输入数据,如图6.50所示。 图6.50 （2）单击Analyze (General linear Model (Univariate，打开Univariate主对话框。如图6.51所示： 图6.51 选择要分析的变量“产量”进入Dependent Variable 框中，选择因素变量“行、列和因素”进入Fixed Factor框中。 （3）单击Model... 按纽选择分析模型，得到Model对话框。如图6.52 所示，在Specify Model框中，指定模型类型。 图6.52 本例选择Custom选项。从左边Factors Covariates框中选择变量“因素、行、列”进入Model框中，单击Build Term下面的小菜单，选择主效应。最后在Sum of Square 中选择分解平方和的方法，选取默认项TypeⅢ。本例中模型包含截距项。单击Continue 返回在主对话框。