第八章 方分析ppt课件

1、第八章方差分析一）第八章方差分析一）Analysis of Variance Analysis of Variance （ANOVA )ANOVA ) 因素也称为处理因素因素也称为处理因素factor）（名义分类变量），每）（名义分类变量），每一处理因素至少有两个水平一处理因素至少有两个水平(level)（也称（也称“处理组处理组”）。）。 一个因素水平间独立）一个因素水平间独立） 单向方差分析单向方差分析 （第十章）（第十章） 两个因素水平间独立或相关）两个因素水平间独立或相关）双向方差分析双向方差分析 （第十一章）（第十一章） 一个个体多个测量值一个个体多个测量值重复测量资料的方差分析重复

2、测量资料的方差分析 ANOVA与回归分析相结合与回归分析相结合协方差分析协方差分析 目的：用这类资料的样本信息来推断各处理组间多个总目的：用这类资料的样本信息来推断各处理组间多个总体均数的差别有无统计学意义。体均数的差别有无统计学意义。SiS1S2S3S4合计值5.99 4.15 3.78 4.71 6.65 ANOVA ANOVA 由英国统由英国统计学家计学家R.A.FisherR.A.Fisher首首创，为纪念创，为纪念FisherFisher，以以F F命名，故方差分析命名，故方差分析又称又称 F F 检验检验 （F F testtest）。用于推断多）。用于推断多个总体均数有无差异个总

3、体均数有无差异 第一节第一节 方差分析的基本思想方差分析的基本思想 将所有测量值间的总变异按照其变异的来源分解为多个部份，然后进行比较，评价由某种因素所引起的变异是否具有统计学意义。一、离均差平方和的分解一、离均差平方和的分解组间变异组间变异总变异总变异组内变异组内变异对于例对于例8-1完全随机设计）完全随机设计）资料，共有三种不同的变异资料，共有三种不同的变异 总变异总变异Total variation）：全部测量值）：全部测量值Yij与总均数与总均数 间的差异间的差异 组间变异（组间变异（ between group variation ）：各）：各组的均数组的均数 与总均数与总均数 间的

