方差分析怎么操作

方差分析

常用于方差分析的主要过程有ANOVA和GLM（广义线性模

型），对于平衡数据资料（各水平下等重复，数据没有丢失），一般用ANOVA过程，对于非平衡数据，应采用GLM过程. 1、

ANOVA过程格式及使用说明

过程格式：

PROC ANOVA [选项]； CLASS 处理因素； MODEL 因变量=效应表[/选择项]； MEANS 效应表 [/选择项]； 过程说明：

◆ PROC ANOVA 语句的选项主要有：

DATA=数据集名 指明要分析的SAS数据集，缺省时SAS将使用最近建立的数据集.

OUTSTAT=输出数据集 ◆CLASS

指定分析计算结果输出的数据集名.

语句指明分类变量，是ANOVA过程的必需语句，并且

必须出现在MODEL语句之前. 分类变量可以为数值型或字符型，分类变量的个数表示方差分析的因素个数.

◆MODEL

语句定义分析所用的效应模型，即方差分析的因变量

和效应变量. 在方差分析过程中，关键在于定义线性数学模型，常用的模型定义语句有：

MODEL y=a 单因素一元方差分析

一、单因素方差分析

1．完全窗口介绍

单因素方差分析的完全窗口管理通过Analyze菜单中的Compare Means由One-Way ANOVA菜单项调用。 （1）主对话框

按Analyze → Compared Means → One-Way Anova的顺序单击。就可以打开“单因素方差分析”主对话框，如图1所示。

图1 “单因素方差分析”对话框

（2）因变量框

在主对话框中可以看到因变量框（Dependent List），该框中列出主要分析的所有因变量。要从左源变量框中选取变量进入该框，只需选中所要选取的变量，然后按向右的箭头即可。可以有多个因变量。 （3）因素框

在主对话框中可以看到因素框（Factor），该框中列出了因素。要从左边源变量框中选取变量进入该框，只需选中所要选取的变量，然后按向右的箭头即可。因素同样也是分组变量，必须满足只取有限个水平的条件。 （4）Contrast对话框

在主对话框中单击【Contrast】键，即可打开“Contrast”对话框，如图2所示。在该框中指定一种要用t检验来检验的priori对比，可以进行均值的多项式比较。

图2 多项式比较对话框

该框中各项意义如下： ① Polynomial复选框 选

北京大学医学部

第五章多组数值变量比较王洪源

北京大学医学部

假设检验

两组数值变量比较

正态性、等方差假设

t-检验 正态性假设成立、不等方差 调整的t-检验 正态性、等方差假设不成立 Wilcoxon秩和检验 在正态性、等方差假设成立时t-检验的效 率是好的。

北京大学医学部

假设检验

多组数值变量比较

正态性、等方差假设 方差分析 正态性、等方差假设不成立 Kruskal-Wallis秩和检验

北京大学医学部

为研究铅对儿童神经行为的影 响，研究者在某铅矿区对儿童的血铅水平及 神经行为评价指标手指敲击测验进行了测定， 第一年和第二年儿童的血铅水平均大于等于 40 mg/dl的17名，为暴露组(group=2),第一 年儿童的血铅水平均大于等于40mg/dl、第 二年儿童的血铅水平小于40mg/dl的15名， 为既往暴露组(group=3),第一年和第二年儿 童的血铅水平均小于40mg/dl的15名，为对 照组(group=1),神经行为评价指标为第二年 的手指敲击测验得分。

例9.1

北京大学医学部

表 9.1 某铅矿区儿童不同铅表露水平的手指敲击测验结果 对照组 手指敲击 No 1 2 3 4 16 17 18 19 group 1 1

方差分析

5.1.1评价不同行业的服务质量，消费者协会分别在零售业、旅游业、航空公司、家电制造业抽取了不同的企业作为样本，其中零售业7家，旅游业6家，航空公司5家，家电制造业5家，然后统计出近期消费者对这23家企业的投诉次数，试分析4个行业之间的服务质量是否存在显著差异？（基本数据见5-1.sav，资料来源：赖国毅等编著，SPSS17.0常用功能与应用，电子工业出版社）。

