Mplus model55 模型讲解

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

Mplus 中介、调节及调节中介效应模型代码速查

• 理论模型
• 数学模型
• 数学推导
• 代码解读

大家好，欢迎来到本期教程！ 这份PPT展示了我们Mplus模型代码速查工具的结构。 我们将会用Mplus代码来构建中介效应、调节效应以及调节中介效应模型。 本教程并非Mplus的系统讲解，而是为了方便大家快速查找代码。 每个模型都包含四个部分：首先我们会用图示的方式展示理论模型，让大家直观地理解模型结构；然后是相应的数学模型，用公式来表达模型；接着，对于比较复杂的模型，我们会给出详细的数学推导过程，帮助大家理解公式的来源；最后，当然是最重要的部分——Mplus代码解读，我会一步步解释代码的含义以及如何运行。 希望这个速查工具能够帮助大家更高效地使用Mplus进行数据分析。 Data Sense 也提供统计咨询服务，如有需要，欢迎联系我们！

理论模型

我们来看，这张图展示了一个比较复杂的模型。其中，X 是自变量，它会影响中介变量 M 和因变量 Y；M 也会影响 Y。此外，W 和 Z 调节 X 对 M 和 Y 的影响；V 和 Q 则调节 M 对 Y 的影响。整个模型考虑了所有可能的二阶和三阶交互作用，例如 X 和 W 的交互，M 和 V 的交互，以及 X、W、Z 三者的交互，还有 M、V、Q 三者的交互。总的来说，这个模型研究的是 X 如何通过 M 影响 Y，以及这些影响如何被 W、Z、V 和 Q 调节。

数学模型

这张图展示了一个复杂的模型，其中包含了自变量X，中介变量M和因变量Y。同时，还有四个调节变量：W, Z, V 和 Q。X影响M，影响路径系数是a1、a2、a3、a4、a5、a6和a7；M 影响Y，影响路径系数是b1；X直接影响Y，影响路径系数是c1’、c2’、c3’、c5’、c6’、c7’。 调节变量W, Z分别与X形成交互项XW, XZ和WZ，三个交互项也同时影响M和Y。调节变量V和Q，也分别与M形成交互项MV, MQ和VQ，这三个交互项同时也影响Y。另外，还有三个变量的交互项XWZ和MVQ，它们也影响Y。图中箭头上的a, b, c,表示了不同路径的影响效果，这些效应共同决定了X如何影响Y。

数学公式1

模型方程:
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + b6VQ + b7MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ

首先，我们来看模型方程。Y是因变量，M是中介变量，V和Q是调节变量，X、W和Z是自变量，带有撇号的c项代表直接效应，其他项则表示变量之间的相互作用，包括中介和调节作用。

数学公式2

