中介效应检验程序及其应用.docx
《中介效应检验程序及其应用.docx》由会员分享,可在线阅读,更多相关《中介效应检验程序及其应用.docx(13页珍藏版)》请在冰豆网上搜索。
中介效应检验程序及其应用
收稿日期:
2004202209
3
本研究得到全国教育科学”十五”规划教育部重点课题(DBA010169以及香港中文大学和华南师范大学心理应用研究中心(教育部文科基地资助。
中介效应检验程序及其应用3
温忠麟1,2 张 雷2 侯杰泰2 刘红云3
(1华南师范大学教育科学学院,广州510631 (2香港中文大学教育学院,香港
(3北京师范大学心理学院,北京100875
摘 要 讨论了中介变量以及相关概念、中介效应的估计;比较了检验中介效应的主要方法;提出了一个检验程序,它包含了依次检验和Sobel检验。
该程序检验的第一类和第二类错误率之和通常比单一检验方法小,既可以做部分中介检验,也可以做完全中介检验。
作为示范例子,引入中介变量研究学生行为对同伴关系的影响。
关键词 中介变量,检验,第一类错误率,功效,同伴关系分类号 B841.2
1 引言
中介变量(mediator是一个重要的统计概念,国外涉及中介变量的研究很多,依PsycINFO检索结果,仅1998~2003年涉及中介变量且在标题中含
有相关词条(包括mediator、mediating或mediation的就有1100多篇。
国内对中介变量的研究很少,依中国期刊网“文史哲”和“教育与社会科学”专栏目录的检索结果,1998~2003年在标题或关键词中含有“中介变量”或“中介效应”的文章不足20篇。
这些文章中,有些只是做定性分析,说说中介变量而已;有些虽然做了统计分析,但没有中介变量分析。
只有少数几篇做了涉及中介变量的统计分析,其中较好的如文[1]和[2],都使用了结构方程分析,但对中介效应的分析还是有点粗略。
以文[1]为例,标题中出现了“中介变量”,可见中介变量及其效应分析应当是该文的重点。
文中虽然估计了中介效应的大小,但没有报告中介效应的相对大小(即中介效应占全部效应的比例,最大的不足是没有对中介效应进行必要的检验。
国内涉及中介变量的统计研究稀少并且总体质量不高的原因可能是多方面的,缺少方法论的研究是一个重要原因。
事实上,国内还未见到专门讨论如何分析中介效应的文章。
检验中介效应的方法很多,MacKinnon等人[3]通过模拟研究比较了十几种检验方法。
他们倾向于
使用该文主要作者1998年提出的一个高功效
(power的检验方法[4]。
但该方法的主要缺点是在有些情况下,第一类错误率太大,远远高于给定的显著性水平[3]。
本文提出了一个包含两种检验方法的检验程序,其中并没有MacKinnon等人的方法[4]。
新的检验程序很好地控制了第一类错误率,同时又有较高的检验功效。
使用新提出的检验程序,我们在“学生行为对同伴关系的影响”研究中对两个变量(教师喜欢程度和教师管教方式进行了中介效应的检验。
2 中介变量和相关概念
在本文中,假设我们感兴趣的是因变量(Y和自变量(X的关系。
虽然它们之间不一定是因果关系,而可能只是相关关系,但按文献上的习惯而使用“X对Y的影响”、“因果链”的说法。
为了简单明确起见,本文在论述中介效应的检验程序时,只考虑一个自变量、一个中介变量的情形。
但提出的检验程序也适合有多个自变量、多个中介变量的模型(见5.1节的讨论。
2.1 中介变量的定义
考虑自变量X对因变量Y的影响,如果X通过影响变量M来影响Y,则称M为中介变量。
例如,“父亲的社会经济地位”影响“儿子的教育程度”,
进而影响“儿子的社会经济地位”[5]。
又如,“工作
心 理 学 报 2004,36(5:
614~620 ActaPsychologicaSinica
614
环境”(如技术条件通过“工作感觉”
(如挑战性影响“工作满意度”[6]
。
在这两个例子中,“儿子的教育程度”和“工作感觉”是中介变量。
假设所有变量都已经中心化(即均值为零,可用下列方程来描述变量之间的关系(相应的路径图见图1:
Y=cX+e1
(1 M=aX+e2
(2 Y=c’X+bM+e3
(3 假设Y与X的相关显著,意味着回归系数c显著(即H0:
c=0的假设被拒绝,在这个前提下考虑
中介变量M。
如何知道M真正起到了中介变量的
作用,或者说中介效应(mediatoreffect显著呢?
目前有三种不同的做法[3]
。
图1 中介变量示意图
传统的做法是依次检验回归系数[7,8]。
如果下
面两个条件成立,则中介效应显著:
(i自变量显著影响因变量;(ii在因果链中任一个变量,当控制了它前面的变量(包括自变量后,显著影响它的后继变量。
这是Baron和Kenny定义的(部分中介过程[8]。
如果进一步要求:
(iii在控制了中介变量后,自变量对因变量的影响不显著,变成了Judd和Kenny定义的完全中介过程[7]。
在只有一个中介变
量的情形,上述条件相当于(见图1:
(i系数c显著(即H0:
c=0的假设被拒绝;(ii系数a显著(即
H0:
a=0被拒绝,且系数b显著(即H0:
b=0被
拒绝。
完全中介过程还要加上:
(iii系数c’不显著。
第二种做法是检验经过中介变量的路径上的回归系数的乘积ab是否显著,即检验H0:
ab=0,如果拒绝原假设,中介效应显著[4,9,10]
这种做法其
实是将ab作为中介效应。
第三种做法是检验c’与c的差异是否显著,即
检验H0:
c-c’=0,如果拒绝原假设,中介效应显著[11,12]。
2.2 中介效应与间接效应 依据路径分析中的效应分解的术语[13],中介效应属于间接效应(indirecteffect。
在图1中,c是X对Y的总效应,ab是经过中介变量M的间接效应(也就是中介效应,c’是直接效应。
当只有一个自变量、一个中介变量时,效应之间有如下关系
c=c’+ab
(4 当所有的变量都是标准化变量时,公式(4就是相关系数的分解公式[13]。
但公式(4对一般的回归系数也成立(证明见[14]。
由公式(4得c2c’=
ab,即c2c’等于中介效应,因而检验H0:
ab=0与H0:
c2c’=0是等价的。
但由于各自的检验统计量
不同,检验结果可能不一样。
中介效应都是间接效应,但间接效应不一定是中介效应。
实际上,这两个概念是有区别的。
首先,当中介变量不止一个时,中介效应要明确是哪个中介变量的中介效应,而间接效应既可以指经过某个特定中介变量的间接效应(即中介效应,也可以指部分或所有中介效应的和。
其次,在只有一个中介变量的情形,虽然中介效应等于间接效应,但两者还是不等同。
中介效应的大前提是自变量与因变量相关显著,否则不会考虑中介变量。
但即使自变量与因变量相关系数是零,仍然可能有间接效应。
下面的人造例子可以很好地说明这一有趣的现象([15],p128;也可参见[13],p48。
设Y是装配线上工人的出错次数,X是他的智力,M是他的厌倦程度。
又设智力(X对厌倦程度(M的效应是0.707(=a,厌倦程度(M对出错次数(Y的效应也是0.707(=b,而智力对出错次数的直接效应是20.50
(=c′
。
智力对出错次数的总效应(=c是零(即智力与出错次数的相关系数是零。
本例涉及效应(或相关系数的遮盖(suppression问题。
由于实际中比较少见,这里不多讨论。
但从这个例子可以看出中介效应和间接效应是有区别的。
当然,如果修改中介效应的定义,不以自变量与因变量相关为前提,则另当别论。
在实际应用中,当两个变量相关不显著时,通常不再进一步讨论它们的关系了。
3 中介效应分析方法
由于中介效应是间接效应,无论变量是否涉及
潜变量,都可以用结构方程模型分析中介效应(可参考[13],中文可参考[16]。
从路径图(图1可以看
5期温忠麟等:
中介效应检验程序及其应用615
出,模型是递归的(recursive,即在路径图上直线箭头都是单向的,没有反向或循环的直线箭头,且误差之间没有弧线箭头联系。
所以,如果所有变量都是显变量,可以依次做方程(1—(3的回归分析,来替代路径分析。
就是说,如果研究的是显变量,只需要做通常的回归分析就可以估计和检验中介效应了。
无论是回归分析还是结构方程分析,用适当的统计软件都可以得到c的估计^c;a,b,c′的估计^a,^b,^c′,以及相应的标准误。
中介效应的估计是^a^b′或^c-^c′,在显变量情形并且用通常的最小二乘回归估计时,这两个估计相等[14]。
在其他情形,使用^a^b比较直观,并且它等于间接效应的估计。
除了报告中介效应的大小外,还应当报告中介效应与总效应之比(^a^b/(^c’+^a^b,或者中介效应与直接效应之比(^a^b/^c’,它们都可以衡量中介效应的相对大小[14]。
与中介效应的估计相比,中介效应的检验要复杂得多。
下面按检验的原假设分别讨论。
3.1 依次检验回归系数
在三种做法中,依次检验回归系数涉及的原假设最多,但其实是最容易的。
如果H0:
a=0被拒绝且H0:
b=0被拒绝,则中介效应显著,否则不显著。
完全中介效应还要检验H0:
c’=0。
检验统计量t等于回归系数的估计除以相应的标准误。
流行的统计软件分析结果中一般都有回归系数的估计值、标准误和t值,检验结果一目了然。
这种检验的第一类错误率很小,不会超过显著性水平,有时会远远小于显著性水平。
问题在于当中介效应较弱时,检验的功效很低。
这容易理解,如果a很小(检验结果是不显著,而b很大(检验结果是显著,因而依次检验的结果是中介效应不显著,但实际上的ab与零有实质的差异(中介效应存在,此时犯了第二类错误。
做联合检验(原假设是H0:
a=0且b=0,即同时检验a和b的显著性,功效要比依次检验的高[3]。
问题是联合检验的显著性水平与通常的不一样,做起来有点麻烦。
3.2 检验H0:
ab=0
检验H0:
ab=0的关键在于求出^a^b的标准误。
目前至少有5种以上的近似计算公式[3]。
当样本容量比较大时(如大于500,各种检验的功效差别不大。
值得在此介绍的是Sobel根据一阶Taylor展式得到的近似公式[9,10]
sab=^a2s2b+^b2s2a(5
3.3 检验H0:
c2c’=0
同样,检验H0:
c2c’=0的关键在于如何计算^c-^c’的标准误。
目前也有多种近似公式。
MacKin2non等人[3]比较的结果是其中有两个公式得到的检验有较高的功效,在总体参数a=0且b=0时的第一类错误率与0.05很接近。
一个是Clogg等人给出的公式[11]
sc2c’=rXMsc’(6其中rXM是X和M的相关系数。
另一个是Freed2man等人推出的公式[12]
sc2c’=s2c+s2c’22scsc’12r2XM(7 当a=0但b≠0时(此时ab=0,即中介效应为零,这两种公式对应的检验(即t=(^c-^c’/sc2c’作为检验统计量的第一类错误率都很高。
特别是公式(6,对应的第一类错误率有可能高达100%。
事实上,由公式(6得到的检验与H0:
b=0的检验等价[11]。
就是说,即使中介效应不存在(ab=0,只要b显著,检验结果就是中介效应显著(犯了第一类错误。
616
心 理 学 报36卷
3.4 中介效应检验方法小结
虽然检验的三种做法(即三种检验类型的原假设不一样,但它们的本质是相同的:
中介效应为零对应于下面三种情形之一:
(1a=b=0,(2仅a=0,(3仅b=0;而中介效应不为零则对应于ab≠0。
为了方便比较,表1列出了上面提到的在同类检验中较好的检验方法,其中第一类错误率和功效的比较是归纳了MacKinnon等人[3]的模拟结果。
表1 中介效应检验方法的比较
检验方法检验统计量
第一类错误率功效a=b=0仅a=0仅b=0ab≠0
依次检验
Judd和KennytN-2=^a/sa,tN-3=^b/sb,tN-3=^c/sc’
完全中介[7]
等人[4]
同上,但临界值与标准正态分布的不同约0.05较大较大很高H0:
c-c’=0
Freedman等人[12]tN2=
^c-^c’
s2c+s2c’-2scsc’1-r2XM
注:
N是样本容量。
”ν0.05”表示远远小于0.05。
考虑c’=0的情形,仅a=0时的第一类错误率与b,N有关;仅b=0时的第一类错误率与a,N有关;功效与a,b,N有关。
3.5 一个实用的中介效应检验程序
为了使一个中介效应检验的第一类错误率和第二类错误率都比较小,既可以检验部分中介效应,又可以检验完全中介效应,而且还比较容易实施,综合表1的结果,我们提出如下检验程序。
1.检验回归系数c,如果显著,继续下面的第2步。
否则停止分析。
2.做Baron和Kenny[8]部分中介检验,即依次检验系数a,b,如果都显著,意味着X对Y的影响至少有一部分是通过了中介变量M实现的,第一类错误率小于或等于0.05,继续下面第3步。
如果至少有一个不显著,由于该检验的功效较低(即第二类错误率较大,所以还不能下结论,转到第4步。
3.做Judd和Kenny[7]完全中介检验中的第三个检验(因为前两个在上一步已经完成,即检验系数c’,如果不显著,说明是完全中介过程,即X对Y的影响都是通过中介变量M实现的;如果显著,说明只是部分中介过程,即X对Y的影响只有一部分是通过中介变量M实现的。
检验结束。
4.做Sobel[9]检验,如果显著,意味着M的中介效应显著,否则中介效应不显著。
检验结束。
整个检验程序见图2。
这个程序有可能只需要依次检验,即使需要Sobel检验,用公式(5直接计算sab和检验统计量z=^a^b/sab都不算难。
如果使用LISREL进行分析,输出结果中可以找到本检验程序所需的全部检验统计量的值和检验结果。
下面看一个实际例子。
图2 中介效应检验程序
4 学生行为对同伴关系影响的中介效应分析
要研究的是初中学生行为(X对同伴关系(Y
5期温忠麟等:
中介效应检验程序及其应用617
的影响。
变量及其数据来自香港中文大学张雷教授主持的儿童同伴关系研究,本文只用到部分变量和数据。
有关的研究背景和量表及其施测方法等说明参见[18,19],这里只简单地介绍有关变量的含义和符号。
学生行为(X是被试的违纪捣乱行为,包括9个题目(如挑起争斗、欺负同学、说脏话等,同伴关系(Y是被试受同学欢迎的程度,具体地说,就是同班同学有多少人将其列入喜欢的名单(每人所列的喜欢名录没有名额限制。
老师的管教方式(U是被试对班主任老师的管教方式的评价,也有9个题目(如班主任愿意听我们的意见,班主任的期望和要求明确清晰,等等。
老师对学生的喜欢程度(W由班主任为被试打分(从“一点都不喜欢”到“非常喜欢”5级记分。
被试人数N=595。
由于潜变量和显变量的中介效应检验方法是一样的,为简单起见,这里将上述变量都作为显变量处理(即用该变量包含的题目得分的平均值作为变量值。
所有变量都已经中心化,数据分析中只需要下面的协方差矩阵:
Y 18.87
W 1.13 0.45
X -9.78 -2.20 94.25
U 0.63 0.09 -0.22 0.56
使用广义最小二乘估计方法进行分析,由于样本容量大,广义最小二乘估计与极大似然估计的结果非常接近。
4.1 教师喜欢程度的中介效应分析
假设我们认为学生行为会影响老师对他的喜欢程度,而同伴关系会受到老师喜欢程度的影响,则喜欢程度是中介变量。
喜欢程度(W的中介效应分析结果见表2,其中的结果是标准化解,用小写字母代表相应变量的标准化变量。
由于依次检验(指前面3个t检验都是显著的,所以喜欢程度的中介效应显著。
由于第四个t检验也是显著,所以是部分中介效应,中介效应占总效应的比例为01338×01349/01232=5018%。
表2 喜欢程度(W的中介效应依次检验
标准化回归方程回归系数检验
第一步y=201232xSE=0.040,t=25.833
第二步w=201338xSE=0.039,t=28.733
第三步y=0.349wSE=0.040,t=8.733
20.114xSE=0.040,t=22.833
注:
SE表示标准误。
33p<0.01。
上述包含了中介变量W的模型分析结果表明:
一方面,学生行为对同伴关系有直接负效应,即违纪捣乱行为多的同学,受同学欢迎的程度往往会低一点。
另一方面,学生行为通过教师喜欢程度对同伴关系有间接负效应,即违纪捣乱行为多的同学,老师往往比较不喜欢,而老师的态度会影响同学,使同学也比较不喜欢。
4.2 教师管教方式的中介效应分析
假设我们认为学生的行为会影响老师的管教方式,而管教方式会影响同伴关系,则管教方式是中介变量。
管教方式(U的中介效应分析结果(标准化解见表3。
由于依次检验中的第二步检验不显著(即u对x的回归系数不显著,t=20.72,p>0.05,根据我们提出的检验程序,需要做Sobel检验,检验统计量是z=^a^b/^a2s2b+^b2s2a,此处
^a=20.030,sa=0.041,^b=0.187,sb=0.039
计算得z=20.72,p>0.05。
所以管教方式(U的中介效应不显著。
表3 管教方式(U的中介效应依次检验
标准化回归方程回归系数检验
第一步y=201232xSE=01040,t=2518133
第二步u=201030xSE=01041,t=20172
第三步y=01187uSE=01039,t=417933
201226xSE=01039,t=2517933
注:
SE表示标准误。
33表示在0101水平上显著。
5 讨论和结语
511 讨论
虽然为了行文简便,本文前面只考虑一个自变量和一个中介变量的模型,但所提出的检验程序也适合有多个自变量或(和多个中介变量的模型,只是此时“完全中介”的概念没有多大意义,即不用考虑做“完全中介”检验。
下面以两个自变量X1、X2和两个中介变量M1、M2,并且都含交互效应项[20](分别是X1X2和M1M2的模型为例加以说明。
假设所有变量都已经中心化(即均值为零,类似于方程(1—(3并加以推广,可用下列方程来描述变量之间的关系:
Y=c1X1+c2X2+c3X1X2+e1
M1=a11X1+a12X2+a13X1X2+e12
M2=a21X1+a22X2+a23X1X2+e22
M1M2=a31X1+a32X2+a33X1X2+e32
618
心 理 学 报36卷
期5温忠麟等:
中介效应检验程序及其应用 619Y=c’X1+c’X2+c’X1X2+b1M1+b2M2+123b3M1M2+e3 在自变量和中介变量不止一个时,3 MacKinnonDP,LockwoodCM,HoffmanJM,WestSG,SheetsV.AComparisonofmethodstotestmediationandotherin2SocietyforPreventionResearch,ParkCity,UT.1998,Junemodels.In:
JSLong(Ed.Commonproblems/propersolutions.Statistics,1992,17(1:
51~74Wiley,1989研究者首先要明确,感兴趣的是哪个自变量经过哪个中介变量的中介效应。
然后找出哪些系数分别对应于本文提出的中介效应检验程序(图2中的系数c,a,b,c’就可以检验了(但不用检验c’,的显著性。
例如,要检验X1经过M2的中介效应是否显著,则与c,a,b,c’对应的分别是c1,a21,b2,c’;要1检验X1X2经过M1的中介效应是否显著,则与c,a,b,c’对应的分别是c3,a13,b1,c’;要检验X2经3terveningvariableeffects.PsychologicalMethods,2002,7(1:
83~1044 MacKinnonDP,LockwoodCM,HoffmanJM.Anewmethodtotestformediation.Paperpresentedattheannualmeetingofthe过M1M2的中介效应是否显著,则与c,a,b,c’对应的分别是c2,a32,b3,c’。
在实际应用中,一般不2会考虑经过M1M2的中介效应,因为解释起来比较困难。
如果画出了路径图,则容易从图上找出与c,a,b,c’对应的路径(系数并用本文提出的检验程序加以检验。
还要注意的一点是,当自变量或中介变量不止一个时,公式(4不再成立。
例如,要检验X1经过M2的中介效应是否显著,与c,a,b,c’对应的分别10 SobelME.DirectandindirecteffectsinlinearstructuralequationBeverlyHills,CA:
Sage,1988.46~64是c1,a21,b2,c’。
此时c1=c’+a21b2一般来说11是不成立的。
5.2 结语 在多变量分析中,除了考虑自变量对因变量的影响外,经常还会涉及中介变量。
例如,有间接效应的路径分析,其实已经涉及中介变量,但研究者如果不知道相应的概念和分析方法,自然不可能进行真正的中介效应分析(特别是中介效应的检验。
本文提出的中介效应检验程序,可以做部分中介效应和完全中介效应的检验。
由于同时考虑了两类错误率,该程序比单一的检验方法要好。
而且,该程序简单可行,计算量少。
该程序可以让读者避免在繁多的检验方法中无所适从,能够按部就班地进行中介效应的检验。
参考文献1 ShenJ,WangD,PengH,TangD.TheeffectsofmediatorsonicaSinica,2003,35(6:
802~809theagingofprimarymentalabilities(inChinese.ActaPsycholog2 (申继亮,王大华,彭华茂,唐丹.基本心理能力老化的中介变量.心理学报,2003,35(6:
802~809dren’academicachievements2 FangP,XiongD,GuoC.Theeffectofparentingstylesonchil2ence,2003,(1:
78~81 (方平,熊端琴,郭春彦.父母教养方式对子女学业成就影响的研究.心理科学,2003,(1:
78~81(inChinese