统计学实验SPSS和R软件应用与实例（第二版） 费宇第5章方差分析-SPSS

2023/1/5《统计学实验》第5章方差分析5-1统计学实验

—SPSS和R软件应用与实例

主编：费宇（第二版）2023/1/5《统计学实验》第5章方差分析5-2第5章

方差分析2023/1/5《统计学实验》第5章方差分析5-3一、实验目的熟悉分析菜单里比较均值子菜单中单因素方差分析；熟悉一般线性模型子菜单中单变量多因素方差分析模块的主要功能，

掌握运用SPSS进行方差分析的基本操作过程，并能读懂SPSS输出的结果。2023/1/5《统计学实验》第5章方差分析5-4二、实验环境系统软件Windows2000或WindowsXP或Windows7；统计软件SPSS20.0或更高版本。2023/1/5《统计学实验》第5章方差分析5-5三、实验内容单因素方差分析

多因素方差分析

2023/1/5《统计学实验》第5章方差分析5-6第5章

方差分析5.1单因素方差分析5.2双因素方差分析2023/1/5《统计学实验》第5章方差分析5-75.1单因素方差分析【例5.1】(数据文件为li5.1.sav)为研究咖啡因对人体的影响，进行如下的试验：咖啡因剂量取三个水平：0mg,100mg,200mg。挑选同一年龄，体质大致相同的30名健康的男大学生进行手指叩击训练。训练结束后，对每个水平随机的选定其中10个人，在服用咖啡因2小时后，请每个人做手指叩击，记录下每分钟叩击的次数。该试验进行双盲试验，即试验者和生物学家均不知道他们接受的是哪一种剂量的咖啡因，只有统计人员知道。试验数据如下表：

2023/1/5《统计学实验》第5章方差分析5-8表5.1咖啡因试验数据咖啡因剂量

重复数据0mg242245244248247248242244246242100mg248246245247248250247246243244200mg246248250252248250246248245250（数据来源：费宇等,《统计学》第5章,高等教育出版社,2010）给定显著性水平

,比较试验中咖啡因用量是否对人体神经功能有显著影响？(2)如果有显著差异，在0.05的显著性水平下，说明试验中咖啡因用量在哪些水平上有显著差异?2023/1/5《统计学实验》第5章方差分析5-9【统计理论】在一个实验中，因素A有r个水平

，

，

，

，在每个水平下作了

次实验。那么，在水平

下的第

次试验的观测值为

描述方差分析的统计模型为2023/1/5《统计学实验》第5章方差分析5-10

其中，

表示观测指标值的总平均，

表示因素

第

水平上的主效应，

表示随机误差；主效应参数

满足约束条件

研究因素的影响是否显著可以归结为比较这

个总体的均值，即检验如下假设【统计理论】该假设的检验可以通过平方和分解得到。2023/1/5《统计学实验》第5章方差分析5-11【统计理论】总平方和:组间平方和：组内平方和：2023/1/5《统计学实验》第5章方差分析5-12三种“

平方和”之间的关系平方和分解：【统计理论】2023/1/5《统计学实验》第5章方差分析5-13

由于上述几种平方和的数值受到样本量和水平数的影响，一种更为科学的方法是将各部分平方和除以相应自由度，其比值称为均方和，简称均方(meansquare，MS)，即

上式中分母的数值为对应平方和的自由度

。【统计理论】2023/1/5《统计学实验》第5章方差分析5-14

，则拒绝零假设，即认为因素对指标值（或试验结果）有显著影响，否则认统计量的数值大于对给定显著性水平

，没有显著影响。为因素为了检验

，定义F统计量【统计理论】2023/1/5《统计学实验》第5章方差分析5-15

将上述主要结果总结成一个表格，称为方差分析表，可以直观反映方差分析的计算及检验过程。表5.2单因素方差分析表【统计理论】2023/1/5《统计学实验》第5章方差分析5-16【统计理论】如果

检验的结论是拒绝

，则说明因素

的个

水平效应有显著差异，也就是说

个均值之间有显著差异。但这仅仅说明在

中至少有两个存在显著差异，这时我们还需要对每一对和作一对一的比较，即多重比较。

具体的说就是要比较第

组和第

组平均数，即检验

2023/1/5《统计学实验》第5章方差分析5-17注意到与是等价的。因此

该假设中共有个不同的成对比较。

【统计理论】

