1、华中科技大学 同济医学院 宇传华制作，2019，101方差分析方差分析Analysis of Variance（ANOVA)因素也称为因素也称为处理处理，每，每一处理因素至少有两个一处理因素至少有两个水平水平(level)（也称（也称“处理组处理组”）。）。一个一个因素因素（水平水平间独立）间独立）单向方差分析单向方差分析 （第十章）（第十章）两个两个因素因素（水平水平间独立或相关）间独立或相关）双向方差分析双向方差分析 （第十一章）（第十一章）一个个体多个测量值一个个体多个测量值重复测量资料的方差分析重复测量资料的方差分析 ANOVA与回归分析相结合与回归分析相结合协方差分析协方差分析 目的
2、：目的：用这类资料的样本信息来推断各处理组间用这类资料的样本信息来推断各处理组间多个总多个总体均数体均数的差别有无统计学意义。的差别有无统计学意义。华中科技大学 同济医学院 宇传华制作，2019，102SiS1S2S3S4合计值5.99 4.15 3.78 4.71 6.65 华中科技大学 同济医学院 宇传华制作，2019，103华中科技大学 同济医学院 宇传华制作，2019，104 ANOVA ANOVA 由英国统由英国统计学家计学家R.A.FisherR.A.Fisher首首创，为纪念创，为纪念FisherFisher，以以F F命名，故方差分析命名，故方差分析又称又称 F F 检验检验
3、（F F testtest）。用于推断）。用于推断多多个总体均数个总体均数有无差异有无差异 华中科技大学 同济医学院 宇传华制作，2019，105第十章第十章 单向方差分析单向方差分析One-way analysis of variance第一节第一节 方差分析的基本思想方差分析的基本思想 将所有测量值间的总变异总变异按照其变异的来源分解为多个部份分解为多个部份，然后进行比较，评价由某种因素某种因素所引起的变异是否具有统计学意义。华中科技大学 同济医学院 宇传华制作，2019，106一、离均差平方和的分解一、离均差平方和的分解组间变异组间变异总变异总变异组内变异组内变异华中科技大学 同济医学院
4、 宇传华制作，2019，107对于例对于例8-1（完全随机设计）（完全随机设计）资料，共有三种不同的变异资料，共有三种不同的变异 总变异总变异（Total variation）：全部测量值）：全部测量值Yij与与总均数总均数 间的差异间的差异 组间变异组间变异（between group variation）：各）：各组的均数组的均数 与总均数与总均数 间的差异间的差异组内变异组内变异（within group variation)：每组的：每组的每个测量值每个测量值Yij与该组均数与该组均数 的差异的差异下面用下面用离均差平方和离均差平方和(sum of squares of(sum of
5、squares of deviations from meandeviations from mean，SSSS)反映变异的大小反映变异的大小 20.0Y YiYiY8 1.1.总变异总变异:所有测量值之间总所有测量值之间总的变异程度，的变异程度，计算公式计算公式22111122,1)iinnaaijijijijNiji jSSYYYCYCNS 总(2211,()()inaNijijiji jYYCNN校正系数校正系数：1N总9 2 2组间变异：组间变异：各组均数与总均数的各组均数与总均数的离均差平方和，离均差平方和，计算公式为计算公式为21211()()inijjaaiiiiiYSSn YY
6、Cn组间1a组间SS组间反映了各组均数 的变异程度组间变异组间变异随机误差随机误差+处理因素效应处理因素效应 iY1021121()(1)inaijiijaiiiSSYYnS 组 内Na组内 3组内变异：在同一处理组内，虽然每个受试对象接受的处理相同，但测量值仍各不相同，这种变异称为组内变异，也称SS误差。用各组内各测量值Yij与其所在组的均数差值的平方和来表示，反映随机误差的影响。计算公式为11三种三种“变异变异”之间的关系之间的关系离均差平方和离均差平方和分解分解：12One-Factor ANOVA Partitions of Total VariationVariation Due t
7、o Treatment SSBVariation Due to Random Sampling SSWTotal Variation SSTw Commonly referred to as:oSum of Squares Within,oroSum of Squares Error,oroWithin Groups Variationw Commonly referred to as:oSum of Squares Among,oroSum of Squares Between,oroSum of Squares Model,oroAmong Groups Variation=+13 均方差
8、，均方均方差，均方(mean square，MS)14 二、二、F 值与值与F分布分布，15F 分布曲线分布曲线10,10215,1215,52122121122/22/12121121)(222)(FFFf16F 界值表界值表附表附表5 5 F F界值表（方差分析用，单侧界值）界值表（方差分析用，单侧界值）上行：上行：P P=0.05 =0.05 下行：下行：P P=0.01=0.01分母自由度分母自由度2 2分子的自由度，分子的自由度，1 11 12 23 34 45 56 6 1 1161161200200216216225225230230234234 4052405249994999
9、54035403562556255764576458595859 2 218.5118.5119.0019.0019.1619.1619.2519.2519.3019.3019.3319.33 98.4998.4999.0099.0099.1799.1799.2599.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华中科技大学 同济医学
10、院 宇传华制作，2019，1019第二节第二节 实例实例8.18.1的方差分析的方差分析华中科技大学 同济医学院 宇传华制作，2019，1020H0：即即4个试验组总体均数相等个试验组总体均数相等 H1：4个试验组总体均数个试验组总体均数不全相等不全相等 检验水准检验水准 12340.05一、一、建立检验假设建立检验假设华中科技大学 同济医学院 宇传华制作，2019，1021SiS1S2S3S4合计值5.99 4.15 3.78 4.71 6.65 华中科技大学 同济医学院 宇传华制作，2019，1022二、二、计算离均差平方、自由度、均方计算离均差平方、自由度、均方华中科技大学 同济医学院
11、宇传华制作，2019，1023三、计算三、计算F值值华中科技大学 同济医学院 宇传华制作，2019，1024四、下结论四、下结论 注意：当组数为注意：当组数为2时，完全随机设计的方时，完全随机设计的方差分析结果与两样本均数比较的差分析结果与两样本均数比较的t检验结果等检验结果等价，对同一资料价，对同一资料,有：有：tF华中科技大学 同济医学院 宇传华制作，2019，1025第三节第三节 平均值之间的多重比较平均值之间的多重比较不拒绝不拒绝H0，表示拒绝总体均数相等的证据，表示拒绝总体均数相等的证据不足不足 分析终止。分析终止。拒绝拒绝H0，接受，接受H1,表示总体均数不全相等表示总体均数不全相
12、等哪两两均数之间相等？哪两两均数之间相等？哪两两均数之间不等？哪两两均数之间不等？需要进一步作多重比较。需要进一步作多重比较。华中科技大学 同济医学院 宇传华制作，2019，1026控制累积控制累积类错误概率增大的方法类错误概率增大的方法采用采用Bonferroni法、法、SNK法和法和Tukey法等方法法等方法华中科技大学 同济医学院 宇传华制作，2019，1027累积累积类错误的概率为类错误的概率为 当有当有k个均数需作两两比较时，比较的次数共有个均数需作两两比较时，比较的次数共有c=k!/(2!(k-2)!)=k(k-1)/2设每次检验所用设每次检验所用类错误的概率水准为类错误的概率水准
13、为，累积，累积类错误的概率为类错误的概率为，则在对同一实验资料进行，则在对同一实验资料进行c次检次检验时，在样本彼此独立的条件下，根据概率乘法原理，验时，在样本彼此独立的条件下，根据概率乘法原理，其累积其累积类错误概率类错误概率与与c有下列关系：有下列关系：1(1)c (8.6)例如，设例如，设0.05，c=3(即即k=3)，其累积，其累积类错误类错误的概率为的概率为1(1-0.05)3=1-(0.95)3=0.1432k 华中科技大学 同济医学院 宇传华制作，2019，1028一、一、BonferroniBonferroni法法方法：采用方法：采用/c作为下结论时所采用的作为下结论时所采用的
14、检验水准。检验水准。c为两两比较次数，为两两比较次数，为累积为累积I类错误的概率。类错误的概率。12,11ihiheYYYYtNaSMSnn组内组内（）华中科技大学 同济医学院 宇传华制作，2019，1029例例8-18-1四个均值的四个均值的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.432
15、271.9271t A Bt At A Dt B Ct BtBDC D 同理只有有统计学意义，其他与其他各无统计组间差异学意义。华中科技大学 同济医学院 宇传华制作，2019，1030BonferroniBonferroni法的适用性法的适用性 当当比较次数不多时比较次数不多时，Bonferroni法的效法的效果较好。果较好。但当但当比较次数较多比较次数较多(例如在例如在10次以上次以上)时，时，则由于其检验水准选择得过低，结论偏于保则由于其检验水准选择得过低，结论偏于保守。守。华中科技大学 同济医学院 宇传华制作，2019，1031二、二、SNKSNK法法 SNK(student-Newma
16、n-Keuls)法又称q检验，是根据q值的抽样分布作出统计推论（例8-1）。1将各组的平均值按由大到小的顺序排列由大到小的顺序排列：顺序顺序(1)(2)(3)(4)平均值平均值28.018.718.514.8 原组号原组号BCAD2.计算两个平均值之间的差值及组间跨度差值及组间跨度k，见表8-3第(2)、(3)两列。3.计算统计量计算统计量q值值4.根据计算的q值及查附表6得到的q界值（p286），作出统计统计推断推断。华中科技大学 同济医学院 宇传华制作，2019，1032附表附表6华中科技大学 同济医学院 宇传华制作，2019，1033华中科技大学 同济医学院 宇传华制作，2019，103
17、4第四节第四节 方差分析的假定条件和数据转换方差分析的假定条件和数据转换 一、方差分析的假定条件一、方差分析的假定条件（上述条件与两均数比较的上述条件与两均数比较的t检检验的应用条件相同。）验的应用条件相同。）1.各处理组样本来自随机、独立的正态总体(D法、W法、卡方检验)；2.各处理组样本的总体方差相等（不等会增加I型错误的概率，影响方差分析结果的判断）二、方差齐性检验二、方差齐性检验1.Bartlett检验法2.Levene等3.最大方差与最小方差之比3,初步认为方差齐同。华中科技大学 同济医学院 宇传华制作，2019，10351.Bartlett 检验法华中科技大学 同济医学院 宇传华制
18、作，2019，10362.Levene 检验法 将原样本观察值作离均差变换，或离均差平方变换，然后执行完全随机设计的方差分析，其检验结果用于判断方差是否齐性。因为levene检验对原数据是否为正态不灵敏，所以比较稳健。目前均推荐采用LEVENE方差齐性检验华中科技大学 同济医学院 宇传华制作，2019，1037 三、数据变换三、数据变换 改善数据的正态性或方差齐性。使之满足方差分析的假定条件。平方根反正弦变换适用于二项分布率（比例）数据。平方根变换适用于泊松分布的计数资料 对数变换适用于对数正态分布资料XY11sinsin180XYXY或10log()XY华中科技大学 同济医学院 宇传华制作，
19、2019，1038第五节第五节 完全随机设计方法简介完全随机设计方法简介将将120120名高血脂患者完全随机分成名高血脂患者完全随机分成4 4个例数相等的组个例数相等的组 1.1.编号：编号：120120名高血脂患者从名高血脂患者从1 1开始到开始到120120，见下面表第，见下面表第1 1行；行；2.2.取随机数字：取随机数字：从附表从附表1515中的任一行任中的任一行任一列开始，如第一列开始，如第5 5行第行第7 7列开始，依次列开始，依次读取三位数作为一个随机数录于编号读取三位数作为一个随机数录于编号下，见下面表的第下，见下面表的第2 2行；行；华中科技大学 同济医学院 宇传华制作，2019，10393.3.排序：排序：按随机数字从小到大按随机数字从小到大 (数据相同则按先后数据相同则按先后顺序）编序号，见下面表的第顺序）编序号，见下面表的第3 3行。行。4.4.事先规定：事先规定：序号序号1-301-30为甲组，序号为甲组，序号31-6031-60为乙组，为乙组，序号序号61-9061-90为丙组，序号为丙组，序号91-12091-120为丁组，见下面表的为丁组，见下面表的第第4 4行。行。华中科技大学 同济医学院 宇传华制作，2019，1040宇传华 制作statdtedm.6to23
