实验八-方差分析2.pptVIP

Foobar 实验八 随机区组设计资料 的方差分析 目的与要求： 1.掌握利用anova过程或glm过程进行随机区组设计方差分析， 能对SAS程序的输出结果作出合理解释。 2.熟悉循环do语句的用法。 一、方差分析的基本思想 根据资料的设计类型和研究目的，把全部观察值的变异（总变异）分解为两个或多个部分，每部分可以用某因素的作用来解释，将某因素解释的变异和误差变异进行比较，作出某因素是否有统计学意义的结论。 二、随机区组设计的方差分析表 三、常用两两比较方法： SNK法 Bonfferoni法 Dunnett法 四、方差分析所用的过程 ANOVA过程(Analysis Of Variance) GLM过程(General Linear Model) PROC ANOVA的过程格式 PROC ANOVA； CLASS 变量表； MODEL 因变量表＝效应； MEANS 效应〈/选择项〉； PROC GLM的过程格式 PROC GLM； CLASS 变量表； MODEL 因变量表＝效应； MEANS 效应〈/选择项〉； 五、循环语句 DO 指针变量＝初值 to 终值 ； END; 七、SAS应用 Sas简明教程P40 例5.3 病情类似的男性哮喘静息期患者20例，条件齐同配为一组。区组中的每一例随机分配一种处理，服药4小时测定各例的最大呼气中期流速MMF(升/秒)，结果如下，问不同药物及不同配伍组之间MMF值的差异有无显著性？ 1.建立检验假设，选定检验水准 2.Sas过程 随机区组设计的方差分析表 3.2结果解释 结果1为对A、B两因素作总的方差分析得F＝111.01，P＜0.001，按＝0.05的水准拒绝H0，可认为不同处理组和不同配伍组总的来说患者的MMF有差异。 进一步对不同处理（药物）组比较（即A因素）分析得F＝175.58，P＜0.001，按＝0.05的水准拒绝H0，可认为不同处理组患者间的MMF有差异；对不同配伍组比较（即B因素）分析得F＝62.58，P＜0.001，按＝0.05的水准拒绝H0，可认为不同配伍组患者间的MMF有差异。 SNK法两两比较 注：为了便于理解，SAS将两两比较的结果直接用英文字母的形式标示出来。两两比较结果的最右侧是处理因素变量的取值，最左侧标以字母A、B、C等等，用以表示该处理组和其它组有无差异。如果两组有相同的字母（如2、3两组），则两者之间无差异；而如果两组间只有不同的字母，则表示两组间的差异有统计学意义。 4.1.1结果解释 经snk两两比较，结果除氨茶碱0.2g组与氨茶碱0.2g、氨茶碱0.2g组和舒喘灵2.4mg组总体均数间的差异无统计学意义外，其余各组总体均数间的差异均有统计学意义，即除对照组MMF较低外，其余各组MMF相等。 4.2.1结果解释 采用SNK方法对配伍组（即B因素）进行多重比较，结果得除第2与第3个配伍组总体均数间的差异无统计学意义外，其余各配伍组总体均数间的差异均有统计学意义，即第2 和第3配伍组MMF相同，配伍组MMF由高到低顺序为第4、5、1、2和3组。 小结 过程歩： Proc anova;proc glm; 语句： Class group block; Model x= group block; Means group/snk; Means block/snk; 作业 课本P153第四题 作业要求： 1.要求机算； 2.内容：方差分析（要求写出方差分析表） 两两比较 SS误差/v误差 SS区组/ b-1 SS处理/ k-1 MS N-k- b+1 SS误差 误差 N-1 SS总 总 MS区组/ MS误差 b-1 SS区组 区组 MS处理/ MS误差 k-1 SS处理 处理组 F df SS 变异来源 表5.3 20例男性哮喘静息期患者治疗后MMF的测定值 2.22 2.24 2.14 1.52 2.32 2.35 2.25 1.74 1.89 1.86 1.78 1.05 1.93 1.92 1.80 1.25 2.01 1.98 1.85 1.42 氨茶碱0.1g及舒喘灵2.4mg 氨茶碱0.2g 氨茶碱0.1g 对照 数据步 data SAS5_3; do a=1 to 4; do b=1 to 5; input x @@; output; end; end; cards; 1.42 1.25 1.05 1.74 1.52 1.85 1.80 1.78 2.25 2.14 1.98 1