计算机sas方差分析

SAS单向分组资料的方差分析周美兰 教授湖南农业大学农学院单向分组资料的方差分析单向分组资料是指观察值仅按一个方向分组的资料示。所用的试验设计为完全随机试验设计。一、组内观察值数目相等的单向分组资料的方差分析二、组内观察值数目不等的单向分组资料的方差分析三、组内又分亚组的单向分组资料的方差分析分类一、组内观察值数目相等的单向分组资料的方差分析这是在 k组处理中，每处理皆含有 n个供试单位的资料如 表 6.1。在作方差分析时，其任一观察值的线性模型皆由 表示，方差分析如表 6.10。表 6.10 组 内 观 察 值 数目相等的 单 向分 组资 料的方差分析变 异来 源自由度DF平方和SS均 方MS F期望均方 EMS固定模型 随机模型处 理 间 k 1 MSt MSt/MSe误 差 k(n 1) MSe总变 异 nk 1例 作一水稻施肥的盆栽试验，设 5个处理， A和 B系分别施用两种不同工艺流程的氨水， C施碳酸氢铵， D施尿素， E不施氮肥。每处理 4盆 (施肥处理的施肥量每盆皆为折合纯氮 1.2克 )，共 54=20盆，随机放置于同一网室中，其稻谷产量 (克 /盆 )列于表 6.11，试测验各处理平均数的差异显著性。表 6.11 水稻施肥盆栽 试验 的 产 量 结 果处 理 观 察 值 (yij)(克 /盆 )A (氨水 1) 24 30 28 26 108 27.0B (氨水 2) 27 24 21 26 98 24.5C (碳酸 氢铵) 31 28 25 30 114 28.5D (尿素 ) 32 33 33 28 126 31.5E (不施 ) 21 22 16 21 80 20.0526 26.3(1) 自由度和平方和的分解总变异自由度 DFT=nk 1=54 1=19处理间自由度 DFt=k 1=5 1=4误差 (处理内 )自由度 DFe=k(n 1)=5(4 1)=15矫正数 (2) F测验将上述结果录入表 6.12表 6.12 表 6.11资 料的方差分析变 异来源 DF SS MS F F0.05 F0.01处 理 间 4 301.2 75.30 11.19* 3.06 4.89处 理内 (试验误 差) 15 101.0 6.73总 变 异 19 402.2假设 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)。p SSR0.05 SSR0.01 LSR0.05 LSR0.012 3.01 4.17 3.90 5.413 3.16 4.37 4.10 5.674 3.25 4.50 4.22 5.845 3.31 4.58 4.29 5.94表 6.13 多重比较时的 值计算表 6.14 施肥效果的 显 著性 (SSR测验 )处 理 平均 产 量(克 /盆 )差异 显 著性5% 1%尿 素 31.5 a A碳酸 氢铵 28.5 ab AB氨水 1 27.0 b AB氨水 2 24.0 b BC不 施 20.0 c C推断：根据表 6.14多重比较结果可知，施用氮肥 (A、 B、 C和 D)与不施氮肥有显著差异，且施用尿素、碳酸氢铵、氨水 1与不施氮肥均有极显著差异；尿素与碳酸氢铵、碳酸氢铵与氨水 1、氨水 1与氨水 2处理间均无显著差异。输入数据保存数据运行程序调用 SulutionnsAnalysis Analyst也可直接输入进入数据窗口进入分析窗口 直接输入数据调用数据Fileopen data By SAS Name （ Work SASUSER map）点击 Work SASUSER map选择数据库名StatisticsANOVAOne -Way ANOVE将处理变量点击到 dependent 框将数据变量点击到 independent 框确定多重比较方法方差分析结果多重比较结果调用 SulutionnsAnalysis Analyst