4-运用SAS进行方差分析.ppt

1、,SAS在方差分析中的应用,计算机在生物工程中的应用,上海应用技术学院 香料香精技术与工程学院,授课老师：王 一 非QQ:46478797,方差分析知识要点回顾,方差分析：利用试验数据，分析各个因素对某事物某指标的影响是否显著,单因素试验的方差分析,双因素等重复试验的方差分析,多因素等重复试验的方差分析,ANOVA过程和GLM过程,单因素方差分析也叫一维方差分析（One-Way ANOVA），或叫完全随机设计资料方差分析，它用于对单因素多个独立样本均值进行比较,多因素方差分析（MANOVA ， Multivariate Analysis Of Variance）是利用

2、方差分析方法处理多个自变量对因变量的影响效果的方法,如施肥量和播种量对产量的影响。,ANOVA过程和GLM过程,方差分析,ANOVA过程和GLM过程,平衡试验设计（ANOVA）,非平衡试验设计（GLM）,效率高,单因素方差分析实例,例 用二氧化硅50mg对大鼠染尘后，不同时期全肺湿重的变化。试比较染尘后1月、3月、6月，三个时期的全肺湿重有无差别(单因素三个水平，完全随机平衡设计的方差分析)? 1 3.3 3 4.4 6 3.6 1 3.6 3 4.4 6 4.4 1 4.3 3 3.4 6 5.1 1 4.1 3 4.2 6 5.0 1 4.2 3 4.7 6 5.5 1 3.3 3 4.2

3、 6 4.7,单因素方差分析实例,file-open-调入模版，写入数据,单因素方差分析实例,运行，生成SAS数据集,单因素方差分析实例,solutionsassistdata analysisANOVAanalysis of variance,table选用oneway1数据集,单因素方差分析实例,在columns中选用dependent（因变量） 选用weight，单击ok,单因素方差分析实例,在columns中选用clasification（分组变量） 选用time，单击ok,单因素方差分析实例,单击additional options（添加选项）,单因素方差分析实例,选用output

4、statistics（输出统计量）,单因素方差分析实例,选择means（输出统计量）,单因素方差分析实例,选择time,单因素方差分析实例,单因素方差分析实例,选择options for means,选择comparison tests，得到14种比较方法,单因素方差分析实例,本例选用1，2，11，12，检验法，单击ok,单因素方差分析实例,选用 specify comparison options（指定比较选择）,单因素方差分析实例,单因素方差分析实例,单因素方差分析实例,比较选择 comparison options包括 （1）Type1/Type2 seriousness ratio (

5、Water-Duncan) 瓦尔-邓肯检验第一类错误/第二类严肃比率 Type1/Type2 seriousness ratio的比值 50 0.1 100 0.05 500 0.001的显著水平 SAS预设的比值为100即预设的显著性水平为0.05,单因素方差分析实例,比较选择 comparison options包括 （2）Level of significance for comparison of means 平均值多重比较的显著性水平，由用户自行选定 Confidence intervals for differences between means 均值间差异的可信区间 Confi

6、dence intervals for means 各均值的可信区间 Do not sort means into descending order -不必以递减顺序将均值排序显示 Indicate nonsignficant subsets of means -指出均值未到显著差异的各分组,单因素方差分析实例,单因素方差分析实例,运行，得到结果,单因素方差分析实例,单因素方差分析实例,结果分析 经4种方法比较检验，只有1月与6月的全肺湿重（weight）有显著差异（p0.05）,单因素方差分析实例,以上计算，可以直接用编程方法进行，得到结果相同：,单因素方差分析实例,输出说明 CLASS L

7、EVEL INFORMATION分类水平。其中包括： CLASS CLASS语句中列出的对应名 LEVELS 因素对应的水平数 VALUES 因素对应各水平的值或标记 SOURCE 变异来源 SUM OF SQUARES(SS) 平方和 MEAN SQUARE(MS) 均方 F VALUE F值,单因素方差分析实例,输出说明 PRF 显著水平 MODEL模型 它的平方和等于各因素效应的平方和之和，其均方等于它的平方和除以自由度 ERROR 误差 CORRECTED TOTAL 校正总变异 R-SQUARE R2，其值为模型的平方和除以校正总平方和。一般来说，R2值越大，模型拟合数据越好 C.V

8、 变异系数。样本的变异系数为该样本的标准差除以均值，表示单位量的变异。 ROOT MSE 误差均方根，是观察变量的标准差的估计值,几何均数的单因素方差分析,科学研究中，对于多组成等比数列的资料，通常要对原始数据进行对数转换后才能进行方差分析 例 某市对四个年龄段的健康男子接种某疫苗后，测定抗体滴度，问各组的抗体平均滴度（即几何均值）之间有无显著性差异？,几何均数的单因素方差分析,fileopen调入oneway2.sas文件 submit，生成SAS数据集：work.oneway2 solutions-assistdata analysisANOVA -analysis of variance

9、 往下的选项，类似于例1,几何均数的单因素方差分析,几何均数的单因素方差分析,几何均数的单因素方差分析,几何均数的单因素方差分析,几何均数的单因素方差分析,几何均数的单因素方差分析,几何均数的单因素方差分析,多因素方差分析,多因素方差分析（MANOVA, Multivariate Analysis Of Variance）是利用方差分析方法处理多个自变量对因变量的影响效果的方法。,多因素方差分析,例 已知3组贫血病患者，其血红蛋白浓度（x1，%）及红细胞计数（x2，104/mm3）的数据如下，试做多因素方差分析： x1 x2 x1 x2 x1 x2 a 3.9 210 b 4.8 270 c

10、4.4 250 a 4.2 190 b 4.7 180 c 3.7 305 a 3.7 240 b 5.4 230 c 2.9 240 a 4.0 170 b 4.5 245 c 4.5 330 a 4.4 220 b 4.6 270 c 3.3 230 a 5.2 230 b 4.4 220 c 4.5 195 a 2.7 160 b 5.9 290 c 3.8 275 a 2.4 260 b 5.5 220 c 3.7 310 a 3.6 240 b 4.3 290 a 5.5 180 b 5.1 310 a 2.9 200 a 3.3 300,多因素方差分析,建立SAS语言程序命令：gl

11、m1.sas,多因素方差分析,submit, get results,多因素方差分析,多因素方差分析,多因素方差分析,多因素方差分析,多因素方差分析,Homeworks,下次课交，由课代表收齐，按学号排好先后次序 总体要求：A4打印，题目不打印，正反面打印 排版规范、清楚,Homeworks,1、一批由同种原料织成的布，用不同的染整工艺处理，每台进行缩水率试验，目的是考察不同的工艺对布的缩水率是否有显著影响.这里采用了五种不同的染整工艺，每种工艺处理四块不养，测得缩水率百分表如下。请运用SAS软件的单因素分析方法（Duncans multiple-range test和Fishers Least-significant-difference test）确认染整工艺对缩水率影响是否显著（=0.05）. 要求：（1）需把SAS运行结果粘帖上（SAS每次运行都会显示时间，所以每个人的运行结果时间如果一样，就是抄袭），编辑合理的字体大小和行间距，节省纸张。 （2）运用Origin软件绘制柱状图，做好标注。 （3）根据第一次的作业讲解，更加合理的运用专业术语描述试验结果，如果采用英文回答，作业加分。,Homewo