SPSS方差分析.docx

1、SPSS方差分析第六章 方差分析第一节 Simple Factorial过程6.1.1 主要功能6.1.2 实例操作第二节 General Factorial过程6.2.1 主要功能6.2.2 实例操作第三节 Multivarite过程6.3.1 主要功能6.3.2 实例操作 方差分析是R.A.Fister发明的，用于两个及两个以上样本均数差别的显著性检验。由于各种因素的影响，研究所得的数据呈现波动状，造成波动的原因可分成两类，一是不可控的随机因素，另一是研究中施加的对结果形成影响的可控因素。方差分析的基本思想是：通过分析研究中不同来源的变异对总变异的贡献大小，从而确定可控因素对研究结果影响力

2、的大小。 方差分析主要用于：1、均数差别的显著性检验，2、分离各有关因素并估计其对总变异的作用，3、分析因素间的交互作用，4、方差齐性检验。第一节 Simple Factorial过程6.1.1 主要功能 调用此过程可对资料进行方差分析或协方差分析。在方差分析中可按用户需要作单因素方差分析（其结果将与第五章第四节相同）或多因素方差分析（包括医学中常用的配伍组方差分析）；当观察因素中存在有很难或无法人为控制的因素时，则可对之加以指定以便进行协方差分析。6.1.2 实例操作 例6-1下表为运动员与大学生的身高（cm）与肺活量（cm3）的数据，考虑到身高与肺活量有关，而一般运动员的身高高于大学生，为

3、进一步分析肺活量的差异是否由于体育锻炼所致，试作控制身高变量的协方差分析。运 动 员大 学 生身高肺活量身高肺活量184.9167.9171.0171.0188.0179.0177.0179.5187.0187.0169.0188.0176.7179.0183.0180.5179.0178.0164.0174.043003850410043004800400054004000480048004500478037005250425048005000370036004050168.7170.8165.0169.7171.5166.5165.0165.0173.0169.0173.8174.0170

4、.5176.0169.5176.3163.0172.5177.0173.0345041003800330034503250360032003950400041503450325041003650395035003900345038506.1.2.1 数据准备激活数据管理窗口，定义变量名：组变量为group（运动员=1，大学生=2），身高为x，肺活量为y，按顺序输入相应数值，建立数据库，结果见图6.1。图6.1 原始数据的输入6.1.2.2 统计分析 激活 Statistics 菜单选ANOVA Models中的Simple Factorial.项，弹出Simple Factorial ANOV

5、A对话框（图6.2）。在变量列表中选变量y，点击钮使之进入Dependent框；选分组变量group，点击钮使之进入Factor(s)框中, 并点击Define Range.钮在弹出的Simple Factorial ANOVA:Define Range框中确定分组变量group的起止值（1,2）；选协变量x，点击钮使之进入Covariate(s)框中。图6.2 协方差分析对话框 点击Options.框，弹出Simple Factorial ANOVA:Options对话框。系统在协方差分析的方法（Method）上有三种选项： 1、Unique：同时评价所有的效应； 2、Hierarchica

6、l：除主效应外，逐一评价各因素的效应； 3、Experimental：评价因素干预之前的主效应。 本例选Unique方法，之后点击Continue钮返回Simple Factorial ANOVA对话框，再点击OK钮即可。6.1.2.3 结果解释 在结果输出窗口中可见如下统计数据： 先输出肺活量总均数和两组的肺活量均数，总均数为4033.25，运用员组均数为4399.00，大学生组为3667.50。 接着协方差分析表明，混杂因素X（身高）两组间是有差异的（F=10.679，P=0.002），控制其影响后，两组间肺活量的差别依然存在（F=9.220，P=0.004），故可以认为两组间肺活量的均数

7、在消除了身高因素的影响之后仍有差别，运动员的肺活量大于大学生，即体育锻炼会提高肺活量。 最后系统输出公共回归系数，= 36.002，该值可用于求修正均数： = - ( - ) 本例为= 4399.00 - 36.002（178.175 - 174.3325）= 4260.6623 = 3667.50 - 36.002（170.49 - 174.3325）= 3805.8377Y by GROUPTotal Population 4033.25 ( 40)GROUP 1 2 4399.00 3667.50 ( 20) ( 20)Y by GROUP with X UNIQUE sums of s

8、quares All effects entered simultaneously Sum of Mean SigSource of Variation Squares DF Square F of FCovariates 1630763 1 1630762.635 10.679 .002 X 1630763 1 1630762.635 10.679 .002Main Effects 1407847 1 1407847.095 9.220 .004 GROUP 1407847 1 1407847.095 9.220 .004Explained 6981685 2 3490842.568 22.

9、860 .000Residual 5649992 37 152702.496Total 12631678 39 323889.16740 cases were processed.0 cases (.0 pct) were missing.Covariate Raw Regression CoefficientX 36.002返回目录 返回全书目录第二节 General Factorial过程6.2.1 主要功能 调用此过程可对完全随机设计资料、配伍设计资料、析因设计资料、正交设计资料等等进行多因素方差分析或协方差分析。返回目录 返回全书目录6.2.2 实例操作 例6-2下表为三因素析因实验的

