用SPSS进行单因素方差分析和多重比较

1、SPSS 单因素方差分析单因素方差分析单因素方差分析 也称作一维方差分析。它检验由单一因素影响的一个 或几个相互独立的因变量由因素各水平分组的均值之间的差异是否具有统计 意义。还可以对该因素的假设干水平分组中哪一组与其他各组均值间具有显著性差异进行分析，即进行均值的多重比拟。On e-Way ANOV过程要求因变量属于 正态分布总体。如果因变量的分布明显的是非正态，不能使用该过程，而应该使 用非参数分析过程。如果几个因变量之间彼此不独立，应该用Repeated Measure过程。例子调查不同水稻品种百丛中稻纵卷叶螟幼虫的数量，数据如表1-1所示。表1-1不同水稻品种百丛中稻纵卷叶螟幼虫数重复

2、水稻品种12345141333837312393735393434035353834数据保存在“文件中，变量格式如图1-1。图1-1分析水稻品种对稻纵卷叶螟幼虫抗虫性是否存在显著性差异。2启动分析过程点击主菜单“ Analyze项，在下拉菜单中点击“ Compare MeanS项, 在右拉式菜单中点击“ One-Way ANOVA项，系统翻开单因素方差分析设置窗口如图1-2。图1-2 单因素方差分析窗口3设置分析变量因变量：选择一个或多个因子变量进入 “Dependent List 框中。本例 选择“幼虫。因素变量：选择一个因素变量进入“Factor 框中。本例选择“品种。4设置多项式比拟单击

3、“ Contrasts 按钮，将翻开如图1-3所示的对话框。该对话框用于 设置均值的多项式比拟。图 1-3“ Contrasts 对话框定义多项式的步骤为：均值的多项式比拟是包括两个或更多个均值的比拟。例如图 1-3中显示 的是要求计算“x mean1-1x mean2的值，检验的假设H0第一组均 值的1.1 倍与第二组的均值相等。单因素方差分析的“ On e-Way ANOVA过程允许进行高 达5次的均值多项式比拟。多项式的系数需要由读者自己根据研究的需要输入。具体的操作步骤如下： 选中“ Polynomial 复选项，该操作激活其右面的“ Degree参数框。 单击Degree参数框右面的

4、向下箭头展开阶次菜单，可以选择“ Linear 线 性、“Quadratic 二次、“ Cubic 三次、“ 4th 四次、“ 5th 五次多项式。 为多项式指定各组均值的系数。方法是在“Coefficients 框中输入一个系数，单击Add按钮，“ Coefficients 框中的系数进入下面 的方框中。依次输 入各组均值的系数，在方形显示框中形成一列数值。因素变量分为几组，输入几 个系数，多出的无意义。如果多项式中只包括第一组与第四组的均值的系数，必须把第二个、第三个系数输入为0值。如果只包括第一组与第二组的均值，那么 只需要输入前两个系数，第三、四个系数可以不输入。可以同时建立多个多项式

5、。一个多项式的一组系数输入结束，激话“Next 按钮，单击该按钮后“ Coefficients 框中清空，准备接受下一组系数数 据。如果认为输入的几组系数中有错误，可以分别单击“ Previous 或“Next按钮 前后翻找出错的一组数据。单击出错的系数，该系数显示在编辑框中，可以在此进行修改，修改后单击“ Changd'按钮在系数显示框中出现正确的系数值。当 在系数显示框中选中一个系数时，同时激话“Remove按钮，单击该按钮将选中的系数去除。 单击“ Previous 或“ Next按钮显示输入的各组系数检查无误后，按“ Continue 按钮确认输入的系数并返回到主对话框。要取消

6、刚刚的输入，单击“Cancel 按钮；需要查看系统的帮助信息，单击“ Help按钮。本例子不做多项式比拟的选择，选择缺省值。5设置多重比拟在主对话框里单击“ Post Hoc按钮，将翻开如图5-4所示的多重比拟 对话框。该对话框用于设置多重比拟和配比照较。方差分析一旦确定各组均值间 存在差异显著，多重比拟检测可以求出均值相等的组；配比照较可找出和其它组均值有差异的组，并输出显著性水平为0.95的均值比拟矩阵，在矩阵中用星号表示有差异的组。*1 *ada>ftce p a 厂 QariilrilHil r sidik r sirtB r ejegw r r FH>GWQr蛊k广r T

