Mplus VERSION 8.8
MUTHEN & MUTHEN
04/19/2022 10:24 PM
INPUT INSTRUCTIONS
title: this is an example of a discrete-time
survival analysis
montecarlo:
names = u1-u4 x;
generate = u1-u4(1);
categorical = u1-u4;
missing = u2-u4;
nobs = 500;
nreps = 1;
save = ex6.19.dat;
model population:
[x@0]; x@1;
[u1$1*2 u2$1*1.5 u3$1*1 u4$1*1];
f by u1-u4@1;
f on x*.5;
f@0;
model missing:
[u2-u4@-15]; ! Low logit intercepts give P(missing) = 0
u2-u4 on u1@30; ! if u1=1, later u's are missing with certainty
u3-u4 on u2@30; ! if u2=1, ...
u4 on u3@30; ! if u3=1, ...
analysis:
estimator = mlr;
model:
[u1$1*2 u2$1*1.5 u3$1*1 u4$1*1];
f by u1-u4@1;
f on x*.5;
f@0;
output:
tech8 tech9;
INPUT READING TERMINATED NORMALLY
this is an example of a discrete-time
survival analysis
SUMMARY OF ANALYSIS
Number of groups 1
Number of observations 500
Number of replications
Requested 1
Completed 1
Value of seed 0
Number of dependent variables 4
Number of independent variables 1
Number of continuous latent variables 1
Observed dependent variables
Binary and ordered categorical (ordinal)
U1 U2 U3 U4
Observed independent variables
X
Continuous latent variables
F
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-02
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 STANDARD
Number of integration points 15
Dimensions of numerical integration 0
Adaptive quadrature ON
Link LOGIT
Cholesky ON
SUMMARY OF DATA FOR THE FIRST REPLICATION
Number of missing data patterns 4
Number of y missing data patterns 0
Number of u missing data patterns 4
SUMMARY OF MISSING DATA PATTERNS FOR THE FIRST REPLICATION
MISSING DATA PATTERNS FOR U (x = not missing)
1 2 3 4
U1 x x x x
U2 x x x
U3 x x
U4 x
MISSING DATA PATTERN FREQUENCIES FOR U
Pattern Frequency Pattern Frequency Pattern Frequency
1 262 3 87
2 65 4 86
COVARIANCE COVERAGE OF DATA FOR THE FIRST REPLICATION
Minimum covariance coverage value 0.100
PROPORTION OF DATA PRESENT FOR U
Covariance Coverage
U1 U2 U3 U4
________ ________ ________ ________
U1 1.000
U2 0.870 0.870
U3 0.696 0.696 0.696
U4 0.524 0.524 0.524 0.524
SAMPLE STATISTICS FOR THE FIRST REPLICATION
SAMPLE STATISTICS
Means
X
________
0.032
Covariances
X
________
X 0.991
Correlations
X
________
X 1.000
MODEL FIT INFORMATION
Number of Free Parameters 5
Loglikelihood
H0 Value
Mean -726.559
Std Dev 0.000
Number of successful computations 1
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.000 -726.559 -726.559
0.980 0.000 -726.559 -726.559
0.950 0.000 -726.559 -726.559
0.900 0.000 -726.559 -726.559
0.800 0.000 -726.559 -726.559
0.700 0.000 -726.559 -726.559
0.500 0.000 -726.559 -726.559
0.300 0.000 -726.559 -726.559
0.200 0.000 -726.559 -726.559
0.100 0.000 -726.559 -726.559
0.050 0.000 -726.559 -726.559
0.020 0.000 -726.559 -726.559
0.010 0.000 -726.559 -726.559
Information Criteria
Akaike (AIC)
Mean 1463.117
Std Dev 0.000
Number of successful computations 1
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.000 1463.117 1463.117
0.980 0.000 1463.117 1463.117
0.950 0.000 1463.117 1463.117
0.900 0.000 1463.117 1463.117
0.800 0.000 1463.117 1463.117
0.700 0.000 1463.117 1463.117
0.500 0.000 1463.117 1463.117
0.300 0.000 1463.117 1463.117
0.200 0.000 1463.117 1463.117
0.100 0.000 1463.117 1463.117
0.050 0.000 1463.117 1463.117
0.020 0.000 1463.117 1463.117
0.010 0.000 1463.117 1463.117
Bayesian (BIC)
Mean 1484.190
Std Dev 0.000
Number of successful computations 1
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.000 1484.190 1484.190
0.980 0.000 1484.190 1484.190
0.950 0.000 1484.190 1484.190
0.900 0.000 1484.190 1484.190
0.800 0.000 1484.190 1484.190
0.700 0.000 1484.190 1484.190
0.500 0.000 1484.190 1484.190
0.300 0.000 1484.190 1484.190
0.200 0.000 1484.190 1484.190
0.100 0.000 1484.190 1484.190
0.050 0.000 1484.190 1484.190
0.020 0.000 1484.190 1484.190
0.010 0.000 1484.190 1484.190
Sample-Size Adjusted BIC (n* = (n + 2) / 24)
Mean 1468.320
Std Dev 0.000
Number of successful computations 1
Proportions Percentiles
Expected Observed Expected Observed
0.990 0.000 1468.320 1468.320
0.980 0.000 1468.320 1468.320
0.950 0.000 1468.320 1468.320
0.900 0.000 1468.320 1468.320
0.800 0.000 1468.320 1468.320
0.700 0.000 1468.320 1468.320
0.500 0.000 1468.320 1468.320
0.300 0.000 1468.320 1468.320
0.200 0.000 1468.320 1468.320
0.100 0.000 1468.320 1468.320
0.050 0.000 1468.320 1468.320
0.020 0.000 1468.320 1468.320
0.010 0.000 1468.320 1468.320
MODEL RESULTS
ESTIMATES S. E. M. S. E. 95% % Sig
Population Average Std. Dev. Average Cover Coeff
F BY
U1 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
U2 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
U3 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
U4 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
F ON
X 0.500 0.5971 0.0000 0.0702 0.0094 1.000 1.000
Thresholds
U1$1 2.000 2.0466 0.0000 0.1383 0.0022 1.000 1.000
U2$1 1.500 1.4553 0.0000 0.1255 0.0020 1.000 1.000
U3$1 1.000 1.1059 0.0000 0.1292 0.0112 1.000 1.000
U4$1 1.000 0.7618 0.0000 0.1371 0.0567 1.000 1.000
Residual Variances
F 0.000 0.0000 0.0000 0.0000 0.0000 1.000 0.000
QUALITY OF NUMERICAL RESULTS
Average Condition Number for the Information Matrix 0.218E+00
(ratio of smallest to largest eigenvalue)
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
TAU
U1$1 U2$1 U3$1 U4$1
________ ________ ________ ________
2 3 4 5
NU
U1 U2 U3 U4 X
________ ________ ________ ________ ________
0 0 0 0 0
LAMBDA
F X
________ ________
U1 0 0
U2 0 0
U3 0 0
U4 0 0
X 0 0
THETA
U1 U2 U3 U4 X
________ ________ ________ ________ ________
U1 0
U2 0 0
U3 0 0 0
U4 0 0 0 0
X 0 0 0 0 0
ALPHA
F X
________ ________
0 0
BETA
F X
________ ________
F 0 1
X 0 0
PSI
F X
________ ________
F 0
X 0 0
STARTING VALUES
TAU
U1$1 U2$1 U3$1 U4$1
________ ________ ________ ________
2.000 1.500 1.000 1.000
NU
U1 U2 U3 U4 X
________ ________ ________ ________ ________
0.000 0.000 0.000 0.000 0.000
LAMBDA
F X
________ ________
U1 1.000 0.000
U2 1.000 0.000
U3 1.000 0.000
U4 1.000 0.000
X 0.000 1.000
THETA
U1 U2 U3 U4 X
________ ________ ________ ________ ________
U1 1.000
U2 0.000 1.000
U3 0.000 0.000 1.000
U4 0.000 0.000 0.000 1.000
X 0.000 0.000 0.000 0.000 0.000
ALPHA
F X
________ ________
0.000 0.000
BETA
F X
________ ________
F 0.000 0.500
X 0.000 0.000
PSI
F X
________ ________
F 0.000
X 0.000 0.500
POPULATION VALUES
TAU
U1$1 U2$1 U3$1 U4$1
________ ________ ________ ________
2.000 1.500 1.000 1.000
NU
U1 U2 U3 U4 X
________ ________ ________ ________ ________
0.000 0.000 0.000 0.000 0.000
LAMBDA
F X
________ ________
U1 1.000 0.000
U2 1.000 0.000
U3 1.000 0.000
U4 1.000 0.000
X 0.000 1.000
THETA
U1 U2 U3 U4 X
________ ________ ________ ________ ________
U1 0.000
U2 0.000 0.000
U3 0.000 0.000 0.000
U4 0.000 0.000 0.000 0.000
X 0.000 0.000 0.000 0.000 0.000
ALPHA
F X
________ ________
0.000 0.000
BETA
F X
________ ________
F 0.000 0.500
X 0.000 0.000
PSI
F X
________ ________
F 0.000
X 0.000 1.000
TECHNICAL 8 OUTPUT
TECHNICAL 8 OUTPUT FOR REPLICATION 1
E STEP ITER LOGLIKELIHOOD ABS CHANGE REL CHANGE ALGORITHM
1 -0.72931501D+03 0.0000000 0.0000000 EM
2 -0.72656249D+03 2.7525241 0.0037741 EM
3 -0.72655862D+03 0.0038652 0.0000053 EM
4 -0.72655862D+03 0.0000000 0.0000000 EM
TECHNICAL 9 OUTPUT
Error messages for each replication (if any)
SAVEDATA INFORMATION
Order of variables
U1
U2
U3
U4
X
Save file
ex6.19.dat
Save file format Free
Save file record length 10000
Missing designated by 999
Beginning Time: 22:24:18
Ending Time: 22:24:18
Elapsed Time: 00:00:00
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-2022 Muthen & Muthen
Back to examples