Mplus VERSION 7.4
MUTHEN & MUTHEN
06/01/2016 7:05 PM
INPUT INSTRUCTIONS
TITLE: Example: MAR simulation
DATA:
FILE = n500mix00xreplist.dat;
TYPE = MONTECARLO;
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 500
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
Maximum number of iterations for H1 2000
Convergence criterion for H1 0.100D-03
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
SUMMARY OF MISSING DATA PATTERNS FOR THE FIRST REPLICATION
MISSING DATA PATTERNS (x = not missing)
1 2
Y x
M x x
TX x x
X x x
MISSING DATA PATTERN FREQUENCIES
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
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
SAMPLE STATISTICS
NOTE: These are average results over 500 data sets.
ESTIMATED 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 9
Loglikelihood
H0 Value
Mean -950.329
Std Dev 22.870
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.988 -1003.531 -1005.435
0.980 0.986 -997.297 -995.047
0.950 0.952 -987.948 -987.906
0.900 0.900 -979.639 -980.898
0.800 0.784 -969.576 -970.172
0.700 0.680 -962.322 -962.970
0.500 0.516 -950.329 -949.698
0.300 0.294 -938.336 -938.665
0.200 0.212 -931.082 -929.737
0.100 0.116 -921.019 -920.144
0.050 0.052 -912.710 -912.400
0.020 0.022 -903.361 -903.085
0.010 0.004 -897.127 -900.815
H1 Value
Mean -950.329
Std Dev 22.870
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.988 -1003.531 -1005.435
0.980 0.986 -997.297 -995.047
0.950 0.952 -987.948 -987.906
0.900 0.900 -979.639 -980.898
0.800 0.784 -969.576 -970.172
0.700 0.680 -962.322 -962.970
0.500 0.516 -950.329 -949.698
0.300 0.294 -938.336 -938.665
0.200 0.212 -931.082 -929.737
0.100 0.116 -921.019 -920.144
0.050 0.052 -912.710 -912.400
0.020 0.022 -903.361 -903.085
0.010 0.004 -897.127 -900.815
Information Criteria
Akaike (AIC)
Mean 1918.658
Std Dev 45.740
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.996 1812.254 1816.891
0.980 0.978 1824.722 1823.844
0.950 0.948 1843.421 1840.384
0.900 0.884 1860.038 1857.773
0.800 0.788 1880.163 1876.765
0.700 0.706 1894.672 1895.083
0.500 0.484 1918.658 1917.392
0.300 0.320 1942.644 1943.929
0.200 0.216 1957.152 1957.856
0.100 0.100 1977.278 1976.789
0.050 0.048 1993.895 1992.245
0.020 0.014 2012.594 2007.567
0.010 0.012 2025.062 2027.424
Bayesian (BIC)
Mean 1956.589
Std Dev 45.740
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.996 1850.185 1854.823
0.980 0.978 1862.654 1861.775
0.950 0.948 1881.352 1878.316
0.900 0.884 1897.969 1895.704
0.800 0.788 1918.095 1914.697
0.700 0.706 1932.604 1933.015
0.500 0.484 1956.589 1955.323
0.300 0.320 1980.575 1981.860
0.200 0.216 1995.084 1995.787
0.100 0.100 2015.209 2014.721
0.050 0.048 2031.827 2030.177
0.020 0.014 2050.525 2045.498
0.010 0.012 2062.994 2065.355
Sample-Size Adjusted BIC (n* = (n + 2) / 24)
Mean 1928.023
Std Dev 45.740
Number of successful computations 500
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.996 1821.619 1826.256
0.980 0.978 1834.087 1833.209
0.950 0.948 1852.786 1849.749
0.900 0.884 1869.403 1867.138
0.800 0.788 1889.528 1886.130
0.700 0.706 1904.037 1904.448
0.500 0.484 1928.023 1926.756
0.300 0.320 1952.009 1953.294
0.200 0.216 1966.517 1967.221
0.100 0.100 1986.643 1986.154
0.050 0.048 2003.260 2001.610
0.020 0.014 2021.958 2016.932
0.010 0.012 2034.427 2036.789
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.986 1.000 1.000
0.980 0.962 1.000 1.000
0.950 0.940 1.000 1.000
0.900 0.904 1.000 1.000
0.800 0.802 1.000 1.000
0.700 0.696 1.000 1.000
0.500 0.532 1.000 1.000
0.300 0.306 1.000 1.000
0.200 0.198 1.000 1.000
0.100 0.090 1.000 1.000
0.050 0.046 1.000 1.000
0.020 0.012 1.000 1.000
0.010 0.006 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 1.000 0.000 0.000
0.980 1.000 0.000 0.000
0.950 0.996 0.000 0.000
0.900 0.970 0.000 0.000
0.800 0.812 0.000 0.000
0.700 0.650 0.000 0.000
0.500 0.414 0.000 0.000
0.300 0.256 0.000 0.000
0.200 0.176 0.000 0.000
0.100 0.100 0.000 0.000
0.050 0.080 0.000 0.000
0.020 0.048 0.000 0.000
0.010 0.034 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.0759 0.0756 0.0058 0.946 1.000
TX 0.000 -0.0082 0.0977 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.5001 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
Intercepts
Y 0.000 0.0036 0.0753 0.0769 0.0057 0.958 0.042
M 0.000 0.0010 0.0462 0.0447 0.0021 0.944 0.056
Residual Variances
Y 0.750 0.7437 0.0556 0.0579 0.0031 0.952 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.665E-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.516 0.053
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.250
X 0.000 0.000 0.034 1.002
DIAGRAM INFORMATION
Mplus diagrams are currently not available for Monte Carlo analysis.
No diagram output was produced.
Beginning Time: 19:05:57
Ending Time: 19:06:01
Elapsed Time: 00:00:04
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