最新变异数方差分析精品课件.ppt

《最新变异数方差分析精品课件.ppt》由会员分享，可在线阅读，更多相关《最新变异数方差分析精品课件.ppt（41页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、华中科技大学 同济医学院 宇传华制作， 2004，102SiS1S2S3S4合计值5.99 4.15 3.78 4.71 6.65 9 2 2组间变异：组间变异：各组均数与总均数的各组均数与总均数的离均差平方和，离均差平方和，计算公式为计算公式为21211()()inijjaaiiiiiYSSn YYCn组间1a组间SS组间反映了各组均数 的变异程度组间变异组间变异随机误差随机误差+ +处理因素效应处理因素效应 iY1021121()(1)inaijiijaiiiSSYYnS 组 内Na组内 3组内变异：在同一处理组内，虽然每个受试对象接受的处理相同，但测量值仍各不相同，这种变异称为组内变异，

2、也称SS误差。 用各组内各测量值Yij与其所在组的均数差值的平方和来表示，反映随机误差的影响。计算公式为11三种三种“变异变异”之间的关系之间的关系离均差平方和离均差平方和分解分解：12One-Factor ANOVA Partitions of Total VariationVariation Due to Treatment SSBVariation Due to Random Sampling SSWTotal Variation SSTw Commonly referred to as:oSum of Squares Within, oroSum of Squares Error, o

