Mplus model71 模型讲解

来自图书《MPlus中介调节模型》

复杂中介调节模型Mplus教程结构

• 理论模型：阐述研究模型背后的理论基础，明确变量之间的关系假设。
• 数学模型：将理论模型转化为具体的数学公式，描述变量之间的线性关系和交互作用。
• 数学推导：详细展示如何从数学模型推导出直接效应、间接效应、总效应以及条件效应的计算方法。
• 代码解读：逐行解释Mplus代码的含义，指导如何构建模型并进行结果分析。

大家好，今天我们要介绍的是一个复杂中介调节模型Mplus教程。这个教程旨在帮助大家理解并应用包含多个中介变量和调节变量的复杂模型。我们的教程结构清晰明了，分为以下四个部分：首先，我们将从理论模型入手，讲解该模型背后的逻辑和变量之间的关系。接着，我们会将理论模型转化为数学模型，用公式来描述变量之间的具体关系。然后，我们会深入探讨数学推导过程，展示如何计算直接效应、间接效应和总效应，特别是当调节变量影响这些效应时。最后，我们将逐行解读Mplus代码，让大家能够自己构建模型并分析结果。希望通过这个教程，大家可以掌握复杂中介调节模型的分析方法。

理论模型

我们来看这张图，它展示了一个比较复杂的模型。在这个模型里，X会影响Y，而它影响Y的方式是通过M。同时，W和V这两个变量，会影响X如何影响M，也会影响M如何影响Y。此外，W和V之间还存在相互影响，并且它们与X的相互作用也会影响M和Y。这个模型考虑了各种变量之间的复杂关系，帮助我们全面理解它们是如何相互作用的。

数学模型

这张图展示了一个统计模型，X、W和V这三个变量会影响变量M。同时，X、W、V以及它们两两组合的XW、XV、WV，还有三个变量组合的XWV，都会直接影响变量Y。此外，M本身，以及M与W、V、W和V组合形成的MW、MV、MWV，也会影响Y。图中的箭头代表变量之间的影响关系，例如a1、a2和a3表示X、W和XW对M的影响程度，b1、b2、b3和b4表示M、MW、MV和MWV对Y的影响程度，而c1’到c7’则表示X、W、V、XW、XV、WV和XWV对Y的直接影响程度。

数学公式1

Y = b0 + b1M + b2MW + b3MV + b4MWV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
M = a0 + a1X + a2W + a3XW

首先，我们有两个关键的等式。第一个等式表示因变量Y受到自变量X，中介变量M以及调节变量W和V的影响。第二个等式表示中介变量M受到自变量X和调节变量W的影响。我们的目标是将M的表达式代入Y的表达式，从而分析X对Y的直接和间接影响。

数学公式2

Y = b0 + b1M + b2MW + b3MV + b4MWV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
M = a0 + a1X + a2W + a3XW
Y = b0 + b1(a0 + a1X + a2W + a3XW) + b2(a0 + a1X + a2W + a3XW)W + b3(a0 + a1X + a2W + a3XW)V + b4(a0 + a1X + a2W + a3XW)WV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV

接下来，我们将M的表达式，即a0 + a1X + a2W + a3XW，代入到Y的表达式中。这样，我们就得到了一个只包含X, W, V的Y的表达式，为后续分析奠定基础。

数学公式3

Y = b0 + b1M + b2MW + b3MV + b4MWV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
M = a0 + a1X + a2W + a3XW
Y = b0 + b1(a0 + a1X + a2W + a3XW) + b2(a0 + a1X + a2W + a3XW)W + b3(a0 + a1X + a2W + a3XW)V + b4(a0 + a1X + a2W + a3XW)WV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1XW + a0b2W + a1b2XW + a2b2WW + a3b2XWW + a0b3V + a1b3XV + a2b3WV + a3b3XWV + a0b4WV + a1b4XWV + a2b4WWV + a3b4XWWV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV

现在，我们将上一步得到的表达式中的括号全部展开。这一步的目标是将所有项都明确地列出来，方便后续进行整理和分组。

数学公式4

