1、实验七实验七完全随机设计资料完全随机设计资料的方差分析的方差分析目的与要求：目的与要求：（一）掌握内容（一）掌握内容 1.anova和和glm过程的格式过程的格式;2.利用利用anova或或glm过程进行完全随机设计的过程进行完全随机设计的 方差分析，能对方差分析，能对SAS程序的输出结果作出合理程序的输出结果作出合理 的解释的解释;（二）熟悉内容（二）熟悉内容 snk、Dunnett、Bonfferoni等多重比较方法在等多重比较方法在SAS中的实现。中的实现。一、方差分析的基本思想一、方差分析的基本思想 方差分析的基本思想是把方差分析的基本思想是把全部全部观察值间观察值间的变异按设计和需要
2、的变异按设计和需要分解分解成成两个两个或或多个多个组成部分，然后将组成部分，然后将各部分各部分的变异的变异与与随机误差随机误差进行进行比较比较，以判断，以判断各部分的变异各部分的变异是否具有统计学意义是否具有统计学意义。二、完全随机设计的方差分析表二、完全随机设计的方差分析表变异来源SSdfMSFP组间(处理组间)SS组间k-1SS组间/v组间MS组间/MS组内组内(误差)SS组内N-kSS组内/v组内总变异SS总N-1三、三、应用条件应用条件1.1.独立性独立性:各样本是相互独立的随机各样本是相互独立的随机 样本样本；2.2.正态性正态性:各样本来自正态总体各样本来自正态总体；3.3.方差齐
3、性方差齐性:各总体方差相等。各总体方差相等。四四、方差齐性检验、方差齐性检验1.F检验检验：适用于正态分布资料：适用于正态分布资料两个总体两个总体方方 差齐性检验差齐性检验2.Bartlett 检验检验：资料服从正态分布的：资料服从正态分布的多个总体多个总体 方差齐性检验法方差齐性检验法 3.Levene检验检验：适用于适用于任何分布资料任何分布资料,可用于可用于两个两个 或多个总体或多个总体方差齐性检验方差齐性检验不满足条件的处理不满足条件的处理（1）轻微轻微l允许允许应用应用t检验、方差分析来作分析检验、方差分析来作分析。（2）严重严重l数据转换数据转换 (对数变换、平方根变换、平方根反正
4、弦变换对数变换、平方根变换、平方根反正弦变换)l非参数统计非参数统计五五、常用两两比较方法、常用两两比较方法lDunnett法法：多个实验组与：多个实验组与1个对照组比较个对照组比较lSNK法法：多组之间两两比较：多组之间两两比较lBonfferoni法法：多组之间两两比较：多组之间两两比较六、方差分析所用的六、方差分析所用的SAS过程过程ANOVA过程过程(Analysis Of Variance)GLM过程过程(General Linear Model)七、七、SAS应用应用Sas简明教程简明教程P34页页例例5.1 盐酸二甲双胍具有降低血糖作用，为了解盐酸二甲双胍具有降低血糖作用，为了解
5、某种新药是否具有较好的降低血糖作用，某种新药是否具有较好的降低血糖作用，18名名病情相近的糖尿病患者被病情相近的糖尿病患者被随机随机分为三组，治疗分为三组，治疗一个月后，测定患者治疗前后餐后两小时一个月后，测定患者治疗前后餐后两小时血糖血糖的下降值的下降值数据见表数据见表5-1。问该新药。问该新药是否具有较是否具有较好的好的降低血糖的作用？降低血糖的作用？表表5-1 18名患者治疗前后餐后两小时血糖的下降值名患者治疗前后餐后两小时血糖的下降值（mmol/L）安慰剂组安慰剂组盐酸二甲双胍组盐酸二甲双胍组新药组新药组1.03.24.31.54.13.9-1.32.14.1-1.03.53.40.5
6、3.12.90.95.13.6血糖的下血糖的下降值降值分组分组变量变量1.建立检验假设建立检验假设,确定检验水准确定检验水准H0：，即接受不同处理患者血糖下降值，即接受不同处理患者血糖下降值总体均数相等总体均数相等H1：、不全相等，即接受不同处理患者血不全相等，即接受不同处理患者血糖下降值总体均数不全相等糖下降值总体均数不全相等=0.05321123 2.建立数据集建立数据集data aa1；input x g；Cards;1.0 1 3.2 24.3 31.5 1 4.1 23.9 3-1.3 1 2.1 24.1 3-1.0 1 3.523.4 30.5 1 3.122.9 30.9 1
7、5.123.6 3;3.统计描述（正态性检验）统计描述（正态性检验）建立检验假设，确定检验水准建立检验假设，确定检验水准proc univariate normal;class g;var x;run;The SAS System The UNIVARIATE ProcedureVariable:x Tests for Normalityg=1 Test -Statistic-p Value-Shapiro-Wilk W 0.874111 Pr W 0.2431g=2 Test -Statistic-p Value-Shapiro-Wilk W 0.975005 Pr W 0.9242 g=3
8、 Test -Statistic-p Value-Shapiro-Wilk W 0.970331 Pr F g 2 2.4751 1.2376 1.76 0.2056 Error 15 10.5421 0.7028经方差齐性检验得，经方差齐性检验得，F=1.76，P0.20560.10,故暂不拒绝故暂不拒绝H0，可认为三组总体的方差齐。可认为三组总体的方差齐。离均差平方和均方The SAS SystemThe ANOVA ProcedureDependent Variable:x变异来源 离均差平方和 均方 Sum ofSource DF Squares Mean Square F Value
9、 Pr FModel 2 44.76777778 22.38388889 25.82 F g 2 44.76777778 22.38388889 25.82 .0001F=25.82，P0.0001，拒绝，拒绝H0,差别有统计学差别有统计学意义，可认为三组患者血糖下降值总体均数意义，可认为三组患者血糖下降值总体均数不相等不相等。完全随机设计的方差分析表完全随机设计的方差分析表变异来源SSdfMSFP组间(处理组间)44.76772 22.383925.82，方差齐；方差齐；(3)方差分析方差分析：F=25.82，P0.0001，拒绝，拒绝H0,差别有统差别有统计计 学意义，学意义，可认为三组患
10、者血糖下降值总体均数不可认为三组患者血糖下降值总体均数不相等相等。(4)两两比较两两比较：即盐酸二甲双胍组和新药组均与第安慰：即盐酸二甲双胍组和新药组均与第安慰剂组比较的剂组比较的P值均值均0.05，在，在 0.05的水准上不拒的水准上不拒绝绝H0，可认为盐酸二甲双胍和新药均有降低血糖作用。，可认为盐酸二甲双胍和新药均有降低血糖作用。盐酸二甲双胍组与新药组比较的盐酸二甲双胍组与新药组比较的P值值0.05，按，按 0.05的水准不拒绝的水准不拒绝H0，可认为两药降低血糖作用相同。，可认为两药降低血糖作用相同。小结小结:1.过程步过程步 Proc univariate proc anova2.选择项选择项 hovtest dunnett(“1”)snk bon 作业作业lP153，计算题第，计算题第2题题要求机算：要求机算：1.正态性正态性2.方差齐性方差齐性3.方差分析方差分析4.两两比较两两比较
