两因子方差分析.docx
- 文档编号:9974961
- 上传时间:2023-02-07
- 格式:DOCX
- 页数:18
- 大小:92.06KB
两因子方差分析.docx
《两因子方差分析.docx》由会员分享,可在线阅读,更多相关《两因子方差分析.docx(18页珍藏版)》请在冰豆网上搜索。
两因子方差分析
两因子析因设计的方差分析
1.概述
本文介绍如何分析多因子方差分析。
先从两因子析因设计开始。
在析因设计方差分析中,一个处理的每个水平出现在另一个处理的所有水平下。
2x3析因设计如下。
2x3因子有6个单元
图1.2x3全因子方差分析示意图
因子B
B1
B2
B3
因子A
A1
A2
水平A1出现在B的所有水平下,水平A2也出现在B的所有水平下。
本文介绍如何用GLM:
Univariate过程分析这种类型的数据。
2.数据
数据来源于表演的行为研究。
参与者是24只猴子。
它们的任务是完成一个“单数”的问题。
给猴子出示三个物品,例如两个圆环和一个方块,有一个奖励食物放置在单数物品下的小坑里。
因变量(score)是猴子选择单数物品和得到奖励的实验总数。
有两个自变量,奖励(reward)有三个水平(1个葡萄,3个葡萄,和5个葡萄),驱赶(drive)有两个水平(剥夺1小时食物和剥夺24小时食物)。
数据保存在文件glm-2way.sav中。
变量如表1所示。
表1.glm-2way.sav数据文件中的变量
变量
变量标签/值标签/缺失值
reward
奖励数量/
1="1个葡萄"
2="3个葡萄"
3="5个葡萄"
drive
动物的驱赶水平(剥夺时数)/
1="1小时"
2="24小时"
score
20次试验中的正确数
3.假设检验
全因子受试者间方差分析的假设检验如表2所示。
表2.方差分析假设检验
1.每个单元中观测对象是独立的。
2.因变量至少是间隔以上的测量尺度。
3.每个单元的分布形状对称。
4.每个单元的分布同质。
假设1(独立性).每个设计单元中有不同的参与者,所以数据是独立的。
假设2(测量尺度).得分的测量尺度是比率。
假设3(正态性). 假设数据在6个设计单元中服从正态分布。
如果分布对称,那么方差分析是稳健的。
这个假设能用探索过程检验每个设计单元的正态性。
然而,探索过程默认的分析是在每个选择因子的主效应上计算所需要的统计量。
对这个数据,当在FactorList框中放置reward和drive时,探索过程将在奖励的三个水平和驱赶的两个水平内分别进行正态性检验。
它不会在6个设计单元内按照假设要求计算所需要的统计量。
所以,为了生成需要的检验,必须使用语法命令改变这一切。
打开探索过程对话框
Analyze
DescriptiveStatistics
Explore...
将变量score移到DependentList窗口,变量drive和reward移到FactorList窗口。
选择茎叶图,因子水平在一起的箱图,带检验的正态图,带Levene检验的展布和水平幂估计和描述统计。
运行探索过程,看到它没有给出方差分析假设的正确检验。
注意,输出显示了每个因子的统计量。
更确切地说,输出描述奖励和驱赶的主效应,没有给出六个设计单元的统计量。
回到探索过程对话框,点击Paste出现语法窗口。
语法窗口的内容如下。
表3.探索过程语法
EXAMINE
VARIABLES=scoreBYdrivereward
/PLOTBOXPLOTSTEMLEAFNPPLOTSPREADLEVEL
/COMPAREGROUP
/STATISTICSDESCRIPTIVES
/CINTERVAL95
/MISSINGLISTWISE
/NOTOTAL.
在2x3设计中,要得到每个单元的统计量,仅仅需要在命令
VARIABLES=scoreBYdrivereward中的驱赶和奖励因子间插入关键字BY。
正确的语法文件显示在表4中。
插入的BY用蓝色显示。
表4修正的探索过程语法
EXAMINE
VARIABLES=scoreBYdriveBYreward
/PLOTBOXPLOTSTEMLEAFNPPLOTSPREADLEVEL
/COMPAREGROUP
/STATISTICSDESCRIPTIVES
/CINTERVAL95
/MISSINGLISTWISE
/NOTOTAL.
然后执行语法命令。
注意:
方差一致性假设能由GLM或探索过程来检验。
探索过程建议使用幂变换简化一致性问题。
GLM做不到。
表5给出偏度和峰度统计量的概要。
指出每个单元中的数据是正态分布。
在24小时和5个葡萄的单元中峰度有些问题,但是这不是我们所关心的,因为当分布对称时,方差分析是稳健的。
表5.Skewness(SESkewness)andKurtosis(SEKurtosis)Statistics
Reward
1grape
3grapes
5grapes
DriveLevel
1hour
skewness
0.63(1.01)
0.00(1.01)
-0.60(1.01)
kurtosis
-1.70(2.62)
-4.34(2.62)
-0.77(2.62)
24-hours
skewness
-0.60(1.01)
0.00(1.01)
0.00(1.01)
kurtosis
-0.77(2.62)
-3.30(2.62)
-5.41(2.62)
科尔莫果诺夫—沙米尔诺夫正态性检验没有提供显著性水平是因为每个单元中观测数太少。
假设4(一致性).Levene一致性检验显示方差是一致的,见表6。
方差一致性假设能由GLM或探索过程来实现。
探索过程建议使用幂变换简化一致性问题。
GLM做不到。
表6.基于均值的方差一致性检验
LeveneStatistic
df1
df2
Sig.
Numbercorrectonthe20trainingtrials
1.000
5
18
.446
无论何时,进行Levene检验,都应该检查自由度。
Df1是设计中组间单元数减1。
在这个设计中,组间单元有6个,所以df1是5。
如果在探索过程的语法中忘记增加BY项,会有几个Levene检验,在设计中有几个因子就有几个Levene检验。
本例中,有两个Levene检验,自由度为1的驱赶因子有一个Levene检验,和自由度为2的奖励因子有一个Levene检验。
箱图从视觉上提供了研究中发生了什么,见图1。
首先,方差看起来相当一致。
在剥夺24小时食物的高驱赶水平中,动物的表演没有受到奖励数量的影响。
然而,在剥夺1小时食物的低驱赶水平中,奖励越多,动物的表演得分越高。
箱图提示方差分析应该显示奖励和驱赶的显著性交互作用。
图1.箱图显示奖励与驱赶的正确响应数。
你的箱图可能与图1不同。
箱图是按照因子进入FactorList窗口的次序组织的。
当先移动驱赶水平后移动奖励数量到FactorList窗口时,箱图的次序(从左到右)如表7所示。
表7.当因子次序是驱赶和奖励时箱图的次序
Boxorder
(lefttoright)
DriveLevel
(2levels)
RewardMagnitude
(3levels)
1stbox
1
1
2ndbox
1
2
3rdbox
1
3
4thbox
2
1
5thbox
2
2
6thbox
2
3
规律是:
后进先变。
最后进入因子的指针最先变动。
在这种情况下,驱赶固定在水平1,最后进入因子(奖励)的指针在奖励的水平间逐渐增加。
然后,驱赶增加到水平2,最后进入因子的指针再重新在奖励的水平间逐渐增加。
当先移动奖励数量后移动驱赶水平到FactorList窗口时,箱图的次序(从左到右)如表8所示。
确定箱图从左到右次序的规则相同。
表8..当因子次序是奖励和驱赶时箱图次序
Boxorder
(lefttoright)
RewardMagnitude
(3levels)
DriveLevel
(2levels)
1stbox
1
1
2ndbox
1
2
3rdbox
2
1
4thbox
2
2
5thbox
3
1
6thbox
3
2
4.选择GLM过程
选择GLM过程
Analyze
GeneralLinearModel
Univariate...
移动变量score到DependentVariable窗口。
Reward和drive都是固定因子,移动它们到FixedFactor(s)窗口。
在这个研究中没有随机因子也没有协变量。
事实上,我不记得在心理学的本科生或研究生论文中用过随机因子(除了对象外)。
变量能在加权最小二乘法(WLSWeight)中指定为加权变量。
在心理学研究中这些是不常用的选项。
考虑这个问题有两种不同的方法。
在对话框中检查所有的选项按钮,选择每一种可能的选项,然后非常吃力地读完GLM过程产生的所有输出。
或用我推荐的方法,查看由GLM产生的基本输出,然后决定增加什么附加信息。
所以,确定因变量,自变量和描述统计量后,点击OK按钮,看看得到什么。
5.基本输出
GLM基本输出有两张表,第一张表是组间因子和因子水平观测数表,见表9。
表9.Between-SubjectsFactors
ValueLabel
N
Drivelevelofanimals(hoursofdeprivation)
1
1hourdeprived
12
2
24hoursdeprived
12
Magnitudeofreward
1
1grape
8
2
3grapes
8
3
5grapes
8
第二张是描述统计表,在下一部分讲解。
即,6.显著性效应解释:
显示均值
第三张表包含方差分析的结果,见表10。
包括平方和,F值,显著性水平。
奖励和驱赶的交互作用显著,F(2,18)=3.927,p=.038。
两个主效应都不显著,驱赶的F(1,18)=1.309,p=.268,奖励的F(2,18)=3.055,p=.072.
表10.TestsofBetween-SubjectsEffects
DependentVariable:
Numbercorrectonthe20trainingtrials
Source
TypeIIISumofSquares
df
MeanSquare
F
Sig.
CorrectedModel
280.000(b)
5
56.000
3.055
.036
Intercept
2400.000
1
2400.000
130.909
.000
DRIVE
24.000
1
24.000
1.309
.268
REWARD
112.000
2
56.000
3.055
.072
DRIVE*REWARD
144.000
2
72.000
3.927
.038
Error
330.000
18
18.333
Total
3010.000
24
CorrectedTotal
610.000
23
aRSquared=.459(AdjustedRSquared=.309)
Intercept.方差分析中的截距项检验总均值是否为零。
因为所有因变量得分为正,所以总均值不为零。
因此截距的检验我们不感兴趣。
Correctedmodel.带5个自由度的校正模型是总体模型。
包括两个主效应和交互作用的变异,因此自由度为5。
当设计平衡时,即当每个单元中观测数相等,没有缺失单元时,校正模型的平方和是主效应和交互效应的平方和之和。
在这种情况下,SScorrectedmodel=SSdrive+SSreward+SSdrive*reward
280.00=24.00+112.00+144.00
有一个RSquared与校正模型相关联,见注释a。
R方是校正模型占因变量变异的比例。
本例中,两个主效应和交互效应占得分变异的46%。
AdjustedRSquare是总体中预测模型的估计。
它总是小于R方。
小于R方多少取决于模型中变量数和样本量。
变量固定时,样本量越小,相关性越大。
模型中样本量固定时,变量越多,相关性越大。
在这种情况下,调整R方占因变量方差的31%。
调整R方有时被叫做缩水R方。
correctedtotal的平方和是校正模型和误差项的平方和之和。
在心理学中,校正模型的统计量在报告不常使用。
我们感兴趣的是主效应和交互效应而不是总体模型。
total平方和是截距,主效应,交互效应和误差项的平方和之和。
Error行给出单元内误差项的统计量。
均方误差用于检验主效应和交互效应。
6.显著性效应解释:
显示均值
方差分析表中,奖励和驱赶交互作用显著。
下一步是设法理解交互作用。
这需要查看交互作用的均值。
对2x3设计,有三种方法查看交互作用的均值。
Descriptives
第一选择是查看描述统计表。
描述统计表给出设计单元的均值,标准差和观测数。
因为是两因子方差分析,所以描述统计表也有两因子交互作用的描述统计量。
如果是三因子或更多因子的设计,描述统计表将给出更高阶的交互作用均值,但不会给出低阶交互作用的均值。
描述统计的复选框位于Options部分。
在左下部找Display,点击DescriptiveStatistics。
描述统计显示在表11中。
看看均值如何在表11中组织的。
当浏览表时,显示的第一个变量是驱赶水平,第二个变量是奖励数量。
因为这种组织,奖励数量的每个水平的均值显示在两个驱赶水平内。
表的底部total汇总的均值是奖励数量的主效应均值。
这个次序是由对话框中初始选择因子的次序决定的。
在这个例子中,驱赶因子列在前面。
如果奖励数量列在前面,那么表11中第一个变量是奖励数量。
驱赶的每个水平的均值显示在奖励的每个水平内。
底部汇总是驱赶水平主效应的主效应均值。
表11.DescriptiveStatistics
Drivelevelofanimals(hoursofdeprivation)
Magnitudeofreward
Mean
Std.Deviation
N
Numbercorrectonthe20trainingtrials
1hourdeprived
1grape
3.00
3.16
4
3grapes
10.00
4.76
4
5grapes
14.00
3.92
4
Total
9.00
5.97
12
24hoursdeprived
1grape
11.00
3.92
4
3grapes
12.00
5.48
4
5grapes
10.00
4.08
4
Total
11.00
4.20
12
Total
1grape
7.00
5.40
8
3grapes
11.00
4.87
8
5grapes
12.00
4.28
8
Total
10.00
5.15
24
主效应均值在Total行中。
例如,1个葡萄,3个葡萄和5个葡萄的奖励主效应分别是7.00,11.0和12.0。
驱赶1小时和24小时的主效应分别是9.0和11.0。
注释:
如果设计单元中观测数相等,描述统计表中主效应均值可以使用DescriptiveStatistics选项。
在这个例子中,每个设计中单元有4个动物,所以主效应均值是正确的。
如果单元数不等,那么报告中用"EstimatedMarginalMeans"的均值而不是DescriptiveStatistics选项的均值。
EstimatedMarginalMeans
还有一种方法呈现显著性效应的均值。
Options对话框的顶部有一个EstimatedMarginalMeans.对话框左边有Factors(s)andFactorInteractions。
移动主效应或交互效应到DisplayMeansfor。
对话框会出现主效应和交互效应。
移动显著的drive*reward交互作用到DisplayMeans框。
drive*reward的估计边际均值显示在表12中。
除了均值外,还有标准误和均值的95%置信区间。
注意:
单元数和标准差在边际估计均值中没有显示。
表12.奖励和驱赶的估计边际均值
DependentVariable:
Numbercorrectonthe20trainingtrials
Mean
Std.Error
95%ConfidenceInterval
Drivelevelofanimals(hoursofdeprivation)
Magnitudeofreward
LowerBound
UpperBound
1hourdeprived
1grape
3.000
2.141
-1.498
7.498
3grapes
10.000
2.141
5.502
14.498
5grapes
14.000
2.141
9.502
18.498
24hoursdeprived
1grape
11.000
2.141
6.502
15.498
3grapes
12.000
2.141
7.502
16.498
5grapes
10.000
2.141
5.502
14.498
ProfilePlots
选择Plots,均值能用图示方式表示。
这个选项会产生特征图。
特征图的纵轴总是表示因变量(本例中是得分)。
可以选择因子作为水平轴。
对这个研究,选择奖励因子作为水平轴。
接下来,选择一个因子作为独立折线或独立图。
如果选择驱赶水平作为独立折线,那么在特征图中每个驱赶水平会以各自的线条来显示。
如果选择驱赶水平为独立图,那么驱赶的两个水平分两个图来显示。
现在选择驱赶水平作为独立折线。
最后按下Add按钮完成定义特征图的处理。
特征图显示在图2中。
图2.奖励与驱赶交互作用的特征图
通过绘制不同驱赶水平的折线,特征图着重于驱赶水平内奖励的效应。
仔细查看特征图发现两条线不平行。
事实上,不平行是表示存在交互作用。
平行线表示没有交互作用。
高驱赶水平(24小时)的折线相当平坦。
即,如果动物有高驱赶水平,奖励数量对表演没有影响。
低驱赶水平(1小时)的折线有正斜率,奖励越多,动物表演越好。
总结如下,特征图提示在高驱赶水平下,奖励对动物的表演没有影响,对低驱赶水平下,奖励肯定与的动物表演有关。
如果选择驱赶因子作为水平轴,特征图的结果着重于奖励水平内驱赶的效应。
看了特征图后,如果认为特征图对描述数据有帮助,就要返回并选择另一个因子作为水平轴。
对这个数据,试着用每一个因子作水平轴。
思考一下那一个因子对描述数据更有帮助。
通过只选择水平轴的一个因子建立主效应特征图。
即,不选择独立折线或独立图作为因子。
试着建立奖励主效应的特征图。
7.显著性效应解释:
事后比较
MainEffects
在这个分析中主效应不显著。
如果显著,GLM中有两种方式可实现显著主效应的事后比较:
(a)进入PostHoc对话框。
选择待选的主效应,选择需要进行的检验类型,然后按继续按钮.事后比较检验与单因子方差分析中描述的相同。
(b)进入Options对话框的估计边际均值部分。
从因子和因子交互作用窗口选择待选的主效应,移动它们到显示均值窗口。
钩选比较主效应选项,按继续按钮。
使用这个选项有三种方式调整置信区间:
(a)none,用未校正的t检验进行个案比较,这种方法也叫做LSD检验;
(b)Bonferroni,个案的alpha水平是Bonferronialpha=alpha/C
这里C是主效应检验可能的比较对数;
(c)Sidak,个案的alpha水平是Sidakalpha=1-(1-alpha)1/C
这里C是主效应检验可能的比较对数。
在奖励主效应上练习TukeyHSD事后比较。
因为奖励主效应不显著,输出显示均值间没有差别。
InteractionEffects
检查交互作用均值暗示高驱赶水平内奖励没有影响,但在低驱赶水平内奖励有影响,奖励越多,表演越好。
统计上检验这些观测的一个可行方法是进行事后比较。
事后比较检验包含在PostHoc选项中。
左上部的Factor(s)显示进行事后比较检验的因子。
移动因子到PostHocTestsfor窗口进行所选因子的事后比较。
如你所见,存在一个问题。
我们想检验交互作用均值,但是这里只列出主均值。
似乎没有办法解决这个问题。
如果试图编辑事后比较检验的GLM语法,会得到一个错误消息,星号是非法符号的,只有主效应能定义事后比较。
GLM只能进行主效应因子上的事后比较。
这部分讨论显著性交互作用成对比较检验的两种选择:
(a)使用Tukey'sHSD用手算检验;
(b)在GLM中建立检验因子交互作用的单因子检验。
交互作用效应:
手算运行Tukey'sHSD检验
交互作用间检验差别的一个方案是通过手算进行事后比较检验。
ψ(HSD)是等观测数n设计的临界值,
这里,MSerror是方差分析中的均方误差;n是一个单元的观测数;qα,p,v是由给定显著性水平为α,设计单元数为p,均方误差项自由度为v的StudentizedRangeStatistic表的百分点。
本例中,MSerror=18.333,n=4,α=.05,p=6个单元,和v=18。
StudentizedRangeStatistic表的qα,p,v值是4.49.替代HSD公式得到
=4.49*2.14
=9.61
大于9.61的均值差彼此间有明显的差别。
研究中所有配对均值差显示在表13中。
表13. MeanDifferencesforallPossiblePairedComparisons
Diff(i-j)=rowi-columnj
1hourdep
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 因子 方差分析