模型方程:
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + b6VQ + b7MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2V + b3Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b5(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b6VQ + b7(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ

接下来，我们将中介变量M的方程代入到Y的方程中，以建立自变量X如何通过M影响Y的整体关系。

数学公式3

模型方程:
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + b6VQ + b7MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2V + b3Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b5(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b6VQ + b7(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
展开括号:
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1Z + a4b1XW + a5b1XZ + a6b1WZ + a7b1XWZ + b2V + b3Q + a0b4V + a1b4XV + a2b4WV + a3b4ZV + a4b4XWV + a5b4XZV + a6b4WZV + a7b4XWZV + a0b5Q + a1b5XQ + a2b5WQ + a3b5ZQ + a4b5XWQ + a5b5XZQ + a6b5WZQ + a7b5XWZQ + b6VQ + a0b7VQ + a1b7XVQ + a2b7WVQ + a3b7ZVQ + a4b7XWVQ + a5b7XZVQ + a6b7WZVQ + a7b7XWZVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ

现在，我们展开所有括号，以方便后续的同类项合并。

数学公式4

模型方程:
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + b6VQ + b7MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2V + b3Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b5(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b6VQ + b7(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
展开括号:
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1Z + a4b1XW + a5b1XZ + a6b1WZ + a7b1XWZ + b2V + b3Q + a0b4V + a1b4XV + a2b4WV + a3b4ZV + a4b4XWV + a5b4XZV + a6b4WZV + a7b4XWZV + a0b5Q + a1b5XQ + a2b5WQ + a3b5ZQ + a4b5XWQ + a5b5XZQ + a6b5WZQ + a7b5XWZQ + b6VQ + a0b7VQ + a1b7XVQ + a2b7WVQ + a3b7ZVQ + a4b7XWVQ + a5b7XZVQ + a6b7WZVQ + a7b7XWZVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
合并同类项,整理成 Y = a + bX 的形式:
Y = (b0 + a0b1 + a2b1W + a3b1Z + a6b1WZ + b2V + b3Q + a0b4V + a2b4WV + a3b4ZV + a6b4WZV + a0b5Q + a2b5WQ + a3b5ZQ + a6b5WZQ + b6VQ + a0b7VQ + a2b7WVQ + a3b7ZVQ + a6b7WZVQ + c2'W + c3'Z + c6'WZ) + (a1b1 + a4b1W + a5b1Z + a7b1WZ + a1b4V + a4b4WV + a5b4ZV + a7b4WZV + a1b5Q + a4b5WQ + a5b5ZQ + a7b5WZQ + a1b7VQ+ a4b7WVQ + a5b7ZVQ + a7b7WZVQ + c1' + c4'W + c5'Z + c7'WZ)X

合并同类项，并将公式整理成 Y = a + bX 的形式，其中a表示截距，b表示X的总效应，包括直接效应和间接效应。

数学公式5

模型方程:
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + b6VQ + b7MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2V + b3Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b5(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b6VQ + b7(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
展开括号:
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1Z + a4b1XW + a5b1XZ + a6b1WZ + a7b1XWZ + b2V + b3Q + a0b4V + a1b4XV + a2b4WV + a3b4ZV + a4b4XWV + a5b4XZV + a6b4WZV + a7b4XWZV + a0b5Q + a1b5XQ + a2b5WQ + a3b5ZQ + a4b5XWQ + a5b5XZQ + a6b5WZQ + a7b5XWZQ + b6VQ + a0b7VQ + a1b7XVQ + a2b7WVQ + a3b7ZVQ + a4b7XWVQ + a5b7XZVQ + a6b7WZVQ + a7b7XWZVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
合并同类项,整理成 Y = a + bX 的形式:
Y = (b0 + a0b1 + a2b1W + a3b1Z + a6b1WZ + b2V + b3Q + a0b4V + a2b4WV + a3b4ZV + a6b4WZV + a0b5Q + a2b5WQ + a3b5ZQ + a6b5WZQ + b6VQ + a0b7VQ + a2b7WVQ + a3b7ZVQ + a6b7WZVQ + c2'W + c3'Z + c6'WZ) + (a1b1 + a4b1W + a5b1Z + a7b1WZ + a1b4V + a4b4WV + a5b4ZV + a7b4WZV + a1b5Q + a4b5WQ + a5b5ZQ + a7b5WZQ + a1b7VQ+ a4b7WVQ + a5b7ZVQ + a7b7WZVQ + c1' + c4'W + c5'Z + c7'WZ)X
X对Y的间接效应（在W, Z, V, Q条件下）:
a1b1 + a4b1W + a5b1Z + a7b1WZ + a1b4V + a4b4WV + a5b4ZV + a7b4WZV + a1b5Q + a4b5WQ + a5b5ZQ + a7b5WZQ + a1b7VQ+ a4b7WVQ + a5b7ZVQ + a7b7WZVQ = (a1 + a4W + a5Z + a7WZ)(b1 + b4V + b5Q + b7VQ)

从整理后的公式中，我们可以提取出X对Y的间接效应，即X通过中介变量M对Y的影响。这个效应受到W、Z、V和Q的调节作用。

数学公式6

模型方程:
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + b6VQ + b7MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2V + b3Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b5(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b6VQ + b7(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
展开括号:
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1Z + a4b1XW + a5b1XZ + a6b1WZ + a7b1XWZ + b2V + b3Q + a0b4V + a1b4XV + a2b4WV + a3b4ZV + a4b4XWV + a5b4XZV + a6b4WZV + a7b4XWZV + a0b5Q + a1b5XQ + a2b5WQ + a3b5ZQ + a4b5XWQ + a5b5XZQ + a6b5WZQ + a7b5XWZQ + b6VQ + a0b7VQ + a1b7XVQ + a2b7WVQ + a3b7ZVQ + a4b7XWVQ + a5b7XZVQ + a6b7WZVQ + a7b7XWZVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ
合并同类项,整理成 Y = a + bX 的形式:
Y = (b0 + a0b1 + a2b1W + a3b1Z + a6b1WZ + b2V + b3Q + a0b4V + a2b4WV + a3b4ZV + a6b4WZV + a0b5Q + a2b5WQ + a3b5ZQ + a6b5WZQ + b6VQ + a0b7VQ + a2b7WVQ + a3b7ZVQ + a6b7WZVQ + c2'W + c3'Z + c6'WZ) + (a1b1 + a4b1W + a5b1Z + a7b1WZ + a1b4V + a4b4WV + a5b4ZV + a7b4WZV + a1b5Q + a4b5WQ + a5b5ZQ + a7b5WZQ + a1b7VQ+ a4b7WVQ + a5b7ZVQ + a7b7WZVQ + c1' + c4'W + c5'Z + c7'WZ)X
X对Y的间接效应（在W, Z, V, Q条件下）:
a1b1 + a4b1W + a5b1Z + a7b1WZ + a1b4V + a4b4WV + a5b4ZV + a7b4WZV + a1b5Q + a4b5WQ + a5b5ZQ + a7b5WZQ + a1b7VQ+ a4b7WVQ + a5b7ZVQ + a7b7WZVQ = (a1 + a4W + a5Z + a7WZ)(b1 + b4V + b5Q + b7VQ)
X对Y的直接效应（在W, Z条件下）:
c1' + c4'W + c5'Z + c7'WZ

最后，我们可以提取出X对Y的直接效应，这个效应受到W和Z的调节作用。

代码解读1

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;

第一步，我们声明分析中需要用到的变量。X是自变量，M是中介变量，W、Z、V、Q是四个调节变量，Y是因变量。后面这些XW、XZ等则是我们通过定义语句生成的交互项，代表变量之间的乘积。

代码解读2

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;

第二步，我们定义交互项。例如MQ等于M乘以Q，也就是中介变量M和调节变量Q的交互作用。其他交互项也类似，包括二阶和三阶的交互作用。这些交互项会在后续的分析中用到。

代码解读3

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;

第三步，我们设置分析的参数。TYPE=GENERAL表示使用一般线性模型进行分析。ESTIMATOR=ML表示使用最大似然估计方法。BOOTSTRAP=10000表示进行一万次自助抽样，用来计算置信区间。

代码解读4

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);

第四步，我们定义模型结构。Y的截距项是b0。Y对M的回归系数是b1，Y对V的回归系数是b2，以此类推。这部分定义了因变量Y和中介变量M所有相关路径的系数，包括自变量X、调节变量及其交互项的影响。例如M对X的回归系数是a1。

代码解读5

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);
MODEL CONSTRAINT:
 NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
 ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
 ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
 ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
 ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
 ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
 ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
 ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
 ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
 ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
 DLOW_LOZ DMEW_LOZ DHIW_LOZ DLOW_MEZ DMEW_MEZ DHIW_MEZ
 DLOW_HIZ DMEW_HIZ DHIW_HIZ
 TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
 TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
 TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
 TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
 TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
 TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
 TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
 TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
 TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);

第五步，也是最关键的一步，我们定义模型约束。这部分代码定义了不同调节变量取值组合下的条件间接效应、条件直接效应和条件总效应。我们先定义了调节变量的低、中、高三个水平的值，需要用户根据自己的数据进行替换。然后定义各种效应的组合形式，以便后续的公式计算。

代码解读6

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);
MODEL CONSTRAINT:
 NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
 ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
 ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
 ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
 ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
 ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
 ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
 ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
 ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
 ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
 DLOW_LOZ DMEW_LOZ DHIW_LOZ DLOW_MEZ DMEW_MEZ DHIW_MEZ
 DLOW_HIZ DMEW_HIZ DHIW_HIZ
 TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
 TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
 TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
 TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
 TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
 TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
 TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
 TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
 TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
 MED_W = #MEDW;
 HIGH_W = #HIGHW;
 LOW_Z = #LOWZ;
 MED_Z = #MEDZ;
 HIGH_Z = #HIGHZ;
 LOW_V = #LOWV;
 MED_V = #MEDV;
 HIGH_V = #HIGHV;
 LOW_Q = #LOWQ;
 MED_Q = #MEDQ;
 HIGH_Q = #HIGHQ;

接下来，我们需要用具体数值替换先前定义的调节变量的低、中、高水平的占位符。这需要根据你自己的数据来设定这些数值。

代码解读7

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);
MODEL CONSTRAINT:
 NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
 ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
 ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
 ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
 ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
 ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
 ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
 ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
 ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
 ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
 DLOW_LOZ DMEW_LOZ DHIW_LOZ DLOW_MEZ DMEW_MEZ DHIW_MEZ
 DLOW_HIZ DMEW_HIZ DHIW_HIZ
 TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
 TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
 TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
 TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
 TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
 TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
 TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
 TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
 TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
 MED_W = #MEDW;
 HIGH_W = #HIGHW;
 LOW_Z = #LOWZ;
 MED_Z = #MEDZ;
 HIGH_Z = #HIGHZ;
 LOW_V = #LOWV;
 MED_V = #MEDV;
 HIGH_V = #HIGHV;
 LOW_Q = #LOWQ;
 MED_Q = #MEDQ;
 HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z +
 a1*b4*LOW_V + a4*b4*LOW_W*LOW_V + a5*b4*LOW_Z*LOW_V +
 a7*b4*LOW_W*LOW_Z*LOW_V + a1*b5*LOW_Q + a4*b5*LOW_W*LOW_Q +
 a5*b5*LOW_Z*LOW_Q + a7*b5*LOW_W*LOW_Z*LOW_Q + a1*b7*LOW_V*LOW_Q +
 a4*b7*LOW_W*LOW_V*LOW_Q + a5*b7*LOW_Z*LOW_V*LOW_Q +
 a7*b7*LOW_W*LOW_Z*LOW_V*LOW_Q;

这部分是一个例子，展示了如何计算当四个调节变量都处于低水平时，自变量X对因变量Y的间接效应ILLLL。公式使用了之前定义的路径系数a1、b1等，和调节变量的低水平值。之后的很多公式和这个类似，只是调节变量的取值不同。

代码解读8

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);
MODEL CONSTRAINT:
 NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
 ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
 ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
 ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
 ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
 ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
 ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
 ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
 ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
 ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
 DLOW_LOZ DMEW_LOZ DHIW_LOZ DLOW_MEZ DMEW_MEZ DHIW_MEZ
 DLOW_HIZ DMEW_HIZ DHIW_HIZ
 TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
 TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
 TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
 TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
 TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
 TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
 TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
 TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
 TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
 MED_W = #MEDW;
 HIGH_W = #HIGHW;
 LOW_Z = #LOWZ;
 MED_Z = #MEDZ;
 HIGH_Z = #HIGHZ;
 LOW_V = #LOWV;
 MED_V = #MEDV;
 HIGH_V = #HIGHV;
 LOW_Q = #LOWQ;
 MED_Q = #MEDQ;
 HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z +
 a1*b4*LOW_V + a4*b4*LOW_W*LOW_V + a5*b4*LOW_Z*LOW_V +
 a7*b4*LOW_W*LOW_Z*LOW_V + a1*b5*LOW_Q + a4*b5*LOW_W*LOW_Q +
 a5*b5*LOW_Z*LOW_Q + a7*b5*LOW_W*LOW_Z*LOW_Q + a1*b7*LOW_V*LOW_Q +
 a4*b7*LOW_W*LOW_V*LOW_Q + a5*b7*LOW_Z*LOW_V*LOW_Q +
 a7*b7*LOW_W*LOW_Z*LOW_V*LOW_Q;
... (其余的MODEL CONSTRAINT部分省略，因为篇幅过长) ...

由于篇幅限制，我们省略了剩余的MODEL CONSTRAINT部分的公式定义。但是它们和前面ILLLL的计算原理是一样的，只是针对不同的调节变量组合和效应类型进行计算。

代码解读9

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);
MODEL CONSTRAINT:
 NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
 ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
 ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
 ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
 ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
 ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
 ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
 ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
 ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
 ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
 DLOW_LOZ DMEW_LOZ DHIW_LOZ DLOW_MEZ DMEW_MEZ DHIW_MEZ
 DLOW_HIZ DMEW_HIZ DHIW_HIZ
 TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
 TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
 TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
 TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
 TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
 TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
 TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
 TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
 TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
 MED_W = #MEDW;
 HIGH_W = #HIGHW;
 LOW_Z = #LOWZ;
 MED_Z = #MEDZ;
 HIGH_Z = #HIGHZ;
 LOW_V = #LOWV;
 MED_V = #MEDV;
 HIGH_V = #HIGHV;
 LOW_Q = #LOWQ;
 MED_Q = #MEDQ;
 HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z +
 a1*b4*LOW_V + a4*b4*LOW_W*LOW_V + a5*b4*LOW_Z*LOW_V +
 a7*b4*LOW_W*LOW_Z*LOW_V + a1*b5*LOW_Q + a4*b5*LOW_W*LOW_Q +
 a5*b5*LOW_Z*LOW_Q + a7*b5*LOW_W*LOW_Z*LOW_Q + a1*b7*LOW_V*LOW_Q +
 a4*b7*LOW_W*LOW_V*LOW_Q + a5*b7*LOW_Z*LOW_V*LOW_Q +
 a7*b7*LOW_W*LOW_Z*LOW_V*LOW_Q;
... (其余的MODEL CONSTRAINT部分省略，因为篇幅过长) ...
PLOT(PLLLL PMLLL PHLLL PLMLL PMMLL PHMLL PLHLL PMHLL PHHLL
 PLLML PMLML PHLML PLMML PMMML PHMML PLHML PMHML PHHML
 PLLHL PMLHL PHLHL PLMHL PMMHL PHMHL PLHHL PMHHL PHHHL
 PLLLM PMLLM PHLLM PLMLM PMMLM PHMLM PLHLM PMHLM PHHLM
 PLLMM PMLMM PHLMM PLMMM PMMMM PHMMM PLHMM PMHMM PHHMM
 PLLHM PMLHM PHLHM PLMHM PMMHM PHMHM PLHHM PMHHM PHHHM
 PLLLH PMLLH PHLLH PLMLH PMMLH PHMLH PLHLH PMHLH PHHLH
 PLLMH PMLMH PHLMH PLMMH PMMMH PHMMH PLHMH PMHMH PHHMH
 PLLHH PMLHH PHLHH PLMHH PMMHH PHMHH PLHHH PMHHH PHHHH);

这部分代码列出了需要绘制的各种效应的组合。每个名称都代表在特定调节变量组合下的间接效应，例如PLLLL代表当所有调节变量都为低值时的间接效应。

代码解读10

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);
MODEL CONSTRAINT:
 NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
 ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
 ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
 ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
 ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
 ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
 ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
 ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
 ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
 ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
 DLOW_LOZ DMEW_LOZ DHIW_LOZ DLOW_MEZ DMEW_MEZ DHIW_MEZ
 DLOW_HIZ DMEW_HIZ DHIW_HIZ
 TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
 TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
 TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
 TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
 TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
 TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
 TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
 TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
 TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
 MED_W = #MEDW;
 HIGH_W = #HIGHW;
 LOW_Z = #LOWZ;
 MED_Z = #MEDZ;
 HIGH_Z = #HIGHZ;
 LOW_V = #LOWV;
 MED_V = #MEDV;
 HIGH_V = #HIGHV;
 LOW_Q = #LOWQ;
 MED_Q = #MEDQ;
 HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z +
 a1*b4*LOW_V + a4*b4*LOW_W*LOW_V + a5*b4*LOW_Z*LOW_V +
 a7*b4*LOW_W*LOW_Z*LOW_V + a1*b5*LOW_Q + a4*b5*LOW_W*LOW_Q +
 a5*b5*LOW_Z*LOW_Q + a7*b5*LOW_W*LOW_Z*LOW_Q + a1*b7*LOW_V*LOW_Q +
 a4*b7*LOW_W*LOW_V*LOW_Q + a5*b7*LOW_Z*LOW_V*LOW_Q +
 a7*b7*LOW_W*LOW_Z*LOW_V*LOW_Q;
... (其余的MODEL CONSTRAINT部分省略，因为篇幅过长) ...
PLOT(PLLLL PMLLL PHLLL PLMLL PMMLL PHMLL PLHLL PMHLL PHHLL
 PLLML PMLML PHLML PLMML PMMML PHMML PLHML PMHML PHHML
 PLLHL PMLHL PHLHL PLMHL PMMHL PHMHL PLHHL PMHHL PHHHL
 PLLLM PMLLM PHLLM PLMLM PMMLM PHMLM PLHLM PMHLM PHHLM
 PLLMM PMLMM PHLMM PLMMM PMMMM PHMMM PLHMM PMHMM PHHMM
 PLLHM PMLHM PHLHM PLMHM PMMHM PHMHM PLHHM PMHHM PHHHM
 PLLLH PMLLH PHLLH PLMLH PMMLH PHMLH PLHLH PMHLH PHHLH
 PLLMH PMLMH PHLMH PLMMH PMMMH PHMMH PLHMH PMHMH PHHMH
 PLLHH PMLHH PHLHH PLMHH PMMHH PHMHH PLHHH PMHHH PHHHH);
LOOP(XVAL,1,5,0.1);
PLLLL = ILLLL*XVAL;
 ... (其余的LOOP部分省略，因为篇幅过长) ...

这里我们使用循环计算不同X值下的条件间接效应。XVAL从1到5以0.1的步长递增，ILLLL乘以XVAL可以得到不同X取值下的间接效应。类似地，后面的代码也计算了其他条件间接效应在不同X值下的取值。篇幅限制省略了重复的部分。

代码解读11

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);
MODEL CONSTRAINT:
 NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
 ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
 ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
 ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
 ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
 ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
 ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
 ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
 ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
 ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
 DLOW_LOZ DMEW_LOZ DHIW_LOZ DLOW_MEZ DMEW_MEZ DHIW_MEZ
 DLOW_HIZ DMEW_HIZ DHIW_HIZ
 TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
 TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
 TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
 TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
 TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
 TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
 TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
 TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
 TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
 MED_W = #MEDW;
 HIGH_W = #HIGHW;
 LOW_Z = #LOWZ;
 MED_Z = #MEDZ;
 HIGH_Z = #HIGHZ;
 LOW_V = #LOWV;
 MED_V = #MEDV;
 HIGH_V = #HIGHV;
 LOW_Q = #LOWQ;
 MED_Q = #MEDQ;
 HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z +
 a1*b4*LOW_V + a4*b4*LOW_W*LOW_V + a5*b4*LOW_Z*LOW_V +
 a7*b4*LOW_W*LOW_Z*LOW_V + a1*b5*LOW_Q + a4*b5*LOW_W*LOW_Q +
 a5*b5*LOW_Z*LOW_Q + a7*b5*LOW_W*LOW_Z*LOW_Q + a1*b7*LOW_V*LOW_Q +
 a4*b7*LOW_W*LOW_V*LOW_Q + a5*b7*LOW_Z*LOW_V*LOW_Q +
 a7*b7*LOW_W*LOW_Z*LOW_V*LOW_Q;
... (其余的MODEL CONSTRAINT部分省略，因为篇幅过长) ...
PLOT(PLLLL PMLLL PHLLL PLMLL PMMLL PHMLL PLHLL PMHLL PHHLL
 PLLML PMLML PHLML PLMML PMMML PHMML PLHML PMHML PHHML
 PLLHL PMLHL PHLHL PLMHL PMMHL PHMHL PLHHL PMHHL PHHHL
 PLLLM PMLLM PHLLM PLMLM PMMLM PHMLM PLHLM PMHLM PHHLM
 PLLMM PMLMM PHLMM PLMMM PMMMM PHMMM PLHMM PMHMM PHHMM
 PLLHM PMLHM PHLHM PLMHM PMMHM PHMHM PLHHM PMHHM PHHHM
 PLLLH PMLLH PHLLH PLMLH PMMLH PHMLH PLHLH PMHLH PHHLH
 PLLMH PMLMH PHLMH PLMMH PMMMH PHMMH PLHMH PMHMH PHHMH
 PLLHH PMLHH PHLHH PLMHH PMMHH PHMHH PLHHH PMHHH PHHHH);
LOOP(XVAL,1,5,0.1);
PLLLL = ILLLL*XVAL;
 ... (其余的LOOP部分省略，因为篇幅过长) ...
PLOT:
 TYPE = plot2;

最后，我们使用PLOT命令生成图形，TYPE=plot2 表示使用第二种绘图类型。这部分代码会生成间接效应随自变量X变化的图像。

代码解读12

USEVARIABLES = X M W Z V Q Y XW XZ WZ VQ MV MQ XWZ MVQ;
DEFINE:
 MQ = M*Q;
 MV = M*V;
 XW = X*W;
 XZ = X*Z;
 WZ = W*Z;
 VQ = V*Q;
 MVQ = M*V*Q;
 XWZ = X*W*Z;
ANALYSIS:
 TYPE = GENERAL;
 ESTIMATOR = ML;
 BOOTSTRAP = 10000;
MODEL:
 [Y] (b0);
 Y ON M (b1);
 Y ON V (b2);
 Y ON Q (b3);
 Y ON MV (b4);
 Y ON MQ (b5);
 Y ON VQ (b6);
 Y ON MVQ (b7);
 Y ON X (cdash1);
 Y ON W (cdash2);
 Y ON Z (cdash3);
 Y ON XW (cdash4);
 Y ON XZ (cdash5);
 Y ON WZ (cdash6);
 Y ON XWZ (cdash7);
 [M] (a0);
 M ON X (a1);
 M ON W (a2);
 M ON Z (a3);
 M ON XW (a4);
 M ON XZ (a5);
 M ON WZ (a6);
 M ON XWZ (a7);
MODEL CONSTRAINT:
 NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
 ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
 ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
 ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
 ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
 ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
 ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
 ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
 ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
 ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
 DLOW_LOZ DMEW_LOZ DHIW_LOZ DLOW_MEZ DMEW_MEZ DHIW_MEZ
 DLOW_HIZ DMEW_HIZ DHIW_HIZ
 TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
 TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
 TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
 TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
 TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
 TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
 TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
 TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
 TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
 MED_W = #MEDW;
 HIGH_W = #HIGHW;
 LOW_Z = #LOWZ;
 MED_Z = #MEDZ;
 HIGH_Z = #HIGHZ;
 LOW_V = #LOWV;
 MED_V = #MEDV;
 HIGH_V = #HIGHV;
 LOW_Q = #LOWQ;
 MED_Q = #MEDQ;
 HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z +
 a1*b4*LOW_V + a4*b4*LOW_W*LOW_V + a5*b4*LOW_Z*LOW_V +
 a7*b4*LOW_W*LOW_Z*LOW_V + a1*b5*LOW_Q + a4*b5*LOW_W*LOW_Q +
 a5*b5*LOW_Z*LOW_Q + a7*b5*LOW_W*LOW_Z*LOW_Q + a1*b7*LOW_V*LOW_Q +
 a4*b7*LOW_W*LOW_V*LOW_Q + a5*b7*LOW_Z*LOW_V*LOW_Q +
 a7*b7*LOW_W*LOW_Z*LOW_V*LOW_Q;
... (其余的MODEL CONSTRAINT部分省略，因为篇幅过长) ...
PLOT(PLLLL PMLLL PHLLL PLMLL PMMLL PHMLL PLHLL PMHLL PHHLL
 PLLML PMLML PHLML PLMML PMMML PHMML PLHML PMHML PHHML
 PLLHL PMLHL PHLHL PLMHL PMMHL PHMHL PLHHL PMHHL PHHHL
 PLLLM PMLLM PHLLM PLMLM PMMLM PHMLM PLHLM PMHLM PHHLM
 PLLMM PMLMM PHLMM PLMMM PMMMM PHMMM PLHMM PMHMM PHHMM
 PLLHM PMLHM PHLHM PLMHM PMMHM PHMHM PLHHM PMHHM PHHHM
 PLLLH PMLLH PHLLH PLMLH PMMLH PHMLH PLHLH PMHLH PHHLH
 PLLMH PMLMH PHLMH PLMMH PMMMH PHMMH PLHMH PMHMH PHHMH
 PLLHH PMLHH PHLHH PLMHH PMMHH PHMHH PLHHH PMHHH PHHHH);
LOOP(XVAL,1,5,0.1);
PLLLL = ILLLL*XVAL;
 ... (其余的LOOP部分省略，因为篇幅过长) ...
PLOT:
 TYPE = plot2;
OUTPUT:
 STAND CINT(bcbootstrap);

最后，我们指定输出内容。STAND表示输出标准化参数估计值。CINT(bcbootstrap) 表示输出基于bootstrap方法的置信区间。

资源汇总

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

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