3、roWithin Groups Variationw Commonly referred to as:oSum of Squares Among, oroSum of Squares Between, oroSum of Squares Model, oroAmong Groups Variation=+13 均方差，均方均方差，均方( (mean square，MS) ) 14 二、二、F 值与值与F分布分布，15F 分布曲线分布曲线10,10215, 1215, 52122121122/22/12121121)(222)(FFFf16F 界值表界值表附表附表5 5 F F界值表（方差分析用

4、，单侧界值）界值表（方差分析用，单侧界值）上行：上行：P P=0.05 =0.05 下行：下行：P P=0.01=0.01分母自由度分母自由度2 2分子的自由度，分子的自由度，1 11 12 23 34 45 56 6 1 1161161200200216216225225230230234234 405240524999499954035403562556255764576458595859 2 218.5118.5119.0019.0019.1619.1619.2519.2519.3019.3019.3319.33 98.4998.4999.0099.0099.1799.1799.2599

5、.2599.3099.3099.3399.33 25254.244.243.393.392.992.992.762.762.602.602.492.49 7.777.775.575.574.684.684.184.183.853.853.633.63 517F F 分布曲线下面积与概率分布曲线下面积与概率18华中科技大学 同济医学院 宇传华制作， 2004，1019第二节第二节 实例实例8.18.1的方差分析的方差分析华中科技大学 同济医学院 宇传华制作， 2004，1020H0： 即即4个试验组总体均数相等个试验组总体均数相等 H1：4个试验组总体均数个试验组总体均数不全相等不全相等 检验水

6、准检验水准 12340.05一、一、 建立检验假设建立检验假设华中科技大学 同济医学院 宇传华制作， 2004，1021SiS1S2S3S4合计值5.99 4.15 3.78 4.71 6.65 华中科技大学 同济医学院 宇传华制作， 2004，1022二、二、 计算离均差平方、自由度、均方计算离均差平方、自由度、均方华中科技大学 同济医学院 宇传华制作， 2004，1023三、计算三、计算F值值华中科技大学 同济医学院 宇传华制作， 2004，1024四、下结论四、下结论 注意：当组数为注意：当组数为2时，完全随机设计的方时，完全随机设计的方差分析结果与两样本均数比较的差分析结果与两样本均数

7、比较的t检验结果等检验结果等价，对同一资料价，对同一资料,有：有：tF华中科技大学 同济医学院 宇传华制作， 2004，1025第三节第三节 平均值之间的多重比较平均值之间的多重比较不拒绝不拒绝H0，表示拒绝总体均数相等的证据，表示拒绝总体均数相等的证据不足不足 分析终止。分析终止。拒绝拒绝H0，接受，接受H1, 表示总体均数不全相等表示总体均数不全相等哪两两均数之间相等？哪两两均数之间相等？哪两两均数之间不等？哪两两均数之间不等？ 需要进一步作多重比较。需要进一步作多重比较。华中科技大学 同济医学院 宇传华制作， 2004，1026控制累积控制累积类错误概率增大的方法类错误概率增大的方法采用

8、采用Bonferroni法、法、SNK法和法和Tukey法等方法法等方法华中科技大学 同济医学院 宇传华制作， 2004，1027累积累积类错误的概率为类错误的概率为 当有当有k个均数需作两两比较时，比较的次数共有个均数需作两两比较时，比较的次数共有c= k!/(2!(k-2)!)=k(k-1)/2设每次检验所用设每次检验所用类错误的概率水准为类错误的概率水准为，累积，累积类错误的概率为类错误的概率为，则在对同一实验资料进行，则在对同一实验资料进行c次检次检验时，在样本彼此独立的条件下，根据概率乘法原理，验时，在样本彼此独立的条件下，根据概率乘法原理，其累积其累积类错误概率类错误概率与与c有下

9、列关系：有下列关系：1(1)c (8.6)例如，设例如，设0.05，c=3(即即k=3)，其累积，其累积类错误类错误的概率为的概率为1(1-0.05)3 =1-(0.95)3 = 0.1432k 华中科技大学 同济医学院 宇传华制作， 2004，1028一、一、BonferroniBonferroni法法方法：采用方法：采用/c作为下结论时所采用的作为下结论时所采用的检验水准。检验水准。c为两两比较次数，为两两比较次数， 为累积为累积I类错误的概率。类错误的概率。12,11ihiheYYYYtNaSMSnn组内组内（）华中科技大学 同济医学院 宇传华制作， 2004，1029例例8-18-1四

10、个均值的四个均值的BonferroniBonferroni法比较法比较 设设/c0.05/6=0.0083,由此由此t的临的临界值为界值为t（0.0083/2,20）=2.927118.528.0(:),244201122.3866(:C)0.072.9271,(:)1.3523.482.92713.402.92714.832.9.9271(:), (:),(:)1.432271.9271t A Bt At A Dt B Ct BtBDC D 同理只有有统计学意义，其他与其他各无统计组间差异学意义。华中科技大学 同济医学院 宇传华制作， 2004，1030BonferroniBonferron

11、i法的适用性法的适用性 当当比较次数不多时比较次数不多时，Bonferroni法的效法的效果较好。果较好。 但当但当比较次数较多比较次数较多(例如在例如在10次以上次以上)时，时，则由于其检验水准选择得过低，结论偏于保则由于其检验水准选择得过低，结论偏于保守。守。华中科技大学 同济医学院 宇传华制作， 2004，1031二、二、SNKSNK法法 SNK(student-Newman-Keuls)法又称q检验，是根据q值的抽样分布作出统计推论（例8-1）。1将各组的平均值按由大到小的顺序排列由大到小的顺序排列： 顺序顺序(1)(2)(3)(4) 平均值平均值28.018.718.514.8 原组

12、号原组号BCAD2. 计算两个平均值之间的差值及组间跨度差值及组间跨度k，见表8-3第(2)、 (3)两列。3. 计算统计量计算统计量q值值4. 根据计算的q值及查附表6得到的q界值（p286），作出统计统计推断推断。华中科技大学 同济医学院 宇传华制作， 2004，1032附表附表6华中科技大学 同济医学院 宇传华制作， 2004，1033华中科技大学 同济医学院 宇传华制作， 2004，1034第四节第四节 方差分析的假定条件和数据转换方差分析的假定条件和数据转换 一、方差分析的假定条件一、方差分析的假定条件（上述条件与两均数比较的上述条件与两均数比较的t检检验的应用条件相同。）验的应用条

13、件相同。）1.各处理组样本来自随机、独立的正态总体(D法、W法、卡方检验)；2.各处理组样本的总体方差相等（不等会增加I型错误的概率，影响方差分析结果的判断） 二、方差齐性检验二、方差齐性检验1. Bartlett检验法2. Levene等3. 最大方差与最小方差之比3,初步认为方差齐同。华中科技大学 同济医学院 宇传华制作， 2004，10351. Bartlett 检验法华中科技大学 同济医学院 宇传华制作， 2004，10362. Levene 检验法 将原样本观察值作离均差变换，或离均差平方变换，然后执行完全随机设计的方差分析，其检验结果用于判断方差是否齐性。 因为levene检验对原

14、数据是否为正态不灵敏，所以比较稳健。目前均推荐采用LEVENE方差齐性检验华中科技大学 同济医学院 宇传华制作， 2004，1037 三、数据变换三、数据变换 改善数据的正态性或方差齐性。使之满足方差分析的假定条件。 平方根反正弦变换适用于二项分布率（比例）数据。 平方根变换适用于泊松分布的计数资料 对数变换适用于对数正态分布资料XY11sinsin180XYXY或10log ( )XY华中科技大学 同济医学院 宇传华制作， 2004，1038第五节第五节 完全随机设计方法简介完全随机设计方法简介将将120120名高血脂患者完全随机分成名高血脂患者完全随机分成4 4个例数相等的组个例数相等的组

15、 1. 1. 编号：编号：120120名高血脂患者从名高血脂患者从1 1开始到开始到120120，见下面表第，见下面表第1 1行；行；2. 2. 取随机数字：取随机数字：从附表从附表1515中的任一行任中的任一行任一列开始，如第一列开始，如第5 5行第行第7 7列开始，依次列开始，依次读取三位数作为一个随机数录于编号读取三位数作为一个随机数录于编号下，见下面表的第下，见下面表的第2 2行；行；华中科技大学 同济医学院 宇传华制作， 2004，10393. 3. 排序：排序：按随机数字从小到大按随机数字从小到大 ( (数据相同则按先后数据相同则按先后顺序）编序号，见下面表的第顺序）编序号，见下面表的第3 3行。行。4. 4. 事先规定：事先规定：序号序号1-301-30为甲组，序号为甲组，序号31-6031-60为乙组，为乙组，序号序号61-9061-90为丙组，序号为丙组，序号91-12091-120为丁组，见下面表的为丁组，见下面表的第第4 4行。行。华中科技大学 同济医学院 宇传华制作， 2004，1040宇传华 制作http:/41 结束语结束语