Y = b0 + b1M + b2MW + b3MV + b4MWV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
M = a0 + a1X + a2W + a3XW
Y = b0 + b1(a0 + a1X + a2W + a3XW) + b2(a0 + a1X + a2W + a3XW)W + b3(a0 + a1X + a2W + a3XW)V + b4(a0 + a1X + a2W + a3XW)WV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1XW + a0b2W + a1b2XW + a2b2WW + a3b2XWW + a0b3V + a1b3XV + a2b3WV + a3b3XWV + a0b4WV + a1b4XWV + a2b4WWV + a3b4XWWV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
Y = (b0 + a0b1 + a2b1W + a0b2W + a2b2WW + a0b3V + a2b3WV + a0b4WV + a2b4WWV + c2'W + c3'V + c6'WV) + (a1b1 + a3b1W + a1b2W + a3b2WW + a1b3V + a3b3WV + a1b4WV + a3b4WWV + c1' + c4'W + c5'V + c7'WV)X

我们将上一步得到的表达式按照是否包含X进行分组。这样，我们就得到了一个类似 Y = a + bX 的形式，其中 a 是截距项， b 是 X 的系数。这使得我们能够清晰地看到X对Y的影响。

数学公式5

Y = b0 + b1M + b2MW + b3MV + b4MWV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
M = a0 + a1X + a2W + a3XW
Y = b0 + b1(a0 + a1X + a2W + a3XW) + b2(a0 + a1X + a2W + a3XW)W + b3(a0 + a1X + a2W + a3XW)V + b4(a0 + a1X + a2W + a3XW)WV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1XW + a0b2W + a1b2XW + a2b2WW + a3b2XWW + a0b3V + a1b3XV + a2b3WV + a3b3XWV + a0b4WV + a1b4XWV + a2b4WWV + a3b4XWWV + c1'X + c2'W + c3'V + c4'XW + c5'XV + c6'WV + c7'XWV
Y = (b0 + a0b1 + a2b1W + a0b2W + a2b2WW + a0b3V + a2b3WV + a0b4WV + a2b4WWV + c2'W + c3'V + c6'WV) + (a1b1 + a3b1W + a1b2W + a3b2WW + a1b3V + a3b3WV + a1b4WV + a3b4WWV + c1' + c4'W + c5'V + c7'WV)X
间接效应: (a1 + a3W)(b1 + b2W + b3V + b4WV)
直接效应: c1' + c4'W + c5'V + c7'WV

最后，我们从整理后的表达式中提取出X对Y的间接效应和直接效应。间接效应表示X通过M影响Y的路径，而直接效应表示X直接影响Y的路径。这些公式可以帮助我们理解在不同的W和V取值下，X是如何影响Y的。

代码解读1

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;

第一步，我们使用USEVARIABLES命令来声明模型中要用到的所有变量。X是预测变量，M是中介变量，W和V是调节变量，Y是结果变量。剩下的XW, XV等等则是通过原始变量相乘得到的交互项。

代码解读2

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above

这里是注释部分，提示我们接下来要创建交互项，并且这些交互项的定义必须放在USEVARIABLES命令之后。

代码解读3

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:

第三步，使用DEFINE命令，表示接下来要定义新的变量或进行计算。

代码解读4

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;

第四步，我们定义了所有需要的交互项。例如，MW是中介变量M和调节变量W的乘积。这些交互项在模型中用来检验调节效应。

代码解读5

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:

第五步，使用ANALYSIS命令，开始定义分析的设置。

代码解读6

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;

第六步，设定分析类型为一般模型。

代码解读7

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;

第七步，设定使用最大似然估计方法进行参数估计。

代码解读8

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;

第八步，设定使用10000次自助抽样，以获得更稳健的标准误和置信区间。

代码解读9

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses

这里是一个注释，提醒我们接下来在模型语句中，要使用括号给每个路径和截距命名。

代码解读10

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:

第十步，使用MODEL命令，开始定义模型本身。

代码解读11

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);

第十一步，定义结果变量Y的回归方程。Y表示Y的截距，命名为b0。Y ON M (b1)表示Y对M的回归系数，命名为b1。其他类似，这些系数代表了在考虑其他变量的情况下，每个预测变量对Y的效应。

代码解读12

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);

第十二步，定义中介变量M的回归方程。M表示M的截距，命名为a0。M ON X (a1)表示M对X的回归系数，命名为a1。其他类似。

