来自图书《MPlus中介调节模型》
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + c'X
M = a0 + a1W + a2Z + a3X + a4XW + a5XZ
Y = b0 + b1(a0 + a1W + a2Z + a3X + a4XW + a5XZ) + b2V + b3Q + b4(a0 + a1W + a2Z + a3X + a4XW + a5XZ)V + b5(a0 + a1W + a2Z + a3X + a4XW + a5XZ)Q + c'X
Y = b0 + a0b1 + a1b1W + a2b1Z + a3b1X + a4b1XW + a5b1XZ + b2V + b3Q + a0b4V + a1b4WV + a2b4ZV + a3b4XV + a4b4XWV + a5b4XZV + a0b5Q + a1b5WQ + a2b5ZQ + a3b5XQ + a4b5XWQ + a5b5XZQ + c'X
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + c'X
M = a0 + a1W + a2Z + a3X + a4XW + a5XZ
Y = b0 + b1(a0 + a1W + a2Z + a3X + a4XW + a5XZ) + b2V + b3Q + b4(a0 + a1W + a2Z + a3X + a4XW + a5XZ)V + b5(a0 + a1W + a2Z + a3X + a4XW + a5XZ)Q + c'X
Y = b0 + a0b1 + a1b1W + a2b1Z + a3b1X + a4b1XW + a5b1XZ + b2V + b3Q + a0b4V + a1b4WV + a2b4ZV + a3b4XV + a4b4XWV + a5b4XZV + a0b5Q + a1b5WQ + a2b5ZQ + a3b5XQ + a4b5XWQ + a5b5XZQ + c'X
Y = (b0 + a0b1 + a1b1W + a2b1Z + b2V + b3Q + a0b4V + a1b4WV + a2b4ZV + a0b5Q + a1b5WQ + a2b5ZQ) + (a3b1 + a4b1W + a5b1Z + a3b4V + a4b4WV + a5b4ZV + a3b5Q + a4b5WQ + a5b5ZQ + c')X
Y = b0 + b1M + b2V + b3Q + b4MV + b5MQ + c'X
M = a0 + a1W + a2Z + a3X + a4XW + a5XZ
Y = b0 + b1(a0 + a1W + a2Z + a3X + a4XW + a5XZ) + b2V + b3Q + b4(a0 + a1W + a2Z + a3X + a4XW + a5XZ)V + b5(a0 + a1W + a2Z + a3X + a4XW + a5XZ)Q + c'X
Y = b0 + a0b1 + a1b1W + a2b1Z + a3b1X + a4b1XW + a5b1XZ + b2V + b3Q + a0b4V + a1b4WV + a2b4ZV + a3b4XV + a4b4XWV + a5b4XZV + a0b5Q + a1b5WQ + a2b5ZQ + a3b5XQ + a4b5XWQ + a5b5XZQ + c'X
Y = (b0 + a0b1 + a1b1W + a2b1Z + b2V + b3Q + a0b4V + a1b4WV + a2b4ZV + a0b5Q + a1b5WQ + a2b5ZQ) + (a3b1 + a4b1W + a5b1Z + a3b4V + a4b4WV + a5b4ZV + a3b5Q + a4b5WQ + a5b5ZQ + c')X
X 对 Y 的间接效应 (条件依赖于 W, Z, V, Q):(a3 + a4W + a5Z)(b1 + b4V + b5Q)
X 对 Y 的直接效应:c'
USEVARIABLES = X M W Z V Q Y XW XZ MV MQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*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 X (cdash);
[M] (a0);
M ON W (a1);
M ON Z (a2);
M ON X (a3);
M ON XW (a4);
M ON XZ (a5);
USEVARIABLES = X M W Z V Q Y XW XZ MV MQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*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 X (cdash);
[M] (a0);
M ON W (a1);
M ON Z (a2);
M ON X (a3);
M ON XW (a4);
M ON XZ (a5);
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 ...);
LOW_W = #LOWW; ...
ILLLL = a3*b1 + a4*b1*LOW_W + ...;
USEVARIABLES = X M W Z V Q Y XW XZ MV MQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*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 X (cdash);
[M] (a0);
M ON W (a1);
M ON Z (a2);
M ON X (a3);
M ON XW (a4);
M ON XZ (a5);
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 ...);
LOW_W = #LOWW; ...
ILLLL = a3*b1 + a4*b1*LOW_W + ...;
TLLLL = ILLLL + cdash; ...
USEVARIABLES = X M W Z V Q Y XW XZ MV MQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*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 X (cdash);
[M] (a0);
M ON W (a1);
M ON Z (a2);
M ON X (a3);
M ON XW (a4);
M ON XZ (a5);
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 ...);
LOW_W = #LOWW; ...
ILLLL = a3*b1 + a4*b1*LOW_W + ...;
TLLLL = ILLLL + cdash; ...
PLOT(...); LOOP(XVAL,1,5,0.1); PLLLL = ILLLL*XVAL; ...
USEVARIABLES = X M W Z V Q Y XW XZ MV MQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*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 X (cdash);
[M] (a0);
M ON W (a1);
M ON Z (a2);
M ON X (a3);
M ON XW (a4);
M ON XZ (a5);
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 ...);
LOW_W = #LOWW; ...
ILLLL = a3*b1 + a4*b1*LOW_W + ...;
TLLLL = ILLLL + cdash; ...
PLOT(...); LOOP(XVAL,1,5,0.1); PLLLL = ILLLL*XVAL; ...
PLOT: TYPE = plot2;
OUTPUT: STAND CINT(bcbootstrap);