卫生统计学专题九：方差分析.doc

专题九 方差分析方差分析（anova），又称f检验，是一种以分析数据的变异为基础，以f值为统计量的计量资料的假设检验方法。用于多个均数的比较。方差分析的基本思想和应用条件 根据资料的设计类型（即变异的不同来源），将全部观察值之间的总变异按设计类型分解为两个或多个组成部分，通过比较不同变异来源的均方，借助f分布做出统计推断。方差分析的应用条件：各观察值相互独立，且每一水平下的观察值均服从正态分布。各总体方差相等，即具有方差齐性。对于不满足方差分析应用条件的资料，可采用以下两种方式进行处理：变量变换，是指满足方差分析的基本假定。应用非参数统计分析方法。完全随机设计的anova 各组样本均数各不相等，这种差异可能由两种原因引起：随机误差；处理因素。完全随机设计的变异分解：完全随机设计（成组设计）资料的总变异可以分解为组间变异和组内变异。总变异：ss总，用所有观察值与均数的离均差平方和表示。ss总=，总=n-1（n为总例数）组间变异：ss组间，用各组均数与总均数的离均差平方和表示。ss组间=，组间=k-1（k为处理组数）组内变异：ss组内，用各组内每个测量值xij与该组均数得离均差平方和表示，仅反映随机误差，又称误差变异。ss组内=，组内=n-k三种变异及相关自由度的关系为：ss总=ss组间+ss组内，总=组间+组内均方ms组间=ss组间组间；ms组内= ss组内组内方差分析的统计量f：f=ms组间ms组内f界值表：纵标目为组间自由度1，横标目为组内自由度2，表中给出了=0.05和=0.01时供方差分析用的单侧f界值，用f，（1，2）表示。 若ff，（1，2），则p，按水准拒绝h0，接受h1，差别有统计学意义，可以认为总体均数不等或不全等（处理因素有效应）；反之，则差别无统计学意义，尚不能认为总体均数不等或不全等（尚不能认为处理因素有效应）。完全随机设计方差分析计算公式变异来源ssmsf组间变异（处理）k-1组内变异（误差）n-k总变异-n-1【小结】完全随即设计资料的总变异可分解为组间变异和组内变异两部分。方差分析常用于三个及以上均数的比较，当用于两个均数的比较时，同一资料所得结果与t检验等价，即f=t2方差分析的结果若拒绝h0，而接受h1，只说明多个总体均数中至少有两个均数不同，不能说明任意两个总体均数都有差别；若要比较那两组均数有差别，需进一步进行多个均数间的两两比较。随机区组设计的anova 随机区组设计的变异分解：ss总=ss处理+ss区组+ss误差，总=处理+区组+误差用k表示处理组数，b表示区组数，实验观察值xij下标i（i=1，2，k）表示组别，下标j（j=1，2，b）表示区组序号。完全随机设计方差分析计算公式变异来源ssmsf处理组k-1区组b-1误差ss总-ss处理-ss区组n-k-b-1或（k-1）（b-1）总变异-n-1随机区组设计资料的anova的基本步骤：建立假设检验，确定检验水准： 对于处理组：h0：1=2=k，即相同；h1：1、2、k不等或不全等，即不同或不全同。 对于区组：h0：k个区组的总体均数相等；h1：k个区组的总体均数不等或不全等。计算检验统计量：f= ms处理ms误差；f= ms区组ms误差确定p值，做出统计推断。 对处理组：以1=k-1，2=n-k-b-1或（k-1）（b-1），查f界值表，确定p值，做出统计推断。对区组：以1=b-1，2= n-k-b-1或（k-1）（b-1），查f界值表，确定p值，做出统计推断。【注意】随机区组设计资料的总变异可以分解为处理组间变异、区组间变异和误差变异三个部分。随机区组设计与完全随机设计相比，由于利用区组控制了可能的混杂因素，并在进行方差分析时，将区组间变异从原组内变异中分离出来，当区组间变异有统计学意义时，由于减少了误差均方使处理间的f值更容易出现显著性，从而提高了实验效率。随机区组设计是配对设计的扩展，当k=2时，随机区组设计的方差分析与配对资料的t检验完全等价，即f=t2。（平方关系）多个样本均数的两两比较 多个样本均数经方差分析后，若有统计学意义，需用两两比较的方法进一步确定哪些均数相等，哪些均数不等。多组均数假的两两比较不能用t检验，会增大型错误。两两比较的方法较多，常用的是多重比较，多重比较法有两种情况：探索性研究，即事先没有计划的任意两个均数的事后比较：snk-q检验、bonfferoni t检验。证实性研究，即事先有明确假设，用于多组中某一对或某几对在专业上有特殊意义的均数间的比较：dunnett-t检验、lsd-t检验、bonfferoni t检验等。snk-q检验与lsd-t检验的比较snk-q检验lsd-t检验应用探索性研究证实性研究h0a=b，即任两对比组的总体均数相等。a=b，即所研究的两对比组的总体均数相等。检验统计量q=t=误差=误差【注意】统计软件可以同时进行十多种多重比较检验，因此在应用时，应根据统计设计和专业知识来确定采用哪一种方法，不能多种方法一起使用，然后选取“有利”的结果。交叉设计的anova 两阶段交叉设计的变异分解： 二阶段交叉设计方差分析的总变异可以分解为处理间变异、阶段间变异、个体变异和误差变异四个部分。ss总=ss处理+ss阶段+ss个体+ss误差，总=处理+阶段+个体+误差两阶段交叉设计的方差分析基本步骤：建立建设检验，确定检验水准。计算检验统计量f值：ms阶段ms误差、ms个体ms误差、ms处理ms误差；确定p值，做出统计推断。【说明】交叉设计可以采用完全随机设计或配对设计方法来安排其受试对象。交叉实验的处理是单因素的，但影响实验结果的因素还有非人为控制的受试者之间的个体差异和试验阶段这两个因素。因此交叉设计试验实际上是一个试验因素和两个重要的非实验因素的多因素实验。析因设计的anova 析因设计资料的方差分析包括主效应分析、交互效应分析和单独效应分析三个层次。单独效应：是指其他因素的水平固定时，同一因素不同水平间的差别。主效应：是指某一因素各水平间的平均差别。交互效应：是指当某一因素的各单独效应随某一因素水平的变化而变化时则称这两个因素间存在交互效应。两因素析因设计即ij析因设计，表示有两种处理因素，第一种处理因素有i个水平，第二种因素有j个水平，共有ij种处理组合。析因设计的变异分解：ss总=ss处理+ ss误差=（ssa+ssb+ssab）+sse，总=处理+误差=（a+b+ab）+e若处理因素a有a个水平，处理因素b有b个水平，每个处理因素组含有n个受试对象，则全部收拾对象的总数n=abn。用xijk表示每个受试对象的观察值，其中i（i=1，2，a）表示a因素的水平。j（j=1，2，b）表示b因素有b个水平；k（k=1，2，n）表示a、b两因素不同水平组合下受试对象的序号。析因设计方差分析计算表变异来源ssmsf处理ab-1aa-1bb-1abss处理-ssa-ssb（a-1）（b-1）误差ss误差=ss总-ss处理n-ab或ab（n-1）总变异n-1或（abn-1）两因素析因设计资料方差分析的基本步骤：对于因素a、b、交互作用ab分别建立假设检验，确定检验水准。 检验假设：对于因素a：h0：a因素各水平总体均数相等；h1：a因素各水平总体均数不相等； 对于因素b：h0：b因素各水平总体均数相等；h1：b因素各水平总体均数不相等； 对于交互作用ab：h0：a、b无交互作用；h1：a、b有交互作用；计算检验统计量：列表计算各组均数及方差分析表。确定p值，做出统计推断。【说明】析因设计是将两个或多个实验因素的各水平进行交叉分组、全面实验的方法。因此析因设计的方差分析被广泛用于需要分析交互效应和选择最佳组合的实验研究中。析因设计不但可以分析主效应和交互效应，也可以分析单独效应，故效率较高，但当因素太多时，所需的样本含量会很大。对析因设计资料，应先分析交互效应，若交互效应有统计学意义，方进一步分析个因素的单独效应。反之，若交互效应无统计学意义，则因素之间的作用相互独立，分析某一因素的作用只需考虑因素的主效应。重复测量设计的anova 重复测量设计是指给予一种或多种处理后，在多个时间点上从同一受试对象重复获得指标的观察值。重复测量数据的两因素多水平设计，两因素是指处理因素和处理时间；多水平是指处理有k（1）个水平；测量时间有i（i2）个水平（时间点），即每个观察对象有i个重复测量数据。重复测量设计的变异分解：两因素重复测量设计的总变异包括两部分：受试对象间变异、受试对象内变异。ss总=ss受试对象间+ss受试对象内=（ss处理+ss个体间误差）+ss时间+ss处理与时间交互作用+ss个体内误差总=受试对象间+受试对象内=（处理+个体间误差）+时间+处理与时间交互作用+个体内误差重复测量设计方差分析的基本步骤：对于处理因素、时间因素、交互作用分别建立检验假设，确定检验水准；计算检验统计量：spss软件确定p值，做出统计推断；重复测量设计方差分析的前提条件：正态性和方差齐性；协方差阵地球形性或复合对称性。球对称性：mauchly检验。【注意】若资料不满足球对称性，则方差分析的f值有偏，通常会增大型错误的概率。资料满足球对称性时，可采用重复测量设计资料的单变量方差分析方法；若资料不满足球对称性，可用greenhouse-geisser、huynh-feldt或lower-bound三种“球对称”系数对自由度进行校正或采用多变量方差分析的方法。【小结】完全随机设计资料的方差分析，若处理无作用，则f值理论上应等于1。重复测量设计资料和随