1、方差分析（方差分析（ANOVAANOVA）2023-1-91 1 1Y1n 2Y2n3Y3n4Y4n2023-1-92 2例子：例子：某研究者在某单位工作人员中进行了体重指数（BMI）抽样调查，随机抽取不同年龄组男性受试者各16名，测量了被调查者的身高和体重值，由此按照BMI=体重/身高2公式计算了体重指数，请问，不同年龄组的体重指数有无差异。项目项目18岁岁30岁岁4560岁岁21.6527.1520.2820.6628.5822.8818.8223.9326.49样本量样本量161616平均值平均值22.0725.9425.49标准差标准差8.978.117.192023-1-93 3一、
2、方差分析的基本思想一、方差分析的基本思想2023-1-94 45 5组间变异组间变异总变异总变异组内变异组内变异思想来源：思想来源：观察值总变异可以分解为组间变异和组内变异观察值总变异可以分解为组间变异和组内变异 2023-1-96 61.1.总变异总变异(Total variation):(Total variation):全部测量值全部测量值X Xij ij与总均与总均数数 间的差异间的差异 2.2.组间变异组间变异(between group variation):(between group variation):各组的均各组的均数数 与总均数与总均数 间的差异间的差异3.3.组内变异
3、组内变异（within group variation)within group variation)：每组的每：每组的每个测量值个测量值 与该组均数与该组均数 的差异的差异可用可用离均差平方和离均差平方和反映变异的大小反映变异的大小 XiXXijXiX2023-1-91.1.总变异总变异:所有测量值之间总的变异所有测量值之间总的变异程度程度，SS总1N总kinjijTiXXSS112)(2023-1-97 7 2 2组间变异：组间变异：各组均数与总均数的离各组均数与总均数的离均差平方和，均差平方和，SS组间1a组间SSSS组间组间反映了各组均数反映了各组均数 的变异程度的变异程度组间变异随机
4、误差组间变异随机误差+处理因素效应处理因素效应 21)(XXnSSikiiTRiX2023-1-98 8Na组内3 3组内变异：组内变异：用各组内各测量值用各组内各测量值X Xij ij与与其所在组的均数差值的平方和来表示，其所在组的均数差值的平方和来表示，SSSS组内组内SSSS组内组内反映反映随机误差随机误差的影响（个体差异和测量误差）。的影响（个体差异和测量误差）。211)(ikinjijeXXSSi9 9 均方差，均方均方差，均方(mean square(mean square，MSMS)变异程度除与离均差平方和的大小有关外，变异程度除与离均差平方和的大小有关外，还与其还与其自由度自由
5、度有关，由于各部分自由度不相等，有关，由于各部分自由度不相等，因此各部分离均差平方和不能直接比较，须将因此各部分离均差平方和不能直接比较，须将 各部分离均差平方和除以相应自由度，其比值各部分离均差平方和除以相应自由度，其比值 称为均方差，简称称为均方差，简称均方均方(mean square(mean square，MSMS)。组。组 间均方和组内均方的计算公式为间均方和组内均方的计算公式为:SSMS组间组间组间 SSMS组内组内组内 10101n2n3nVSVS1n2n3n1Y2Y3Y2Y3Y1Y11111n2n3n1n2nVSVS3Y1Y2Y2Y3Y1Y1212组间均方与组内均方比值越小，样
6、本越可能来组间均方与组内均方比值越小，样本越可能来源于同一个总体，比值越大，样本越可能不是源于同一个总体，比值越大，样本越可能不是来源于一个总体来源于一个总体 1313 二、二、F F 值与值与F F分布分布，如果各组样本的总体均数相等（H0成立），即各处理组的样本来自相同总体，无处理因素的作用，则组间变异同组内变异一样，只反映随机误差作用的大小。组间均方与组内均方的比值称为F统计量 F值接近于1，就没有理由拒绝H0；反之，F值越大，拒绝H0的理由越充分。数理统计的理论证明，当H0成立时，F统计量服从F分布。MSFMS组间组内1组间2组内1414F F 分布曲线分布曲线10,10215,121
7、5,52122121122/22/12121121)(222)(FFFf2023-1-91515回忆回忆t t分布和分布和t t检验检验2023-1-9161617172023-1-91818F F 界值表界值表2023-1-9二、完全随机设计方差分析二、完全随机设计方差分析(单因素方差分析单因素方差分析)2023-1-91919关于因素与水平关于因素与水平因素也称为处理因素（因素也称为处理因素（factorfactor）每一处理因素至少有两个水平每一处理因素至少有两个水平(level)(level)（也称（也称“处理组处理组”）。）。2023-1-92020完全随机设计：完全随机设计：将实验
8、对象随机分配到不同处理组的单因素设计方法。针对一个处理因素，通过比较该因素不同水平组均值，推断该处理因素不同水平组的均值是否存在统计学差异。2023-1-92121例例 在评价某药物耐受性及安全性的I期临床试验中，对符合纳入标准的30名健康自愿者随机分为3组每组10名，各组注射剂量分别为0.5U、1U、2U，观察48小时部分凝血活酶时间（s）试问不同剂量的部分凝血活酶时间有无不同？2023-1-922222023-1-92323方差分析步骤方差分析步骤 ：（1 1）提出检验假设，确定检验水准）提出检验假设，确定检验水准 H0:1=2=3 H1:1，2，3不全相同 a=0.052023-1-92
9、424 （2 2）计算检验统计量）计算检验统计量F F 值值 2023-1-92525 （3 3）确定）确定P P值，做出推断结论值，做出推断结论 F0.05(2，26)=2.52，FF0.05(2，26)，P0.05，拒绝 H0。三种不同剂量48小时部分凝血活酶时间 不全相同。2023-1-92626例子：某研究者在某单位工作人员中进行了体重指数（BMI）抽样调查，随机抽取不同年龄组男性受试者各16名，测量了被调查者的身高和体重值，由此按照BMI=体重/身高2公式计算了体重指数，请问，不同年龄组的体重指数有无差异。项目项目18岁岁30岁岁4560岁岁21.6527.1520.2820.662
10、8.5822.8818.8223.9326.49样本量样本量161616平均值平均值22.0725.9425.49标准差标准差8.978.117.19方差分析适合于任何多组独立均衡可比的数据方差分析适合于任何多组独立均衡可比的数据2023-1-92727基本步骤基本步骤（1 1）建立假设，确定检验水准）建立假设，确定检验水准H0：三个总体均数相等，即三组工作人员的体重指数总体均数相等H1：三个总体均数不等或不全相等a=0.052023-1-92828（2 2）计算检验统计量）计算检验统计量F F值值变异来源变异来源SSSS自由度（自由度（dfdf）MSMSF F组间143.406271.703
11、8.87组内363.86458.09总变异507.36472023-1-92929（3 3）确定）确定p p值，作出统计推断值，作出统计推断P2,45=3.20-3.218.87，本次F值处于F界值之外，说明组间均方组内均方比值属于小概率事件，因此拒绝H0，接受H1，三个总体均数不等或不全相等2023-1-93030方差分析的关键条件方差分析的关键条件第一、各组服从正态分布！第一、各组服从正态分布！第二、各组符合方差齐性！第二、各组符合方差齐性！第三、独立性第三、独立性2023-1-93131方差齐性检验方差齐性检验Bartlett检验法Levene F 检验最大方差与最小方差之比3，初步认为
12、方差齐同。2023-1-93232问题：问题：不符合条件怎么办？不符合条件怎么办？第一招：数据转换第一招：数据转换 方差齐性转换；正态性转换方差齐性转换；正态性转换第二招：特别分析方法第二招：特别分析方法 非参数检验非参数检验2023-1-93333三、多个样本均数的两两比较三、多个样本均数的两两比较2023-1-93434方差分析能说明什么问题？方差分析能说明什么问题？不拒绝不拒绝H H0 0，表示拒绝总体均数相等的证据不，表示拒绝总体均数相等的证据不足足 分析终止分析终止 拒绝拒绝H H0 0，接受，接受H H1 1,表示总体均数不全相等表示总体均数不全相等 哪两两均数之间相等？哪两哪两两
13、均数之间相等？哪两两均数之间不等？两均数之间不等？需要进一步作多重比较需要进一步作多重比较2023-1-93535能否用T检验呢当有k个均数需作两两比较时，比较的次数共有c=k!/(2!(k-2)!)=k(k-1)/2设每次检验所用类错误的概率水准为，累积类错误的概率为，则在对同一实验资料进行c次检验时，在样本彼此独立的条件下，根据概率乘法原理，其累积类错误概率与c有下列关系：1(1)c 例如，设0.05，c=3(即k=3)，其累积类错误的概率为1(1-0.05)3=1-(0.95)3=0.1432023-1-93636多重比较的方法：多重比较的方法：SNK检验（q 检验）：探索性研究，进行两
14、两比较。LSD-t 检验：证实性检验，可认为LSD法是最灵敏的Turkey 检验方法，探索性研究，要求样本量相同。Duncan 检验方法，探索性研究Dunnet 检验方法，证实性检验，常用于多个试验组与一个对照组间的比较。2023-1-93737 例例1 1 在肾缺血再灌注过程的研究中，将在肾缺血再灌注过程的研究中，将3636只雄性大鼠随机等只雄性大鼠随机等分成三组，分别为正常对照组、肾缺血分成三组，分别为正常对照组、肾缺血6060分组和肾缺血分组和肾缺血6060分分再灌注组，测得各个体的再灌注组，测得各个体的NONO数据见数据文件数据见数据文件no.savno.sav，试问各，试问各组的组的
15、NONO平均水平是否相同？平均水平是否相同？单因素方差分析单因素方差分析2023-1-93838分析：分析：对于单因素方差分析，其资料在对于单因素方差分析，其资料在SPSSSPSS中的数据结构应当由两中的数据结构应当由两列数据构成，其中一列是观察指标的变量值，另一列是用以列数据构成，其中一列是观察指标的变量值，另一列是用以表示分组变量。实际上，几乎所有的统计分析软件，包括表示分组变量。实际上，几乎所有的统计分析软件，包括SASSAS，STATASTATA等，都要求方差分析采用这种数据输入形式，等，都要求方差分析采用这种数据输入形式，这一点也暗示了方差分析与线性模型间千丝万缕的联系。这一点也暗示
16、了方差分析与线性模型间千丝万缕的联系。单因素方差分析单因素方差分析2023-1-93939 预分析（重要）：检验其应用条件预分析（重要）：检验其应用条件单因素方差分析单因素方差分析选择选择data data 中的中的split filesplit file，出现如下对话框：，出现如下对话框：2023-1-94040单因素方差分析单因素方差分析2023-1-94141单因素方差分析单因素方差分析2023-1-94242单因素方差分析单因素方差分析v 这里仅取其中一组结果，表明该资料符合这里仅取其中一组结果，表明该资料符合分组正态性的条件。分组正态性的条件。2023-1-94343单因素方差分析单
17、因素方差分析注意分组检验正态性后，要先回到注意分组检验正态性后，要先回到datadata菜单下的菜单下的split file split file，如下操作取消拆分后才能进行后续的方差分析：如下操作取消拆分后才能进行后续的方差分析：2023-1-94444单因素方差分析单因素方差分析2023-1-94545单因素方差分析单因素方差分析选入分组变量选入分组变量选入因变量选入因变量2023-1-94646给出各组间样本给出各组间样本均数的折线图均数的折线图指定进行方差指定进行方差齐性检验齐性检验单因素方差分析单因素方差分析2023-1-94747结果分析结果分析单因素方差分析单因素方差分析（1 1
18、）方差齐性检验方差齐性检验v Levene Levene方法检验统计量为方法检验统计量为3.2163.216，其，其P P值为值为0.0530.053，可认，可认为样本所来自的总体满足方差齐性的要求。为样本所来自的总体满足方差齐性的要求。Test of Homogeneity of Variancesno3.216233.053Levene Statisticdf1df2Sig.2023-1-94848单因素方差分析单因素方差分析结果分析结果分析（2 2）方差分析表方差分析表v 第第1 1列为变异来源，第列为变异来源，第2 2、3 3、4 4列分别为离均差平方和、自列分别为离均差平方和、自由度
19、、均方，检验统计量由度、均方，检验统计量F F值为值为5.5645.564，P P0.0080.008，组间均数差，组间均数差别统计学意义，可认为各组的别统计学意义，可认为各组的NONO不同。不同。ANOVAno46925.950223462.9755.564.008139157.6334216.898186083.635Between GroupsWithin GroupsTotalSum ofSquaresdfMean SquareFSig.变变异异来来源源2023-1-94949单因素方差分析单因素方差分析结果分析结果分析（3 3）各组样本均数折线图各组样本均数折线图Means plot
20、s Means plots 选项给出，更直观。选项给出，更直观。注意注意：当分组变量体现出顺序的趋势时，绘制这种折线图可以提示：当分组变量体现出顺序的趋势时，绘制这种折线图可以提示我们选择正确的趋势分析模型。我们选择正确的趋势分析模型。2023-1-95050通过以上分析得到了拒绝通过以上分析得到了拒绝H H0 0的结论，但实际上单因素方差分的结论，但实际上单因素方差分析并不这样简单。在解决实际问题时，往往仍需要回答多个析并不这样简单。在解决实际问题时，往往仍需要回答多个均数间到底是哪些存在差异。虽然结论提示不同组别个体的均数间到底是哪些存在差异。虽然结论提示不同组别个体的NONO量不同，但研
21、究者并不知道到底是三者之间均有差别，量不同，但研究者并不知道到底是三者之间均有差别，还是某一组与其他两组有差别。这就应当通过两两比较（多还是某一组与其他两组有差别。这就应当通过两两比较（多重比较）进行考察。重比较）进行考察。均数两两比较方法均数两两比较方法2023-1-95151直接校正检验水准直接校正检验水准（相对粗糙）（相对粗糙）专用的两两比较方法：专用的两两比较方法：计划好的多重比较（计划好的多重比较（Planned ComparisonsPlanned Comparisons）非计划的多重比较（非计划的多重比较（PostPostHoc ComparisonsHoc Comparison
22、s）均数两两比较方法均数两两比较方法ContrastsContrasts按钮按钮Post HocPost Hoc按钮按钮2023-1-95252点击单因素方差分析主对话框中的点击单因素方差分析主对话框中的Post HocPost Hoc按钮，总共有按钮，总共有1414种两两比较的方法，如下：种两两比较的方法，如下：均数两两比较方法均数两两比较方法2023-1-95353LSDLSD法：法：最灵敏最灵敏，会犯假阳性错误；，会犯假阳性错误；SidakSidak法：比法：比LSDLSD法保守；法保守；BonferroniBonferroni法：比法：比SidakSidak法更为保守一些；法更为保守一
23、些；ScheffeScheffe法：多用于进行比较的两组间样本含量不等时；法：多用于进行比较的两组间样本含量不等时；DunnetDunnet法：常用于多个试验组与一个对照组的比较；法：常用于多个试验组与一个对照组的比较；S-N-KS-N-K法：寻找同质亚组的方法；法：寻找同质亚组的方法；TurkeyTurkey法：最迟钝，要求各组样本含量相同；法：最迟钝，要求各组样本含量相同；DuncanDuncan法：与法：与SidakSidak法类似。法类似。均数两两比较方法均数两两比较方法2023-1-95454仍以例仍以例1 1为例，为例，LSDLSD法的输出格式：法的输出格式：均数两两比较方法均数两
24、两比较方法结果分析结果分析2023-1-95555仍以例仍以例1 1为例，为例，SNKSNK法的输出格式：法的输出格式：结果分析结果分析均数两两比较方法均数两两比较方法v 该方法的目的是寻找同质子集，故各组在表格的纵向上，均该方法的目的是寻找同质子集，故各组在表格的纵向上，均数按大小排序，然后根据多重比较的结果将所有的组分为若干数按大小排序，然后根据多重比较的结果将所有的组分为若干个子集，子集间有差别，子集内均数无差别。个子集，子集间有差别，子集内均数无差别。2023-1-95656 当各组样本含量不同，选择当各组样本含量不同，选择S Scheffecheffe法，得结果：法，得结果：均数两两
25、比较方法均数两两比较方法结果分析结果分析Multiple ComparisonsDependent Variable:noScheffe13.6125026.51068.877-54.338981.563982.48167*26.51068.01414.5303150.4330-13.6125026.51068.877-81.563954.338968.86917*26.51068.046.9178136.8205-82.48167*26.51068.014-150.4330-14.5303-68.86917*26.51068.046-136.8205-.9178(J)group231312(
26、I)group123MeanDifference(I-J)Std.ErrorSig.Lower BoundUpper Bound95%Confidence IntervalThe mean difference is significant at the.05 level.*.2023-1-95757 假设在调查的设计阶段，就计划好了第二组和第一组，以及第三组假设在调查的设计阶段，就计划好了第二组和第一组，以及第三组和第一组的比较，可以使用主对话框中的和第一组的比较，可以使用主对话框中的contrast contrast 按钮按钮实现。实现。v 在在coefficientscoefficien
27、ts后面的框中输入后面的框中输入1 1，-1-1，0 0，每次输入后点击，每次输入后点击addadd，就可，就可以比较第一组和第二组的以比较第一组和第二组的NONO；再点击；再点击nextnext按钮，继续输入下一个组合，按钮，继续输入下一个组合，即即0 0，-1-1，1 1。均数两两比较方法均数两两比较方法2023-1-95858均数两两比较方法均数两两比较方法结果分析结果分析可见，第一个组合无统计学意义，而第二个组合有显著性差异。可见，第一个组合无统计学意义，而第二个组合有显著性差异。Contrast Coefficients1-1010-1Contrast12123groupContra
28、st Tests13.6125 26.51068.51333.61182.4817 26.510683.11133.00413.6125 22.16579.61421.981.54582.4817 28.306862.91418.900.009Contrast1212Assume equal variancesDoes not assume equalvariancesnoValue ofContrast Std.ErrortdfSig.(2-tailed)2023-1-95959四、多因素方差分析四、多因素方差分析2023-1-96060多因素方差分析多因素方差分析一个因素（水平间独立）单因
29、素方差分析两个因素（水平间独立或相关）多(两)因素方差分析一个个体多个测量值重复测量资料的方差分析目的：用这类资料的样本信息来推断各处理组间多个总体均数的差别有无统计学意义。2023-1-96161（一）随机区组方差分析（水平间独立两因素）（一）随机区组方差分析（水平间独立两因素）随机区组设计：随机区组设计：先将受试对象按条件相同或相近分成m个区组（或配伍组），每个区组有k个受试对象，再将其随机地分配到k个处理组中，称之为随机区组设计，属于两因素方差分析。区组（区组（B B）处理因素（处理因素（A A）X11X21Xk1区组1X12X22Xk2区组2X13X23Xk3区组mX1mX2mXkm2
30、023-1-96262总变异分解 SS总=SS组间+SS组内（完全随机）SS总=SS组间+（SS区组+SS组内）（随机区组）在存在着区组的情况下，随机区组设计的方差分析效率要高于单因素方差分析2023-1-96363计算公式计算公式总变异：，自由度N-1。处理组变异：，自由度k-1区组变异：，自由度ni-1组内：SS总-SS处理-SS区组，自由度N-k-ni-1 Nxx22)(kiiixxn12)(injjjxxn12)(2023-1-96464举例 为探讨Rgl对镉诱导大鼠睾丸损伤的保护作用，某研究者将同一窝别的3只大鼠随机地分到T1、T2、T3三组，进行不同处理,共观察了10个窝别大鼠的睾
31、丸MT含量（g/g）。试问不同处理对大鼠MT含量有无影响？2023-1-965652023-1-96666方差分析方差分析 =3.55，P0.05，三组大鼠 MT含量的总体均值不全相同。0.05(2,18)F0.05(2,18)FF2023-1-96767随机区组设计方差分析前提条件：随机区组设计方差分析前提条件：数据服从正态分布数据要求具有独立性数据要求具有方差齐性2023-1-96868举例举例对三组病人进行降血压药物治疗，一组病人为安慰剂，一组病人为常规药物，另外一组为新研制药物。分别在服药5、10、20天测量血糖这种例子看似随机区组设计这种例子看似随机区组设计实则不然，数据之间存在相关性实则不然，数据之间存在相关性需采用重复测量方差分析需采用重复测量方差分析2023-1-96969小结小结完全随机设计t检验：单因素两样本分析完全随机设计方差分析：单因素多样本分析随机区组设计方差分析：两因素多样本分析2023-1-97070
