多个样本均数比较的方差分析

第九章方差分析讲授内容方差分析的基本思想及应用条件完全随机设计资料的方差分析随机区组设计资料的方差分析重复测量资料的方差分析多个样本均数间的多重比较基本思想及应用条件1、基本思想：根据资料设计的类型及研究目的，可将总变异分解为两个或多个部分，每个部分的变异可由某因素的作用来解释。通过比较可能由某因素所至的变异与随机误差，即可了解该因素对测定结果有无影响。2、应用条件：1）各样本是相互独立的随机样本2）各样本来自正态总体3）各处理组总体方差相等，即方差齐基本思想及应用条件完全随机设计资料的方差分析一、完全随机设计完全随机设计：(completelyrandomdesign)是采用完全随机化的分组方法，将全部试验对象分配到g个处理组（水平组），各组分别接受不同的处理，试验结束后比较各组均数之间的差别有无统计学意义，推论处理因素的效应。完全随机设计资料的方差分析例1某医生为研究一种四类降糖新药的疗效，以统一的纳入标准和排除标准选择了60名2型糖尿病患者，按完全随机设计方案将患者分为三组进行双盲临床试验。其中，降糖新药高剂量组21人、低剂量组19人、对照组20人。对照组服用公认的降糖药物，治疗4周后测得其餐后2小时血糖的下降值(mmol/L)，结果如表9-1所示。问治疗4周后，餐后2小时血糖下降值的三组总体平均水平是否不同？完全随机设计资料的方差分析分组方法：先将60名糖尿病患者从1开始到60编号；从随机数字表（附表15）中的任一行任一列开始，依次读取三位数作为一个随机数录于编号下；然后将全部随机数从小到大编序号(数据相同的按先后顺序编序号)，将每个随机数对应的序号记录；规定序号1-21为甲组，序号22-40为乙组，序号41-60为丙组。完全随机设计资料的方差分析二、变异分类与表达表4-1g个处理组的试验结果处理分组测量值统计量1水平X11X12…X1j…1nX1n11XS12水平X21X22…X2j…2nX2n22XS2…………………………g水平Xg1Xg2…Xgj…ggnXnggXSg完全随机设计资料的方差分析完全随机设计资料在进行统计分析时，需根据数据的分布特征选择方法，对于正态分布且方差齐同的资料，常采用完全随机设计的单因素方差分析(one-wayANOVA)或成组资料的t检验（g=2）；对于非正态分布或方差不齐的资料，可进行数据变换或采用Wilcoxon秩和检验。完全随机设计资料的方差分析完全随机设计资料的方差分析记总均数为，各处理组均数为，总例数为N＝nl+n2+…+ng，g为处理组数。11/ingijijXXN1/iniijijXXn完全随机设计资料的方差分析完全随机设计资料的方差分析1.总变异60名2型糖尿病患者的餐后2小时血糖Xij大小各不相同，与它们的总均数(overallmean)也不相同，这种变异称为总变异(totalvariation)。该变异既包含了随机误差(即2型糖尿病患者的个体差异和测量误差)，又包含了三组用药即处理的不同，其大小用所有数据(N=60)的方差即均方MS(meansquare)来表示。完全随机设计资料的方差分析完全随机设计资料的方差分析1．总变异全部测量值大小不同，这种变异称为总变异。总变异的大小可以用离均差平方和(sumofsquaresofdeviationsfrommean，SS)表示，即各测量值Xij与总均数差值的平方和，记为SS总。总变异SS总反映了所有测量值之间总的变异程度完全随机设计资料的方差分析完全随机设计资料的方差分析221111iinnggijijijijSSXXXC总211()/ingijijCXN1N总完全随机设计资料的方差分析完全随机设计资料的方差分析2．组内变异在同一处理组中，虽然每个受试对象接受的处理相同，但测量值仍各不相同，各组内2型糖尿病患者的餐后2小时血糖Xij大小各不相同，与本组的样本均数也不相同，这种变异称为组内变异（误差）。组内变异可用组内各测量值Xij与其所在组的均数的差值的平方和表示，记为SS组内,表示随机误差(含个体差异和测量误差)的影响。又称误差变异完全随机设计资料的方差分析完全随机设计资料的方差分析211()ingijiijSSXX组内Ng组内完全随机设计资料的方差分析完全随机设计资料的方差分析3．组间变异各处理组由于接受处理的水平不同，各组的样本均数(i＝1，2，…，g)也大小不等，三组2型糖尿病患者餐后2小时血糖的样本均数各不相同，它与总均数也不相同，这种变异称为组间变异。其大小可用各组均数与总均数的离均差平方和表示，记为SS组间，它反映了三组用药不同的影响(如处理确实有作用)，同时也包括了随机误差完全随机设计资料的方差分析完全随机设计资料的方差分析21211()()inijjggiiiiiXSSnXXCn组间1g组间完全随机设计资料的方差分析完全随机设计资料的方差分析各组均数之间相差越悬殊，它们与总均数的差值越大，SS组间就越大,反之SS组间越小。SS组间反映了各间的变异程度。存在组间变异的原因有：①随机误差(包括个体变异和测量误差)②处理的不同水平可能对试验结果的影响。完全随机设计资料的方差分析完全随机设计资料的方差分析基本思想：SSSSSS总组间组内总组间组内完全随机设计资料的方差分析完全随机设计资料的方差分析变异程度除与离均差平方和的大小有关外，还与其自由度有关，将各部分离均差平方和除以相应自由度，其比值称为均方差，简称均方(meansquare，MS)。SSMS组间组间组间SSMS组内组内组内完全随机设计资料的方差分析完全随机设计资料的方差分析如果各组样本的总体均数相等（H0：…），即各处理组的样本来自相同总体，无处理因素的作用（处理效应），则组间变异同组内变异一样，只反映随机误差作用的大小。组间均方与组内均方的比值称为F统计量MSFMS组间组内1组间2组内完全随机设计资料的方差分析完全随机设计资料的方差分析MSFMS组间组内误差误差处理组内变异组间变异完全随机设计资料的方差分析完全随机设计资料的方差分析F值（Fisher）接近于l，就没有理由拒绝H0；反之，F值越大，拒绝H0的理由越充分。数理统计的理论证明，当H0成立时，F统计量服从F分布。F分布有两个自由度，分子自由度为1，分母自由度为2，记为F~。12,F完全随机设计资料的方差分析完全随机设计资料的方差分析由F界值表（附表3），可查出按水准（一般取=0.05）F分布的单尾界值12(,)F，作为判断统计量F大小的标准。若根据试验结果计算的F值偏大，如120.05(,)FF时，则0.05P，拒绝H0，接受H1：i不全相等(i＝1，2，…，g)，说明各样本来自不同总体，即认为各样本的总体均数不等。反之，当F120.05(,)F时，P0.05，则不拒绝H0，还不能下各样本的总体均数不等的结论。完全随机设计资料的方差分析完全随机设计资料的方差分析通过上述变异的分解，可以看出，方差分析的基本思想就是根据试验设计的类型，将全部测量值总的离均差平方和及其自由度分解为两个或多个部分，除随机误差作用外，每个部分的变异可由某个因素的作用(或某几个因素的交互作用)加以解释，如组间变异SS组间可由处理因素的作用加以解释。通过比较不同变异来源的均方，借助F分布做出统计推断，从而推论各种研究因素对试验结果有无影响。完全随机设计资料的方差分析完全随机设计资料的方差分析应用条件：1）各样本是相互独立的随机样本2）各样本来自正态总体3）各处理组总体方差相等，即方差齐完全随机设计资料的方差分析完全随机设计资料的方差分析分析步骤(1)建立检验假设，确定检验水准H0：三个总体均数全相等，即μ1=μ2=μ3H1：三个总体均数不全相等，亦即至少有两个总体均数不等。即μ1≠μ2≠μ3或μ1=μ2≠μ3或μ1=μ3≠μ2或μ2=μ3≠μ1α=0.05完全随机设计资料的方差分析完全随机设计资料的方差分析完全随机设计资料的方差分析(2)计算检验统计量可根据下表的公式和前面表9-1下半部分数据来计算。也可用统计软件包如SAS或SPSS等进行计算，直接获得表9-4的方差分析表。完全随机设计资料的方差分析本例的资料是按完全随机设计方法获得的试验结果,可将总变异分解成组间变异和组内变异,并列方差分析表.表完全随机设计资料的方差分析表变异来源自由度SSMSF总变异N－1211ingijijXC组间g－1211()inijjgiiXCnSS组间组间MSMS组间组内组内N－gSSSS总组间SS组内组内完全随机设计资料的方差分析完全随机设计资料的方差分析按表中的公式计算各离均差平方和SS、自由度、均方MS和F值。完全随机设计资料的方差分析表9-4例9-1的方差分析表变异来源dfSSMSFP总变异591086.63组间2176.7688.385.5370.01组内（误差）57909.8715.96完全随机设计资料的方差分析(3)确定P值，作出推断结论以求F值时分子的自由度ν1=ν组间、分母的自由度ν2=ν组内查附表3的F界值表得P值。若F≥Fα(ν1,ν2)，则P≤α，按α水准，拒绝H0，接受H1，有统计学意义。本例：ν1=3−1=2，ν2=60−3=57。因附表3中ν2无57，故取最接近者ν2=60,得P0.01。按α=0.05水准，拒绝H0，接受H1，有统计学意义。可以认为2型糖尿病患者经药物(新药和标准药物)治疗4周，其餐后2小时血糖的总体平均水平不全相同，即三个总体均数中至少有两个不同。完全随机设计资料的方差分析完全随机设计资料的方差分析注意：方差分析的结果若拒绝H0，接受H1，不能说明各组总体均数两两间都有差别。如果要分析哪些两组间有差别，要进行多个均数间的多重比较（见本章第六节）。当g=2时，方差分析的结果与两样本均数比较的t检验等价，有。tF完全随机设计资料的方差分析随机区组设计资料的方差分析随机区组设计随机区组设计(randomizedblockdesign)又称为配伍组设计，是配对设计的扩展。具体做法是：先按影响试验结果的非处理因素（如性别、体重、年龄、职业、病情、病程等）将受试对象配成区组(block)，再分别将各区组内的受试对象随机分配到各处理或对照组。随机区组设计资料的方差分析与完全随机设计相比，随机区组设计的特点是随机分配的次数要重复多次，每次随机分配都对同一个区组内的受试对象进行，且各个处理组受试对象数量相同，区组内均衡。在进行统计分析时，将区组变异离均差平方和从完全随机设计的组内离均差平和中分离出来，从而减小组内平方和（误差平方和），提高了统计检验效率。若将区组作为另一处理因素的不同水平，随机区组设计等同于无重复观察的两因素设计。随机区组设计资料的方差分析随机区组设计资料的方差分析例9-2为探索丹参对肢体缺血再灌注损伤的影响，将30只纯种新西兰实验用大白兔，按窝别相同、体重相近划分为10个区组。每个区组3只大白兔随机采用A、B、C三种处理方案，即在松止血带前分别给予丹参2ml/kg、丹参1ml/kg、生理盐水2ml/kg，在松止血带前及松后1小时分别测定血中白蛋白含量(g/L)，算出白蛋白减少量如下表9-6所示，问A、B两方案分别与C方案的处理效果是否不同？随机区组设计资料的方差分析方法:先将小白鼠的体重从轻到重编号，体重相近的3只小白鼠配成一个区组，共10个区组。在随机数字表（附表15）中任选一行一列开始的2位数作为1个随机数；在每个区组内将随机数按大小排序；各区组中内序号为1的接受A方案、序号为2的接受B方案、序号为3的接受C方案.随机区组设计资料的方差分析随机区组设计资料的方差分析符合随机区组设计的资料在进行数据统计分析时也需根据数据的分布特征选择方法，对于正态分布且方差齐同的资料，应采用两因素方差分析(two-wayANOVA)或配