7、hlcy-e P ；Efraw r UMhlERB-B QTIGHiri祕MfaHtf DuricanI丹.屮VP中I IF DumrcHC<iIi-Tj itl L fHrIIE.i*b4 "Vhri ARcei-kiMn益dr r*mih«iwlB T3 厂 Dwin-it4! I F G«<i>ri： Hciwll L Dwif'HrW'a CCancElioiCnrwtiniie对话框图 1-4“Post Hoc Multiple Comparis ons多重比拟的选择项: 方差具有齐次性时(Equal Variances

8、 Assumed)，该矩形框中有如下方法供选 择：LSD (Least-significant differenee)最小显著差数法，用 t检验完成各组均值间的配比照较。对多重比拟误差率不进行调整。Bonferroni(LSDMOD)用t检验完成各组间均值的配比照较，但通过设置每个检验的误差率来控制整个误差率。Sidak 计算t统计量进行多重配比照较。可以调整显著性水平，比Bofferro ni 方法的界限要小。Scheffe对所有可能的组合进行同步进入的配比照较。这些选择项可以同时选择假设干个。以便比拟各种均值比拟方法的结果。R-E-G-WF (Ryan-Einot-Gabriel-Wels

9、ch F)用 F 检验进行多重比较检验。R-E-G-WQ (Ryan-Einot-Gabriel-Welsch range test)正态分布范围进行多重配比照较。S-N-K (Student-Newmnan-Keuls)用 Student Range 分布进行所有各组均值间的配比照较。如果各组样本含量相等或者选择了“Harm onic average of all groups即用所有各组样本含量的调和平均数进行样本量估计时还用逐步过程进行齐次子集(差异较小的子集)的均值配比照较。在该比拟过程中， 各组均值从大到小按顺序排列，最先比拟最末端的差异。Tukey (Tukey's ，ho

10、nestly signicant differenee)用 Student-Range统计量进行所有组间均值的配比照较，用所有配比照较误差率作为实验误差率。Tukey's-b 用“ stn de nt Ra nge 分布进行组间均值的配比照较。其精确值为前两种检验相应值的平均值。Duncan (Duncan's multiple range test)新复极差法SSR，指定一系列的“ Rangd'值，逐步进行计算比拟得出结论。Hochberg's GT2 用正态最大系数进行多重比拟。Gabriel用正态标准系数进行配比照较，在单元数较大时，这种方法较自由。Wal

11、ler-Dunca 用t统计量进行多重比拟检验,使用贝叶斯逼近。Dunn ett指定此选择项，进行各组与对照组的均值比拟。默认的对 照组是最后一组。选择了该项就激活下面的“ Co ntrol Category 参数框。展开 下拉列表，可以重新选择对照组。“ Test 框中列出了三种区间分别为：“2-sides双边检验；“ <Control 左边检验 “ >Conbo1 “右边检验。 方差不具有齐次性时(Equal Varance not assumed)，检验各均数间是否有差 异的方祛有四种可供选择：Tamha ne's T2, t检验进行配比照较。Dunnett'

12、s T3，采用基于学生氏最大模的成比照较法。Games-Howell, Games-Howell 比拟，该方法较灵活。Dunnett's C，采用基于学生氏极值的成比照较法。 Significanee选择项，各种检验的显著性概率临界值，默认值为0.05，可由用户重新设定。本例选择“ LSD和“Dun can比拟，检验的显著性概率临界值 0.05。6)设置输出统计量单击“ Options 按钮，翻开“ Options 对话框，如图1-5所示。选择 要求输出的统计量。并按要求的方式显示这些统计量。 在该对话框中还可以选择 对缺失值的处理要求。各组选择项的含义如下：One-lay AHOVA

13、： Options图1-5输出统计量的设置“Statistics 栏中选择输出统计量：Descriptive ,要求输出描述统计量。选择此项输出观测量数目、 均值、标准差、标准误、最小值、最大值、各组中每个因变量的95%置信区间。Fixed and ran dom effects,固定和随机描述统计量Homoge neity-of-varia nee ，要求进行方差齐次性检验，并输出检验 结果。用“ Levene lest 检验，即计算每个观测量与其组均值之差，然后对这 些差值进行一维方差分析。Brow n-Forsythe布朗检验Welch,韦尔奇检验Means plot ,即均数分布图，根

14、据各组均数描绘出因变量的分布情 况。“ Missing Values 栏中，选择缺失值处理方法。Exclude cases an alysis by an alysis选项，被选择参与分析的变量含缺失值的观测量，从分析中剔除。Exclude cases listwise选项，对含有缺失值的观测量，从所有分析中剔除。以上选择项选择完成后，按“ Continue 按钮确认选择并返回上一级对 话框；单击“ Can cel按钮作废本次选择；单击“ Help 按钮，显示有关的帮助 信息。本例子选择要求输出描述统计量和进行方差齐次性检验，缺失值处理方 法选系统缺省设置。6提交执行设置完成后，在单因素方差分

