Mplus VERSION 7.4
MUTHEN & MUTHEN
06/03/2016 2:59 PM
INPUT INSTRUCTIONS
title:
Simulating Z moderation of X to M using a random slope,
saving the data for external Monte Carlo analysis
montecarlo:
names = y m x z;
nobs = 400;
nreps = 500;
repsave = all;
save = xzrep*.dat;
cutpoints = x(0);
model population:
x-z@1;
x with z@0.5;
y on m*.5 x*.2 z*.1;
y*.5; [y*0];
gamma1 | m on x;
[gamma1*.3];
gamma1 on z*.2;
gamma1@0;
m on z*.3;
m*.5; [m*0];
analysis:
type = random;
model:
y on m*.5 (b)
x*.2 z*.1;
y*.5; [y*0];
gamma1 | m on x;
[gamma1*.3] (gamma1);
gamma1 on z*.2 (gamma3);
gamma1@0;
m on z*.3;
m*.5; [m*0];
model constraint:
new(indavg*.15 indlow*.05 indhigh*.25);
indavg = b*gamma1;
indlow = b*(gamma1-gamma3);
indhigh = b*(gamma1+gamma3);
INPUT READING TERMINATED NORMALLY
Simulating Z moderation of X to M using a random slope,
saving the data for external Monte Carlo analysis
SUMMARY OF ANALYSIS
Number of groups 1
Number of observations 400
Number of replications
Requested 500
Completed 500
Value of seed 0
Number of dependent variables 2
Number of independent variables 2
Number of continuous latent variables 1
Observed dependent variables
Continuous
Y M
Observed independent variables
X Z
Continuous latent variables
GAMMA1
Estimator MLR
Information matrix OBSERVED
Maximum number of iterations 100
Convergence criterion 0.100D-05
Maximum number of EM iterations 500
Convergence criteria for the EM algorithm
Loglikelihood change 0.100D-02
Relative loglikelihood change 0.100D-05
Derivative 0.100D-03
Minimum variance 0.100D-03
Maximum number of steepest descent iterations 20
Optimization algorithm EMA
SAMPLE STATISTICS FOR THE FIRST REPLICATION
ESTIMATED SAMPLE STATISTICS
Means
Y M X Z
________ ________ ________ ________
1 0.195 0.241 0.538 0.007
Covariances
Y M X Z
________ ________ ________ ________
Y 0.764
M 0.448 0.705
X 0.159 0.156 0.249
Z 0.361 0.425 0.234 1.054
Correlations
Y M X Z
________ ________ ________ ________
Y 1.000
M 0.610 1.000
X 0.365 0.373 1.000
Z 0.402 0.493 0.456 1.000
MODEL FIT INFORMATION
Number of Free Parameters 10
Loglikelihood
H0 Value
Mean -852.092
Std Dev 19.695
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.994 -897.908 -897.070
0.980 0.982 -892.539 -891.917
0.950 0.956 -884.488 -883.800
0.900 0.914 -877.333 -876.734
0.800 0.788 -868.667 -869.650
0.700 0.694 -862.420 -862.758
0.500 0.484 -852.092 -853.131
0.300 0.294 -841.764 -842.160
0.200 0.194 -835.517 -835.999
0.100 0.110 -826.851 -826.193
0.050 0.056 -819.696 -819.522
0.020 0.016 -811.645 -814.453
0.010 0.008 -806.276 -807.064
Information Criteria
Akaike (AIC)
Mean 1724.184
Std Dev 39.390
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.992 1632.552 1633.888
0.980 0.984 1643.290 1647.604
0.950 0.944 1659.393 1659.013
0.900 0.890 1673.703 1671.096
0.800 0.806 1691.034 1691.902
0.700 0.706 1703.529 1704.235
0.500 0.516 1724.184 1726.097
0.300 0.306 1744.840 1745.513
0.200 0.212 1757.335 1758.908
0.100 0.086 1774.666 1772.914
0.050 0.044 1788.976 1787.428
0.020 0.018 1805.079 1803.239
0.010 0.006 1815.816 1813.613
Bayesian (BIC)
Mean 1764.099
Std Dev 39.390
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.992 1672.467 1673.803
0.980 0.984 1683.205 1687.518
0.950 0.944 1699.307 1698.927
0.900 0.890 1713.617 1711.010
0.800 0.806 1730.949 1731.816
0.700 0.706 1743.443 1744.150
0.500 0.516 1764.099 1766.011
0.300 0.306 1784.755 1785.428
0.200 0.212 1797.249 1798.823
0.100 0.086 1814.581 1812.828
0.050 0.044 1828.891 1827.343
0.020 0.018 1844.993 1843.154
0.010 0.006 1855.731 1853.527
Sample-Size Adjusted BIC (n* = (n + 2) / 24)
Mean 1732.368
Std Dev 39.390
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.992 1640.736 1642.072
0.980 0.984 1651.474 1655.788
0.950 0.944 1667.576 1667.197
0.900 0.890 1681.887 1679.280
0.800 0.806 1699.218 1700.086
0.700 0.706 1711.713 1712.419
0.500 0.516 1732.368 1734.281
0.300 0.306 1753.024 1753.697
0.200 0.212 1765.519 1767.092
0.100 0.086 1782.850 1781.098
0.050 0.044 1797.160 1795.612
0.020 0.018 1813.263 1811.423
0.010 0.006 1824.000 1821.797
MODEL RESULTS
ESTIMATES S. E. M. S. E. 95% % Sig
Population Average Std. Dev. Average Cover Coeff
GAMMA1 ON
Z 0.200 0.2010 0.0775 0.0771 0.0060 0.950 0.744
Y ON
M 0.500 0.5007 0.0524 0.0494 0.0027 0.922 1.000
X 0.200 0.2056 0.0783 0.0784 0.0061 0.938 0.754
Z 0.100 0.0963 0.0470 0.0433 0.0022 0.926 0.604
M ON
Z 0.300 0.2999 0.0531 0.0545 0.0028 0.964 1.000
Intercepts
Y 0.000 -0.0017 0.0527 0.0522 0.0028 0.934 0.066
M 0.000 -0.0008 0.0543 0.0545 0.0029 0.946 0.054
GAMMA1 0.300 0.3010 0.0776 0.0770 0.0060 0.962 0.978
Residual Variances
Y 0.500 0.4938 0.0341 0.0347 0.0012 0.928 1.000
M 0.500 0.4940 0.0331 0.0346 0.0011 0.950 1.000
GAMMA1 0.000 0.0000 0.0000 0.0000 0.0000 1.000 0.000
New/Additional Parameters
INDAVG 0.150 0.1505 0.0417 0.0416 0.0017 0.956 0.974
INDLOW 0.050 0.0497 0.0546 0.0548 0.0030 0.958 0.138
INDHIGH 0.250 0.2514 0.0628 0.0603 0.0039 0.928 0.988
QUALITY OF NUMERICAL RESULTS
Average Condition Number for the Information Matrix 0.587E-03
(ratio of smallest to largest eigenvalue)
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
NU
Y M X Z
________ ________ ________ ________
1 0 0 0 0
LAMBDA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
Y 0 0 0 0 0
M 0 0 0 0 0
X 0 0 0 0 0
Z 0 0 0 0 0
THETA
Y M X Z
________ ________ ________ ________
Y 0
M 0 0
X 0 0 0
Z 0 0 0 0
ALPHA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
1 1 2 3 0 0
BETA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
GAMMA1 0 0 0 0 4
Y 0 0 5 6 7
M 0 0 0 0 8
X 0 0 0 0 0
Z 0 0 0 0 0
PSI
GAMMA1 Y M X Z
________ ________ ________ ________ ________
GAMMA1 0
Y 0 9
M 0 0 10
X 0 0 0 0
Z 0 0 0 0 0
PARAMETER SPECIFICATION FOR THE ADDITIONAL PARAMETERS
NEW/ADDITIONAL PARAMETERS
INDAVG INDLOW INDHIGH
________ ________ ________
1 11 12 13
STARTING VALUES
NU
Y M X Z
________ ________ ________ ________
1 0.000 0.000 0.000 0.000
LAMBDA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
Y 0.000 1.000 0.000 0.000 0.000
M 0.000 0.000 1.000 0.000 0.000
X 0.000 0.000 0.000 1.000 0.000
Z 0.000 0.000 0.000 0.000 1.000
THETA
Y M X Z
________ ________ ________ ________
Y 0.000
M 0.000 0.000
X 0.000 0.000 0.000
Z 0.000 0.000 0.000 0.000
ALPHA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
1 0.300 0.000 0.000 0.000 0.000
BETA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
GAMMA1 0.000 0.000 0.000 0.000 0.200
Y 0.000 0.000 0.500 0.200 0.100
M 0.000 0.000 0.000 0.000 0.300
X 0.000 0.000 0.000 0.000 0.000
Z 0.000 0.000 0.000 0.000 0.000
PSI
GAMMA1 Y M X Z
________ ________ ________ ________ ________
GAMMA1 0.000
Y 0.000 0.500
M 0.000 0.000 0.500
X 0.000 0.000 0.000 0.500
Z 0.000 0.000 0.000 0.000 0.500
STARTING VALUES FOR THE ADDITIONAL PARAMETERS
NEW/ADDITIONAL PARAMETERS
INDAVG INDLOW INDHIGH
________ ________ ________
1 0.150 0.050 0.250
POPULATION VALUES
NU
Y M X Z
________ ________ ________ ________
1 0.000 0.000 0.000 0.000
LAMBDA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
Y 0.000 1.000 0.000 0.000 0.000
M 0.000 0.000 1.000 0.000 0.000
X 0.000 0.000 0.000 1.000 0.000
Z 0.000 0.000 0.000 0.000 1.000
THETA
Y M X Z
________ ________ ________ ________
Y 0.000
M 0.000 0.000
X 0.000 0.000 0.000
Z 0.000 0.000 0.000 0.000
ALPHA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
1 0.300 0.000 0.000 0.000 0.000
BETA
GAMMA1 Y M X Z
________ ________ ________ ________ ________
GAMMA1 0.000 0.000 0.000 0.000 0.200
Y 0.000 0.000 0.500 0.200 0.100
M 0.000 0.000 0.000 0.000 0.300
X 0.000 0.000 0.000 0.000 0.000
Z 0.000 0.000 0.000 0.000 0.000
PSI
GAMMA1 Y M X Z
________ ________ ________ ________ ________
GAMMA1 0.000
Y 0.000 0.500
M 0.000 0.000 0.500
X 0.000 0.000 0.000 1.000
Z 0.000 0.000 0.000 0.500 1.000
SAVEDATA INFORMATION
Order of variables
Y
M
X
Z
Save file
xzrep*.dat
Save file format Free
Save file record length 10000
DIAGRAM INFORMATION
Mplus diagrams are currently not available for Monte Carlo analysis.
No diagram output was produced.
Beginning Time: 14:59:14
Ending Time: 15:00:08
Elapsed Time: 00:00:54
MUTHEN & MUTHEN
3463 Stoner Ave.
Los Angeles, CA 90066
Tel: (310) 391-9971
Fax: (310) 391-8971
Web: www.StatModel.com
Support: Support@StatModel.com
Copyright (c) 1998-2015 Muthen & Muthen