4、差异间的差异组内变异组内变异within group variation )：每组的：每组的每个测量值每个测量值Yij与该组均数与该组均数 的差异的差异下面用离均差平方和下面用离均差平方和(sum of squares of (sum of squares of deviations from meandeviations from mean，SS)SS)反映变异的大小反映变异的大小 20.0Y YiYiY 1. 总变异: 所有测量值之间总的变异程度，计算公式22111122,1)iinnaaijijijijNiji jSSYYYCYCNS 总(2211,()()inaNijijiji jYY

5、CNN校正系数：校正系数：1N总 2组间变异：各组均数与总均数的离均差平方和，计算公式为21211()()inijjaaiiiiiYSSn YYCn组间1a组间SS组间反映了各组均数 的变异程度组间变异随机误差+处理因素效应 iY21121()(1)inaijiijaiiiSSYYnS 组 内Na组内 3组内变异：在同一处理组内，虽然每个受试对象接受的处理相同，但测量值仍各不相同，这种变异称为组内变异，也称SS误差。 用各组内各测量值Yij与其所在组的均数差值的平方和来表示，反映随机误差的影响。计算公式为三种三种“变异之间的关系变异之间的关系离均差平方和分解：离均差平方和分解：One-Fact

6、or ANOVA Partitions of Total VariationVariation Due to Treatment SSBVariation Due to Random Sampling SSWTotal Variation SSTw Commonly referred to as:w Sum of Squares Within, orw Sum of Squares Error, orw Within Groups Variationw Commonly referred to as:w Sum of Squares Among, orw Sum of Squares Betw

7、een, orw Sum of Squares Model, orw Among Groups Variation=+ 均方差，均方均方差，均方(mean square(mean square，MS) MS) 二、二、F 值与值与F分布分布，F 分布曲线分布曲线10,10215, 1215, 52122121122/22/12121121)(222)(FFFfF 界值表界值表附表附表5 F5 F界值表方差分析用，单侧界值）界值表方差分析用，单侧界值）上行：上行：P=0.05 P=0.05 下行：下行：P=0.01P=0.01分母自由度分母自由度2 2分子的自由度，分子的自由度，1 11 12

8、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.2599.3099.3099.3399.33 25254.244.243.393.392.992.992.762.762.602.602.492.49 7.777.775.575.574.6

9、84.684.184.183.853.853.633.63 5F F 分布曲线下面积与概率分布曲线下面积与概率第二节第二节 方差分析的步骤方差分析的步骤H0： 即即4个试验组总体均数相等个试验组总体均数相等 H1：4个试验组总体均数不全相等个试验组总体均数不全相等 检验水准检验水准 12340.05一、一、 建立检验假设建立检验假设SiS1S2S3S4合计值5.99 4.15 3.78 4.71 6.65 二、二、 计算离均差平方、自由度、均方计算离均差平方、自由度、均方三、计算三、计算F值值四、下结论四、下结论 注意：当组数为注意：当组数为2时，完全随机设计的方时，完全随机设计的方差分析结果

10、与两样本均数比较的差分析结果与两样本均数比较的t检验结果等检验结果等价，对同一资料价，对同一资料,有：有：tF第三节第三节 平均值之间的多重比较平均值之间的多重比较不拒绝不拒绝H0，表示拒绝总体均数相等的证据，表示拒绝总体均数相等的证据不足不足 分析终止。分析终止。拒绝拒绝H0，接受，接受H1, 表示总体均数不全相等表示总体均数不全相等哪两两均数之间相等？哪两两均数之间相等？哪两两均数之间不等？哪两两均数之间不等？ 需要进一步作多重比较。需要进一步作多重比较。控制累积控制累积类错误概率增大的方法类错误概率增大的方法采用采用Bonferroni法、法、SNK法和法和Tukey法等方法法等方法累积

11、累积类错误的概率为类错误的概率为当有当有k个均数需作两两比较时，比较的次数共有个均数需作两两比较时，比较的次数共有c= k!/(2!(k-2)!)=k(k-1)/2设每次检验所用设每次检验所用类错误的概率水准为类错误的概率水准为，累积，累积类错误的概率为类错误的概率为，则在对同一实验资料进行，则在对同一实验资料进行c次检次检验时，在样本彼此独立的条件下，根据概率乘法原理，验时，在样本彼此独立的条件下，根据概率乘法原理，其累积其累积类错误概率类错误概率与与c有下列关系：有下列关系：1(1)c (8.6)例如，设例如，设0.05，c=3(即即k=3)，其累积，其累积类错误类错误的概率为的概率为1(

12、1-0.05)3 =1-(0.95)3 = 0.1432k 一、一、BonferroniBonferroni法法方法：采用方法：采用/c作为下结论时所采用的作为下结论时所采用的检验水准。检验水准。c为两两比较次数，为两两比较次数， 为累积为累积I类错误的概率。类错误的概率。12,11ihiheYYYYtNaSMSnn组内组内（）例例8-18-1四个均值的四个均值的BonferroniBonferroni法比较法比较 设设/c0.05/6=0.0083,由此由此t的临的临界值为界值为t0.0083/2,20）=2.927118.528.0(:),244201122.3866(:C)0.072.9

13、271,(:)1.3523.482.92713.402.92714.832.9.9271(:), (:),(:)1.432271.9271t A Bt At A Dt B Ct BtBDC D 同理只有有统计学意义，其他与其他各无统计组间差异学意义。BonferroniBonferroni法的适用性法的适用性 当比较次数不多时，当比较次数不多时，Bonferroni法的效法的效果较好。果较好。 但当比较次数较多但当比较次数较多(例如在例如在10次以上次以上)时，时，则由于其检验水准选择得过低，结论偏于保则由于其检验水准选择得过低，结论偏于保守。守。二、二、SNKSNK法法 SNK(studen

14、t-Newman-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），作出统计推断。附表附表6第四节第四节 方差分析的假定条件和数据转换方差分析的假定条件和数据转换 一、方差分析的假定条件上述条件与两均数比较的t检验的应用条件相同。）1.各处理组样本来自随机、独立的正态总体(D法、W法、卡方检验)；2.各处理组样本的总体方差相等不等会增加I型错误的概率，影响方差分析结果的判断） 二、方差齐性检验1. Bartlett检验法2. Levene等3. 最大方差与最小方差之比3,初步认为方差齐同。1. Bartlett 检验法2. Levene