《Ch7方差分析》PPT课件

《《Ch7方差分析》PPT课件》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、方差分析AnalysisofVariance一、方差分析的数据格式要求统计软件对方差分析的数据格式要求与对两个样本均数比较的t检验的要求是完全相同的：应变量(Y)和不同的处理因素(Factors)均应以单独的变量(Variables)在数据集中体现，某个处理因素的各个水平就是其对应变量的取值。例如：考查3种解毒药A、B、C的效果，同时设立一个空白对照D。每组6只大白鼠，用药一段时间后测定其血中胆碱酯酶含量(Y)。（《医学统计学》余松林主编，p92）此例方差分析的数据结构是怎样的？A组，6个数据B组，6个数据C组，

2、6个数据D组，6个数据如果只有两组，就是t检验的格式二、方差分析的SAS过程方差分析由SAS/STAT模块来完成，常用的有ANOVA过程和GLM过程。前者运算速度较快，但主要运用于平衡数据，即每个组内的非缺失观测值数目应相同；后者运算速度较慢，但功能强大，既可用于平衡数据也可用于非平衡数据。ANOVA与GLM的常规用法几乎完全相同，因此，这里就合并介绍。提示：要深入了解它们，请自行参阅SAS的帮助或相关参考文献。ANOVA、GLM过程的常用语句：PROCANOVA;CLASSvariables;

3、MODELdependents=effects;MEANSeffects;PROCGLM;CLASSvariables;MODELdependents=effects;MEANSeffects;以上两个过程中带下划线的语句是必需的，其余语句是可选的。注：此处仅仅列出了这两个过程的常用语句，其它所有相关语句请参阅SAS帮助及相关文献。1、CLASS语句CLASS语句指明需要在MODEL语句中使用的“考查因素”，主要是处

4、理因素，其类型应是分类变量；CLASS语句应该出现在MODEL语句之前。2、MODEL语句MODEL语句指定分析的模型，是方差分析中至关重要的一条语句，指定模型方式为：“结果变量=考查因素”，注意它们都是变量名。例如：CLASSgrp;MODELy=grp;考查一个因素grp对y的影响CLASSab;两个分类变量（考查因素）MODELy=aba*b;模型中包含三部份：a因素对y的影响；b因素对y的影响；a与b两个因素在y上面的交互作用。3、MEANS语句MEANS语句指定要进行多重比较的考查因素及比较方法多重比较

5、的方法有：BON（Bonferroni检验）、DUNCAN（新复极差法）、DUNNETT、GABRIEL、LSD（最小显著差法）、SCHEFFE、SIDAK、SNK（q检验）、TUKEY、WALLER最为常用的有两个：SNK法（Student-Newman-Keulsmultiplerangetest）DUNNETT：多个试验组与一个共用对照组的比较，使用时需指明对照组是谁。例如：MEANSgrp/SNK;/*grp是要进行多重比较的因素，方法SNK*/MEANSgrp/DUNNETT(‘Control’);/*

6、Control是对照的取值*/MEANSgrp/DUNNETT(‘1’);/*1是对照的取值*/三、不同设计类型资料的方差分析例6-1：考查3种解毒药A、B、C的效果，同时设立一个空白对照D。每组6只大白鼠，用药一段时间后测定其血中胆碱酯酶含量(Y)，结果如下表。（《医学统计学》余松林主编，p92）试问不同解毒药的效果有没有差异？（一）完全随机设计资料的方差分析首先建立SAS数据集建立数据集的方式比较灵活，建议选择熟悉的一种。以下程序供参考dataex61a;inputgrp$@@;doi=1to6;inputy

7、@@;output;end;cards;A231218162814B283123242834C142417191622D81221191415;run;dataex61b;inputy@@;grp=int((_n_-1)/6);cards;23121816281428312324283414241719162281221191415;run;dataex61c;inputygrp$@@;cards;23A12A18A16A28A14A28B31B23B24B28B34B14C24C17C19C16C22C8D12

8、D21D19D14D15D;run;一行数据中有一个分组标识数据中没有分组标识，需要用程序产生每一个数据都有分组标识调用GLM过程进行方差分析PROCGLMDATA=ex61a;CLASSgrp;/*处理因素*/MODELy=grp;/*一个处理因素对Y的影响*/RUN;QUIT;/*退出GLM过程，否则处于运行状态*/输出结果如下：分析：此处需要分析的变量（应变量、结果