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

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

2、-1不同水稻品种百丛中稻纵卷叶螟幼虫数重复水 稻 品 种12345141333837312393735393434035353834 数据保存在“文件中，变量格式如图1-1。图1-1 分析水稻品种对稻纵卷叶螟幼虫抗虫性是否存在显著性差异。c:iknowdocsharedatacur_workwwv.renrenxn.do?ss=10791&rt=1。2启动分析过程 点击主菜单“Analyze项，在下拉菜单中点击“pare Means项，在右拉式菜单中点击“0ne-Way ANOVA项，系统打开单因素方差分析设置窗口如图1-2。图1-2单因素方差分析窗口3设置分析变量 因变量:选择一个或多个因子

3、变量进入“Dependent List框中。本例选择“幼虫。 因素变量:选择一个因素变量进入“Factor框中。本例选择“品种。4设置多项式比拟 单击“Contrasts按钮，将打开如图1-3所示的对话框。该对话框用于设置均值的多项式比拟。图1-3 “Contrasts对话框定义多项式的步骤为： 均值的多项式比拟是包括两个或更多个均值的比拟。例如图1-3中显示的是要求计算“mean1-1mean2的值，检验的假设H0：第一组均 值的1.1倍与第二组的均值相等。单因素方差分析的“0ne-Way ANOVA过程允许进展高达5次的均值多项式比拟。多项式的系数需要由读者自己根据研究的需要输入。具体的操

4、作步骤如下：选中“Polynomial复选项，该操作激活其右面的“Degree参数框。单击Degree参数框右面的向下箭头展开阶次菜单，可以选择“Linear线性、“Quadratic二次、“Cubic三次、“4th四次、“5th五次多项式。为多项式指定各组均值的系数。方法是在“Coefficients框中输入一个系数，单击Add按钮，“Coefficients框中的系数进入下面 的方框中。依次输入各组均值的系数，在方形显示框中形成列数值。因素变量分为几组，输入几个系数，多出的无意义。如果多项式中只包括第一组与第四组的均 值的系数，必须把第二个、第三个系数输入为0值。如果只包括第一组与第二组的

5、均值，如此只需要输入前两个系数，第三、四个系数可以不输入。可以同时建立多个多项式。一个多项式的一组系数输入完毕，激话“Next按钮，单击该按钮后“Coefficients框中清空，准备承受下一组系数数据。如果认为输入的几组系数中有错误，可以分别单击“Previous或“Next按钮前后翻找出错的一组数据。单击出错的系数，该系数显示在编辑框中， 可以在此进展修改，修改后单击“Change按钮在系数显示框中出现正确的系数值。当在系数显示框中选中一个系数时，同时激话“Remove按钮，单 击该按钮将选中的系数去除。单击“Previous或“Next按钮显示输入的各组系数检查无误后，按“Continu

6、e按钮确认输入的系数并返回到主对话框。要取消刚刚的输入，单击“Cancel按钮；需要查看系统的帮助信息，单击“Help按钮。本例子不做多项式比拟的选择,选择缺省值。5设置多重比拟 在主对话框里单击“Post Hoc按钮，将打开如图5-4所示的多重比拟对话框。该对话框用于设置多重比拟和配比照拟。方差分析一旦确定各组均值间存在差异显著，多重比拟检测可以 求出均值相等的组；配比照拟可找出和其它组均值有差异的组，并输出显著性水平为0.95的均值比拟矩阵，在矩阵中用星号表示有差异的组。图1-4 “Post Hoc Multiple parisons对话框(1)多重比拟的选择项：方差具有齐次性时(Equa

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

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

9、icant difference)用Student-Range统计量进展所有组间均值的配比照拟，用所有配比照拟误差率作为实验误差率。Tukeys-b用“stndent Range分布进展组间均值的配比照拟。其准确值为前两种检验相应值的平均值。 Duncan (Duncans multiple range test)新复极差法SSR，指定一系列的“Range值，逐步进展计算比拟得出结论。Hochbergs GT2用正态最大系数进展多重比拟。Gabriel用正态标准系数进展配比照拟，在单元数较大时，这种方法较自由。Waller-Dunca用t统计量进展多重比拟检验,使用贝叶斯逼近。Dunnett指

10、定此选择项，进展各组与对照组的均值比拟。默认的对照组是最后一组。选择了该项就激活下面的“Control Category参数框。展开下拉列表，可以重新选择对照组。“Test框中列出了三种区间分别为：“2-sides双边检验;“Conbo1“右边检验。方差不具有齐次性时(Equal Varance not assumed)，检验各均数间是否有差异的方祛有四种可供选择：Tamhanes T2, t检验进展配比照拟。Dunnetts T3，采用基于学生氏最大模的成比照拟法。Games-Howell，Games-Howell比拟，该方法较灵活。Dunnetts C，采用基于学生氏极值的成比照拟法。 S

