来自图书《MPlus中介调节模型》
USEVARIABLES = X M W Y XW;
DEFINE: XW = X*W;
ANALYSIS: TYPE = GENERAL; ESTIMATOR = ML; BOOTSTRAP = 10000;
[Y] (b0); Y ON M (b1); Y ON X (cdash);
[M] (a0); M ON X (a1); M ON W (a2); M ON XW (a3);
MODEL CONSTRAINT: NEW(LOW_W MED_W HIGH_W IND_LOWW IND_MEDW IND_HIW IMM TOT_LOWW TOT_MEDW TOT_HIW);
LOW_W = #LOWW; MED_W = #MEDW; HIGH_W = #HIGHW;
USEVARIABLES = X M W Y XW;
DEFINE: XW = X*W;
ANALYSIS: TYPE = GENERAL; ESTIMATOR = ML; BOOTSTRAP = 10000;
[Y] (b0); Y ON M (b1); Y ON X (cdash);
[M] (a0); M ON X (a1); M ON W (a2); M ON XW (a3);
MODEL CONSTRAINT: NEW(LOW_W MED_W HIGH_W IND_LOWW IND_MEDW IND_HIW IMM TOT_LOWW TOT_MEDW TOT_HIW);
LOW_W = #LOWW; MED_W = #MEDW; HIGH_W = #HIGHW;
IND_LOWW = a1*b1 + a3*b1*LOW_W; IND_MEDW = a1*b1 + a3*b1*MED_W; IND_HIW = a1*b1 + a3*b1*HIGH_W;
USEVARIABLES = X M W Y XW;
DEFINE: XW = X*W;
ANALYSIS: TYPE = GENERAL; ESTIMATOR = ML; BOOTSTRAP = 10000;
[Y] (b0); Y ON M (b1); Y ON X (cdash);
[M] (a0); M ON X (a1); M ON W (a2); M ON XW (a3);
MODEL CONSTRAINT: NEW(LOW_W MED_W HIGH_W IND_LOWW IND_MEDW IND_HIW IMM TOT_LOWW TOT_MEDW TOT_HIW);
LOW_W = #LOWW; MED_W = #MEDW; HIGH_W = #HIGHW;
IND_LOWW = a1*b1 + a3*b1*LOW_W; IND_MEDW = a1*b1 + a3*b1*MED_W; IND_HIW = a1*b1 + a3*b1*HIGH_W;
IMM = a3*b1; TOT_LOWW = IND_LOWW + cdash; TOT_MEDW = IND_MEDW + cdash; TOT_HIW = IND_HIW + cdash;
USEVARIABLES = X M W Y XW;
DEFINE: XW = X*W;
ANALYSIS: TYPE = GENERAL; ESTIMATOR = ML; BOOTSTRAP = 10000;
[Y] (b0); Y ON M (b1); Y ON X (cdash);
[M] (a0); M ON X (a1); M ON W (a2); M ON XW (a3);
MODEL CONSTRAINT: NEW(LOW_W MED_W HIGH_W IND_LOWW IND_MEDW IND_HIW IMM TOT_LOWW TOT_MEDW TOT_HIW);
LOW_W = #LOWW; MED_W = #MEDW; HIGH_W = #HIGHW;
IND_LOWW = a1*b1 + a3*b1*LOW_W; IND_MEDW = a1*b1 + a3*b1*MED_W; IND_HIW = a1*b1 + a3*b1*HIGH_W;
IMM = a3*b1; TOT_LOWW = IND_LOWW + cdash; TOT_MEDW = IND_MEDW + cdash; TOT_HIW = IND_HIW + cdash;
PLOT(LOMOD MEDMOD HIMOD); LOOP(XVAL,1,5,0.1); LOMOD = IND_LOWW*XVAL; MEDMOD = IND_MEDW*XVAL; HIMOD = IND_HIW*XVAL;
USEVARIABLES = X M W Y XW;
DEFINE: XW = X*W;
ANALYSIS: TYPE = GENERAL; ESTIMATOR = ML; BOOTSTRAP = 10000;
[Y] (b0); Y ON M (b1); Y ON X (cdash);
[M] (a0); M ON X (a1); M ON W (a2); M ON XW (a3);
MODEL CONSTRAINT: NEW(LOW_W MED_W HIGH_W IND_LOWW IND_MEDW IND_HIW IMM TOT_LOWW TOT_MEDW TOT_HIW);
LOW_W = #LOWW; MED_W = #MEDW; HIGH_W = #HIGHW;
IND_LOWW = a1*b1 + a3*b1*LOW_W; IND_MEDW = a1*b1 + a3*b1*MED_W; IND_HIW = a1*b1 + a3*b1*HIGH_W;
IMM = a3*b1; TOT_LOWW = IND_LOWW + cdash; TOT_MEDW = IND_MEDW + cdash; TOT_HIW = IND_HIW + cdash;
PLOT(LOMOD MEDMOD HIMOD); LOOP(XVAL,1,5,0.1); LOMOD = IND_LOWW*XVAL; MEDMOD = IND_MEDW*XVAL; HIMOD = IND_HIW*XVAL;
PLOT: TYPE = plot2; OUTPUT: STAND CINT(bcbootstrap);