Mplus VERSION 7.4
MUTHEN & MUTHEN
06/01/2016 7:06 PM
INPUT INSTRUCTIONS
TITLE: Example: MAR simulation
DATA:
FILE = n500mix00xreplist.dat;
TYPE = MONTECARLO;
LISTWISE = ON;
VARIABLE:
NAMES = y m tx x;
USEVARIABLES = y-x;
MISSING = ALL(999);
ANALYSIS:
ESTIMATOR = MLR;
MODEL:
y on m*1 tx*0 x*0;
y*.75;
m on tx*-.5 x*.5;
m*.5;
[y-m*0];
INPUT READING TERMINATED NORMALLY
Example: MAR simulation
SUMMARY OF ANALYSIS
Number of groups 1
Average number of observations 327
Number of replications
Requested 500
Completed 500
Number of dependent variables 2
Number of independent variables 2
Number of continuous latent variables 0
Observed dependent variables
Continuous
Y M
Observed independent variables
TX X
Estimator MLR
Information matrix OBSERVED
Maximum number of iterations 1000
Convergence criterion 0.500D-04
Maximum number of steepest descent iterations 20
Input data file(s)
Multiple data files from
n500mix00xreplist.dat
Input data format FREE
SAMPLE STATISTICS
NOTE: These are average results over 500 data sets.
SAMPLE STATISTICS
Means
Y M TX X
________ ________ ________ ________
1 -0.617 -0.616 0.558 -0.225
Covariances
Y M TX X
________ ________ ________ ________
Y 1.324
M 0.574 0.574
TX -0.090 -0.088 0.246
X 0.353 0.355 0.023 0.907
Correlations
Y M TX X
________ ________ ________ ________
Y 1.000
M 0.658 1.000
TX -0.158 -0.235 1.000
X 0.322 0.492 0.048 1.000
MODEL FIT INFORMATION
Number of Free Parameters 9
Loglikelihood
H0 Value
Mean -728.408
Std Dev 29.909
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.988 -797.984 -799.208
0.980 0.982 -789.831 -789.159
0.950 0.946 -777.605 -778.958
0.900 0.906 -766.739 -766.244
0.800 0.816 -753.579 -751.914
0.700 0.694 -744.092 -745.226
0.500 0.490 -728.408 -729.763
0.300 0.322 -712.724 -711.069
0.200 0.204 -703.237 -703.067
0.100 0.114 -690.077 -688.868
0.050 0.048 -679.212 -679.873
0.020 0.012 -666.985 -670.761
0.010 0.010 -658.832 -665.437
H1 Value
Mean -728.408
Std Dev 29.909
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.988 -797.984 -799.208
0.980 0.982 -789.831 -789.159
0.950 0.946 -777.605 -778.958
0.900 0.906 -766.739 -766.244
0.800 0.816 -753.579 -751.914
0.700 0.694 -744.092 -745.226
0.500 0.490 -728.408 -729.763
0.300 0.322 -712.724 -711.069
0.200 0.204 -703.237 -703.067
0.100 0.114 -690.077 -688.868
0.050 0.048 -679.212 -679.873
0.020 0.012 -666.985 -670.761
0.010 0.010 -658.832 -665.437
Information Criteria
Akaike (AIC)
Mean 1474.816
Std Dev 59.817
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.990 1335.664 1334.028
0.980 0.988 1351.970 1359.278
0.950 0.952 1376.423 1376.763
0.900 0.886 1398.155 1394.045
0.800 0.796 1424.474 1423.992
0.700 0.678 1443.448 1439.995
0.500 0.510 1474.816 1476.512
0.300 0.306 1506.184 1507.825
0.200 0.184 1525.158 1521.571
0.100 0.094 1551.478 1550.183
0.050 0.054 1573.209 1573.846
0.020 0.018 1597.663 1593.280
0.010 0.012 1613.969 1616.390
Bayesian (BIC)
Mean 1508.948
Std Dev 60.048
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.990 1369.258 1367.181
0.980 0.988 1385.627 1392.907
0.950 0.952 1410.175 1410.622
0.900 0.886 1431.990 1427.616
0.800 0.794 1458.411 1457.650
0.700 0.680 1477.459 1473.803
0.500 0.510 1508.948 1510.483
0.300 0.308 1540.437 1542.285
0.200 0.182 1559.485 1555.872
0.100 0.094 1585.906 1584.670
0.050 0.054 1607.721 1608.438
0.020 0.018 1632.269 1627.794
0.010 0.012 1648.638 1651.086
Sample-Size Adjusted BIC (n* = (n + 2) / 24)
Mean 1480.400
Std Dev 60.047
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.990 1340.713 1338.639
0.980 0.988 1357.082 1364.362
0.950 0.952 1381.629 1382.076
0.900 0.886 1403.444 1399.072
0.800 0.794 1429.865 1429.106
0.700 0.680 1448.912 1445.255
0.500 0.510 1480.400 1481.936
0.300 0.308 1511.889 1513.736
0.200 0.182 1530.936 1527.323
0.100 0.094 1557.356 1556.120
0.050 0.054 1579.171 1579.887
0.020 0.018 1603.718 1599.244
0.010 0.012 1620.087 1622.535
Chi-Square Test of Model Fit
Degrees of freedom 0
Mean 0.000
Std Dev 0.000
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 1.000 0.000 0.000
0.980 1.000 0.000 0.000
0.950 1.000 0.000 0.000
0.900 1.000 0.000 0.000
0.800 1.000 0.000 0.000
0.700 1.000 0.000 0.000
0.500 1.000 0.000 0.000
0.300 1.000 0.000 0.000
0.200 1.000 0.000 0.000
0.100 1.000 0.000 0.000
0.050 1.000 0.000 0.000
0.020 1.000 0.000 0.000
0.010 1.000 0.000 0.000
RMSEA (Root Mean Square Error Of Approximation)
Mean 0.000
Std Dev 0.000
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.000 0.000 0.000
0.980 0.000 0.000 0.000
0.950 0.000 0.000 0.000
0.900 0.000 0.000 0.000
0.800 0.000 0.000 0.000
0.700 0.000 0.000 0.000
0.500 0.000 0.000 0.000
0.300 0.000 0.000 0.000
0.200 0.000 0.000 0.000
0.100 0.000 0.000 0.000
0.050 0.000 0.000 0.000
0.020 0.000 0.000 0.000
0.010 0.000 0.000 0.000
CFI/TLI
CFI
Mean 1.000
Std Dev 0.000
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.990 1.000 1.000
0.980 0.980 1.000 1.000
0.950 0.964 1.000 1.000
0.900 0.914 1.000 1.000
0.800 0.828 1.000 1.000
0.700 0.718 1.000 1.000
0.500 0.492 1.000 1.000
0.300 0.234 1.000 1.000
0.200 0.144 1.000 1.000
0.100 0.080 1.000 1.000
0.050 0.060 1.000 1.000
0.020 0.050 1.000 1.000
0.010 0.044 1.000 1.000
TLI
Mean 1.000
Std Dev 0.000
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.000 1.000 1.000
0.980 0.000 1.000 1.000
0.950 0.000 1.000 1.000
0.900 0.000 1.000 1.000
0.800 0.000 1.000 1.000
0.700 0.000 1.000 1.000
0.500 0.000 1.000 1.000
0.300 0.000 1.000 1.000
0.200 0.000 1.000 1.000
0.100 0.000 1.000 1.000
0.050 0.000 1.000 1.000
0.020 0.000 1.000 1.000
0.010 0.000 1.000 1.000
SRMR (Standardized Root Mean Square Residual)
Mean 0.000
Std Dev 0.000
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.972 0.000 0.000
0.980 0.966 0.000 0.000
0.950 0.946 0.000 0.000
0.900 0.912 0.000 0.000
0.800 0.860 0.000 0.000
0.700 0.776 0.000 0.000
0.500 0.524 0.000 0.000
0.300 0.268 0.000 0.000
0.200 0.160 0.000 0.000
0.100 0.070 0.000 0.000
0.050 0.036 0.000 0.000
0.020 0.016 0.000 0.000
0.010 0.012 0.000 0.000
MODEL RESULTS
ESTIMATES S. E. M. S. E. 95% % Sig
Population Average Std. Dev. Average Cover Coeff
Y ON
M 1.000 0.9995 0.0760 0.0756 0.0058 0.946 1.000
TX 0.000 -0.0082 0.0978 0.1009 0.0096 0.966 0.034
X 0.000 -0.0009 0.0605 0.0586 0.0037 0.936 0.064
M ON
TX -0.500 -0.3965 0.0691 0.0698 0.0155 0.684 1.000
X 0.500 0.4013 0.0376 0.0367 0.0111 0.256 1.000
Intercepts
Y 0.000 0.0036 0.0753 0.0769 0.0057 0.958 0.042
M 0.000 -0.3046 0.0509 0.0521 0.0954 0.000 1.000
Residual Variances
Y 0.750 0.7437 0.0556 0.0579 0.0031 0.952 1.000
M 0.500 0.3940 0.0330 0.0307 0.0123 0.120 1.000
QUALITY OF NUMERICAL RESULTS
Average Condition Number for the Information Matrix 0.710E-01
(ratio of smallest to largest eigenvalue)
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
NU
Y M TX X
________ ________ ________ ________
1 0 0 0 0
LAMBDA
Y M TX X
________ ________ ________ ________
Y 0 0 0 0
M 0 0 0 0
TX 0 0 0 0
X 0 0 0 0
THETA
Y M TX X
________ ________ ________ ________
Y 0
M 0 0
TX 0 0 0
X 0 0 0 0
ALPHA
Y M TX X
________ ________ ________ ________
1 1 2 0 0
BETA
Y M TX X
________ ________ ________ ________
Y 0 3 4 5
M 0 0 6 7
TX 0 0 0 0
X 0 0 0 0
PSI
Y M TX X
________ ________ ________ ________
Y 8
M 0 9
TX 0 0 0
X 0 0 0 0
STARTING VALUES
NU
Y M TX X
________ ________ ________ ________
1 0.000 0.000 0.000 0.000
LAMBDA
Y M TX X
________ ________ ________ ________
Y 1.000 0.000 0.000 0.000
M 0.000 1.000 0.000 0.000
TX 0.000 0.000 1.000 0.000
X 0.000 0.000 0.000 1.000
THETA
Y M TX X
________ ________ ________ ________
Y 0.000
M 0.000 0.000
TX 0.000 0.000 0.000
X 0.000 0.000 0.000 0.000
ALPHA
Y M TX X
________ ________ ________ ________
1 0.000 0.000 0.542 -0.145
BETA
Y M TX X
________ ________ ________ ________
Y 0.000 1.000 0.000 0.000
M 0.000 0.000 -0.500 0.500
TX 0.000 0.000 0.000 0.000
X 0.000 0.000 0.000 0.000
PSI
Y M TX X
________ ________ ________ ________
Y 0.750
M 0.000 0.500
TX 0.000 0.000 0.248
X 0.000 0.000 0.056 1.005
DIAGRAM INFORMATION
Mplus diagrams are currently not available for Monte Carlo analysis.
No diagram output was produced.
Beginning Time: 19:06:11
Ending Time: 19:06:14
Elapsed Time: 00:00:03
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
Back to examples