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

SPSS——单因素方差分析单因素方差分析单因素方差分析 也称作一维方差分析它检验由单一因素影响的一个(或几个相互独立的)因变量由因素各水平分组的均值之间的差异是否具有统计意义还可以对该因素的若干水 平分组中哪一组与其他各组均值间具有显著性差异进行分析，即进行均值的多重比较One-Way ANOVA 过程要求因变量属于正态分布总体如果因变量的分布明显的是非正态，不能使用该过程，而应该使用非参数分析过程如果几个因变量之间彼此不独立，应该用 Repeated Measure过程[例子]调查不同水稻品种百丛中稻纵卷叶螟幼虫的数量，数据如表 1-1 所示表 1-1 不同水稻品种百丛中稻纵卷叶螟幼虫数水 稻 品 种重复1 2 3 4 51 41 33 38 37 312 39 37 35 39 343 40 35 35 38 34数据保存在“data1.sav”文件中，变量格式如图 1-1图 1-1分析水稻品种对稻纵卷叶螟幼虫抗虫性是否存在显著性差异2）启动分析过程点击主菜单“Analyze”项，在下拉菜单中点击“Compare Means”项，在右拉式菜单中点击“0ne-Way ANOVA”项，系统 打开单因素方差分析设置窗口如图 1-2。

图 1-2 单因素方差分析窗口3）设置分析变量因变量: 选择一个或多个因子变量进入“Dependent List”框中本例选择“幼虫”因素变量: 选择一个因素变量进入“Factor”框中本例选择“品种”4）设置多项式比较单击“Contrasts”按钮，将打开如图 1-3 所示的对话框该对话框用于设置均值的多项式比较图 1-3 “Contrasts”对话框定义多项式的步骤为：均值的多项式比较是包括两个或更多个均值的比较例如图 1-3 中显示的是要求计算“1.1×mean1-1×mean2”的值，检验的假设 H0：第一组均 值的 1.1 倍与第二组的均值相等单因素方差分析的“0ne-Way ANOVA”过程允许进行高达 5 次的均值多项式比较多项式的系数需要由读者自己根据研究的需要输入具体的操作步骤如下：① 选中“Polynomial”复选项，该操作激活其右面的“Degree”参数框② 单击 Degree 参数框右面的向下箭头展开阶次菜单，可以选择“Linear”线性、“Quadratic”二次、“Cubic”三次、“4th”四次、“5th”五次多项式③ 为多项式指定各组均值的系数方法是在“Coefficients”框中输入一个系数，单击 Add 按钮，“Coefficients”框中的系数进入下面 的方框中。

依次输入各组均值的系数，在方形显示框中形成—列数值因素变量分为几组，输入几个系数，多出的无意义如果多项式中只包括第一组与第四组的均 值的系数，必须把第二个、第三个系数输入为 0 值如果只包括第一组与第二组的均值，则只需要输入前两个系数，第三、四个系数可以不输入 可以同时建立多个多项式一个多项式的一组系数输入结束，激话“Next”按钮，单击该按钮后“Coefficients”框中清空，准备接受下一组系数数据如果认为输入的几组系数中有错误，可以分别单击“Previous”或“Next”按钮前后翻找出错的一组数据单击出错的系数，该系数显示在编辑框中， 可以在此进行修改，修改后单击“Change”按钮在系数显示框中出现正确的系数值当在系数显示框中选中一个系数时，同时激话“Remove”按钮，单 击该按钮将选中的系数清除④单击“Previous”或“Next”按钮显示输入的各组系数检查无误后，按“Continue”按钮确认输入的系数并返回到主对话框要取消刚刚的输入，单击“Cancel”按钮；需要查看系统的帮助信息，单击“Help”按钮 本例子不做多项式比较的选择,选择缺省值5）设置多重比较在主对话框里单击“Post Hoc”按钮，将打开如图 5-4 所示的多重比较对话框。

该对话框用于设置多重比较和配对比较方差分析一旦确定各组均值间存在差异显著，多重比较检测可以 求出均值相等的组；配对比较可找出和其它组均值有差异的组，并输出显著性水平为 0.95 的均值比较矩阵，在矩阵中用星号表示有差异的组 图 1-4 “Post Hoc Multiple Comparisons”对话框(1)多重比较的选择项：①方差具有齐次性时(Equal Variances Assumed)，该矩形框中有如下方法供选择：LSD (Least-significant difference) 最小显著差数法，用 t 检验完成各组均值间的配对比较对多重比较误差率不进行调整Bonferroni (LSDMOD) 用 t 检验完成各组间均值的配对比较，但通过设置每个检验的误差率来控制整个误差率Sidak 计算 t 统计量进行多重配对比较可以调整显著性水平，比 Bofferroni 方法的界限要小Scheffe 对所有可能的组合进行同步进入的配对比较这些选择项可以同时选择若干个以便比较各种均值比较方法的结果R-E-G-WF (Ryan-Einot-Gabriel-Welsch F) 用 F 检验进行多重比较检验。