代码解读13

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects

这里是一个注释，说明接下来我们将使用MODEL CONSTRAINT命令来检验条件间接效应。

代码解读14

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations

这里是注释，提醒我们需要为调节变量W和V选择低，中，高三个值，比如平均值加减一个标准差。两个调节变量，每个取三个值，一共会产生9种组合。

代码解读15

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.

这里是注释，解释下面要使用的条件间接效应和总效应的命名规则，比如MEV_LOQ表示V为中值而Q为低值。

代码解读16

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:

第十六步，使用MODEL CONSTRAINT命令，表示接下来定义模型约束，用于计算条件效应。

代码解读17

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);

第十七步，使用NEW命令，定义模型中需要用到的新变量，包括W和V的低中高值，以及各种条件间接效应（I开头），条件直接效应（D开头）和条件总效应（T开头）的变量名。

代码解读18

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V

第十八步，我们将实际的低，中，高值分配给变量LOW_W, MED_W, HIGH_W, LOW_V, MED_V, HIGH_V，这些值需要用户根据实际数据替换#LOWW, #MEDW等等。

代码解读19

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;

第十九步，根据前面定义的回归系数(a1, a3, b1, b2, b3, b4)和调节变量的值(LOW_W, MED_W, HIGH_W, LOW_V, MED_V, HIGH_V)，计算不同调节变量值组合下的条件间接效应。例如，ILOW_LOV表示W和V都取低值时的间接效应。

代码解读20

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;

第二十步，根据回归系数(cdash1, cdash4, cdash5, cdash7)和调节变量的值，计算不同调节变量值组合下的条件直接效应。例如，DLOW_LOV表示W和V都取低值时的直接效应。

代码解读21

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;

第二十一步，计算不同调节变量值组合下的条件总效应，通过将相应的条件间接效应和条件直接效应相加得到。例如，TLOW_LOV表示W和V都取低值时的总效应。

代码解读22

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;
! Use loop plot to plot conditional indirect effect of X on Y for each combination of low, med, high moderator values
! Could be edited to show conditional direct or conditional total effects instead
! NOTE - values of 1,5 in LOOP() statement need to be replaced by
! logical min and max limits of predictor X used in analysis

这里是注释，说明接下来这部分代码用于绘制条件间接效应图，并且可以修改为绘制条件直接或总效应图。同时提醒用户，需要根据实际数据，替换LOOP函数中的最小值和最大值。

代码解读23

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;
! Use loop plot to plot conditional indirect effect of X on Y for each combination of low, med, high moderator values
! Could be edited to show conditional direct or conditional total effects instead
! NOTE - values of 1,5 in LOOP() statement need to be replaced by
! logical min and max limits of predictor X used in analysis
PLOT(PLOW_LOV PMEW_LOV PHIW_LOV PLOW_MEV PMEW_MEV PHIW_MEV
PLOW_HIV PMEW_HIV PHIW_HIV);

第二十三步，定义要绘制的变量，这些变量将在接下来的LOOP命令中计算。

代码解读24

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;
! Use loop plot to plot conditional indirect effect of X on Y for each combination of low, med, high moderator values
! Could be edited to show conditional direct or conditional total effects instead
! NOTE - values of 1,5 in LOOP() statement need to be replaced by
! logical min and max limits of predictor X used in analysis
PLOT(PLOW_LOV PMEW_LOV PHIW_LOV PLOW_MEV PMEW_MEV PHIW_MEV
PLOW_HIV PMEW_HIV PHIW_HIV);
LOOP(XVAL,1,5,0.1);

第二十四步，使用LOOP命令创建一个名为XVAL的变量，它从1增加到5，步长为0.1。该循环用于在预测变量X的不同值上计算条件效应。

