Mplus VERSION 7.4
MUTHEN & MUTHEN
06/01/2016 7:26 PM
INPUT INSTRUCTIONS
data:
file = n500mix00xreplist.dat;
type = montecarlo;
Variable:
names = y m tx x;
missing = all(999);
usev = y m tx x missing;
categorical = missing;
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 500
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
n500mix00xreplist.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 319 2 181
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.638
M 0.638 1.000
TX 0.638 1.000 1.000
X 0.638 1.000 1.000 1.000
UNIVARIATE PROPORTIONS FOR CATEGORICAL VARIABLES FOR THE FIRST REPLICATION
MISSING
Category 1 0.638
Category 2 0.362
SAMPLE STATISTICS
NOTE: These are average results over 500 data sets.
SAMPLE STATISTICS
Means
Y M TX X
________ ________ ________ ________
1 -0.249 -0.249 0.500 0.002
Covariances
Y M TX X
________ ________ ________ ________
Y 1.563
M 0.812 0.812
TX -0.126 -0.124 0.249
X 0.500 0.501 0.001 1.001
Correlations
Y M TX X
________ ________ ________ ________
Y 1.000
M 0.721 1.000
TX -0.202 -0.277 1.000
X 0.399 0.556 0.001 1.000
MODEL FIT INFORMATION
Number of Free Parameters 11
Loglikelihood
H0 Value
Mean -1180.425
Std Dev 23.173
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.994 -1234.332 -1233.706
0.980 0.982 -1228.015 -1226.821
0.950 0.948 -1218.542 -1219.052
0.900 0.900 -1210.123 -1210.605
0.800 0.790 -1199.927 -1200.597
0.700 0.690 -1192.577 -1193.410
0.500 0.486 -1180.425 -1181.610
0.300 0.304 -1168.273 -1167.798
0.200 0.206 -1160.922 -1160.401
0.100 0.110 -1150.726 -1149.800
0.050 0.050 -1142.307 -1142.436
0.020 0.022 -1132.834 -1132.451
0.010 0.008 -1126.517 -1127.971
Information Criteria
Akaike (AIC)
Mean 2382.849
Std Dev 46.346
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.992 2275.034 2276.427
0.980 0.978 2287.668 2284.023
0.950 0.950 2306.614 2306.515
0.900 0.890 2323.452 2318.850
0.800 0.794 2343.844 2340.682
0.700 0.696 2358.545 2357.567
0.500 0.514 2382.849 2384.913
0.300 0.310 2407.153 2408.645
0.200 0.210 2421.854 2422.991
0.100 0.100 2442.246 2442.188
0.050 0.052 2459.084 2460.048
0.020 0.018 2478.030 2474.993
0.010 0.006 2490.664 2487.779
Bayesian (BIC)
Mean 2429.210
Std Dev 46.346
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.992 2321.395 2322.788
0.980 0.978 2334.029 2330.384
0.950 0.950 2352.975 2352.876
0.900 0.890 2369.812 2365.211
0.800 0.794 2390.205 2387.042
0.700 0.696 2404.906 2403.928
0.500 0.514 2429.210 2431.274
0.300 0.310 2453.514 2455.006
0.200 0.210 2468.215 2469.352
0.100 0.100 2488.607 2488.548
0.050 0.052 2505.445 2506.408
0.020 0.018 2524.391 2521.353
0.010 0.006 2537.025 2534.140
Sample-Size Adjusted BIC (n* = (n + 2) / 24)
Mean 2394.295
Std Dev 46.346
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.992 2286.480 2287.873
0.980 0.978 2299.114 2295.469
0.950 0.950 2318.060 2317.961
0.900 0.890 2334.898 2330.296
0.800 0.794 2355.290 2352.128
0.700 0.696 2369.991 2369.013
0.500 0.514 2394.295 2396.359
0.300 0.310 2418.599 2420.091
0.200 0.210 2433.300 2434.437
0.100 0.100 2453.692 2453.634
0.050 0.052 2470.530 2471.494
0.020 0.018 2489.476 2486.439
0.010 0.006 2502.110 2499.225
MODEL RESULTS
ESTIMATES S. E. M. S. E. 95% % Sig
Population Average Std. Dev. Average Cover Coeff
Y ON
M 1.000 1.3233 0.0783 0.0763 0.1106 0.004 1.000
TX 0.000 -0.0125 0.0994 0.1048 0.0100 0.964 0.036
X 0.000 0.0020 0.0617 0.0606 0.0038 0.950 0.050
M ON
TX -0.500 -0.5000 0.0653 0.0631 0.0043 0.936 1.000
X 0.500 0.5008 0.0315 0.0315 0.0010 0.952 1.000
MISSING ON
Y 2.000 1.9290 0.2311 0.2287 0.0583 0.908 1.000
Intercepts
Y 0.000 0.4810 0.0823 0.0860 0.2381 0.000 1.000
M 0.000 0.0010 0.0462 0.0447 0.0021 0.944 0.056
Thresholds
MISSING$1 0.750 1.6704 0.2718 0.2643 0.9210 0.010 1.000
Residual Variances
Y 0.750 0.9234 0.0809 0.0826 0.0366 0.448 1.000
M 0.500 0.4966 0.0316 0.0311 0.0010 0.946 1.000
QUALITY OF NUMERICAL RESULTS
Average Condition Number for the Information Matrix 0.315E-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.6742 0.1284 0.1306 0.0468 0.712 1.000
Tot indirect -0.500 -0.6617 0.0946 0.0912 0.0351 0.560 1.000
Specific indirect
Y
M
TX -0.500 -0.6617 0.0946 0.0912 0.0351 0.560 1.000
Direct
Y
TX 0.000 -0.0125 0.0994 0.1048 0.0100 0.964 0.036
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:26:29
Ending Time: 19:49:36
Elapsed Time: 00:23:07
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