10、资料，请用方差分析说明不同基础液与不同血清种类对钩端螺旋体的培养计数的影响。基础液（A）血清种类（B）兔血清浓度（C）胎盘血清浓度（C）5858缓冲液64812461398909114418771671184583085344110305786696431002蒸馏水1763124113812421144718831896192692070984857493310241092742自来水580 10261026830178912151434165111261176128012126855465955666.2.2.1 数据准备 激活数据管理窗口，定义变量名：基础液为base，血清种类为sero

11、，血清浓度为pct，钩端螺旋体的培养计数为X，按顺序输入相应数值，建立数据库。6.2.2.2 统计分析 激活Statistics菜单选ANOVA Models中的General Factorial.项，弹出General Factorial ANOVA对话框（图6.3）。在对话框左侧的变量列表中选变量x，点击钮使之进入Dependent Variable框；选要控制的分组变量base、sero和pct，点钮使之进入Factor(s)框中，并分别点击Define Range钮，在弹出的General Factorial ANOVA:Define Range对话框中确定各变量的起止值，本例变量ba

12、se的起止值为1、3，变量sero的起止值为1、2，变量pct的起止值为1、2。之后点击OK钮即可。图6.3 析因方差分析对话框6.2.2.3 结果解释 在结果输出窗口中，系统显示48个观察值进入统计，三个因素按其各自水平共产生12种组合。 分析表明，模型总效应的F值为10.55，P值 0.001，说明三因素间存在有交互作用。单因素效应和交互效应导致的组间差别比较结果是： 单因素组间比较： A：基础液（BASE） F = 4.98，P = 0.012，说明三种培养基培养钩体的计数有差别； B：血清种类（SERO） F = 61.265，P 0.001，说明两种血清培养钩体的计数有差别； C：血

13、清浓度（PCT） F = 3.49，P = 0.070，说明两种血清浓度培养钩体的计数无差别。 两因素构成的一级交互作用： AB：基础液（BASE）血清种类（SERO） F = 5.16，P = 0.011，交互作用明显； BC：血清种类（SERO）血清浓度（PCT） F = 15.96，P 0.001，交互作用明显； AC：基础液（BASE）血清浓度（PCT） F = 0.78，P = 0.465，交互作用不明显。 三因素构成的二级交互作用： ABC：基础液（BASE）血清种类（SERO）血清浓度（PCT） F = 6.75，P = 0.003，交互作用明显。48 cases accepte

14、d. 0 cases rejected because of out-of-range factor values. 0 cases rejected because of missing data.12 non-empty cells. 1 design will be processed. - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Univariate Homogeneity of Variance Tests Variable . X Cochrans C(3,12) = .3400

15、4, P = .036 (approx.) Bartlett-Box F(11,897) = 1.69822, P = .069 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -* * * * * * A n a l y s i s o f V a r i a n c e - design 1 * * * * * * Tests of Significance for X using UNIQUE sums of squares Source of Variation SS DF MS F Sig

16、 of F WITHIN+RESIDUAL 2459233.75 36 68312.05 BASE 679967.38 2 339983.69 4.98 .012 PCT 238713.02 1 238713.02 3.49 .070 SERO 4184873.52 1 4184873.5 61.26 .000 BASE BY PCT 107005.54 2 53502.77 .78 .465 BASE BY SERO 705473.04 2 352736.52 5.16 .011 PCT BY SERO 1089922.69 1 1089922.7 15.96 .000 BASE BY PC

17、T BY SERO 922307.37 2 461153.69 6.75 .003 (Model) 7928262.56 11 720751.14 10.55 .000 (Total) 10387496.31 47 221010.56 R-Squared = .763 Adjusted R-Squared = .691返回目录 返回全书目录第三节 Multivarite过程6.3.1 主要功能 调用此过程可进行多元方差分析。此外，对于一元设计，如涉及混合模型的设计、分割设计（又称列区设计）、重复测量设计、嵌套设计、因子与协变量交互效应设计等，此过程均能适用。返回目录 返回全书目录6.3.2 实

18、例操作 例6-3甲地区为大城市，乙地区为县城，丙地区为农村。某地分别调查了上述三类地区8岁男生三项身体生长发育指标：身高、体重和胸围，数据见下表，问：三类地区之间男生三项身体生长发育指标的差异有无显著性？学生编号甲地区乙地区丙地区身高体重胸围身高体重胸围身高体重胸围123456789101112131415161718192021222324252627282930119.80121.70121.40124.40120.00117.00118.10118.80124.20124.90124.70123.00125.30124.20127.40128.20126.10128.70129.5012

19、6.90126.50128.20131.40130.80133.90130.40131.30130.20136.00141.0022.6021.5019.1021.8021.4020.1018.8022.0021.3024.0023.3022.5022.9019.5022.9022.3022.7023.5024.5025.5025.0026.1027.9026.8027.2024.4024.4023.0026.3031.9060.5055.5056.5060.5057.7057.0057.1061.7058.4060.8060.0060.0065.2053.8059.5060.0057.406