R-E-G-WQ (Ryan-Einot-Gabriel-Welsch range test) 正态分布范围进行多重配对比较S-N-K (Student-Newmnan-Keuls) 用 Student Range 分布进行所有各组均值间的配对比较如果各组样本含量相等或者选择了“Harmonic average of all groups”即用所有各组样本含量的调和平均数进行样本量估计时还用逐步过程进行齐次子集(差异较小的子集)的均值配对比较在该比较过程中，各组均值从大到小按顺序排列，最先比较最末端的差异Tukey (Tukey's，honestly signicant difference) 用 Student-Range统计量进行所有组间均值的配对比较，用所有配对比较误差率作为实验误差率Tukey's-b 用“stndent Range”分布进行组间均值的配对比较其精确值为前两种检验相应值的平均值 Duncan (Duncan's multiple range test) 新复极差法（SSR），指定一系列的“Range”值，逐步进行计算比较得出结论Hochberg's GT2 用正态最大系数进行多重比较。

Gabriel 用正态标准系数进行配对比较，在单元数较大时，这种方法较自由Waller-Dunca 用 t 统计量进行多重比较检验,使用贝叶斯逼近Dunnett 指定此选择项，进行各组与对照组的均值比较默认的对照组是最后一组选择了该项就激活下面的“Control Category”参数框展开下拉列表，可以重新选择对照组Test”框中列出了三种区间分别为： “2-sides” 双边检验 ; “Conbo1”“右边检验②方差不具有齐次性时(Equal Varance not assumed)，检验各均数间是否有差异的方祛有四种可供选择：Tamhane's T2, t 检验进行配对比较Dunnett's T3，采用基于学生氏最大模的成对比较法Games-Howell，Games-Howell 比较，该方法较灵活Dunnett's C，采用基于学生氏极值的成对比较法③ Significance 选择项，各种检验的显著性概率临界值，默认值为 0.05，可由用户重新设定本例选择“LSD”和“Duncan”比较，检验的显著性概率临界值 0.056) 设置输出统计量单击“Options”按钮，打开“Options”对话框，如图 1-5 所示。

选择要求输出的统计量并按要求的方式显示这些统计量在该对话框中还可以选择对缺失值的处理要求各组选择项的含义如下：图 1-5 输出统计量的设置“Statistics”栏中选择输出统计量:Descriptive，要求输出描述统计量选择此项输出观测量数目、均值、标准差、标准误、最小值、最大值、各组中每个因变量的 95％置信区间Fixed and random effects, 固定和随机描述统计量Homogeneity-of-variance，要求进行方差齐次性检验，并输出检验结果用“Levene lest ”检验，即计算每个观测量与其组均值之差，然后对这些差值进行一维方差分析Brown-Forsythe 布朗检验Welch,韦尔奇检验Means plot，即均数分布图，根据各组均数描绘出因变量的分布情况Missing Values”栏中，选择缺失值处理方法Exclude cases analysis by analysis 选项，被选择参与分析的变量含缺失值的观测量，从分析中剔除Exclude cases listwise 选项，对含有缺失值的观测量，从所有分析中剔除以上选择项选择完成后，按“Continue”按钮确认选择并返回上一级对话框；单击“Cancel”按钮作废本次选择；单击“Help”按钮，显示有关的帮助信息。

本例子选择要求输出描述统计量和进行方差齐次性检验，缺失值处理方法选系统缺省设置6）提交执行 设置完成后，在单因素方差分析窗口框中点击“OK”按钮，SPSS 就会根据设置进行运算，并将结算结果输出到 SPSS 结果输出窗口中7) 结果与分析输出结果：表 5-2 描述统计量，给出了水稻品种分组的样本含量 N、平均数 Mean、标准差 Std.Deviation、标准误 Std.Error、95%的置信区间、最小值和最大值表 5-3 为方差齐次性检验结果，从显著性慨率看，p>0.05，说明各组的方差在 a=0.05 水平上没有显著性差异，即方差具有齐次性这个结论在选择多重比较方法时作为一个条件表 5-4 方差分析表：第 1 栏是方差来源，包括组间变差“Between Groups”；组内变差“Within Groups”和总变差“Total”第 2 栏是离差平方和“Sum of Squares”，组间离差平方和 87.600，组内离差平方和为 24.000，总离差平方和为 111.600，是组间离差平方和与组内离差平方和相加之 和第 3 栏是自由度 df，组间自由度为 4，组内自由度为 10；总自由度为 14。

第 4 栏是均方“Mean Square”，是第 2 栏与第 3 栏之比；组间均方为 21.900，组内均方为 2.400第 5 栏是 F 值 9.125（组间均方与组内均方之比）第 6 栏：F 值对应的概率值，针对假设 H0：组间均值无显著性差异(即 5 种品种虫数的平均值无显著性差异)计算的 F 值 9.125，对应的概率值为 0.002表 5-5 LSD 法进行多重比较表，从表 5-4 结论已知该例子的方差具有其次性，因此 LSD 方法适用第 1 栏的第 1 列“[i]品种”为比较基准品种，第 2 列“[j] 品种”是比较品种第 2 栏是比较基准品种平均数减去比较品种平均数的差值（Mean Difference），均值之间具有 0.05 水平（可图 5-4 对话框里设置）上有显著性差异，在平均数差值上用“*”号表明第 3 栏是差值的标准误第 4 栏是差值检验的显著性水平第 5 栏是差值的 95%置信范围的下限和上限表 5-6 是多重比较的 Duncan 法进行比较的结果第 1 栏为品种，按均数由小到大排列第 2 栏列出计算均数用的样本数第 3 栏列出了在显著水平 0.05 上的比较结果，表的最后一行是均数方差齐次性检验慨率水平，p>0.05 说明各组方差具有齐次性。

多重比较比较表显著性差异差异的判读：在 同一列的平均数表示没有显著性差异，反之则具有显著性的差异例如，品种 3 横向看，平均数显。