方差分析过程.ppt

分析数据集中变量个数(2)例9-2为探索丹参对肢体缺血再灌注损伤的影响，将30只纯种新西兰实验用大白兔，按窝别相同、体重相近划分为10个区组。每个区组3只大白兔随机采用A、B、C三种处理方案，即在松止血带前分别给予丹参2ml/kg、丹参1ml/kg、生理盐水2ml/kg，在松止血带前及松后1小时分别测定血中白蛋白含量(g/L)，算出白蛋白减少量如下表9-6所示，问A、B两方案分别与C方案的处理效果是否不同？第30页,共60页，星期六，2024年，5月第31页,共60页，星期六，2024年，5月两因素方差分析（two-wayanova）分析：数据集中应该有3个变量，不妨命名为albumin、group、block。albumin记录白蛋白的减少量；group取值1、2、3，分别代表A,B,C三种方案；block取值为1-10，代表10个区组。第32页,共60页，星期六，2024年，5月方差分析数据集的建立技巧我们可以一个一个的输入变量的值。但是方差分析的数据集其变量取值有一定的规律，因此可以利用循环语句和判断语句来简化输入。第33页,共60页，星期六，2024年，5月例9-1数据的整理一个一个的输入变量的值，不仅麻烦，还容易出错。下面数据的输入是我们书上附带光盘的程序，录入有个小错误，你能找出么？第34页,共60页，星期六，2024年，5月DATAex9_1;INPUTxgrp@@;CARDS;5.619.516.018.719.215.013.515.818.0115.5111.8116.3111.8114.614.918.113.816.1113.2116.519.21-0.625.7212.824.12-1.82-0.126.3212.729.8212.622.025.627.027.924.326.427.025.423.1212.430.937.033.931.636.433.033.932.231.132.737.836.931.539.433.837.538.4312.236.03;第35页,共60页，星期六，2024年，5月例9-1数据的整理看不出吧！很多同学录入错误了，自己发现不了！举手问老师，哪里错了，你帮我找找？晕……………！加个print过程，很容易看出了！第36页,共60页，星期六，2024年，5月DATAex9_1;INPUTxgrp@@;CARDS;5.619.516.018.719.215.013.515.818.0115.5111.8116.3111.8114.614.918.113.816.1113.2116.519.21-0.625.7212.824.12-1.82-0.126.3212.729.8212.622.025.627.027.924.326.427.025.423.1212.430.937.033.931.636.433.033.932.231.132.737.836.931.539.433.837.538.4312.236.03;procprint;Run;第37页,共60页，星期六，2024年，5月例9_1:问治疗4周后，餐后2小时血糖下降值的三组总体平均水平是否不同？

解：由于是三组做比较，因此要用方差分析法。在建立数据集后可用ANOVA过程或GLM过程进行分析，同时给出选项进行各组均数的两两比较(此处用SNK法),两过程的输出格式基本一致。对于单因素的方差分析，在比较之前一般要考虑方差齐性的问题。这里可以用GLM(anova)过程中的HOVTEST选项实现，实际上这也是在SAS中做方差齐性检验的唯一途径。第38页,共60页，星期六，2024年，5月procanova;classgroup;modelbsvalue=group;meansgroup/hovtestsnk;run;procglm;