20、0.4051.0061.5063.9063.0063.1061.5065.8062.6059.5062.6060.0063.70125.10127.00125.70114.90124.90117.60124.20117.90120.40115.00126.20125.10114.90121.50114.00118.70120.60122.90119.60112.30121.30121.20120.20120.30120.00123.30122.10123.30109.90125.6023.0021.5023.4017.5023.5018.9020.8020.3020.0019.7021.202

21、2.1019.7022.0019.0019.1020.0018.5019.5020.0020.0021.2023.1021.0022.2020.1021.0021.5017.8023.3062.0059.0061.5052.5058.5057.0058.5061.0056.0056.5056.5058.5056.0057.0054.5054.5055.5056.0059.5058.0058.0059.0059.5059.5059.5056.5057.5061.0056.5060.50118.30121.30121.80124.20123.50123.00134.90123.70105.2011

22、2.20118.60112.00121.50124.50119.50122.50115.50122.50124.50125.00117.50127.30122.30121.30120.50116.00120.50114.50131.00122.5020.4020.0026.6022.1023.2022.9032.3022.7020.2020.8021.0023.2024.0021.5020.5023.0019.0022.5025.0025.5023.0022.5022.0021.0022.0019.0020.0019.0025.5024.5054.4054.3061.1058.6060.205

23、8.2064.8059.9054.5057.5057.6058.2060.3055.6055.5056.7054.2057.6057.9060.3059.0058.9058.2055.6055.1053.5054.4053.4058.3058.706.3.2.1 数据准备 激活数据管理窗口，定义变量名：地区为G，身高为X1，体重为X2，胸围为X3，按顺序输入相应数值，变量G的数值是：甲地区为1，乙地区为2，丙地区为3。6.3.2.2 统计分析 激活Statistics菜单选ANOVA Models中的Multivarite.项，弹出Multivarite ANOVA 对话框（图6.8）。首先指

24、定供分析用的变量x1、x2、x3，故在对话框左侧的变量列表中选变量x1、x2、x3，点击钮使之进入Dependent Variable框；然后选变量g（分组变量）点击钮使之进入Factor(s)框中，并点击Define Range钮，确定g的起始值和终止值。图6.4 多元方差分析对话框 点击Options.钮，弹出Multivarite ANOVA:Options对话框，选择需要计算的指标。在Factor(s)栏内选变量g，点击钮使之进入Display Means for框，要求计算平均值指标；在Matriced Within Cell栏内选Correlation、Covariance、SSC

25、P项，要求计算单元内的相关矩阵、方差协方差矩阵和离均差平方和交叉乘积矩阵；在Error Matrices栏内也选上述三项，要求计算误差的相关矩阵、方差协方差矩阵和离均差平方和交叉乘积矩阵；在Diagnostics栏内选Homogeneity test项，要求作变量的方差齐性检验。之后点击Continue钮返回Multivarite ANOVA对话框，最后点击OK钮即可。6.3.2.3 结果解释 在结果输出窗口中将看到如下分析结果： 系统首先显示共90个观察值进入统计分析，因分组变量g为三个地区，故分析的单元数为3。然后输出3个应变量（x1、x2、x3）的方差齐性检验结果，分别输出了Cochra

26、n C检验值及其显著性水平P值、Bartlett-Box F检验值及其显著性水平P值。其中 身高：C = 0.39825，P = 0.540；F = 1.01272，P = 0.363； 体重：C = 0.43787，P = 0.227；F = 4.48624, P = 0.011； 胸围：C = 0.47239, P = 0.089；F = 2.06585, P = 0.127； 可见3项指标的方差基本整齐（P值均大于0.05）。90 cases accepted. 0 cases rejected because of out-of-range factor values. 0 cases

27、 rejected because of missing data. 3 non-empty cells. 1 design will be processed. CELL NUMBER 1 2 3 Variable G 1 2 3 Univariate Homogeneity of Variance Tests Variable . X1 Cochrans C(29,3) = .39825, P = .540 (approx.) Bartlett-Box F(2,17030) = 1.01272, P = .363 Variable . X2 Cochrans C(29,3) = .4378

28、7, P = .227 (approx.) Bartlett-Box F(2,17030) = 4.48624, P = .011 Variable . X3 Cochrans C(29,3) = .47239, P = .089 (approx.) Bartlett-Box F(2,17030) = 2.06585, P = .127 Cochran C检验和Bartlett-Box F检验对考查协方差矩阵的相等性比较方便，但还不够。于是系统接着分别输出了三类地区（即各个单元）各生长发育指标的离均差平方和交叉乘积矩阵和方差协方差矩阵。之后作Box M检验，Box M检验提供矩阵一致性的多元测试，本例Boxs M = 36.93910，在基于方差分析的显著性检验中F = 2.92393；在基于2的显著性检验中2 = 35.09922, 两者P 0