ch5-方差分析I-完全随机化析因设计

《ch5-方差分析I-完全随机化析因设计》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、上海财经大学统计与管理学院方差分析I-平衡设计问题和数据描述研究高血压治疗方案效果比较方法介绍方差分析、多重比较方差齐性检验参考资料汪仁官，陈荣昭（译），蒙哥马利（著），1998，实验设计与分析，中国统计出版社。SAS分析tabulate过程anova过程glm过程gplot过程版权所有：吴纯杰制表方差分析（平衡设计）广义线性模型(不平衡数据的方差分析)作效应图数据描述数据描述几种高血压治疗方案效果研究。数据集：hyper变量：drug-X,Y,Z三种治疗药物biofeed-心理护理Present有，Absent无diet-特殊膳食Y有，N无

2、问题提出几种治疗方案是否有差异？治疗、护理和膳食三因子是否存在交互效应？上海财经大学统管院数据描述上海财经大学统管院图：hyper数据集数据录入和整理datahyper;/*Chapter5*/infile'd:hypertension.dat';inputn1-n12;if_n_<4thenbiofeed='P';elsebiofeed='A';if_n_in(1,4)thendrug='X';if_n_in(2,5)thendrug='Y';if_n_in(3,6)thendrug='Z';arraynall{12}n1-n12;doi

3、=1to12;ifi>6thendiet='Y';elsediet='N';bp=nall{i};cell=drug

4、

5、biofeed

6、

7、diet;output;end;dropin1-n12;run;上海财经大学统管院注：原始数据为6行×12列条数据，即共有6条行观测；其中1-3行对应着进行心理护理，而4-6行则无；而三种药物分别为X在1和4行，Y在2和5行，Z在3和6行。每行的12条治疗后的血压数据使用数组nall存储。通过循环将每行12条数据转换成12行血压数据，并且前6条对应辅助特殊膳食方案，后6条则无。单元cell存储不同水平的组合

8、，共12种。注：_n_代表观测值的序号；

9、

10、为字符串连接命令；output将读入的数据和产生的新变量数据录入到当前数据集；drop舍弃录入数据所产生的辅助变量i和数组n1-n12。带交互效应的方差分析模型上海财经大学统管院模型检验(F检验)独立性正态性方差齐次性制表-tabulate过程procprintdata=hyper;run;proctabulatedata=hyper;classdrugdietbiofeed;varbp;tabledrug*diet*biofeed,bp*(meanstdn);run;注：对血压bp按照drug,d

11、iet,Biofeed水平分类，并在表格中显示各分类的均值、标准差和数据个数。思考：过程语句添加选项order=data，结果有何不同？上海财经大学统管院方差齐性检验procanovadata=hyper;classcell;modelbp=cell;meanscell/hovtest;run;注：对12个水平组合的方差通过/hovtest选项进行齐性检验；此处方差分析仅对不同水平处理的差异性进行检验。上海财经大学统管院Cell方差分析结果上海财经大学统管院模型显著单元显著差异cell方差齐性检验结果上海财经大学统管院p值>0.05,未通过方

12、差齐性检验，说明各水平组合观察值的方差相等hyper方差分析procanovadata=hyper;classdietdrugbiofeed;modelbp=diet

13、drug

14、biofeed;meansdiet*drug*biofeed;odsoutputmeans=outmeans;run;procprintdata=outmeans;run;procsortdata=outmeans;bydrug;run;注：model语句等价于modelbp=dietdrugdiet*drugbiofeeddiet*biofeeddrug*biofe

15、eddiet*drug*biofeed上海财经大学统管院方差分析模型：做三因子主效应以及所有可能的交互效应模型（包括一阶或二阶三水平交互）。运算符

16、为逻辑或操作。means语句生成三个水平不同组合即cell的均值和方差等描述性结果。odsoutput为通过输出传递系统将means语句所得结果输出到outmeans数据集。方差分析结果上海财经大学统管院方差分析结果上海财经大学统管院结果显示：三因子diet,drug和biofeed的主效应显著（水平=0.05)，其中一阶两水平交互效应并不显著，但三因子交互效应显著，可能是由于一因子固定前提下，另

17、外两因子交互效应存在很大差异造成。outmean输出结果上海财经大学统管院按drug排序输出绘制交互效应图procsortdata=outmeans;bydrug;