代码解读25

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;
! Use loop plot to plot conditional indirect effect of X on Y for each combination of low, med, high moderator values
! Could be edited to show conditional direct or conditional total effects instead
! NOTE - values of 1,5 in LOOP() statement need to be replaced by
! logical min and max limits of predictor X used in analysis
PLOT(PLOW_LOV PMEW_LOV PHIW_LOV PLOW_MEV PMEW_MEV PHIW_MEV
PLOW_HIV PMEW_HIV PHIW_HIV);
LOOP(XVAL,1,5,0.1);
PLOW_LOV = ILOW_LOV*XVAL;
PMEW_LOV = IMEW_LOV*XVAL;
PHIW_LOV = IHIW_LOV*XVAL;
PLOW_MEV = ILOW_MEV*XVAL;
PMEW_MEV = IMEW_MEV*XVAL;
PHIW_MEV = IHIW_MEV*XVAL;
PLOW_HIV = ILOW_HIV*XVAL;
PMEW_HIV = IMEW_HIV*XVAL;
PHIW_HIV = IHIW_HIV*XVAL;

第二十五步，这些代码计算在循环中的XVAL的每个值下，不同调节变量组合的条件间接效应。这些值将用于绘图。

代码解读26

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;
! Use loop plot to plot conditional indirect effect of X on Y for each combination of low, med, high moderator values
! Could be edited to show conditional direct or conditional total effects instead
! NOTE - values of 1,5 in LOOP() statement need to be replaced by
! logical min and max limits of predictor X used in analysis
PLOT(PLOW_LOV PMEW_LOV PHIW_LOV PLOW_MEV PMEW_MEV PHIW_MEV
PLOW_HIV PMEW_HIV PHIW_HIV);
LOOP(XVAL,1,5,0.1);
PLOW_LOV = ILOW_LOV*XVAL;
PMEW_LOV = IMEW_LOV*XVAL;
PHIW_LOV = IHIW_LOV*XVAL;
PLOW_MEV = ILOW_MEV*XVAL;
PMEW_MEV = IMEW_MEV*XVAL;
PHIW_MEV = IHIW_MEV*XVAL;
PLOW_HIV = ILOW_HIV*XVAL;
PMEW_HIV = IMEW_HIV*XVAL;
PHIW_HIV = IHIW_HIV*XVAL;
PLOT:

第二十六步，使用PLOT命令开始定义绘图设置。

代码解读27

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;
! Use loop plot to plot conditional indirect effect of X on Y for each combination of low, med, high moderator values
! Could be edited to show conditional direct or conditional total effects instead
! NOTE - values of 1,5 in LOOP() statement need to be replaced by
! logical min and max limits of predictor X used in analysis
PLOT(PLOW_LOV PMEW_LOV PHIW_LOV PLOW_MEV PMEW_MEV PHIW_MEV
PLOW_HIV PMEW_HIV PHIW_HIV);
LOOP(XVAL,1,5,0.1);
PLOW_LOV = ILOW_LOV*XVAL;
PMEW_LOV = IMEW_LOV*XVAL;
PHIW_LOV = IHIW_LOV*XVAL;
PLOW_MEV = ILOW_MEV*XVAL;
PMEW_MEV = IMEW_MEV*XVAL;
PHIW_MEV = IHIW_MEV*XVAL;
PLOW_HIV = ILOW_HIV*XVAL;
PMEW_HIV = IMEW_HIV*XVAL;
PHIW_HIV = IHIW_HIV*XVAL;
PLOT:
TYPE = plot2;

第二十七步，指定绘制二维图。

代码解读28

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;
! Use loop plot to plot conditional indirect effect of X on Y for each combination of low, med, high moderator values
! Could be edited to show conditional direct or conditional total effects instead
! NOTE - values of 1,5 in LOOP() statement need to be replaced by
! logical min and max limits of predictor X used in analysis
PLOT(PLOW_LOV PMEW_LOV PHIW_LOV PLOW_MEV PMEW_MEV PHIW_MEV
PLOW_HIV PMEW_HIV PHIW_HIV);
LOOP(XVAL,1,5,0.1);
PLOW_LOV = ILOW_LOV*XVAL;
PMEW_LOV = IMEW_LOV*XVAL;
PHIW_LOV = IHIW_LOV*XVAL;
PLOW_MEV = ILOW_MEV*XVAL;
PMEW_MEV = IMEW_MEV*XVAL;
PHIW_MEV = IHIW_MEV*XVAL;
PLOW_HIV = ILOW_HIV*XVAL;
PMEW_HIV = IMEW_HIV*XVAL;
PHIW_HIV = IHIW_HIV*XVAL;
PLOT:
TYPE = plot2;
OUTPUT:

第二十八步，使用OUTPUT命令开始定义输出设置。

代码解读29

USEVARIABLES = X M W V Y XW XV WV MW MV XWV MWV;
! Create interaction terms
! Note that they have to be placed at end of USEVARIABLES subcommand above
DEFINE:
MW = M*W;
MV = M*V;
XW = X*W;
XV = X*V;
WV = W*V;
MWV = M*W*V;
XWV = X*W*V;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
! In model statement name each path and intercept using parentheses
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MW (b2);
Y ON MV (b3);
Y ON MWV (b4);
Y ON X(cdash1);
Y ON W (cdash2);
Y ON V (cdash3);
Y ON XW (cdash4);
Y ON XV (cdash5);
Y ON WV (cdash6);
Y ON XWV (cdash7);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON XW (a3);
! Use model constraint subcommand to test conditional indirect effects
! You need to pick low, medium and high moderator values for W, V
! for example, of 1 SD below mean, mean, 1 SD above mean
! 2 moderators, 3 values for each, gives 9 combinations
! arbitrary naming convention for conditional indirect and total effects used below:
! MEV_LOQ = medium value of V and low value of Q, etc.
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_V MED_V HIGH_V
ILOW_LOV IMEW_LOV IHIW_LOV ILOW_MEV IMEW_MEV IHIW_MEV
ILOW_HIV IMEW_HIV IHIW_HIV
DLOW_LOV DMEW_LOV DHIW_LOV DLOW_MEV DMEW_MEV DHIW_MEV
DLOW_HIV DMEW_HIV DHIW_HIV
TLOW_LOV TMEW_LOV THIW_LOV TLOW_MEV TMEW_MEV THIW_MEV
TLOW_HIV TMEW_HIV THIW_HIV);
LOW_W = #LOWW;
! replace #LOWW in the code with your chosen low value of W

MED_W = #MEDW;
! replace #MEDW in the code with your chosen medium value of W

HIGH_W = #HIGHW;
! replace #HIGHW in the code with your chosen high value of W
LOW_V = #LOWV;
! replace #LOWV in the code with your chosen low value of V

MED_V = #MEDV;
! replace #MEDV in the code with your chosen medium value of V