多重比较的特点是它同时对多个成对假设进行比较。多种比较的思想有两种，一是寻找每一个成对假设的检验统计量，给出检验临界值，通过比较界定显著程度；二是使用同时置信区间(simultaneousconfidenceinterval)的概念。

多重比较有许多种方法，使用比较多的包括Fisher的LSD方法，Turkey方法，Bonferroni方法等。2023/1/5《统计学实验》第5章方差分析5-18【菜单方式】打开数据文件li5.1.sav，选择“分析”→“比较均值”→“单因素方差分析”,将因变量y选入到“因变量列表”中，将因素KFTJL选入到“因素”中,点击进入“选项”对话框，

在“统计量”下选择“方差同质性检验”，点击“继续”，最后点击“确定”。(1)这是单因素方差分析问题。2023/1/5《统计学实验》第5章方差分析5-19表5.3咖啡因用量实验的方差一致性检验表5.4咖啡因用量实验的方差分析表输出结果方差齐性检验Levene

统计量df1df2显著性.292227.749单因素方差分析平方和df均方F显著性组间61.400230.7006.181.006组内134.100274.967

总数195.500292023/1/5《统计学实验》第5章方差分析5-20【程序方式】ONEWAYCSHBYKFYJL/STATISTICSHOMOGENEITY/MISSINGANALYSIS.2023/1/5《统计学实验》第5章方差分析5-21【菜单方式】打开数据文件li5.1.sav，选择“分析”→“比较均值”→“单因素ANOVA”,将因变量y选入到“因变量列表”中，将因素KFTJL选入到“因素”中,(２)这是一个多重比较的问题。

2023/1/5《统计学实验》第5章方差分析5-22在单因素方差分析的对话框中点击“两两比较”，进入后在“假设方差齐性”下选定，Turkey，LSD和Bonferroni选项，点击“继续”，最后点击“确定”。

【菜单方式】2023/1/5《统计学实验》第5章方差分析5-23表5.5咖啡因用量实验的多重比输出结果多重比较(I)

咖啡因剂量(J)咖啡因剂量均值差

(I-J)标准误显著性95%置信区间下限上限TukeyHSD0100-1.60000.99666.261-4.0711.8711200-3.50000*.99666.004-5.9711-1.028910001.60000.99666.261-.87114.0711200-1.90000.99666.156-4.3711.571120003.50000*.99666.0041.02895.97111001.90000.99666.156-.57114.3711LSD0100-1.60000.99666.120-3.6450.4450200-3.50000*.99666.002-5.5450-1.455010001.60000.99666.120-.44503.6450200-1.90000.99666.067-3.9450.145020003.50000*.99666.0021.45505.54501001.90000.99666.067-.14503.9450Bonferroni0100-1.60000.99666.360-4.1439.9439200-3.50000*.99666.005-6.0439-.956110001.60000.99666.360-.94394.1439200-1.90000.99666.202-4.4439.643920003.50000*.99666.005.95616.04391001.90000.99666.202-.64394.4439*.均值差的显著性水平为

0.05。2023/1/5《统计学实验》第5章方差分析5-24【程序方式】ONEWAYyBYKFYJL/MISSINGANALYSIS/POSTHOC=TUKEYLSDBONFERRONIALPHA(0.05).2023/1/5《统计学实验》第5章方差分析5-255.2双因素方差分析

5.2.1有可加效应的双因素方差分析【例5.2】（数据文件为li5.2.sav）有四种品牌(brand)的饮料在五个地区(district)销售，在每一个地区对每一种品牌的饮料销售量观测两次（上半年一次，下半年一次）得到数据如表5.6所示。在0.05的显著性水平下，问品牌及地区对饮料的销售量是否有显著影响？2023/1/5《统计学实验》第5章方差分析5-26因素A因素B123451360375356366367380336356375390228330130332531233332435529531033573703083253063303663863173304302312294310288301322340355370表5.6四种饮料品牌在五个地区的销售数据2023/1/5《统计学实验》第5章方差分析5-27【统计理论】在一个实验中，设有A和B两个因素，因素A有

个水平

，

，

，

，因素B有

个水平

，

，

，

，这样因素A和B就有

个水平组合。在因素A，B的每一种水平组合下均有

个样本观测值(次试验)，第

个样本值记为

2023/1/5《统计学实验》第5章方差分析5-28常数方差

。有可加效应的双因素方差分析模型其中，表示观

