方差分析-华北煤炭医学院卫生统计学

,第九章方差分析（analysisofvariance）,无论是在试验研究还是在调查研究中，如比较组数为k个，当k2时，两组总体均数比较是否相等的假设可采用前面介绍的t检验或z检验（当然也可采用今天所介绍的方差分析）；当k2时，即检验两组以上的总体均数是否相等时，如采用前面所学两样本t检验进行两两比较，而得出结论，会使犯一类错误的概率增大。,如从已知正态总体N（10，52）进行随机抽样，共抽取k10组样本，每组样本含量均20，每组样本均数和标准差结果见表92。,表92从已知正态总体N（10，52）进行随机抽取10个样本的(ni=20)的结果,如要进行两两比较的t检验，即10个样本每两个进行，其比较的次数为：,实验结果表示：若0.05，则在45次比较中，发现有5次有统计学意义，结果见表93。从理论上讲10个样本均来自同一正态总体，应当无差别，但我们用两样本比较的t检验时，规定0.05，其实际犯第一类错误的概率为5/450.11,显然比所要控制的0.05要大。,表93.45次比较中5次有统计学意义的结果,需采用本章介绍的方差分析(ANOVA),方差分析的用途及应用条件进行两个或两个以上样本均数的比较；方差分析对分析数据的要求及条件比较严格，即要求各样本为随机样本，各样本来自正态总体，各样本所代表的总体方差齐性或相等。a.可以同时分析一个、两个或多个因素对试验结果的作用和影响；b.分析多个因素的独立作用及多个因素之间的交互作用；进行回归方程的线性假设检验。,第一节方差分析的基本概念,第二节完全随机设计资料的方差分析,第三节随机区组设计资料的方差分析,第四节多个样本均数间的两两比较,第五节方差分析的前提条件和变量变换,第九章方差分析,一、方差分析的几个名词,离均差,离均差平方和SS,方差（2S2）均方（MS）,标准差：S,自由度：,关系：MS=SS/,二、方差分析的含义方差是描述变异的一种指标，方差分析是一种假设检验的方法。方差分析就是对变异的分析。是对总变异进行分析。看总变异是由哪些部分组成的，这些部分间的关系如何。,三、方差分析的基本思想根据变异的来源，将全部观察值总的离均差平方和及自由度分解为两个或多个部分，除随机误差外，其余每个部分的变异可由某些特定因素的作用加以解释。通过比较不同来源变异的方差（也叫均方MS），得出F值，借助F分布做出统计推断，从而判断某因素对观察指标有无影响。,第一节方差分析的基本概念,第二节完全随机设计资料的方差分析,第三节随机区组设计资料的方差分析,第四节多个样本均数间的两两比较,第五节方差分析的前提条件和变量变换,第九章方差分析,完全随机设计资料的方差分析用于完全随机设计(成组设计)的多个样本均数的比较。该分析中仅涉及一个研究因素，此因素有K（K2）个水平或状态。,无论是实验还是观察，研究目的都是比较不同水平下，各组平均值之间的差别是否具有统计学意义，这种多个样本均数的比较都可用单因素方差分析。,例91某医生为研究一种四类降糖新药的疗效，以统一的纳入标准和排除标准选择了60名2型糖尿病患者，按完全随机设计方案将患者分为三组进行双盲临床试验。其中，降糖新药高剂量组21人，低剂量组19人，对照20人。对照组服用公认的降糖药物，治疗4周后测得其餐后2小时血糖的下降值（mmol/L）,结果如表91所示。问治疗4周后，餐后2小时血糖下降值的总体平均水平是否不同？,表912型糖尿病患者治疗4周后餐后2小时血糖下降值（mmol/L),表94方差分析试验数据示意图,用Xij表示第i个处理组的第j个观察值；i=1、2k;j=1,2,nini第i个处理组的例数总例数N=ni。第i个处理组的均数用表示，全部实验结果的总均数用表示。,方差分析的目的就是通过分析各个处理组均数之间的差别，推断它所代表的k个总体均数间是否存在差别。,一.总变异的分解,总变异（三个组共60人的餐后2小时血糖测定值大小不等；该变异既包含随机误差，又包含了三组用药即处理的不同）,组内变异：个体差异及血糖的随机测定误差。,组间变异：不同用药的影响及随机测量误差。,各种变异的表示方法,SS总总MS总,SS组内组内MS组内,SS组间组间MS组间,三者之间的关系：SS总=SS组内+SS组间总=组内+组间,变异比较：F=F=1：处理因素（三种药物）的作用相同。F1：处理因素（三种药物）的作用不相同F服从自由度为1与2的F分布，见:p435附表3F与比较：F，p处理因素的作用相同F，P处理因素的作用不同,表95成组设计方差分析的计算公式,SS总-SS组间,二、分析计算步骤（1）建立假设并确定检验水准H0：三个总体均数相等；即1=2=3H1：三个总体均数不等或不全相等=0.05（2）计算检验统计量F值,表96例91的方差分析表,3）.确定P值并作出推断结论以分子的自由度组间=2为1，分母的自由度组内=57为2，查附表3.1，方差分析用F界值表，F0.05(2,60)=3.15F0.01(2,60)=4.98,F=5.537F0.01(2,60)=4.98，P0.01。,在=0.05水准上拒绝H0，接受H1可以认为三种人的血糖水平不同。,注意：以上结论表明，总的说三种人的血糖水平有差别，但并不能说明任何两种人的血糖值均有差别。只能说可能至少有两组人的血糖值有差别，可能有的组间没有差别。要了解那些组均数间有差别，那些组均数间没有差别，需要进一步做两两比较。当k=2时，对同一资料，单因素方差分析等价于成组设计的t检验，且F=t2。,第一节方差分析的基本概念,第二节完全随机设计资料的方差分析,第三节随机区组设计资料的方差分析,第四节多个样本均数间的两两比较,第五节方差分析的前提条件和变量变换,第九章方差分析,随机区组设计的多个样本均数的比较可用多个样本均数比较的两因素方差分析。两因素是指主要的研究因素和配伍组（区组）因素，研究因素有k个水平，共有n个区组。,表97A.B.C三种方案处理后大白兔血中白蛋白减少量(g/L),一、总变异的分解,总变异,处理组间变异：A、B、C不同方案的影响及机测量误差。,区组间变异：既反映了十个区组不同的影响同时又包括了随机测量误差。,误差变异：个体差异及血白蛋白的随机测定误差。,SS总总,SS组内组内MS组内,SS组间组间MS组间,变异之间的关系：SS总=SS组内+SS组间+SS区间总=组内+组间+区间,变异间的关系,SS区间区间MS区间,统计量F的计算,F1=MS组间/MS组内F2=MS区间/MS组内,表910随机区组设计方差分析的计算公式,SS总-SS处理-SS配伍,误差,二.分析计算步骤,（1）建立假设并建立检验水准对于处理组：H0：三个总体均数全相等，即1=2=3H1：三个总体均数不等或不全相等对于区组：H0：十个总体均数全相等H1：十个总体均数不等或不全相等=0.05（2）计算检验统计量F值,表99例92的方差分析表,3).确定P值并作出推断结论以分子的自由度处理=2为1，分母的自由度误差=18为2，查附表3，方差分析用F界值表，F0.05(2,18)=2.62，F处理=32.639F0.05(2,18)=2.62，P0.05。在=0.05水准上不拒绝H0，还不能认为十个区组间有差别.,第一节方差分析的基本概念,第二节完全随机设计资料的方差分析,第三节随机区组设计资料的方差分析,第四节多个样本均数间的两两比较,第五节方差分析的前提条件和变量变换,第九章方差分析,在检验多组均数差别的无效假设H0时，常见的有以下两种情况：,1.检验某几个特定的总体均数是否相等，其无效假设称为部分无效假设即某几个组所对应的总体均数相等，H0：i=j（ij），如研究者对实验结果有一个大致的设想，在设计阶段就根据研目的或专业知识决定了某些均数间的两两比较，常见于事先有明确假设的证实性实验研究。,2.检验全部k个总体均数是否相等，其无效假设称为完全无效假设即所有各组所对应的总体均数都相等，H0：1=2=3=k。例如在研究设计阶段对实验结果知之不多的探索性研究，或经数据结果的提示后，才决定的多个均数间的两两比较，这类情况往往涉及到每两个均数的两两比较。,一、多个样本均数间每两个均数的比较：即SNK-q检验。,为两个对比组的样本均数,为比较两组差值的标准误，当各处理组例数相等时，,也相等，,MS误差为方差分析中算得的误差均方（组内均方），nA和nB分别为两对比组的样本例数。,计算的检验统计量为q值，q的分布与两比较组之间的组间跨度a及自由度有关。组间跨度a是指之间涵盖的均数个数（包括自身在内）,每个对比组所包含的组数（组间跨度）a=2，3，K。,根据检验统计量q值，组间跨度a，误差自由度误差及检验水准，查q界值表，确定P值。,例9-5对例9-1资料做两两比较,H0：任两对比组的总体均数相等，即A=BH1：任两对比组的总体均数不等，即AB=0.05,将三个样本均数从小到大排列，并编上组次：组次123均数9.19525.80005.4300组别高剂量组低剂量组对照组列出两两计算表(见表9-18)表9-18例9-1的SNK检验计算表,二、Dunnett检验：它适用于多个实验组与一个对照组的比较,根据算得的t值、误差的自由度、试验组数（K-1）及检验水准查Dunnett-t界值表，作出统计推断结论。,例9-6分析例9-2比较A、B两方案(均为实验组)与方案(对照组的总体均数是否不同？,A（B）方案与C方案相比H0：A（B）方案与C方案的总体均数相等H1：A（B）方案与C方案的总体均数不等=0.052)计算检验统计量表9-19例9-2的Dunnett-t检验计算表,3)确定P值，作出统计推断结论，查Dunnett界值表，得双侧t0.05=2.04,t0.01=2.84,P0.01,在=0.05水准上，拒绝H0，接受H1。故认为A方案与C方案及B方案与C方案的总体均数均有统计学意义。,第一节方差分析的基本概念,第二节完全随机设计资料的方差分析,第三节随机区组设计资料的方差分析,第四节多个样本均数间的两两比较,第五节方差分析的前提条件和变量变换,第九章方差分析,方差分析的前提条件为各样本是相互独立的随机样本，均服从正态分布且各样本的总体方差相等。在进行方差分析时，实际资料有时不能完全满足任何观察值都独立地来自具有等方差正态总体的假定，此时进行方差分析时，可能导致F值偏大，从而有增大第一类错误的危险。在样本例数较多的情况下，数据可看成近似正态分布，此时，方差分析对总体的非正态性并不苛求；当每组样本例数相等时，方差分析对于方差的非齐性并不苛求，故在方差分析时，最好采用每组例数相等的平衡设计方案。,变量变换的目的:1）使各组达到方差齐性；2）使资料转换为正态分布，以满足方差分析和t检验的应用条件。通常情况下，一种适当的函数转换可使上述两个目的同时达到。3）直线化。常用于曲线拟合。,1.对数变换：即将原始数据X的对数值作为新的分析数据。常用于