SPSS学习系列22方差分析.docx
- 文档编号:4718350
- 上传时间:2022-12-07
- 格式:DOCX
- 页数:19
- 大小:34.08KB
SPSS学习系列22方差分析.docx
《SPSS学习系列22方差分析.docx》由会员分享,可在线阅读,更多相关《SPSS学习系列22方差分析.docx(19页珍藏版)》请在冰豆网上搜索。
SPSS学习系列22方差分析
22.方差分析
一、方差分析原理
1.方差分析概述
方差分析可用来研究多个分组的均值有无差异,其中分组是按影响因素的不同水平值组合进行划分的。
方差分析是对总变异进行分析。
看总变异是由哪些部分组成的,这些部分间的关系如何。
方差分析,是用来检验两个或两个以上均值间差别显著性(影响观察结果的因素:
原因变量(列变量)的个数大于2,或分组变量(行变量)的个数大于1)。
一元时常用F检验(也称一元方差分析),多元时用多元方差分析(最常用Wilks'A检验)。
方差分析可用于:
(1)完全随机设计(单因素)、随机区组设计(双因素)、析因设计、
拉丁方设计和正交设计等资料;
(2)可对两因素间交互作用差异进行显著性检验;
(3)进行方差齐性检验。
要比较几组均值时,理论上抽得的几个样本,都假定来白正态总体,
且有一个相同的方差,仅仅均值可以不相同。
还需假定每一个观察值都由若干部分累加而成,也即总的效果可分成若干部分,而每一部分都有一个特定的含义,称之谓效应的可加性。
所谓的方差是离均差平方和除以白由度,在方差分析中常简称为均方(MeanSquare)。
2.基本思想
基本思想是,将所有测量值上的总变异按照其变异的来源分解为多个部份,然后进行比较,评价由某种因素所引起的变异是否具有统计学意义。
根据效应的可加性,将总的离均差平方和分解成若干部分,每一部分都与某一种效应相对应,总白由度也被分成相应的各个部分,各部分的离均差平方除以各白的白由度得出各部分的均方,然后列出方差分析表算出F检验值,作出统计推断。
方差分析的关键是总离均差平方和的分解,分解越细致,各部分的含义就越明确,对各种效应的作用就越了解,统计推断就越准确。
£1方差分析表形式
变异素海
戚姜平方郛自由度均方F♦计室P—率值
SS甫MSFP
效应鬲栽6跖
效应龌
SStMSj-SSj/^;MSi/MSEPL
班VS2=SS2/^再(啪:
,矶户MWMS,P2
rii■A■Ad.Ji*Adka
s%,办的Fe压硕尸MWM%
SSe审宅MSe—ss#就
把变异勺
ssrssi+SSI*-**呢瓯*再―“…,
xiu*MSt=SSH所&如或Ft
+SS^n+SSC+茹)+就
效应项与试验设计或统计分析的目的有关,一般有:
主效应(包括各种因素),交互影响项(因素间的多级交互影响),协变量(来白回归的变异项),等等。
当分析和确定了各个效应项S后,根据原始观察资料可计算出各个离均差平方和SS再根据相应的白由度df,由公式MS=SSdf,求出均方MS,最后由相应的均方,求出各个变异项的F值,F值实际上是两个均方之比值,通常情况下,分母的均方是误差项的均方。
根据F值的分子、分母均方的白由度fl和f2,在确定显著性水平为0C情况下,由F(fl,f2)临界值表查得单侧Fa界限值。
当F
拒绝原假设Ho,也即这个效应项是很可能对总变异有实质影响的。
3.方差分析的实验设计
为了确定方差分析表中各个有关效应项,需要在试验设计阶段就作出安排,再根据设计要求进行试验,得出原始观察值,按原来设计方案算出方差分析表中的各项。
在试验设计阶段通常需要考虑如下4个方面:
(1)研究的主要变量(因变量)
即试验所要观察的主要指标,一次试验时可以有多个观察指标,方差分析时也可以同时对多个因变量进行分析;
(2)因素和水平
试验的因素(factor)可以是品种、人员、方法、时间、地区等等,
因素所处的状态叫水平(level)。
在每一个因素下面可以分成若干水平。
例如,某工厂的原料来白4个不同地区,那么用不同地区的原料生产的产品质量是否一致呢所要比较的地区就是因素,4个地区便是地区这一因素的4个水平。
当某个主要因素的各个水平间的主要因变量的均值呈现统计显著性时,必要时可作两两水平间的比较,称为均值间的两两比较。
(3)因素间的交互影响
多因素的试验设计,有时需要分析因素间的交互影响(interaction),2个因素间的交互影响称为一级交互影响(AXB);3个因素间的交互影响称为二级交互影响(AXBXC)。
当交互影响项呈现统计不显著时,表明各个因素独立,当呈现统计显著时,就需要列出这个交互影响项的效应,以助于作出正确的统计推断。
二、单因素方差分析
1个因变量,1个影响因素:
总差异Yj=平均差异(1+因素差异E+随机差异可
例1比较4种品牌的胶合板的耐磨性,各抽取5个样品,相同转速磨损相同时间测得磨损深度(mm),如下:
wear
brand
1
2.30A
2
232
A
3
2.40A
4
245
A
匚5
2.58A
6
235
「7
230
D
8
242H
g
2.60B
10
235B
11
2.20C
12
200
C
比较4个品牌胶合板的耐磨性有无差异
总差异Yij=平均磨损(1+品牌差异0(1+随机差异肉
1.【分析】一一【一般线性模型】一一【单变量】,打开“单变量”窗口,将变量“wear磨损深度”选入【因变量】框,"brand品牌”
2.点【两两比较】,打开“观测均值的两两比较”子窗口,勾选【假
定方差齐性】下的“LSf3、"S-N-K',点【继续】;
3.点【选项】,打开“选项”子窗口,勾选“描述统计”、“方差齐性
检验”,点【继续】;
3单突晶选改
-估计边阮均值——1
皂三岂因子交互正):
墅
(OVERALL)
U比原王苗成£.)置信区间调矛但)一L$D(无)
辅出
《方差卉it检瑜呵)一一.分布-水平图yyB度差图担)
□蛛芝报台U度检睑心_!
一段估计函教(。
)
W措述统计(功
U幻效怡计(自
□检验散能但)
皿参教估计(!
)-届好繇凝承可
显著理玖平(Y):
瓯置信区间为96.0%
蟠阕[取消(帮助||••■-
点【确定】,得到
描述性统计量
因变量:
磨损深度(mm)
地板品牌
均值
标准偏差
N
A
.11269
5
B
.11760
5
C
.11216
5
D
.03271
5
总计
.21771
20
给出每个品牌的均值、标准差、样本数。
误差方差等同性的Levene检验a
因变量:
磨损深度(mm)
F
df1
df2
Sig.
3
16
.311
检验零假设,即在所有组中因变量的误差方差均相等
a.设计:
截距+brand
方差齐性检验结果,P值=>,故接受原假设H0:
方差齐
主体间效应的检验
因变量:
磨损深度(mm)
源
III型平方和
df
均方
F
Sig.
校正模型
.740a
3
.247
.000
截距
1
.000
brand
.740
3
.247
.000
误差
.161
16
.010
总计
20
校正的总计
.901
19
a.R方=.822(调整R方=.788)
方差分析结果,“校正模型”是整个方差分析模型的检验,原假设
H0:
所有系数(沽宛应)都=0;P值<<,故拒绝原假设。
“截距”检验均值沽原假设H。
:
视(即不考虑品牌时,平均磨损为0);P值<<,故拒绝原假设
“brand”对因素品牌的检验,原假设H0:
按因素水平值的各分组的因变量无差异,即品牌因素对磨损深度无影响;P值<<,故拒绝原假设,即不同品牌的耐磨性有差异。
参数估计
因变量:
磨损深度(mm)
参数
B
标准误差
t
Sig.
95%置信区间
下限
上限
截距
.045
.000
[brand=A]
.063
.021
[brand=B]
.063
.017
[brand=C]
.063
.000
[brand=D]
0a
.
.
.
.
.
a.此参数为冗余参数,将被设为零。
B列为各品牌均值与均值pi(截距)的差
截距
参数
对比
_1
截距
1
[brand=A]
.250
[brand=B]
.250
[brand=C]
.250
[brand=D]
.250
此矩阵的缺省显示是相应的L矩
阵的转置。
基于III型平方和。
估计常数项时使用的L矩阵,均为即总样本的均值是按四种品牌等量混合的情况计算的。
brand
参数
对比
L2
L3
L4
截距
0
0
0
[brand=A]
1
0
0
[brand=B]
0
1
0
[brand=C]
0
0
1
[brand=D]
-1
-1
-1
此矩阵的缺省显示是相应的L矩阵的转置。
基于III型平方和。
对比系数矩阵,默认将最后一组“品牌D”作为对照组,故上上表的截距(均值。
的估计值=品牌D的均值=
L2=[0100-ij,对于L2列,令[(1aio(2心必]XL2=0化简得od=必即前表对淌作的假设检验。
多个比较
因变量:
磨损深度(mm)
I)地板品牌
(J)地板品牌
均值差值(I-J)
标准误差
Sig.
95%置信区间
下限
上限
B
.0060
.06339
.926
.1404
A
C
.3640*
.06339
.000
.2296
.4984
D
*
.06339
.021
A
.06339
.926
.1284
B
C
.3580*
.06339
.000
.2236
.4924
D
*
.06339
.017
_SD
A
*
.06339
.000
C
B
*
.06339
.000
D
*
.06339
.000
A
.1620*
.06339
.021
.0276
.2964
D
B
.1680*
.06339
.017
.0336
.3024
C
.5260*
.06339
.000
.3916
.6604
基于观测到的均值。
误差项为均值方(错误)=.010。
*.均值差值在.05级别上较显著。
LSD法给出的两两比较,将各组均和一个参照水平做比较,未指定默认,则每一个水平都作为参照比较一次。
每两个之间的差异有无
统计学意义,看对应的P值判断(原假设H。
:
无差异)。
磨损深度(mm)
地板品牌
N
子集
1
2
3
C
5
Student-Newman-Keulsa,b
B
5
A
5
D
5
Sig.
.926
已显示同类子集中的组均值。
基于观测到的均值。
误差项为均值方(错误)=.010。
a,使用调和均值样本大小=。
b.Alpha=.05。
LSD法给出的两两比较结果,将各组的值从小到大排序,注意4个品牌共被分成了3个亚组(无差异的作为一组),品牌B和A放在一个亚组,二者的P值=(无差异)。
三、两因素方差分析
1个因变量,2个影响因素:
总差异Yijk=平均差异"因素1差异od+因素2差异&i
+因素1,2交互作用差异下+随机差异邮例2分析超市某商品的销售量在不同的超市规模(小型、中型、大型)、货架位置(A、B、GD)是否有差异部分数据文件如下:
size
position
salps
1
1
A
低0
2
1
A
50.0
1
B
56.0
4
1
B
63.0
5!
1
C
65.。
6
1
C
71.0
1
1
D
阳0
8!
1
D
53.0
g
2A
57.0
匚10
2A
65.0
11
2
Bby0
「12
2
B
78.0
变量size超市规模:
1=小型,2=中型,3=大型。
总差异Yijk=平均差异[1+超市规模差异cd+货架位置差异&+超市规模货架位置交互作用差异吊+随机差异8jk
1.【分析】一一【一般线性模型】一一【单变量】,打开“单变量”窗口,将变量“sale销售量”选入【因变量】框,将变量“size超市
2.点【选项】,打开“选项”子窗口,勾选【输出】下的“描述统计”、
“方差齐性检验”,点【继续】;
点【确定】,得到
主体间因子
值标签
N
1
小型
8
超市规模2
中型
8
3
大型
8
A
6
B
6
摆放位置
C
6
D
6
描述性统计量
因变量:
周销售量
照市规模摆放位置均值
标准偏差
B
小型
C
D
N
2
2
2
2
总计
8
A
2
B
2
中型C
2
D
2
总计
8
A
2
B
2
大型C
2
D
2
总计
8
A
6
B
6
总计C
6
D
6
总计
24
误差方差等同性的Levene检验a
因变量:
周销售量
F
df1
df2
Sig.
.
11
12
.
检验零假设,即在所有组中因变量的误差方差均相等
a.设计:
截足巨+size+position+size*position
超市规模3个水平,货架位置4个水平,共将样本分成3X4=12
组,由于有单组样本数<3个,故无法做方差齐性检验(值缺失)。
主体间效应的检验
因变量:
周销售量
源
II型平方和
df
均方
F
Sig.
校正模型
11
.000
截距
1
.000
size
2
.000
position
3
.000
size*position
6
.689
.663
误差
12
总、计
24
校正的总计
23
a.R方=.921(调整R方=.849)
整个方差分析模型的检验结果,交互作用项size*position的P值=>,故接受原假设H0:
该交互作用无差异。
下面去掉交互因子继续做两因素方差分析。
3.在第1步的窗口点【模型】,打开“模型”子窗口,选择【指定模型】下的“设定”,将【构建项】下的【类型】设为“主效应”,将变
量“size”、“position”选入【模型】框,点【继续】;
4.原窗口点【两两比较】,打开“观测均值的两两比较”子窗口,将因子“sizd、"position”选入【两两比较检验】框,勾选【假定方差齐性】下的“S-N-K,点【继续】;
注:
若已明确对照组,考察其它组与它的比较,宜采用进行多个均值间的两两比较,且各组人数相等,宜采用
LSD法;若要
Tukey法或
S-N-K法(若比较的组数特别多,不宜用S-N-K法,宜用Scheffe法);对于不平衡设计或含有协变量的模型,应采用LSD法、Bonferroni法、
Sidak法。
点【确定】得到:
误差方差等同性的Levene检验a
因变量:
周销售量
F
df1
df2
Sig.
.171
11
12
.997
检验零假设,即在所有组中因变量的误差方差均相等
a.设计:
截足巨+size+position
方差齐性检验,P值=>,故接受原假设H0,即方差齐
主体间效应的检验
因变量:
周销售量
源
III型平方和
df
均方
F
Sig.
校正模型
5
.000
截距
1
.000
size
position
误差
总计
校正的总计
2
3
18
24
23
.000
.000
a.R方=.894(调整R方=.865)
整个方差模型的检验结果(解释参考例1)
周销售量
Student-Newman-Keuls
超市规模
N
子集
1
2
3
小型
8
中型
8
大型
8
Sig.
已显示同类子集中的组均值。
基于观测到的均值。
误差项为均值方(错误)=。
a.使用调和均值样本大小=
b.Alpha=.05。
周销售量
Student-Newman-Keuls
摆放位置
N
子集
1
2
3
D
A
B
C
Sig.
6
6
6
6
.948
已显示同类子集中的组均值。
基于观测到的均值。
误差项为均值方(错误)=。
a.使用调和均值样本大小=。
b.Alpha=.05。
用S-N-K法进行两两比较,可见超市规模越大,销售量越大;货架位置对销售量也有影响,位置AD在同一亚组,销售量最小,位置B销售量居中,位置C销售量最大,三个亚组之间有统计学差异;另
外,由于交互作用被合理地剔除,故上述差异不受另一因素(超市规模)取值的影响
5.若要绘制轮廓图。
原窗口点【绘制】,打开“轮廓图”子窗口,将因子"sizd'、aposition”分别选入【水平轴】点【添加】,点【继续】;
注:
若要得到两变量的联合轮廓图,将另一变量选入【单图】框即可点【确定】,得到单变量的轮廓图:
边际均值,是基于现有模型,控制了其它因素作用后,根据样本情况计算某因素各水平的均值估计值(若模型中有协变量,会按协变量均值加以修正)。
轮廓图,即以边际均值为纵轴,以考察因素为横轴的折线图。
用以比较该因素取不同水平值时,样本均值的变化情况。
另外,轮廓图也可用来检验两因素是否存在交互作用:
对于单因素模
型或包含全部交互项的全模型,边际均值就是各分组的样本均值,其轮廓图就呈现一组平行线;若剔除某交互作用后各曲线明显不平行,则说明两因素存在交互作用。
另外,【选项】子窗口也提供了“缺乏拟合优度检验”,勾选它,运行得到
失拟检验
因变量:
周销售量
源
平方和
df
均方
F
Sig.
失拟
纯误差
6
12
.689
.663
用来检验当前模型(剔除交互项)与全模型(包括全部交互项)的比
较,原假设板两模型无差别;本例的P值=>,接受原假设,即两因素超市规模、货架位置的交互作用可以忽略。
6.若要绘制残差图。
原窗口点【选项】,勾选【输出】下的“残差图”,运行得到
匹I变量:
同销忡筮
。
箫
C°°
Q0
容廿fto(5*
OGO0
&冲°
◎°%。
心砖
0旨0,0g88
吨生-祗市+size■+pasfion
残差图给出了因变量的实测值、预测值、标准化残差的散点图,若预
测值与实测值有明显的相关性(接近直线趋势),标准化残差在0附近随机分布,则表明拟合结果较好。
7.除两两比较外,也可以白定义比较。
下面只说明原理,具体操作需要借助代码实现。
例如,前文比较货架位置A与D时,L矩阵=[100-1]T,有
[ABCD卜[100-1]T=0等价于A=D
前面分析发现位置A与D的销售量基本无差异,现在想将A与D合并再与B比较有无差异,则可以指定L矩阵=[1-201]T,贝U
[ABCD卜[1-201]T=0等价于(A+D)2=B
注意:
是从(A+D)2=B倒推L矩阵,该式即A-2B+0C+D=0.
四、含随机因素的方差分析
随机因素设为固定因素作为分析,可能得到错误的结果。
例3研究4种广告方式(店内展示、发放传单、推销员展示、广播广告)有无差异。
该地区有几百个销售网点,经费有限只随机选取了18个网点,记录了固定时间段内使用某种广告方式的销售额(为减小误差,各网点重复测量两次):
area
adstyp白
salps
\1
1
1
52
匚2
1
1
41
3
2
1
61
4
2
1
44
5
3
1
61
6
3
1
86
7
4
1
76
8
4
1
75
「9
5
1
57
10
5
1
75
11
6
1
52
12
6
1
变量area表示网点;adstype表示广告类型:
1=店内展示,2=发放传单,3=推销员展示,4=广播广告;sales表示销售额。
由于网点是随机选取的,若重复研究重新抽取的网点可能完全不同,故变量area属于随机因素。
注:
若对区域进行细分归类,每类区域选代表网点,则不是随机因素
口,将变量“sales销售额”选入【因变量】框,将“adstype广告类
型”选入【固定因子】框,将“area网点”选入【随机因子】框;
点【确定】得到
主体间效应的检验
因变量:
销售额
源
III型平方和
df
均方
F
Sig.
假设
1
.000
截距
误差
17
假设
3
.000
adstype
误差
51
假设
17
.000
area
误差
51
假设
51
.286
adstype*area
误差
72
a.MS(area)
b.MS(adstype*area)
c.MS(错误)
整个方差分析模型的检验结果,注意当模型含有随机因素时,不再进行总模型的检验,而是分别对每个因素做单独检验,并给出单独的误差项。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- SPSS 学习 系列 22 方差分析