实验七方差分析精品文稿.ppt

《实验七方差分析精品文稿.ppt》由会员分享，可在线阅读，更多相关《实验七方差分析精品文稿.ppt（32页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、实验七方差分析实验七方差分析第1页，本讲稿共32页目的与要求：目的与要求：（一）掌握内容（一）掌握内容 1.anova和和glm过程的格式过程的格式;2.利用利用anova或或glm过程进行完全随机设计的过程进行完全随机设计的 方方差分析，能对差分析，能对SAS程序的输出结果作出合理程序的输出结果作出合理 的解释的解释;（二）熟悉内容（二）熟悉内容 snk、Dunnett、Bonfferoni等多重比较方法在等多重比较方法在SAS中中的实现。的实现。第2页，本讲稿共32页一、方差分析的基本思想一、方差分析的基本思想 方差分析的基本思想是把方差分析的基本思想是把全部观察值全部观察值间间的变异按设

2、计和需要的变异按设计和需要分解分解成成两个两个或或多多个个组成部分，然后将组成部分，然后将各部分的变异各部分的变异与与随随机误差机误差进行进行比较比较，以判断各部分的变异，以判断各部分的变异是否具有统计学意义是否具有统计学意义。第3页，本讲稿共32页二、完全随机设计的方差分析表二、完全随机设计的方差分析表变异来源SSdfMSFP组间(处理组间)SS组间k-1SS组间/v组间MS组间/MS组内组内(误差)SS组内N-kSS组内/v组内总变异SS总N-1第4页，本讲稿共32页三、三、应用条件应用条件1.1.独立性独立性:各样本是相互独立的随机各样本是相互独立的随机 样本；样本；2.2.正态性正态性

3、:各样本来自正态总体；各样本来自正态总体；3.3.方差齐性方差齐性:各总体方差相等。各总体方差相等。第5页，本讲稿共32页四四、方差齐性检验、方差齐性检验1.F检验检验：适用于正态分布资料：适用于正态分布资料两个总体两个总体方方 差齐性检验差齐性检验2.Bartlett 检验检验：资料服从正态分布的：资料服从正态分布的多个总体多个总体 方差齐性检验法方差齐性检验法 3.Levene检验检验：适用于适用于任何分布资料任何分布资料,可用于可用于两个两个 或多个总体或多个总体方差齐性检验方差齐性检验第6页，本讲稿共32页不满足条件的处理不满足条件的处理（1）轻微）轻微l允许允许应用应用t检验、方差分

4、析来作分析。检验、方差分析来作分析。（2）严重）严重l数据转换数据转换 (对数变换、平方根变换、平方根反正弦变换对数变换、平方根变换、平方根反正弦变换)l非参数统计非参数统计第7页，本讲稿共32页五五、常用两两比较方法、常用两两比较方法lDunnett法法：多个实验组与：多个实验组与1个对照组比较个对照组比较lSNK法法：多组之间两两比较：多组之间两两比较lBonfferoni法法：多组之间两两比较：多组之间两两比较第8页，本讲稿共32页六、方差分析所用的六、方差分析所用的SAS过程过程ANOVA过程过程(Analysis Of Variance)GLM过程过程(General Linear

5、Model)第9页，本讲稿共32页七、七、SAS应用应用Sas简明教程简明教程P34页页例例5.1 盐酸二甲双胍具有降低血糖作用，为了解盐酸二甲双胍具有降低血糖作用，为了解某种新药是否具有较好的降低血糖作用，某种新药是否具有较好的降低血糖作用，18名名病情相近的糖尿病患者被病情相近的糖尿病患者被随机随机分为三组，治疗分为三组，治疗一个月后，测定患者治疗前后餐后两小时一个月后，测定患者治疗前后餐后两小时血糖血糖的下降值的下降值数据见表数据见表5-1。问该新药。问该新药是否具有较是否具有较好的好的降低血糖的作用？降低血糖的作用？第10页，本讲稿共32页表表5-1 18名患者治疗前后餐后两小时血糖的

6、下降值名患者治疗前后餐后两小时血糖的下降值（mmol/L）安慰剂组安慰剂组盐酸二甲双胍酸二甲双胍组新药组新药组1.03.24.31.54.13.9-1.32.14.1-1.03.53.40.53.12.90.95.13.6第11页，本讲稿共32页血糖的下降血糖的下降值值分组变分组变量量第12页，本讲稿共32页1.建立检验假设建立检验假设,确定检验水准确定检验水准H0：，即接受不同处理患者血糖下降值，即接受不同处理患者血糖下降值总体均数相等总体均数相等H1：、不全相等，即接受不同处理患者血不全相等，即接受不同处理患者血糖下降值总体均数不全相等糖下降值总体均数不全相等=0.05第13页，本讲稿共3

7、2页 2.建立数据集建立数据集data aa1；input x g；Cards;1.0 1 3.2 24.331.5 1 4.1 23.93-1.3 1 2.1 24.13-1.0 1 3.523.430.5 1 3.122.930.9 1 5.123.63;第14页，本讲稿共32页3.统计描述（正态性检验）统计描述（正态性检验）建立检验假设，确定检验水准建立检验假设，确定检验水准proc univariate normal;class g;var x;run;第15页，本讲稿共32页The SAS System The UNIVARIATE ProcedureVariable:x Tests

8、 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 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

9、，可认，可认为三组总体的方差齐。为三组总体的方差齐。离均差平方和均方第19页，本讲稿共32页The SAS SystemThe ANOVA ProcedureDependent Variable:x变异来源 离均差平方和 均方 Sum ofSource DF Squares Mean Square F Value Pr FModel 2 44.76777778 22.38388889 25.82 F g 2 44.76777778 22.38388889 25.82 .0001F=25.82，P0.0001，拒绝，拒绝H0,差别有统计学意义，差别有统计学意义，可认为三组患者血糖下降值总体均数不

10、相等。可认为三组患者血糖下降值总体均数不相等。第20页，本讲稿共32页完全随机设计的方差分析表完全随机设计的方差分析表变异来源SSdfMSFP组间(处理组间)44.76772 22.383925.82，方差齐；，方差齐；(3)方差分析方差分析：F=25.82，P0.0001，拒绝，拒绝H0,差别有统计差别有统计 学意义，可认为三组患者血糖下降值总体均数不相等。学意义，可认为三组患者血糖下降值总体均数不相等。(4)两两比较两两比较：即盐酸二甲双胍组和新药组均与第安慰剂组比较：即盐酸二甲双胍组和新药组均与第安慰剂组比较的的P值均值均0.05，在，在 0.05的水准上不拒绝的水准上不拒绝H0，可认为

11、盐，可认为盐酸二甲双胍和新药均有降低血糖作用。盐酸二甲双胍组与新酸二甲双胍和新药均有降低血糖作用。盐酸二甲双胍组与新药组比较的药组比较的P值值0.05，按，按 0.05的水准不拒绝的水准不拒绝H0，可认，可认为两药降低血糖作用相同。为两药降低血糖作用相同。第30页，本讲稿共32页小结小结:1.过程步过程步 Proc univariate proc anova2.选择项选择项 hovtest dunnett(“1”)snk bon 第31页，本讲稿共32页作业作业lP153，计算题第，计算题第2题题要求机算：要求机算：1.正态性正态性2.方差齐性方差齐性3.方差分析方差分析4.两两比较两两比较第32页，本讲稿共32页