5.1.2.某企业有4条生产线生产同一中型号的产品，对每条生产线观测其一周的日产量，要求判断不同生产线的日产量是否有显著的差异(基本数据见5-2.sav)。

5.2.1.某商家有商品销售的数据资料，分析销售额是否受到促销方式和售后服务的影响。用变量“促销”对促销方式进行区分，取值为0表示无促销，取值为1表示被动促销，取值为2表示主动促销。变量“售后”对所采取的售后服务进行刻画，取值为0表示没有售后服务，取值为l表示有售后服务(基本数据见5-4.sav，资料来源：徐秋艳等，SPSS统计分析方法与应用实验教程，中国水利水电出版社，2011)。

5.3.1. 政府为了帮助年轻人提高工作技能，进行了一系列有针对性的就业能力和工作技能培训项目，为检验培训工作的成效，对1000

方差分析

5.1.1评价不同行业的服务质量，消费者协会分别在零售业、旅游业、航空公司、家电制造业抽取了不同的企业作为样本，其中零售业7家，旅游业6家，航空公司5家，家电制造业5家，然后统计出近期消费者对这23家企业的投诉次数，试分析4个行业之间的服务质量是否存在显著差异？（基本数据见5-1.sav，资料来源：赖国毅等编著，SPSS17.0常用功能与应用，电子工业出版社）。

5.1.2.某企业有4条生产线生产同一中型号的产品，对每条生产线观测其一周的日产量，要求判断不同生产线的日产量是否有显著的差异(基本数据见5-2.sav)。

5.2.1.某商家有商品销售的数据资料，分析销售额是否受到促销方式和售后服务的影响。用变量“促销”对促销方式进行区分，取值为0表示无促销，取值为1表示被动促销，取值为2表示主动促销。变量“售后”对所采取的售后服务进行刻画，取值为0表示没有售后服务，取值为l表示有售后服务(基本数据见5-4.sav，资料来源：徐秋艳等，SPSS统计分析方法与应用实验教程，中国水利水电出版社，2011)。

5.3.1. 政府为了帮助年轻人提高工作技能，进行了一系列有针对性的就业能力和工作技能培训项目，为检验培训工作的成效，对1000

1、方差齐性检验

由于方差分析的前提是各水平下的总体服从正态分布并且方差相等，因此有必要对方差齐性进行检验，即对控制变量不同水平下各观测变量不同总体方差是否相等进行分析。

SPSS单因素方差分析中，方差齐性检验采用了方差同质性（Homogeneity of Variance）的检验方法，其零假设是各水平下观测变量总体方差无显著性差异，实现思路同SPSS两独立样本t检验中的方差齐性检验。

2、多重比较检验

上面的基本分析可以判断控制变量是否对观测变量产生了显著影响。如果控制变量确实对观测变量产生了显著影响，进一步还应确定，控制变量的不同水平对观测变量的影响程度如何，其中哪个水平的作用明显大于其它水平，哪些水平的作用是不显著的。例如已经确定不同施肥量会对农作物的产量产生显著影响，便希望进一步了解究竟是10公斤、20公斤还是30公斤施肥量最有利于提高产量，哪种施肥量对农作物产量没有显著影响。掌握了这些信息，我们就能够制定合理的施肥方案。

多重比较检验就是分别对每个水平下的观测变量均值进行逐对比较，判断两均值之间是否存在显著差异。其零假设是相应组的均值之间无显著差异。

第五节 方差分析的SPSS操作

一、完全随机设计的单因素方差分析 1．数据

采用本章第二节所用的例1中的数据，在数据中定义一个group变量来表示五个不同的组，变量math表示学生的数学成绩。数据输入格式如图6－3（为了节省空间，只显示部分数据的输入）：

图 6-3 单因素方差分析数据输入 将上述数据文件保存为“6-6-1.sav”。 2．理论分析

要比较不同组学生成绩平均值之间是否存在显著性差异，从上面数据来看，总共分了5个组，也就是说要解决比较多个组（两组以上）的平均数是否有显著的问题。从要分析的数据来看，不同组学生成绩之间可看作相互独立，学生的成绩可以假设从总体上服从正态分布，在各组方差满足齐性的条件下，可以用单因素的方差分析来解决这一问题。单因素方差分析不仅可以检验多组均值之间是否存在差异，同时还可进一步采取多种方法进行多重比较，发现存在差异的究竟是哪些均值。

