计算机SAS方差分析.ppt

SAS单向分组资料的方差分析 周美兰 教授 湖南农业大学农学院 输入数据 保存数据 运行程序 调用Sulutionns→Analysis →Analyst也可直接输入进入数据窗口 进入分析窗口→直接输入数据 调用数据 File→open data By SAS Name →（Work SASUSER map） 点击Work SASUSER map 选择数据库名 Statistics→ANOVA→One-Way ANOVE 将处理变量点击到dependent 框将数据变量点击到independent 框 确定多重比较方法 方差分析结果 多重比较结果 调用Sulutionns→Analysis →Analyst 进入表格输入数据（教材P25表27） Statistics→Descriptive→Summary Statistics… 描述性统计分析 产量点如Analysis,因素点入class选择统计方法 * * 单向分组资料的方差分析 单向分组资料是指观察值仅按一个方向分组的资料示。所用的试验设计为完全随机试验设计。 一、组内观察值数目相等的单向分组资料的方差分析 二、组内观察值数目不等的单向分组资料的方差分析 三、组内又分亚组的单向分组资料的方差分析 分类 一、组内观察值数目相等的单向分组资料的方差分析 这是在k组处理中，每处理皆含有n个供试单位的资料如表6.1。 　 在作方差分析时，其任一观察值的线性模型皆由　　　 　　　　　　　 表示，方差分析如表6.10。 表6.10 组内观察值数目相等的单向分组资料的方差分析 nk－ 1 总变异 MSe k(n－ 1) 误 差 MSt/MSe MSt k－1 处理间 随机模型 固定模型 期望均方EMS F 均 方 MS 平方和 SS 自由度 DF 变 异 来 源 [例] 作一水稻施肥的盆栽试验，设5个处理，A和B系分别施用两种不同工艺流程的氨水，C施碳酸氢铵，D施尿素，E不施氮肥。每处理4盆(施肥处理的施肥量每盆皆为折合纯氮1.2克)，共5×4=20盆，随机放置于同一网室中，其稻谷产量(克/盆)列于表6.11，试测验各处理平均数的差异显著性。 表6.11 水稻施肥盆栽试验的产量结果 26.3 526 20.0 80 21 22 16 21 E (不施) 31.5 126 32 33 33 28 D (尿素) 28.5 114 31 28 25 30 C (碳酸氢铵) 24.5 98 27 24 21 26 B (氨水2) 27.0 108 24 30 28 26 A (氨水1) 观察值(yij)(克/盆) 处 理 (1) 自由度和平方和的分解 总变异自由度DFT=nk－1=5×4－1=19 处理间自由度DFt=k－1=5－1=4 误差(处理内)自由度DFe=k(n－1)=5×(4－1)=15 矫正数 (2) F测验 将上述结果录入表6.12 表6.12 表6.11资料的方差分析 402.2 19 总 变 异 6.73 101.0 15 处理内(试验误差) 4.89 3.06 11.19** 75.30 301.2 4 处 理 间 F0.01 F0.05 F MS SS DF 变异来源 假设H0: ，HA: 不全相等。 为了测验H0，计算处理间均方对误差均方的比率，算得 F =75.3/6.73=11.19 查F表当 v1=4，v2=15时，F0.01=4.89，现实得F=11.19F0.01，故否定H0， 推断这个试验的处理平均数间是有极显著差异的。 (3) 各处理平均数的比较 算得单个平均数的标准误 根据 =15，查SSR表得p=2，3，4，5时的SSR0.05 与SSR0.01值，将 值分别乘以SE值，即得 值，列于表6.13。进而进行多重比较(表6.14)。 5.94 4.29 4.58 3.31 5 5.84 4.22 4.50 3.25 4 5.67 4.10 4.37 3.16 3 5.41 3.90 4.17 3.01 2 LSR0.01 LSR0.05 SSR0.01 SSR0.05 p 表6.13 多重比较时的 值计算 表6.14 施肥效果