有调节的中介模型检验方法_竞争还是替补_温忠麟.pdf

心理学报 2014,Vol.46,No.5 714726 Acta Psychologica Sinica DOI:10.3724/SP.J.1041.2014.00714 收稿日期:2013-05-02*国家自然科学基金(31271116)、教育部人文社会科学重点研究基地项目(11JJD190005)、中国博士后科学基金项目(2013M540535)、教育部人文社会科学研究青年基金项目(13YJC190029)和江西省博士后科研择优资助项目(2013KY08)资助。通讯作者:温忠麟,E-mail: 714 有调节的中介模型检验方法：竞争还是替补?*温忠麟1 叶宝娟2(1华南师范大学心理应用研究中心/心理学院,广州 510631)(2江西师范大学心理学院,南昌 330022)摘 要 在心理和其他社科研究领域,经常遇到中介和调节变量。模型的变量多于3个时,可能同时包含中介和调节变量,一种常见的模型是有调节的中介模型。本文检视文献上各种检验有调节的中介模型的方法,理清方法之间是竞争关系(分清优劣)还是替补关系(分清先后),在此基础上总结出检验有调节的中介模型的步骤,并用一个实例进行演示。文中也讨论了有调节的中介模型与有中介的调节模型的联系与区别。关键词 中介;调节;有调节的中介;有中介的调节 分类号 B841 1 引言 在 心 理 和 行 为 研 究 中,经 常 遇 到 中 介(mediation)和调节(moderation)效应。如果自变量X通过影响第三个变量W来影响因变量Y,此时W是中介变量。例如,青少年的“感恩”影响“日常性学业复原力”,进而影响“学业成就”(叶宝娟,杨强,胡竹菁,2013)。中介效应分析的目的是探究X如何影响Y(Baron&Kenny,1986;MacKinnon,2008;Muller,Judd,&Yzerbyt,2005;Yuan&MacKinnon,2009)。如果自变量X与因变量Y的关系受到第三个变量U的作用,此时U是调节变量,影响X和Y之间关系的方向(正或负)和强弱。例如,“同伴关系”与“学生违纪行为”的关系,受到“教师管教方式”的影响(温忠麟,侯杰泰,张雷,2005)。调节效应分析的目的是探究X何时影响Y或何时影响较大(Baron&Kenny,1986;Muller et al.,2005)。如果一个模型包含的变量多于3个,可能同时包含中介变量和调节变量,这些变量在模型中的位置和作用不同会产生 不 同 的 模 型,有 调 节 的 中 介 模 型(moderated mediation model)就是同时包含中介变量和调节变量的一种常见模型,这种模型意味着自变量通过中介变量对因变量产生影响,而中介过程受到调节变量的调节(Baron&Kenny,1986;温忠麟,刘红云,侯杰泰,2012;温忠麟,张雷,侯杰泰,2006)。文献上出现了多种检验有调节的中介模型的方法,讨论不同方法时建立的有调节的中介模型也不尽相同。不同的有调节的中介模型需要不同的检验方法吗？不同方法之间是竞争关系、互补关系还是替补关系？实际应用中究竟应当从何着手分析？有调 节的 中 介模 型与 有 中介 的调 节 模型(mediated moderation model;Baron&Kenny,1986;温忠麟等,2012)有什么联系和区别？这些都是应用工作者感到疑惑的问题。本文先介绍有调节的中介模型类型,然后用一般的模型讨论文献上出现的检验方法,明晰方法之间的关系,推荐比较好的检验策略,并总结出检验流程,用一个实例演示如何用本文总结的流程检验有调节的中介模型。文中也讨论了有调节的中介模型与有中介的调节模型的联系与区别。2 有调节的中介模型类型 Edwards和Lambert(2007)在讨论既有调节又有中介的模型时,将调节放到中介分析背景中,按中介过程的前半路径、后半路径、直接路径是否受 温忠麟 等:有调节的中介模型检验方法：竞争还是替补 715 到调节,组合出8种模型,但没有进一步分类。这里只考虑有调节的中介模型,将所有路径都没有受到调节的模型(即通常的简单中介模型)作为基准模型看待,也不考虑只有直接路径受到调节的模型(因为不属于有调节的中介模型),并将模型分成两类：只调节间接效应、同时调节间接效应和直接效应。2.1 只调节间接效应 设要分析青少年感恩(X)对其学业成就(Y)的影响。研究发现,日常性学业复原力(W)是感恩影响学业成就的中介变量(图1),而压力性生活事件(U)调节了学业复原力对学业成就的影响(叶宝娟等,2013,图2)。如所知,中介效应等于中介路径上两个路径系数(回归系数)的乘积(MacKinnon,Warsi,&Dwyer,1995;温忠麟等,2012)。对于图1所示模型,中介效应为ab,是不依赖于任何变量的参数。对于图2所示模型,中介效应为a(b1+b2U),是调节变量U的函数,即中介效应受到U的调节。本节有关回归系数的含义请见第2.2节的回归方程。在上述例子中,调节变量调节了中介过程的后半路径(见图2)。但理论上,调节变量也可以调节中介过程的前半路径,如图3,中介效应为(a1+a3U)b,图1 简单中介模型 图2 调节了中介过程后半路径 图3 调节了中介过程前半路径 图4 调节了中介过程前后路径 图5 调节了中介过程后半路径和直接路径 图6 调节了中介过程前半路径和直接路径 5期 716 心 理 学 报 46卷 或者同时调节中介过程的前后两个路径,如图4,中介效应为(a1+a3U)(b1+b2U)(Muller et al.,2005;温忠麟等,2012)。显然,图2和图3所示的模型,都是图4所示模型的特例。2.2 同时调节间接效应和直接效应 前面所述的模型(图2图4),都没有考虑对直接效应的调节。如果在这些模型中调节变量还可以同时调节直接效应(Edwards&Lambert,2007),就可以分别得到如下模型：调节变量同时调节直接路径和中介过程的后半路径(图5);同时调节直接路径和中介过程的前半路径(图6);同时调节直接路径和中介过程的前后两个路径(图7)。显然图7所示的模型,是最一般的有调节的中介模型,其余模型都是其特例。为了比较模型增加了调节直接效应前后的区别,我们比较图7所示的模型和图4所示的模型。为了检验直接效应是否受到调节,在做中介效应分析之前,先建立如下回归方程(即Y对X,U和UX的回归,见图8)。01231Ycc Xc Uc UXe(1)如果c3显著,则应当考虑调节了直接效应的模型(图7)。图7 调节了中介过程前后路径和直接路径 图8 简单调节模型 无论是图7还是图4所示模型,中介效应路径上变量之间关系没有差别,以W为因变量,可以得到如下回归方程(即W对X,U和UX的回归)01232Waa Xa Ua UXe(2)写成02132()Waa Uaa U Xe可知,X对W的效应是a1+a3U。但因为直接路径上变量之间关系不同,以Y为因变量的回归方程有差别。对于图4所示的模型,以Y为因变量,可以得到如下回归方程(即Y对X,U,W和UW的回归)012123Ycc Xc UbWb UWe(3)写成012123()Ycc Xc Ubb U We可知,W对Y的效应是b1+b2U,因而X经过W对Y的中介效应为(a1+a3U)(b1+b2U),这不难将方程(2)中的W代入方程(3)进行验算得到。对于图7所示的模型,以Y为因变量,可以得到如下回归方程(即Y对X,U,UX,W和UW的回归)0123124Ycc Xc Uc UXbWb UWe(4)不难看出,X经过W对Y的中介效应还是(a1+a3U)(b1+b2U)。所以,无论直接效应是否受到调节,中介效应的代数表达式相同。2.3 一个路径图两种解释 对于图7所示的路径图,从有调节的中介模型角度,是与简单中介路径图(图1)比较,无论是X与Y关系的直接路径还是间接路径,都受到U的调节。明确一点说就是,W是X与Y关系的中介变量,中介过程的前后路径以及直接路径,都受到U的调节。但如果将图7所示路径图与简单调节路径图(图8,其中的c3显著)比较,就会发现图7可以有另一种解释：U是X与Y关系的调节变量,调节效应有一部分是经过中介变量W起作用的,间接的调节效应由(a1+a3U)(b1+b2U)中与U有关的系数反映,直接的调节效应由Ucc31中U的系数反映。如果图8中的c3显著而图7中的3c不再显著(相当于变成了图4的模型),则调节效应“完全”经过中介变量W起作用。从这个角度看,分析的是有中介的调节模型了(Baron&Kenny,1986;温忠麟等,2012;叶宝娟,温忠麟,2013)。这样,图2至图7一个路径图可以反映两种模型,与图1比较,可以分析有调节的中介效应;与图8比较(前提是其中的c3显著),可以分析有中介的调节效应。借用Muller等人(2005)的话来说就是“同一个硬币的两面”。而且,在考虑有中介的调节模型时,直接效应不受调节意味着调节效应“完全”经过中介变量W起作用。而直接效应也受到调节意味着调节效应“部分地”经过中介变量W起作用。温忠麟 等:有调节的中介模型检验方法：竞争还是替补 717 但两种模型有不同的研究目的和研究重心,立论和解释也很不同,这一点要从文献综述开始就注意区分。对于有调节的中介模型,重心在于考虑自变量对因变量的作用机制,即中介效应;其次考虑中介过程是否受到调节,即中介作用何时较强、何时较弱。显然,简单中介模型中的中介效应(见图1)相当于U=0时的中介效应,这个中介效应不要求显著。只要中介效应受到U的影响(从而至少在某个U处中介效应是显著的),就说明是有调节的中介了。对于有中介的调节模型,重心在于考虑自变量与因变量之间关系的方向(正或负)和强弱受到的影响,即调节效应;其次考虑调节变量是如何起作用的,即是否通过中介变量而起作用。这个简单调节效应(见图8)不仅要求显著,而且应当报告增加调节项前后的R 2变化,通常认为R2变化要超过2%(甚至3%)才有实质性意义,否则进一步考虑这个调节作用的机制意义不大。重要的是,两种模型分析方法和步骤不同,选择的模型要与研究目的和研究重心配合,对结果的解释则要与模型配合。值得指出的是,Muller等人(2005)以图8所示的模型为出发点,如果直接路径没有受到调节(即c3不显著),则建立有调节的中介模型;如果直接路径受到调节(即c3显著),则建立有中介的调节模型(虽然文末他们有点松口,但整篇文章倾向性很明显)。我们赞同Edwards和Lambert(2007)的观点,无论直接路径是否受到调节,都可以建立有调节的中介模型。但Edwards和Lambert将调节放到中介分析背景中,认为只是当中介过程的前半路径受到调节时(图3和图6),有调节的中介模型与有中介的调节模型从分析角度看是等价的,当中介过程的后半路径也受到调节时(图2、图4、图5和图7),两者不等价,因为在他们看来,此时模型不属于有中介的调节模型了。但其实图5和图7还是可以解释为有中介的调节模型(Muller et al.,2005;叶宝娟,温忠麟,2013)。总结一下,如果方程(1)中的c3显著,则既可以考虑建立有中介的调节模型,也可以考虑建立有调节的中介模型,并且两者的路径图还是一样的。如果方程(1)中的c3不显著,则只能考虑建立有调节的中介模型。从这个意义上说,有调节的中介模型的范围比较大。3 有调节的中介模型检验方法 为了对文献上出现的各种检验有调节的中介效应方法有更本质的了解,我们从中介过程前后路径都受到调节的一般模型出发进行讨论。前面说过,无论直接效应是否受到调节,X经过W对Y的中介效应的表达式都一样。具体来说,图7和图4所示模型的中介效应都是(a1+a3U)(b1+b2U)。显然,只要中介效应(a1+a3U)(b1+b2U)a1b1(a1b2a3b1)Ua3b2U2与U有关,或者说随U变化,则中介效应是有调节的。为简单明确起见,下面以图4所示的模型进行讨论。3.1 系数乘积的检验 为了检验(a1+a3U)(b1+b2U)与U有关,可以检验如下关于系数乘积的假设：a1b2=0,a3b1=0,a3b2=0,只要其中有一个被拒绝,中介效应是有调节的。例如,设a1b20,如果a1b2a3b10,则(a1+a3U)(b1+b2U)与U有关;如果a1b2a3b10,则a3b1=?a1b2 0,从而a3b20,(a1+a3U)(b1+b2U)也与U有关。对于系数乘积检验,文献上多有讨论。一种方法是依次检验,将一个系数乘积的检验变成两个系数的检验。例如,要检验a1b2=0,变成检验a1和b2是否显著(即检验a1=0和b2=0)。现在3个系数乘积涉及4个系数,先检验方程(2)中的a1、a3是否显著,再检验方程(3)中的b1、b2是否显著。如果10a 且20b,或者30a 且10b,或者30a 且20b,至少有一组成立,则中介效应受到调节。虽然所用的模型不尽相同,但许多作者都使用了这种依次检验方法(如Muller et al.,2005;温忠麟等,2012)。显然,如果只考虑中介效应过程的一个路径是否受到调节,则只需要检验一个系数乘积,例如,只考虑前半路径的调节(此时20b),只需要检验a3b1=0,即依次检验a3和b1是否显著。对于系数乘积的检验,已有不少研究发现,依次检验的第一类错误率较低,往往远低于设定的显著性水平(MacKinnon,Lockwood,Hoffman,West,&Sheets,2002;温忠麟,张雷,侯杰泰,刘红云,2004),说明如果依次检验结果显著,已经足够支持所要的结果,即中介效应受到调节。例如,如果检验结果是10a 且20b,足以支持a1b20。但依次检验的检验力(power)也较低,即系数乘积实际上显著,但依次检验比较容易得出不显著的结论(Fritz&MacKinnon,2007;MacKinnon et al.,2002)。直接检验系数乘积,检验力比较高。传统的系数乘积检验方法是Sobel检验(1982),检验力高于依次检验(MacKinnon et al.,2002;温忠麟等,2012),但Sobel检验基于正态性假设,而两个系数5期 718 心 理 学 报 46卷 的乘积(估计量)通常不会是正态的,所以Sobel检验备受质疑(如Hayes,2009;MacKinnon,2008)。新近的研究发现,用偏差校正的非参数百分位Bootstrap法或者有先验信息的马尔科夫链蒙特卡罗(MCMC)法计算系数乘积的置信区间比Sobel法得到的置信区间更精确,有更高的检验力(方杰,张敏强,2012;MacKinnon,2008;Preacher&Hayes,2004;Yuan&MacKinnon,2009)。以后称这种通过构造系数乘积的置信区间的检验为“系数乘积的区间检验”。不熟悉Bootstrap 法和MCMC法的读者,可以将其和最小二乘法、极大似然法那样看待,利用统计软件提供的功能实现计算。3.2 中介效应差异的检验 为了检验中介效应(a1+a3U)(b1+b2U)随U变化,有一种直观的想法是代入不同的U值,看看中介 效 应 是 否 会 变 化,这 就 是 条 件 间 接 效 应(conditional indirect effects)分析法(MacKinnon et al.,2002;Preacher,Rucker,&Hayes,2007)。只要对某个U值,中介效应(a1+a3U)(b1+b2U)显著,则存在 条 件 中 介 效 应(conditional mediation effect,Preacher et al.,2007)。如果还有另一个U值,中介效应不显著了,说明中介效应随U变化,这是所谓的亚组分析法(subgroup approach,Preacher et al.,2007)。亚组分析法可用Bootstrap法或Sobel检验进行检验(Preacher et al.,2007)。亚组分析法有多个明显的缺点(Edward&Lambert,2007)。这里只列出容易想到的两个缺点：一是即使中介效应随U变化,在某两个点处也未必反映出来;二是当调节变量U的不同水平上的条件中介效应都显著或都不显著时,不能确定中介效应是否受到调节。因为有更好的替代方法,本文不推荐使用亚组分析法。针对亚组分析法的不足,Edward和Lambert(2007)提出用Bootstrap 法检验中介效应的差异来判断中介效应是否随U变化,即检验U的不同取值上的中介效应之差是否显著,例如,取U的平均值上下一个标准差的值,分别记为UH和UL,如果(a1+a3UH)(b1+b2UH)与(a1+a3UL)(b1+b2UL)差异显著,则中介效应受到U的调节。但如果差异不显著,仅凭U的两个取值对应的中介效应进行判断,可想而知检验力不高。对连续型调节变量,刘东、张震和汪默(2012)建议,使用调节变量最大和最小观测值来检验中介效应的差异。对于只有一个中介路径受到调节的情形,这个建议是好的,因为此时中介效应是U的线性函数,在U的最大和最小观测值处,中介效应的差异最大。但对于前后路径都受到调节的情形,中介效应是U的二次函数,我们建议从中介效应(a1+a3U)(b1+b2U)出发找最值,然后检验最大 值 和 最 小 值 之 差 是 否 显 著。如 果(a1+a3U)(b1+b2U)的最大值和最小值之差显著,则中介效应受到U的调节。如果不显著,可以认定(a1+a3U)(b1+b2U)不会随U变化,即中介效应不受U的调节。二次函数2yaxbxc的图像是抛物线,在/(2)xba 处达到最小值(当a0)或者最大值(当a0)。不难计算二次函数在区间12,x x上的最大值和最小值。如果/(2)ba在区间12,x x中,二次函数最大值和最小值可以通过计算下面3个点的函数值比较后得到：12,x x和/(2)ba。如果/(2)ba不在区间12,x x中,则二次函数在区间12,x x上是单调的,在区间端点达到最大值和最小值。对于调节变量U,合理的取值区间是其均值上下两个标准差,这样就不难求得(a1+a3U)(b1+b2U)的最大值和最小值了。前面说过,系数乘积的区间检验的检验力高于Sobel检验,而后者高于依次检验。检验力最高的则是中介效应的差异检验,道理很简单,如果(a1+a3U)(b1+b2U)与U有关,则(a1+a3U)(b1+b2U)的最大值和最小值的差异一定显著。看一个简单的例子：ybxc,若系数b显著,则在1x 和0 x 的y值差异(等于b)显著,y值的最大值和最小值之差自然显著。小结一下,检验中介效应是否受到调节,可以通过依次检验或者系数乘积的区间检验,检验中介效应是否与U有关;也可以通过中介效应差异检验,检验中介效应是否随U变化。依次检验与通常的回归系数显著性检验相同,而系数乘积的区间检验和中介效应差异检验都可以使用Bootstrap法,对于系数乘积的区间检验还可以使用MCMC法。就检验力而言,系数乘积的区间检验高于依次检验,而中介效应差异检验是最高的。3.3 检验方法的比较 现在要比较上面小结的3种检验方法：依次检验、系数乘积的区间检验和中介效应差异检验,目的是理清它们之间的关系。如果是竞争关系,应当弄清谁优谁劣;如果是互补关系,应当知道什么情况下用哪种方法较好;如果是替补关系,应当说明谁先谁后。目前文献上还未见从这个角度对有调节的中介模型检验方法进行讨论。前面说过,依次检验的第一类错误率较低,如 温忠麟 等:有调节的中介模型检验方法：竞争还是替补 719 果检验结果是显著,已经可以知道中介效应受到调节。但因为依次检验的检验力也较低,所以当依次检验结果是不显著时(3对系数乘积无一显著),还不能下结论。改用检验力较高的系数乘积的区间检验,还有可能得到显著的结果。无论是依次检验还是系数乘积的区间检验,都是将整个中介效应(a1+a3U)(b1+b2U)分成了几部分进行检验的,检验结果就算是每部分都不显著,还有可能合起来是显著的。所以,即使系数乘积的区间检验结果不显著,检验中介效应的最大值与最小值之差,还有可能得到显著的结果。正如说面说过的,这3种方法的检验力是一个比一个高。既然如此,是不是应当首选检验力最高的中介效应差异检验,一次就解决问题呢？不是的,正好相反,面对有调节的中介模型,我们建议先做依次检验,如果不显著,再做系数乘积的区间检验,如果还不显著,最后做中介效应差异检验。如果前面的检验已经显著,后面的检验不需要做了,除非研究者有特别的用意。下面详述这样建议的理由。首先,从3个检验方法得到的显著性(即中介效应受到调节)结果看,依次检验的最强,中介效应差异检验的最弱。就是说,当用前一种方法检验结果显著的时候,用后一种方法检验结果几乎可以肯定也是显著的。但反之不然。第二,从显著性结果包含的信息看,依次检验的最丰富,中介效应差异检验的最少。依次检验显著的时候,10a 且20b,或者30a 且10b,或者30a 且20b,至少有一组成立,可以知道图4的中介路径上每个系数是否显著,从而知道究竟是前半路径受到调节、后半路径受到调节、还是前后路径都受到调节。写出中介效应(a1+a3U)(b1+b2U)的表达式,可以清晰反映中介效应是如何随U变化的。系数乘积的区间检验显著时,整个中介路径作为一个整体与U有关。根据显著的系数乘积可以判断哪一段中介路径受到调节,如a3b2显著,前后路径都受到调节。此时还是可以写出有意义的中介效应表达式a1b1(a1b2a3b1)Ua3b2U2。中介效应差异检验显著时,不仅不能判断哪一段中介路径受到调节,而且中介效应的表达式a1b1(a1b2a3b1)Ua3b2U2没有统计意义了(如果系数乘积都不显著)。虽然知道中介效应会随U变化,但不知道如何变化了。第三,从显著性结果的解释看,依次检验的最清晰,中介效应差异检验的最含糊。显著性结果的解释依赖显著性包含的信息,信息越多,解释越清晰。如所知,对于简单的调节效应分析(像图8那样的模型),除了报告检验结果外,通常取调节变量的均值上下一个标准差报告对应的效应大小。对于有调节的中介效应,当依次检验或系数乘积的区间检验显著时,也可以取调节变量均值上下一个标准差报告对应的中介效应值。由于此时中介效应的表达式a1b1(a1b2a3b1)Ua3b2U2是有统计意义的,这保证了上面所说的报告是有意义的。对于依次检验显著的,还可以报告调节变量均值上下一个标准差相应的前后路径的效应值。但对于中介效应差异检验,检验的是中介效应最大值和最小值的差异,所以取调节变量均值上下一个标准差报告对应的中介效应值还不知道是否有意义,还应当检验其差异是否显著。其实,差异显著只是说差异在统计上是可以分辨的,结果很弱。当依次检验或系数乘积的区间检验显著时,得到的中介效应与调节变量的关系表达式比上述差异显著性重要得多,此时通常不用关注调节变量的均值上下一个标准差对应的中介效应差异是否显著,就像回归分析中得到显著的回归方程后,通常不用关注自变量的均值上下一个标准差对应的因变量差异是否显著。第四,从检验的难易程度和方便性看,依次检验最简单方便,中介效应之差最麻烦。做依次检验,使用SPSS(显变量情形)或者LISREL或Mplus等软件(潜变量情形)就能解决问题。系数检验与通常的回归系数检验相同。系数乘积的区间检验和中介效应差异检验,用SPSS或者LISREL比较麻烦,不过Mplus有现成的指令计算置信区间(见后面例子)。综上所述,3种检验方法不是竞争关系和互补关系,而是替补关系。研究者当然会希望用比较简便的方法,得到一个信息丰富、容易解释的比较强的检验结果。如果通过依次检验已经得出“中介效应受到调节”的结果,除非特别需要,可以不用后面的检验方法。例如,在需要比较某两个U值对应的中介效应时,才检验中介效应之差是否显著。当需要计算系数乘积的置信区间时,Sobel检验和Bootstrap 法、MCMC法是竞争关系,后两者胜出。如果没有先验信息,使用偏差校正的非参数百分位Bootstrap 法;如果有先验信息,使用MCMC法(方杰,张敏强,2012)。从现有研究结果看,区间检验的Bootstrap法或MCMC法,是两种可选的相互替代方法。在Mplus中只需要添加几个额外参数就很容易用这两种方法得到所要的置信区间5期 720 心 理 学 报 46卷 (见后面例子)。至于条件间接效应分析,中介效应的亚组分析法与差异分析法是竞争关系,上一节已经说过是后者胜出,所以前面我们用差异分析法与依次检验、系数乘积的区间检验进行比较。诚然,各种检验方法都有其优点和缺点,我们在充分了解文献上对已有方法讨论的基础上,作出上述的比较和判断。4 有调节的中介模型分析步骤 我们自然会感兴趣自变量X对因变量Y的影响是否显著,同时考虑直接效应是否受到调节,所以第一步是检验回归方程(1)的系数c1和c3。检验c3可以知道在未考虑中介效应的时候,直接效应是否受到调节。检验系数c1的结果,对将来解释模型有用。但即使c1不显著,都可以继续后面的分析,只是解释结果的时候略有不同,参见温忠麟等(2012)说的广义中介分析。第二步是检验(a1+a3U)(b1+b2U)是否与U有关。做依次检验,先检验方程(2)中的a1、a3是否显著,再检验方程(3)(直接效应没有受到调节)或者方程(4)(直接效应受到调节)中的b1、b2是否显著。如果10a 且20b(调节后半路径),或者30a 且10b(调节前半路径),或者30a 且20b(调节前后路径),至少有一组成立,则中介效应受到调节。将不显著的a3或 b2固定为零,重新估计其他系数。报告结果时,写出中介效应(a1+a3U)(b1+b2U)(不显著的a3或 b2变成零),并报告其在U的均值以及均值上下一个标准差处的中介效应值,这样就对中介效应受到调节的情况有比较清晰的了解。如果是纯粹验证模型,则不论系数是否显著,都保留在模型中,无需重新估计。如果依次检验结果是显著的,已经可以知道中介效应受到调节。如果依次检验结果不显著(3对系数乘积无一显著),还不能下结论。下一步是使用非参数百分位Bootstrap 法或者MCMC法对系数乘积做区间检验。如果至少有一对乘积显著,则中介效应受到调节。如果U2的系数a3b2显著,则前后路径都受到调节,否则将调节不显著的路径改成没有调节,重新估计和检验模型。如果是纯粹验证模型,则不论系数(乘积)是否显著,都保留在模型中,无需重新估计。报告结果时,写出中介效应的表达式a1b1(a1b2a3b1)Ua3b2U2来,并报告其在U的均值以及均值上下一个标准差处的中介效应值。如果系数乘积的区间检验结果还不显著,最后一步是检验中介效应的最大值与最小值之差。中介效应(a1+a3U)(b1+b2U)是U的二次函数,在有了系数估计值后,可以在U值正常区间(例如均值上下两个标准差)内计算二次函数最大值和最小值,使用Bootstrap法检验它们的差异。如果(a1+a3U)(b1+b2U)的最大值和最小值之差显著,则中介效应受到调节,并做出报告。否则,连最大值与最小值之差都不显著,已经可以认定(a1+a3U)(b1+b2U)不随U变化,即中介效应不受U的调节。综上所述,可以给出有调节的中介模型的层次检验流程(图9)。这个流程是针对前半路径和后半路径都受到调节的中介模型,如果理论假设是只有前半或后半路径受到调节,相应的回归方程和检验略有调整。如果只有前半路径受到调节,则方程(3)或(4)中不包含调节变量和中介变量的交互乘积项UW(即没有b2);如果只有后半路径受到调节,则方程(2)中不包含自变量和调节变量的交互乘积项UX(即没有a3)。为了说明使用上面的层次检验流程与已有的检验方法的差别以及层次检验的优点,不妨设想有4个研究者用同一组数据检验同一个有调节的中介模型,甲使用传统的依次检验,乙使用系数乘积的区间检验,丙使用中介效应差异检验,而丁使用本文提出的层次检验。根据3.2节的讨论可知,就甲乙丙而言,甲的检验力最低,乙的次之,而丙的检验力最高。因此,如果甲的结果显著,则乙的结果也显著(不排除出现例外的可能,但现有的模拟和实际数据都未发现例外);如果乙的结果显著,则丙的结果也显著(从3.2节可知,这里不会有例外)。因而,可以推知甲乙丙的检验结果只有表1所列的4种情形,使用层次检验的丁的检验结果可以根据甲乙丙的检验结果而确定。情形1：甲的结果显著,因而乙和丙的结果也显著;丁做了甲的检验就停止了,和甲的结果相同。情形2：甲的结果不显著,但乙的结果显著,因而丙的结果也显著;丁首先得到甲的结果(不显著),然后得到乙的结果(显著)。情形3：甲乙的结果都不显著,但丙的结果显著;丁首先得到甲的结果(不显著),然后得到乙的结果(不显著),最后得到丙的结果(显著)。情形4：甲乙丙的结果都不显著;丁首先得到甲的结果(不显著),然后得到乙的结果(不显著),最后得到丙的结果(不显著)。温忠麟 等:有调节的中介模型检验方法：竞争还是替补 721 图9 有调节的中介模型层次检验流程 丁的检验力和丙的一样,高于甲和乙(情形3);在各种情形,丁获得的信息是最多的,在情形2和3中都多于其他3人;在检验结果解释方面,丁是最强的一个,在情形1-3中,其解释都与甲乙丙中显著性解释最强的一个相同,详见表1。当第一步的检验结果是3c显著,此时可以考虑两种模型,一种是有调节的中介模型,并且调节变量同时调节了直接效应,另一种是有中介的调节模型。采用哪种模型检验并进行解释,取决于所研究的问题背景、学科理论和经验常识。如果是前者,使用本文总结的方法和步骤进行分析,如果是后者,可参考叶宝娟和温忠麟(2013)对有中介的调节模型分析方法和步骤,这里不赘述。5 示例 接下来用一个实际例子演示如何用图9的流程检验有调节的中介模型。本例感兴趣的是感恩(自变量X)通过学业复原力(W)影响学业成就(Y)的中介过程,是否受到了压力性生活事件(U)的调节,即压力性生活事件调节了中介过程的前半路径和后半路径,模型如图4所示。本次研究共有8所初级中学的1148名青少年参与调查,其中女生567人,男生581人,平均年龄15.13(SD=1.65)。有关测量工具、模型假设成立条件的评估等具体细节参见叶宝娟等(2013)的研究。但本实例数据与叶宝娟等(2013)研究中使用的数据不同,这里只作本文的建模示例,不涉及实证研究的问卷信度等其他议题。叶宝娟等(2013)仅检验了压力性生活事件调节了中介过程的后半路径的情况,本研究还感兴趣于压力性生活事件是否调节了中介过程的前半路径。将Y,W,X和U标准化变成Z分数(变量名称不变),然后将相应的Z分数相乘产生交互作用项UX和UW的分数(但UX和UW不是标准化变量)。变量间的协方差矩阵见附录1。用路径分析建模检验。第一步,建立学业成就(Y)与感恩(X)关系的简单调节模型(见图8),检验直接效应是否受到生活事件(U)的调节,结果见图10。感恩(X)对学业成就(Y)的效应显著(1c=0.156,t=5.334,p 0.001),感恩(X)与生活事件(U)的交互项(UX)对学业成就(Y)的效应不显著(c3=0.036,t=1.298,p=0.194)。第二步,建立有调节的中介模型(见图4,直接效应不受调节),检验感恩(X)经过复原力(W)对学 5期 722 心 理 学 报 46卷 表1 层次检验与依次检验、系数乘积的区间检验和中介效应差异检验比较 检验方法 情形1 情形2 情形3 情形4 甲：依次检验3对系数 检验结果显著 检验结果不显著 检验结果不显著 检验结果不显著 信息：知道乙丙丁的结果也显著 信息：不知乙丙丁的 结果 信息：与情形2相同 信息：与情形2相同 解释：明确知道中介路径上每个系数的显著性,以及前、后路径哪个受到调节 解释：中介效应没有受到调节 解释：与情形2相同 解释：与情形2相同 乙：3个系数乘积的区间检验 检验结果显著 检验结果显著 检验结果不显著 检验结果不显著 信息：知道丙丁的结果也显著,但不知甲的结果 信息：与情形1相同 信息：知道甲的结果不显著,但不知丙丁的结果 信息：与情形3相同 解释：明确知道前、后路径哪个受到调节 解释：与情形1相同 解释：中介效应没有受到调节 解释：与情形3相同 丙：中介效应差异检验 检验结果显著 检验结果显著 检验结果显著 检验结果不显著 信息：知道丁的结果也显著,但不知甲乙的结果 信息：与情形1相同 信息：与情形1相同 信息：知道其他3人的结果也不显著 解释：只知道中介过程受到调节,但不知哪个路径受到调节 解释：与情形1相同 解释：与情形1相同 解释：中介效应没有受到调节 丁：层次检验 检验结果显著 检验结果显著 检验结果显著 检验结果不显著 信息：知道甲乙丙的结果也显著 信息：知道甲的结果不显著、乙丙的结果显著信息：知道甲乙的结果不显著、丙的结果显著 信息：知道其他3人的结果也不显著 解释：与甲的情形1 相同 解释：与乙的情形1 相同 解释：与丙的情形1 相同 解释：中介效应没有受到调节 层次检验的优点 丁和甲的信息多于乙和丙,丁和甲的结果解释性也强于乙和丙 丁乙丙的检验力高于甲;丁的信息多于甲乙丙,丁和乙的结果解释性强于丙 丁和丙的检验力高于甲和乙;丁的信息多于甲乙丙 丁和丙的信息多于甲和乙 图10 检验直接效应是否受到调节 业成就(Y)的中介效应是否受到生活事件(U)的调节。先做依次检验,可以通过两个回归分析(2)和(3)完成,但我们用路径分析进行,附录2给出了Mplus程序,结果见图11。感恩(X)对复原力(W)的效应显著(a1=0.529,t=20.937,p 0.001),感恩(X)与生活事件(U)的交互项(UX)对复原力(W)的效应不显著(a3=?0.009,t=?0.386,p=0.699);复原力(W)对学业成就(Y)的效应显著(b1=0.222,t=6.646,p 0.001),生活事件(U)与复原力(W)的交互 项(UW)对 学 业 成 就(Y)的 效 应 不 显 著(b2=?0.048,t=?1.807,p=0.071)。对照图9的流程,依次检验还 图11 检验中介效应的前后路径是否受到调节 不能确定感恩经过复原力对学业成就的中介效应是否受到调节。接下来,使用偏差校正的百分位Bootstrap 法计算1 2a b、3 1a b和3 2a b的置信区间,附录3给出了Mplus程序。系数乘积95%的置信区间为：1 2a b为?0.053,?0.001,3 1a b为?0.016,0.012,3 2a b为?0.002,0.005。因为1 2a b的置信区间不包括0,所以,感恩经过复原力对学业成就的中介效应的后半路径受到生活事件的调节。因为3 1a b和3 2a b的置信 温忠麟 等:有调节的中介模型检验方法：竞争还是替补 723 区间包括0,所以,模型的前半路径没有受到调节,将前半路径设置成没有调节,重新估计和检验模型(如图2,不难对附录3的Mplus程序作相应修改)。因为a3=?0.009几乎等于零,所以设置模型前半路径不受调节后各参数的估计没有什么变化。感恩(X)经过复原力(W)对学业成就(Y)的中介效应为a1b1a1b2U=0.1170.025U,当U取值为1、0和1时,中介效应分别为0.092、0.117和0.142。感恩(X)对学业成就(Y)的总效应是0.181(见附录1),当U取值为1、0和1时,中介效应分别占了总效应的50.8、64.6和78.5。对U取值为1和0检 验 中 介 效 应 之 差,95 置 信 区 间 为0.053,0.001,差异显著。显然,U取值为1和?1对应的中介效应之差更大,可见差异检验显著其实是很弱的结果。上述统计结果说明,感恩(X)经过复原力(W)对学业成就(Y)的中介效应受到生活事件(U)的调节。当个体经历的生活事件属于平均水平时(U=0),感恩(X)经过复原力(W)对学业成就(Y)的间接效应占了总效应的2/3左右。当个体经历较多的生活事件时,这个间接效应变小,直接效应变大;当个体经历较少的生活事件时,这个间接效应变大,直接效应变小。从而,对于经历较多生活事件的个体,提高其感恩水平往往可以直接提高其学业成