11、ignificance选择项，各种检验的显著性概率临界值，默认值为0.05，可由用户重新设定。本例选择“LSD和“Duncan比拟，检验的显著性概率临界值0.05。6)设置输出统计量 单击“Options按钮，打开“Options对话框，如图1-5所示。选择要求输出的统计量。并按要求的方式显示这些统计量。在该对话框中还可以选择对缺失值的处理要求。各组选择项的含义如下：图1-5输出统计量的设置“Statistics栏中选择输出统计量:Descriptive，要求输出描述统计量。选择此项输出观测量数目、均值、标准差、标准误、最小值、最大值、各组中每个因变量的95置信区间。Fixed and ran

12、dom effects,固定和随机描述统计量Homogeneity-of-variance，要求进展方差齐次性检验，并输出检验结果。用“Levene lest 检验，即计算每个观测量与其组均值之差，然后对这些差值进展一维方差分析。Brown-Forsythe布朗检验Welch,韦尔奇检验Means plot，即均数分布图，根据各组均数描绘出因变量的分布情况。“Missing Values栏中，选择缺失值处理方法。Exclude cases analysis by analysis选项，被选择参与分析的变量含缺失值的观测量，从分析中剔除。Exclude cases listwise选项，对含有缺

13、失值的观测量，从所有分析中剔除。 以上选择项选择完成后，按“Continue按钮确认选择并返回上一级对话框；单击“Cancel按钮作废本次选择；单击“Help按钮，显示有关的帮助信息。 本例子选择要求输出描述统计量和进展方差齐次性检验，缺失值处理方法选系统缺省设置。6提交执行 设置完成后，在单因素方差分析窗口框中点击“OK按钮，SPSS就会根据设置进展运算，并将结算结果输出到SPSS结果输出窗口中。7)结果与分析输出结果： 表5-2描述统计量，给出了水稻品种分组的样本含量N、平均数Mean、标准差Std.Deviation、标准误Std.Error、95%的置信区间、最小值和最大值。 表5-3

14、为方差齐次性检验结果，从显著性慨率看，p0.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栏

15、与第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栏是差值的标

16、准误。第4栏是差值检验的显著性水平。第5栏是差值的95%置信X围的下限和上限。 表5-6是多重比拟的Duncan法进展比拟的结果。第1栏为品种，按均数由小到大排列。第2栏列出计算均数用的样本数。第3栏列出了在显著水平0.05上的比拟结果，表的最后一行是均数方差齐次性检验慨率水平，p0.05说明各组方差具有齐次性。 多重比拟比拟表显著性差异差异的判读：在 同一列的平均数表示没有显著性差异，反之如此具有显著性的差异。例如，品种3横向看，平均数显示在第3列“2小列，与它同列显示的有品种2的平均数，说明 与品种2差异不显著0.05水平，再往右看，平均数显示在第3列“3小列，与它同列显示的有品种4的平均

17、数，说明与品种4差异不显著0.05水 平。如此品种3与品种5和品种1具有显著性的差异0.05水平。品种3和品种4都显示有平均数值。结果分析: 根据方差分析表输出的p值为0.002可以看出，无论临界值取0.05，还是取0.01，p值均小于临界值。因此否认Ho假设，水稻品种对稻纵卷叶螟幼虫抗 虫性有显著性意义，结论是稻纵卷叶螟幼虫数量的在不同品种间有明显的不同。根据该结论选择抗稻纵卷叶螟幼虫水稻品种，犯错误的概率几乎为0.008。只有在方差分析中F检验存在差异显著性时，才有比拟的统计意义。LSD法多重比拟明确：品种1与品种2、品种3和品种5之间存在显著性差异；品种2与品种1和品种4之间存在显著性差异；品种3与品种1和品种5之间存在显著性差异；品种4与品种2和品种5之间存在显著性差异；品种5与品种1、品种3和品种4之间存在显著性差异。Duncan法多重比拟明确：品种5与品种3、品种4和品种1之间存在显著性差异。品种2与品种4和品种1之间存在显著性差异；品种3与品种5和品种1之间存在显著性差异；品种4与品种5和品种2之间存在显著性差异；品种1与品种5、品种2和品种3之间存在显著性差异；两种方法比拟结果一致。