第三节拉丁方设计报告.pptVIP

第三节 拉丁方设计与方差分析 用r个拉丁字母排成r行r列的方阵，使每行、每列中每个字母都只出现一次。这样的方阵叫r阶拉丁方或r× r拉丁方（Latin square）。按拉丁方的字母、行和列安排处理及影响因素的试验称为拉丁方试验。 拉丁方设计（Latin square design）是随机区组设计的进一步扩展，可以是考核3个处理因素，也可用于考核1或2个处理因素，而同时试图控制两种非处理因素，其中处理因素、行单位组因素、列单位组因素的水平均为r。 与随机区组设计比较，拉丁方设计控制了两个非处理因素。进一步缩小了实验误差。可以得到比随机区组设计更多一个项目的均衡，因而误差更小，效率更高。拉丁方设计的优点是可以大大减少试验次数，尤其适合于动物实验和实验室研究。缺点是要求处理数必须等于拉丁方的行（列）数，一般的试验不满足此条件，而且数据缺失会增加统计分析的难度。 拉丁方设计的设计要求有： 1.必须是三个因素的试验，且三个因素的水平数相等； 2.行间、列间、处理间均无交互作用； 3.各行、列、处理的方差齐。 一、拉丁方设计分组 例4-5 为比较五种防护服对脉搏数的影响，安排五个受试者在五个不同日期穿五种防护服测量脉搏数，试进行拉丁方设计。 （1）根据处理数在本书后面的附表14中选定5×5基本（标准）拉丁方。 （2）将拉丁方随机化：① 随机排列拉丁方的行。例如，读取4个两位数的随机数，66、05、32、88，则R=3、1、2、4，即先将3、1行对调，然后将2、4行对调。 ② 随机排列拉丁方的列。例如，读取4个两位数的随机数，53、85、39、06，则R=3、4、2、1，即先将3、4列对调，然后1、2列对调。 表4-9 拉丁方设计随机分配结果 （3）随机分配处理。例如，读取5个两位随机数10、28、81、47、20，则R=1、3、5、4、2，于是有A（甲）、B（丙）、C（戊）、D（丁）、E（乙）。将上述最后一个拉丁方的行、列和拉丁字母分别对应于试验日期、受试者和防护服的最终试验方案见表4-9。 二、方差分析步骤 例4-6 为研究注射甲状腺素对甲状腺体的影响，以豚鼠5个种系，每个种系各5只。分养于5个笼子。每笼内置放各种系豚鼠1只，并以甲状腺素的5个不同剂量分别予以注射。以上剂量、种系和笼子三个因素的分组如下，全部数据列于表4-10。试分析不同剂量甲状腺素组均数之间的差别。 表4-10 豚鼠注射不同剂量甲状腺素后的甲状腺 体重量(mg) 剂量：A B C D E，处理数=5 种系：甲 乙 丙 丁 戊，行数=5 笼子：Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ，列数：5 （1）建立假设、确定检验水准α。 剂量间：Ho：剂量组间的总体均数相同。 （处理）H1：剂量组间的总体均数不相等或不全相等， α=0.05。 种系间：Ho：种系间的总体均数相同。 （行间）H1：种系间的总体均数不相等或不全相等， α=0.05。 笼子间：Ho：5个笼子组间的总体均数相同 （列间）H1：5个笼子组间的总体均数不相等或不全相等， α=0.05。 （2）计算检验统计量（F值） 1）离差平方和的分解：根据变异来源，拉丁方设计资料总的离差平方和（SS总）可分解为SS处理、SS列和SS行及SS误差四部分。且 SS总=SS处理+SS列+SS行+SS误差 2）列方差分析表，计算各离差平方和SS，自由度υ，均方MS和F值：其中Xk为第k种处理小计，Xi为第i行小计，Xj为第j列小计，Xij为第i行第j列观察值，校正数C=（∑Xij）2/r2，r为拉丁方的阶，即行数、列数和处理数。 各部分离差平方和、自由度、均方、F值的计算与随机单位组设计的一样，本例： SS总=120719－115736.04=4982.96 （3）查F值表，确定P值，下结论。本例υ处理=4，υ误差=12，查附表10（F界值表）得，F0.05(4，12)=3.26,F0.05(4,12)=5.67。因F处理＞F0.01(3,12)，故P＜0.01。同理，种系间、笼子间P＞0.05。 表4-12 例3.9资料方差分析结果表 结论：对剂量（处理）间差别的检验，在α=0.05水准处，拒绝Ho，接受H1；对种系（行）间差别、笼子（列）间差别的检验，在α=0.05水准处，均不拒绝Ho，接受Ho。故认为不同剂量甲状腺素组的甲状腺体重量总体均数不等。不同种系间无差别，不同笼子间亦无差别。 若想进一步了解哪两组间有差别，可进行多个均数的两两比较。 第四节 配对实验设计与分析