讲稿方差分析.doc

方差分析（ANOVA）一. 方差分析的概念方差分析（analysis of variance, ANOVA）又称变异数分析或F检验，其目的是推断两组或多组资料的总体均数是否相同，检验两个或多个均数的差异是否有统计学意义。我们要学习的主要内容包括单因素方差分析，即完全随机（成组）设计的方差分析和两因素方差分析即随机区组（配伍组）设计的方差分析。方差分析除用于两个或多个均数的比较外，还可分析两个或多个研究因素的交互作用以及回归方程的线性假设检验等。根据资料设计的类型及研究目的，可将总变异分解为两个或多个部分，每个部分的变异可由某因素的作用来解释，通过比较可能由某因素所致的变异与随机误差（如组内变异），即可了解该因素对测定结果有无影响。此外，多个均数间的比较，由于涉及的对比组数大于2，若仍用两组比较的检验，对每两个对比组作比较，会使犯第一类错误的概率（a）增大。例如：有5个均数的比较，可以用10次t检验比较，若每次比较的检验水准a=0.05，则每次比较不犯第一类错误的概率为（10.05）=0.95，那么10次比较均不犯第一类错误的概率为0.9510=0.5987，这时犯第一类错误的概率，也就是总的显著性水准变为10.9510=0.4013，比0.05大多了，大大增加了犯第一类错误的概率，可能把本来无差别的两个总体均数判为有差别。因此，多个均数比较不宜用前述检验分别作两两比较。二方差分析的应用条件方差分析的应用条件如下：1正态分布，要求资料服从正态分布，偏态分布资料不适用方差分析。对偏态分布的资料应考虑使用秩和检验，或用变量变换方法将数据转变为正态或接近正态后再进行方差分析。2方差齐性，即各组的方差要相等，若组间方差不齐则不适用方差分析。故方差分析前，应做多个方差的齐性检验。3各样本是相互独立的随机样本。三完全随机（成组）设计的方差分析下面我们用一个简单的例子来说明其基本思想：三种方案治疗后血红蛋白增加量（g/L）方案观察值nA243625142634237B2018171019246C2011630-1458三种人的载脂蛋白测定结果糖尿病IGT正常人85.796.0144.0105.2124.5117.0109.5105.1110.096.076.4109.0115.295.3103.095.3110.0123.0110.095.2127.0100.099.0121.0125.0120.0159.0111.0115.0106.5从以上资料可以看出，30个观对象的载脂蛋白量各不相同，就其变异情况来分析，其总变异有以下两个来源：组内变异：每组内部各个观察对象的载脂蛋白量各不相同，就其原因主要是由于随机误差所引起的。这种变异，叫做组内变异。组间变异：三组患者的载脂蛋白量的均数不相等，导致不相等的原因，一是不同人群（处理不同）的影响，二是随机误差所引起。这种变异，叫做组间变异。 方差分析就是对数据中所包含的变异进行分析，并找出总变异中的各组成部分，如上述完全随机设计的例子，其总变异与及各组部分的关系为：SS总= SS组间+SS组内 v总= v组间+v组内 完全随机设计的方差分析，其假设如下：H0: 各总体均数相等H1: 各总体均数不等或不全相等a=0.05完全随机设计的方差分析变异分解表变异来源SS自由度vMS（方差，均方）F值P值组间SS组间v组间MS组间= SS组间/ v组间MS组间/ MS组内FFa,n1,n2Pa组内SS组内v组内MS组内= SS组内/ v组内F a合计SS总v总完全随机设计的方差分析，其F值就是用组间均方去除以组内均方（MS组间/MS组内）的商。若Fa，各均数间的差异没有统计学意义；若F Fa,n1,n2 则Pa，各总体均数不等或不全相等，即各组均数间的差异有统计学意义，此时，还需进行后述的多重比较。四. 随机区组(配伍)设计的方差分析对随机区组(配伍)设计的多个均数比较，应采用随机区组(配伍)设计的方差分析，即两因素方差分析。三种营养素喂养小白鼠的增重（g）区组号A营养素B营养素C营养素150.158.264.5247.848.562.4353.153.858.6463.564.272.5571.268.479.3641.445.738.4761.953.051.2842.239.846.2随机区组(配伍)设计的方差分析的基本步骤相同，只是随机区组(配伍)设计方差分析变异的分解方式不同，对随机区组（配伍组）设计的资料，总变异可分解为处理组变异、区组组变异和随机误差，其总变异与及各组部分的关系为：SS总=SS处理+SS区组+SS误差v总=v处理+v区组+v误差随机区组设计的方差分析，可有两个假设，具体如下：各处理间比较的假设H0: 各总体均数相等H1: 各总体均数不等或不全相等a=0.05各区组间比较的假设H0: 各总体均数相等H1: 各总体均数不等或不全相等a=0.05随机区组设计的方差分析变异分解表变异来源SSv自由度MS（方差，均方）F值P值处理SS处理v处理MS处理= SS处理/ v处理MS处理/ MS误差FFa,n1,n2Pa区组SS区组v区组MS区组= SS区组/ v区组MS区组/ MS误差F a合计SS总v总随机区组设计的方差分析，可计算得到两个F值，即F1=MS处理/ MS误差和F2=MS区组/ MS误差，若F1a，各均数间的差异没有统计学意义；若F1Fa,n1,n2 则Pa，各组均数间的差异有统计学意义；若F2a，各区组间的差异没有统计学意义；若F2Fa,n1,n2 则Pa，各区组间的差异有统计学意义。五多个均数的两两比较经过方差分析若拒绝了检验假设，只能说明多个样本所代表的总体均数不相等或不全相等。若要得到各组均数间更详细的信息，应在方差分析的基础上进行多个均数间的两两比较，即多个均数间的多重比较。多个均数间的多重比较（multiple comparison），又称两两比较。由于涉及的对比组数大于2，若仍用前述比较的检验，对每两个对比组作比较，会使犯第一类错误的概率（a）增大。例如：有5个均数的比较，可以用10次t检验比较，若每次比较的检验水准a=0.05，则每次比较不犯第一类错误的概率为（10.05）=0.95，那么10次比较均不犯第一类错误的概率为0.9510=0.5987，这时犯第一类错误的概率，也就是总的显著性水准变为10.9510=0.4013，比0.05大多了，大大增加了犯第一类错误的概率，可能把本来无差别的两个总体均数判为有差别。因此，多重比较不宜用前述检验分别作两两比较。多个均数的比较一般分为二种情况：一是在研究设计阶段未预先考虑或未预料到，经数据结果的提示后，才决定的多个均数间的两两比较，常见于探索性研究（exploratory research），这类情况下往往涉及到每两个均数的比较。二是在设计阶段就根据研究目的或专业知识而决定的某些均数间的两两比较，常见于事先有明确假设的证实性实验研究（con-firmatory research），这种情况又可分两种不同的比较，即多个处理组与对照组的比较，如处理后不同时间与处理前的比较等，或某一对或某几对专业上有特殊意义的均数间的差别比较；现分述这两种情况下适用的统计方法。1SNK检验SNK (Student-Newman-kueuls)检验，亦称q检验，主要用于多个均数两两间的全面比较。其基本步骤为： 建立检验假设样本均数从大到小排序计算q值q = SEAB = 确定a值，查q界值表判断结果。2Dunnett检验Dunnett检验适用于多个试验组均数与一个对照组均数间的比较。q = SEi0= 查Dunnett检验用表，判断结果。3LSD检验最小显著差法(LSD,least significant difference)，主要用于检验某一对或某几对专业上有特殊意义的均数间的差别比较。t = SEAB= 查t检验用表，判断结果。六、方差分析注意事项1方差分析是一大类分析方法，除用于两个或多个均数的比较外，还有交叉设计、析因设计、拉丁方、正交等设计的方差分析，还可分析两个或多个研究因素的交互作用以及回归方程的线性假设检验等。做方差分析时，可根据资料设计的类型及研究目的，可将总变异分解为两个或多个部分，每个部分的变异可由某因素的作用来解释，通过比较可能由某因素所致的变异与误差（或组内）变异，即可了解该因素对测定结果有无影响。2方差分析的应用条件与t检验相似，即要求资料服从正态分布，方差齐性，即各组的方差要相等。故方差分析前，应做多个方差齐性检验，正态