测指标值的总平均，

表示因素第水平上的主效应，

表示因素第水平上的主效应，

表示随机误差。

且

，

。【统计理论】

假设不同水平上观测数据相互独立，同一水平中，

n个重复观测数据也相互独立，

具有

2023/1/5《统计学实验》第5章方差分析5-29

对双因素方差模型，将涉及两个因素主效应的检验。因素的显著性假设为：

而对因素，显著性假设为【统计理论】

双因素方差分析与单因素方差分析的统计原理基本相同，上述两个零假设也可以通过平方和分解进行检验。

2023/1/5《统计学实验》第5章方差分析5-30【统计理论】因素A的平方和：误差平方和：因素B的平方和：2023/1/5《统计学实验》第5章方差分析5-31可以证明，各平方和之间满足：【统计理论】为了给出假设问题和的检验过程，

定义如下两个统计量：

2023/1/5《统计学实验》第5章方差分析5-32可以证明，当成立时，

当成立时，

【统计理论】时，当对给定的显著性水平

，拒绝

，即认为因素A对指标值有显著影响。时，当拒绝

，即认为因素B对指标值有显著影响。2023/1/5《统计学实验》第5章方差分析5-33【菜单方式】打开数据文件li5.2.sav，选择“分析”→“一般线性模型”→“单变量”

，将因变量y选入到“因变量”中，将brand和district两个因素选入到“固定因子”中

，点击“模型”，

然后选择“设定”（这样是模型中不包括交互效应），在“构建项”中选择“主效应”，再把brand和district选入“模型”，选择“在模型中包含截距”，以确定模型中包含常数项，点击“继续”，最后点击“确定”。

2023/1/5《统计学实验》第5章方差分析5-34表5.7

饮料品牌及销售数据的方差分析表主体间效应的检验因变量:销售量源III型平方和df均方FSig.校正模型20158.125a72879.7325.6640.000截距4482972.02514482972.0258817.7780.000brand16674.27535558.09210.9320.000district3483.8504870.9631.7130.171误差16268.85032508.402总计4519399.00040校正的总计36426.97539a.R方

=.553（调整

R方

=.456）2023/1/5《统计学实验》第5章方差分析5-35【程序方式】UNIANOVAyBYbranddistrict/METHOD=SSTYPE(3)/INTERCEPT=INCLUDE/CRITERIA=ALPHA(0.05)/DESIGN=branddistrict.2023/1/5《统计学实验》第5章方差分析5-365.2.2有交互效应的双因素方差分析【例5.3】(数据文件为li5.3.sav)在例5.2中，问品牌和地区是否存在交互效应？2023/1/5《统计学实验》第5章方差分析5-37当考虑有交互效应时，双因素方差分析模型表述为：上式中参数表示交互效应，

它满足约束条件

，【统计理论】2023/1/5《统计学实验》第5章方差分析5-38对具有交互效应的双因素方差模型中效应参数的检验，除5.2.1节中的和之外，与因素的交互效应显著性假设表述为：还涉及交互效应的假设检验。因素【统计理论】仿单因素方差分析的方法，各平方和之间满足2023/1/5《统计学实验》第5章方差分析5-39【统计理论】因素A的离差平方和因素B的离差平方和因素A和B的离差平方和误差平方和2023/1/5《统计学实验》第5章方差分析5-40

对具有交互效应的双因素方差分析问题，对参数的检验分为两个步骤：

第一步：首先检验交互效应，定义检验统计量：【统计理论】对给定的显著性水平

，时，当拒绝

，即认为交互效应显著。可以证明，当成立时，

。2023/1/5《统计学实验》第5章方差分析5-41第二步：如果在第一步中交互效应显著，定义检验

和的两个统计量分别为用类似的方法检验因素与因素的主效应是否显著。【统计理论】如果第一步中交互效应不显著，则因素与因素的主效应的检验采用5.2.1节中可加效应模型的方法检验主效应2023/1/5《统计学实验》第5章方差分析5-42

两因素方差分析的检验过程可以通过双因素方差分析表反映出来：表5.8

双因素方差分析表【统计理论】2023/1/5《统计学实验》第5章方差分析5-43【菜单方式】打开数据文件li5.2.sav，选择“分析”→“一般线性模型”→

“单变量”，将因变量y选入到“因变量”中，将brand和district两个因素选入到“固定因子”中

，2023/