3．单因素方差分析过程 （1）主效应的检验

假如我们现在想检验五组被试的数学成绩（math）的均值差异是否显著性，可依下列操作进行。 ①单击主菜单Analyze/Compare Means/One-Way Anova…，进入主对话框，请把math选入到因变量表列（Depend

1、某研究者观测大白鼠的肝重与体重之比，了解正氟醚对该指标的作用，同时考虑生理盐水和用戊巴比妥作为诱导药对正氟醚毒性作用有无影响，对不同性别大白鼠的作用有何不同？数据如下：

是否用 不用 不用 不用 不用 用 用 用 用 诱导药 生理盐水 生理盐水 戊巴比妥 戊巴比妥 生理盐水 生理盐水 戊巴比妥 戊巴比妥 性别 雌性 雄性 雌性 雄性 雌性 雄性 雌性 雄性 5.68 5.52 5.50 6.46 5.6 7.02 4.6 5.7 肝重与体重之比 5.26 5.00 5.87 6.13 5.42 6.3 4.64 6.02

5.83 5.38 6.20 5.21 5.7 5.9 5.44 5.48

2、在1990 年秋对“亚运会期间收看电视的时间”调查结果如下表所示。问：收看电视的时间比平日减少了（第一组）、与平日无增减（第二组）、比平日增加了（第三组）的三组居民在“对亚运会的总态度得分”上有没有显著的差异？

3、研究5种类型的军装在两种环境、两种活动状态下的散热效果，将100名受试者随机等分20组，观察指标为受试者的主观热感觉，结果如下：

军装类型 活动环境 活动状态 a1 0.3 0.1 a2 0.25 -0.25 1

方差分析举例 一、什么是方差分析

例1：某饮料生产企业研制出一种新型饮料。饮料的颜色共有四种，分别为橘黄色、粉色、绿色和无色透明。这四种饮料的营养含量、味道、价格、包装等可能影响销售量的因素全部相同，先从地理位置相似、经营规模相仿的五家超级市场上收集了前一期该种饮料的销售量情况，见表10-1。

表10-1 该饮料在五家超市的销售情况 单位：箱

超市 1 2 3 4 5 合计 无色 26.5 28.7 25.1 29.1 27.2 136.6 粉色 31.2 28.3 30.8 27.9 29.6 147.8 橘黄色 27.9 25.1 28.5 24.2 26.5 132.2 绿色 30.8 29.6 32.4 31.7 32.8 157.3 问饮料的颜色是否对销售量产生影响。

解：从表10－1中看到，20个数据各不相同，其原因可能有两个方面： 一是销售地点不同的影响。

即使是相同颜色的饮料，在不同超市的销售量也是不同的。但是，由于这五个超市地理位置相似、经营规模相仿，因此，可以把不同地点产品销售量的差异看成是随机因素的影响。

二是饮料颜色不同的影响。

即使在同一个超市里，不同颜色的饮料的销售量也是不同

方差分析

（一）填空题

1、将在实验中或在抽样时发生变化的“量”称为 ；

2、方差分析的目的就是分析 对实验或抽样的结果有显著影响；

3、设因素A有r 个水平，且每一水平下样本容量为n 情况下：组间方差平方和SSA的

自由度为 ，组内方差平方和SSE的自由度为 。

（二）单项选择题

1、利用 “方差分析表”进行方差分析时，该表不包括的项目有（ ） A、方差来源， B、离差平方和及其分解，

C、各离差平方和的自由度， D、原假设的统计判断， 2、以下对方差分析叙述不正确的是（ ）

A、 方差分析可以对若干均值是否相等同时进行检验 B、 进行方差分析要求各水平下的样本容量相同 C、 离差平方和能分解为组内方差与组间方差的和 D、 方差分析方法在社会科学领域也大有用武之地

3、双因素方差分析有两种类型：一个是有交互作用的，一个是无交互作用的。区别的关键是看这对因子（ ）

A、是否独立 B、是否都服从正态分布 C、是否因子的水平相同