HIGH_V = #HIGHV;
! replace #HIGHV in the code with your chosen high value of V
! Calc conditional indirect effects for each combination of moderator values
ILOW_LOV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*LOW_V + a3*b3*LOW_W*LOW_V + a1*b4*LOW_W*LOW_V +
a3*b4*LOW_W*LOW_W*LOW_V;
IMEW_LOV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*LOW_V + a3*b3*MED_W*LOW_V + a1*b4*MED_W*LOW_V +
a3*b4*MED_W*MED_W*LOW_V;
IHIW_LOV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*LOW_V + a3*b3*HIGH_W*LOW_V + a1*b4*HIGH_W*LOW_V +
a3*b4*HIGH_W*HIGH_W*LOW_V;
ILOW_MEV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*MED_V + a3*b3*LOW_W*MED_V + a1*b4*LOW_W*MED_V +
a3*b4*LOW_W*LOW_W*MED_V;
IMEW_MEV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*MED_V + a3*b3*MED_W*MED_V + a1*b4*MED_W*MED_V +
a3*b4*MED_W*MED_W*MED_V;
IHIW_MEV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*MED_V + a3*b3*HIGH_W*MED_V + a1*b4*HIGH_W*MED_V +
a3*b4*HIGH_W*HIGH_W*MED_V;
ILOW_HIV = a1*b1 + a3*b1*LOW_W + a1*b2*LOW_W + a3*b2*LOW_W*LOW_W +
a1*b3*HIGH_V + a3*b3*LOW_W*HIGH_V + a1*b4*LOW_W*HIGH_V +
a3*b4*LOW_W*LOW_W*HIGH_V;
IMEW_HIV = a1*b1 + a3*b1*MED_W + a1*b2*MED_W + a3*b2*MED_W*MED_W +
a1*b3*HIGH_V + a3*b3*MED_W*HIGH_V + a1*b4*MED_W*HIGH_V +
a3*b4*MED_W*MED_W*HIGH_V;
IHIW_HIV = a1*b1 + a3*b1*HIGH_W + a1*b2*HIGH_W + a3*b2*HIGH_W*HIGH_W +
a1*b3*HIGH_V + a3*b3*HIGH_W*HIGH_V + a1*b4*HIGH_W*HIGH_V +
a3*b4*HIGH_W*HIGH_W*HIGH_V;
! Calc conditional direct effects for each combination of moderator values
DLOW_LOV = cdash1 + cdash4*LOW_W + cdash5*LOW_V + cdash7*LOW_W*LOW_V;
DMEW_LOV = cdash1 + cdash4*MED_W + cdash5*LOW_V + cdash7*MED_W*LOW_V;
DHIW_LOV = cdash1 + cdash4*HIGH_W + cdash5*LOW_V + cdash7*HIGH_W*LOW_V;
DLOW_MEV = cdash1 + cdash4*LOW_W + cdash5*MED_V + cdash7*LOW_W*MED_V;
DMEW_MEV = cdash1 + cdash4*MED_W + cdash5*MED_V + cdash7*MED_W*MED_V;
DHIW_MEV = cdash1 + cdash4*HIGH_W + cdash5*MED_V + cdash7*HIGH_W*MED_V;
DLOW_HIV = cdash1 + cdash4*LOW_W + cdash5*HIGH_V + cdash7*LOW_W*HIGH_V;
DMEW_HIV = cdash1 + cdash4*MED_W + cdash5*HIGH_V + cdash7*MED_W*HIGH_V;
DHIW_HIV = cdash1 + cdash4*HIGH_W + cdash5*HIGH_V + cdash7*HIGH_W*HIGH_V;
! Calc conditional total effects for each combination of moderator values
TLOW_LOV = ILOW_LOV + DLOW_LOV;
TMEW_LOV = IMEW_LOV + DMEW_LOV;
THIW_LOV = IHIW_LOV + DHIW_LOV;
TLOW_MEV = ILOW_MEV + DLOW_MEV;
TMEW_MEV = IMEW_MEV + DMEW_MEV;
THIW_MEV = IHIW_MEV + DHIW_MEV;
TLOW_HIV = ILOW_HIV + DLOW_HIV;
TMEW_HIV = IMEW_HIV + DMEW_HIV;
THIW_HIV = IHIW_HIV + DHIW_HIV;
! Use loop plot to plot conditional indirect effect of X on Y for each combination of low, med, high moderator values
! Could be edited to show conditional direct or conditional total effects instead
! NOTE - values of 1,5 in LOOP() statement need to be replaced by
! logical min and max limits of predictor X used in analysis
PLOT(PLOW_LOV PMEW_LOV PHIW_LOV PLOW_MEV PMEW_MEV PHIW_MEV
PLOW_HIV PMEW_HIV PHIW_HIV);
LOOP(XVAL,1,5,0.1);
PLOW_LOV = ILOW_LOV*XVAL;
PMEW_LOV = IMEW_LOV*XVAL;
PHIW_LOV = IHIW_LOV*XVAL;
PLOW_MEV = ILOW_MEV*XVAL;
PMEW_MEV = IMEW_MEV*XVAL;
PHIW_MEV = IHIW_MEV*XVAL;
PLOW_HIV = ILOW_HIV*XVAL;
PMEW_HIV = IMEW_HIV*XVAL;
PHIW_HIV = IHIW_HIV*XVAL;
PLOT:
TYPE = plot2;
OUTPUT:
STAND CINT(bcbootstrap);

第二十九步，输出标准化系数和基于自助法的偏差校正置信区间。

资源汇总

• 本视频讲义地址: https://mlln.cn/mplus-model-templates/model71.html
• 图书《MPlus中介调节模型》打包下载: 点击下载
• 图书《MPlus中介调节模型》在线看: 点击查看
• 视频教程: 点击这里打开视频
• Mplus 模型模板教程列表: https://mlln.cn/mplus-model-templates
• 统计咨询: https://wx.zsxq.com/group/88888188828842

感谢大家的收看！ 本节课的资源已全部提供，请大家点击页面上的链接下载或查看。 希望这份资料对大家的学习有所帮助，也欢迎大家点赞收藏，方便日后查阅！