Mplus VERSION 7.4
MUTHEN & MUTHEN
06/01/2016 7:09 PM
INPUT INSTRUCTIONS
data:
file = n500replist.dat;
type = montecarlo;
Variable:
names = y m tx x missing;
categorical = missing;
missing = all(999);
define:
if(y EQ _MISSING) THEN missing = 1 ELSE missing = 0;
Analysis:
estimator = mlr;
integration = montecarlo;
mconv = 0.00001;
Model:
y on m*1 tx*0 x*0;
y*.75;
m on tx*-.5 x*.5;
m*.5;
[y-m*0];
missing on y*2;
[missing$1*.75];
model indirect:
y IND tx;
INPUT READING TERMINATED NORMALLY
SUMMARY OF ANALYSIS
Number of groups 1
Average number of observations 250
Number of replications
Requested 500
Completed 500
Number of dependent variables 3
Number of independent variables 2
Number of continuous latent variables 0
Observed dependent variables
Continuous
Y M
Binary and ordered categorical (ordinal)
MISSING
Observed independent variables
TX X
Estimator MLR
Information matrix OBSERVED
Optimization Specifications for the Quasi-Newton Algorithm for
Continuous Outcomes
Maximum number of iterations 100
Convergence criterion 0.100D-05
Optimization Specifications for the EM Algorithm
Maximum number of iterations 500
Convergence criteria
Loglikelihood change 0.100D-02
Relative loglikelihood change 0.100D-05
Derivative 0.100D-04
Optimization Specifications for the M step of the EM Algorithm for
Categorical Latent variables
Number of M step iterations 1
M step convergence criterion 0.100D-02
Basis for M step termination ITERATION
Optimization Specifications for the M step of the EM Algorithm for
Censored, Binary or Ordered Categorical (Ordinal), Unordered
Categorical (Nominal) and Count Outcomes
Number of M step iterations 1
M step convergence criterion 0.100D-02
Basis for M step termination ITERATION
Maximum value for logit thresholds 15
Minimum value for logit thresholds -15
Minimum expected cell size for chi-square 0.100D-01
Maximum number of iterations for H1 2000
Convergence criterion for H1 0.100D-03
Optimization algorithm EMA
Integration Specifications
Type MONTECARLO
Number of integration points 250
Dimensions of numerical integration 1
Adaptive quadrature ON
Monte Carlo integration seed 0
Link LOGIT
Cholesky OFF
Input data file(s)
Multiple data files from
n500replist.dat
Input data format FREE
SUMMARY OF DATA FOR THE FIRST REPLICATION
Number of missing data patterns 2
Number of y missing data patterns 2
Number of u missing data patterns 1
SUMMARY OF MISSING DATA PATTERNS FOR THE FIRST REPLICATION
MISSING DATA PATTERNS FOR Y (x = not missing)
1 2
Y x
M x x
TX x x
X x x
MISSING DATA PATTERN FREQUENCIES FOR Y
Pattern Frequency Pattern Frequency
1 154 2 96
COVARIANCE COVERAGE OF DATA FOR THE FIRST REPLICATION
Minimum covariance coverage value 0.100
PROPORTION OF DATA PRESENT FOR Y
Covariance Coverage
Y M TX X
________ ________ ________ ________
Y 0.616
M 0.616 1.000
TX 0.616 1.000 1.000
X 0.616 1.000 1.000 1.000
UNIVARIATE PROPORTIONS FOR CATEGORICAL VARIABLES FOR THE FIRST REPLICATION
MISSING
Category 1 0.616
Category 2 0.384
SAMPLE STATISTICS
NOTE: These are average results over 500 data sets.
SAMPLE STATISTICS
Means
Y M TX X
________ ________ ________ ________
1 -0.585 -0.246 0.501 0.004
Covariances
Y M TX X
________ ________ ________ ________
Y 1.087
M 0.635 0.811
TX -0.098 -0.124 0.249
X 0.394 0.500 0.002 1.002
Correlations
Y M TX X
________ ________ ________ ________
Y 1.000
M 0.676 1.000
TX -0.189 -0.275 1.000
X 0.378 0.555 0.004 1.000
MODEL FIT INFORMATION
Number of Free Parameters 11
Loglikelihood
H0 Value
Mean -582.993
Std Dev 16.122
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.992 -620.497 -618.922
0.980 0.986 -616.102 -614.406
0.950 0.956 -609.512 -609.314
0.900 0.906 -603.655 -603.435
0.800 0.810 -596.561 -596.198
0.700 0.694 -591.447 -591.788
0.500 0.480 -582.993 -584.573
0.300 0.286 -574.539 -575.190
0.200 0.194 -569.425 -570.290
0.100 0.104 -562.332 -561.953
0.050 0.058 -556.475 -555.627
0.020 0.024 -549.884 -549.227
0.010 0.014 -545.489 -543.909
Information Criteria
Akaike (AIC)
Mean 1187.986
Std Dev 32.243
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.986 1112.979 1105.357
0.980 0.976 1121.768 1119.570
0.950 0.942 1134.949 1131.669
0.900 0.896 1146.663 1144.301
0.800 0.806 1160.850 1161.304
0.700 0.714 1171.078 1172.239
0.500 0.520 1187.986 1191.037
0.300 0.306 1204.895 1205.414
0.200 0.190 1215.122 1214.178
0.100 0.094 1229.309 1228.442
0.050 0.044 1241.023 1239.797
0.020 0.014 1254.204 1250.676
0.010 0.008 1262.994 1259.695
Bayesian (BIC)
Mean 1226.722
Std Dev 32.243
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.986 1151.715 1144.093
0.980 0.976 1160.504 1158.306
0.950 0.942 1173.685 1170.405
0.900 0.896 1185.399 1183.037
0.800 0.806 1199.587 1200.040
0.700 0.714 1209.814 1210.975
0.500 0.520 1226.722 1229.773
0.300 0.306 1243.631 1244.150
0.200 0.190 1253.858 1252.914
0.100 0.094 1268.045 1267.178
0.050 0.044 1279.759 1278.533
0.020 0.014 1292.941 1289.412
0.010 0.008 1301.730 1298.431
Sample-Size Adjusted BIC (n* = (n + 2) / 24)
Mean 1191.852
Std Dev 32.243
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.986 1116.844 1109.222
0.980 0.976 1125.633 1123.435
0.950 0.942 1138.815 1135.534
0.900 0.896 1150.529 1148.166
0.800 0.806 1164.716 1165.170
0.700 0.714 1174.943 1176.104
0.500 0.520 1191.852 1194.902
0.300 0.306 1208.760 1209.279
0.200 0.190 1218.987 1218.043
0.100 0.094 1233.175 1232.308
0.050 0.044 1244.889 1243.662
0.020 0.014 1258.070 1254.541
0.010 0.008 1266.859 1263.561
MODEL RESULTS
ESTIMATES S. E. M. S. E. 95% % Sig
Population Average Std. Dev. Average Cover Coeff
Y ON
M 1.000 0.9973 0.0979 0.0961 0.0096 0.946 1.000
TX 0.000 -0.0082 0.1368 0.1330 0.0188 0.928 0.072
X 0.000 0.0043 0.0787 0.0768 0.0062 0.944 0.056
M ON
TX -0.500 -0.5011 0.0882 0.0892 0.0078 0.952 1.000
X 0.500 0.5006 0.0421 0.0445 0.0018 0.964 1.000
MISSING ON
Y 2.000 2.1035 0.5119 0.4737 0.2722 0.946 0.998
Intercepts
Y 0.000 0.0025 0.1160 0.1119 0.0134 0.932 0.068
M 0.000 0.0029 0.0637 0.0631 0.0041 0.948 0.052
Thresholds
MISSING$1 0.750 0.8043 0.3495 0.3228 0.1249 0.946 0.804
Residual Variances
Y 0.750 0.7404 0.1043 0.1061 0.0109 0.940 1.000
M 0.500 0.4949 0.0440 0.0438 0.0020 0.944 1.000
QUALITY OF NUMERICAL RESULTS
Average Condition Number for the Information Matrix 0.476E-02
(ratio of smallest to largest eigenvalue)
TOTAL, TOTAL INDIRECT, SPECIFIC INDIRECT, AND DIRECT EFFECTS
ESTIMATES S. E. M. S. E. 95% % Sig
Population Average Std. Dev. Average Cover Coeff
Effects from TX to Y
Total -0.500 -0.5083 0.1545 0.1557 0.0239 0.946 0.916
Tot indirect -0.500 -0.5001 0.1034 0.1016 0.0107 0.940 1.000
Specific indirect
Y
M
TX -0.500 -0.5001 0.1034 0.1016 0.0107 0.940 1.000
Direct
Y
TX 0.000 -0.0082 0.1368 0.1330 0.0188 0.928 0.072
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
TAU
MISSING$
________
1 11
NU
MISSING Y M TX X
________ ________ ________ ________ ________
1 0 0 0 0 0
LAMBDA
MISSING Y M TX X
________ ________ ________ ________ ________
MISSING 0 0 0 0 0
Y 0 0 0 0 0
M 0 0 0 0 0
TX 0 0 0 0 0
X 0 0 0 0 0
THETA
MISSING Y M TX X
________ ________ ________ ________ ________
MISSING 0
Y 0 0
M 0 0 0
TX 0 0 0 0
X 0 0 0 0 0
ALPHA
MISSING Y M TX X
________ ________ ________ ________ ________
1 0 1 2 0 0
BETA
MISSING Y M TX X
________ ________ ________ ________ ________
MISSING 0 3 0 0 0
Y 0 0 4 5 6
M 0 0 0 7 8
TX 0 0 0 0 0
X 0 0 0 0 0
PSI
MISSING Y M TX X
________ ________ ________ ________ ________
MISSING 0
Y 0 9
M 0 0 10
TX 0 0 0 0
X 0 0 0 0 0
STARTING VALUES
TAU
MISSING$
________
1 0.750
NU
MISSING Y M TX X
________ ________ ________ ________ ________
1 0.000 0.000 0.000 0.000 0.000
LAMBDA
MISSING Y M TX X
________ ________ ________ ________ ________
MISSING 1.000 0.000 0.000 0.000 0.000
Y 0.000 1.000 0.000 0.000 0.000
M 0.000 0.000 1.000 0.000 0.000
TX 0.000 0.000 0.000 1.000 0.000
X 0.000 0.000 0.000 0.000 1.000
THETA
MISSING Y M TX X
________ ________ ________ ________ ________
MISSING 0.000
Y 0.000 0.000
M 0.000 0.000 0.000
TX 0.000 0.000 0.000 0.000
X 0.000 0.000 0.000 0.000 0.000
ALPHA
MISSING Y M TX X
________ ________ ________ ________ ________
1 0.000 0.000 0.000 0.000 0.000
BETA
MISSING Y M TX X
________ ________ ________ ________ ________
MISSING 0.000 2.000 0.000 0.000 0.000
Y 0.000 0.000 1.000 0.000 0.000
M 0.000 0.000 0.000 -0.500 0.500
TX 0.000 0.000 0.000 0.000 0.000
X 0.000 0.000 0.000 0.000 0.000
PSI
MISSING Y M TX X
________ ________ ________ ________ ________
MISSING 1.000
Y 0.000 0.750
M 0.000 0.000 0.500
TX 0.000 0.000 0.000 0.500
X 0.000 0.000 0.000 0.000 0.500
DIAGRAM INFORMATION
Mplus diagrams are currently not available for Monte Carlo analysis.
No diagram output was produced.
Beginning Time: 19:09:13
Ending Time: 19:20:29
Elapsed Time: 00:11:16
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