田间统计实验：第三节 方差分析.ppt

第三节 方差分析 不同的试验设计有其相应的线性数学模型，而方差分析正是根据这一线性数学模型进行的，因此所获数据的试验设计决定了其分析方法（即自由度和平方和的分解以及度量各效应是否显著的尺度）。 正是如此，方差分析的SAS程序中模型的确定是关键。以下结合教材内容顺序说明各种情况下的SAS程序编写方法。 DATA tb611; DO trt=1 TO 5; (或DO trt=”A”,”B”,”C”,”D”,”E”; ) DO r=1 To 4; INPUT y @@; OUTPUT; END; END; CARDS; 24 30 28 26 27 24 21 26 31 28 25 30 32 33 33 28 21 22 16 21 ; 在PROC ANOVA后可加上选项，包括DATA=输入数据集名和OUTSTAT=输出数据集名（用于存储方差分析结果）。 CLASS一句用于指明分类变量，此语句一定要设定，并且应出现在MODEL语句之前。 MODEL语句用于定义分析所用的线性数学模型，通常试验误差项（Error）一项不列。 MEANS语句用于计算效应平均数，并在“/”号后设定多重比较方法、显著水平 (二) 两向分组资料 1．组合内只有单个观察值的两向分组资料 （单因素随机区组或二因素完全随机设计） 本例用ANOVA过程，效应项包括生长素效应和区组效应以及省去的误差。对生长素效应分别用LSD法和DUNNETT氏最小显著差数法（DLSD法，默认对照是处理项中第一项）进行多重比较。 CARDS; 60 62 61 60 65 65 68 65 63 61 61 60 64 67 63 61 62 65 62 64 61 62 62 65 ; 2．组合内有重复的两向分组资料的方差分析（二因素完全随机设计） CARDS; 21.4 21.2 20.1 19.6 18.8 16.4 17.6 16.6 17.5 12.0 14.2 12.1 13.0 13.7 12.0 13.3 14.0 13.9 12.8 13.8 13.7 14.2 13.6 13.3 12.0 14.6 14.0 ; PROC ANOVA; CLASS a b; MODEL y= a | b; MEANS a b a*b / DUNCAN; RUN; (三) 数据转换 方差分析是在可加性，正态性和同质性的假定下进行的，为了满足这三项假定有时必须对原始数据进行数据转换。 以P126例6.15为例。 DATA ex61; DO trt=1 to 4; DO r=1 TO 6; INPUT y @@; yt=ARSIN(SQRT(y/100))*180/3.1415926; OUTPUT; END; END; DATA tb631; DO a=1 TO 3; DO b=1 TO 3; DO r=1 to 3; INPUT y @@; OUTPUT; END; END; END; “MODEL y=a | b;”是“MODEL y=a b a*b;”的简化形式，使用简化形式可使模型语句简洁，以下是一些简化表示法及其所代表的意义： a b a*b c a*c b*c a*b*c a b a*b c a*c b*c a b(a) c a*c b*c(a) a b(a) c a*c a c(b) a*c(b) a(b) c(b) a*c(b) a | b | c a | b | c @ 2 a | b(a) | c a | b(a) | c @ 2 a | c(b) a(b) | c(b) 意 义 简化示法 MEANS语句对肥料、土壤及两者互作的平均数进行多重比较。 [例6.15] 研究华农2号玉米花粉在不同贮藏条件下的生活力： （1）花粉盛于烧杯内，上盖纱布，藏于冰箱中； （2）花粉盛于烧杯内，置于干燥器中，藏于冰箱内； （3）花粉盛于烧杯内，在室温下贮藏。 经贮藏4小时后，在显微镜下检查有生活力花粉的百分数，对照为新鲜花粉。每处理检查了6个视野，其结果如表6.38。试作方差分析。 70 68 66 49 55 64 93 78 75 76 63 71