15、析窗口框中点击“ 0K按钮，SPSS就会根 据设置进行运算，并将结算结果输出到 SPSS结果输出窗口中。7结果与分析输出结果：忖Meanstd.D&vlatlnnstd.Error95% Confidence Interl forMeanMinimumMaKimurriLower BouncUpper Bound12G45Total333331540.Q036.0036.0030.0033.0033.451,002.001 731.001732.92朋1.151.00.581.0D.7337.523103317035.52 2&7034.8442,4639 9740.30404

16、B37.3037.96393335373131413738393441S5-2 描述痂十基Destnatives百梵幼中表5-2描述统计量，给出了水稻品种分组的样本含量N、平均数Mean标准差Std.Deviation 、标准误Std.Error、95%勺置信区间、最小值和最大值。表£刁方歪齐出悝橙验 Test ofHomogerteity ofVariances百梵幼虫Levene Statisticdfldf?Sig.75041C.560表5-3为方差齐次性检验结果，从显著性慨率看，p>0.05，说明各组的方差在a=0.05水平上没有显著性差异，即方差具有齐次性。这个结论在

17、选择多 重比拟方法时作为一个条件。表5讯方着分折表ABOVA百理勃虫Sruii of SquaresdfIII can SquareFSig.Between GroupsTithi 舟 Gro叩 mTotal37.eoo £4, 000 111,6004101421, 90D2.4009 12S.002表5-4方差分析表：第1栏是方差来源，包括组间变差“ Between Grou ps；组内变差“ Within Groups 和总变差“ Total 。第2栏是离差平方和“ S um of Squares ，组间离差平方和87.600，组内离差平方和为24.000，总离差 平方和为11

18、1.600，是组间离差平方和与组内离差平方和相加之和。第3栏是自由度df,组间自由度为4,组内自由度为10;总自由度为14。第4栏是均方 “ Mea n Square，是第2栏与第3栏之比；组间均方为21.900，组内均方为2.400。第5栏是F值9.125组间均方与组内均方之比。第 6栏：F值对应的概率值，针对假设H0组间均值无显著性差异即5种品种虫数的平均值无显著性差异)。计算的F值9.125，对应的概率值为0.002。表多重比拟表 Multiple Compari soneDependent Variable：百直幼虫Mean95M Con£i denct Int erralD

19、ifferenceStd. Err orSic.Upper(I)寫料(I)吴和Li* er EfliutdCL-BoundLSD12s.ao'1.26.0032. 187.9234 (JO"i 26.0101 135. &242 oa1 2E.145-.824.9251 26.0004. 139 SE2 1-s.oo'1.2.003-?.92-2 ie3-1.001 26.w-3 32has4-3.00'1 26.039-5 BZ1852 001 26-324. &231-4.00i1.26.010-5-i ia2I 0D1.26.448-1

20、323. 321-2.001 26.145-4.82.fl253.1 26439.185. S241-2 oa1.26.145-4 B2.fi223, (X)'1.空.039.185. &23s.oai 26.145l 82az55.00'1.2G.0032 187.5251-7.00'1 £6.000-9 3E-162-2. on1.2E.145-4 32.323-3.00'i 26-5.32:IS4-5.0Q11.26.003-7 32£0Th®di ifi s si gxii ficant at th*.05 l

21、71;v«l表5-5 LSD法进行多重比拟表，从表5-4结论该例子的方差具有其 次性，因此LSD方法适用。第1栏的第1列“ i品种为比拟基准品种，第2 列“ j品种是比拟品种。第2栏是比拟基准品种平均数减去比拟品种平均数的差值Mean Differenee，均值之间具有0.05水平可图5-4对话框里设 置上有显著性差异，在平均数差值上用“ *号说明。第3栏是差值的标准误。第4栏是差值检验的显著性水平。第5栏是差值的95汕信范围的下限和上限Hunt m 祛进行比較叢HcwiogeiLeaiiE Substts百融虫品种可SubsetQpha = .051234DiiZLe an'S333 002335.0035.003336.003 & 004339. 003S. 0Qt340 00Si g-.145448145.145Jkan5 for groups in homg也电qwwaie displ ayed l 均数相似的显示在同一列Usas Jiarm cm t H电血£血口1_电Size - 3. OOC <计算均數用的溜调和平均数差值表5-6 是多重比拟的Duncan法进行比拟的结果。第1栏为品种，按均 数由小到大排列。第2栏列出