Mplus VERSION 8.8
MUTHEN & MUTHEN
04/19/2022 10:22 PM
INPUT INSTRUCTIONS
title:
this is an example of multiple imputation of
plausible values generated from a multiple
indicator linear growth model for categorical
outcomes using Bayesian estimation
The Monte Carlo setup is the same as mcex6.15,
except Bayes estimation is used instead of WLSMV
montecarlo:
names = u11 u21 u31 u12 u22 u32 u13 u23 u33;
generate = u11-u33(1);
categorical = u11-u33;
nobs = 500;
nreps = 1;
save = ex11.7.dat;
model montecarlo:
f1 by u11@1
u21*1 (1)
u31*1 (2);
f2 by u12@1
u22*1 (1)
u32*1 (2);
f3 by u13@1
u23*1 (1)
u33*1 (2);
[u11$1*-.5 u12$1*-.5 u13$1*-.5] (3);
[u21$1*0 u22$1*0 u23$1*0] (4);
[u31$1*.5 u32$1*.5 u33$1*.5] (5);
i s | f1@0 f2@1 f3@2;
[i@0 s*1];
i*.4 s*0;
f1-f3*.4;
analysis:
estimator = bayes;
process = 2;
model:
f1 by u11@1
u21*1 (1)
u31*1 (2);
f2 by u12@1
u22*1 (1)
u32*1 (2);
f3 by u13@1
u23*1 (1)
u33*1 (2);
[u11$1*-.5 u12$1*-.5 u13$1*-.5] (3);
[u21$1*0 u22$1*0 u23$1*0] (4);
[u31$1*.5 u32$1*.5 u33$1*.5] (5);
i s | f1@0 f2@1 f3@2;
[i@0 s*1];
i*.4 s*0;
f1-f3*.4;
output:
tech9;
INPUT READING TERMINATED NORMALLY
this is an example of multiple imputation of
plausible values generated from a multiple
indicator linear growth model for categorical
outcomes using Bayesian estimation
The Monte Carlo setup is the same as mcex6.15,
except Bayes estimation is used instead of WLSMV
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 9
Number of independent variables 0
Number of continuous latent variables 5
Observed dependent variables
Binary and ordered categorical (ordinal)
U11 U21 U31 U12 U22 U32
U13 U23 U33
Continuous latent variables
F1 F2 F3 I S
Estimator BAYES
Specifications for Bayesian Estimation
Point estimate MEDIAN
Number of Markov chain Monte Carlo (MCMC) chains 2
Random seed for the first chain 0
Starting value information UNPERTURBED
Algorithm used for Markov chain Monte Carlo GIBBS(PX1)
Convergence criterion 0.500D-01
Maximum number of iterations 50000
K-th iteration used for thinning 1
Link PROBIT
MODEL FIT INFORMATION
Number of Free Parameters 12
Bayesian Posterior Predictive Checking using Chi-Square
Posterior Predictive P-Value
Mean 0.517
Std Dev 0.000
Number of successful computations 1
Cumulative Distribution Function
Value Function Value
0.990 1.000
0.980 1.000
0.950 1.000
0.900 1.000
0.800 1.000
0.700 1.000
0.500 0.000
0.300 0.000
0.200 0.000
0.100 0.000
0.050 0.000
0.020 0.000
0.010 0.000
MODEL RESULTS
ESTIMATES S. E. M. S. E. 95% % Sig
Population Average Std. Dev. Average Cover Coeff
I |
F1 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
F2 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
F3 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
S |
F1 0.000 0.0000 0.0000 0.0000 0.0000 1.000 0.000
F2 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
F3 2.000 2.0000 0.0000 0.0000 0.0000 1.000 0.000
F1 BY
U11 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
U21 1.000 1.0857 0.0000 0.1190 0.0073 1.000 1.000
U31 1.000 1.1478 0.0000 0.1346 0.0218 1.000 1.000
F2 BY
U12 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
U22 1.000 1.0857 0.0000 0.1190 0.0073 1.000 1.000
U32 1.000 1.1478 0.0000 0.1346 0.0218 1.000 1.000
F3 BY
U13 1.000 1.0000 0.0000 0.0000 0.0000 1.000 0.000
U23 1.000 1.0857 0.0000 0.1190 0.0073 1.000 1.000
U33 1.000 1.1478 0.0000 0.1346 0.0218 1.000 1.000
S WITH
I 0.000 -0.1163 0.0000 0.0725 0.0135 0.000 1.000
Means
I 0.000 0.0000 0.0000 0.0000 0.0000 1.000 0.000
S 1.000 0.9328 0.0000 0.0614 0.0045 1.000 1.000
Intercepts
F1 0.000 0.0000 0.0000 0.0000 0.0000 1.000 0.000
F2 0.000 0.0000 0.0000 0.0000 0.0000 1.000 0.000
F3 0.000 0.0000 0.0000 0.0000 0.0000 1.000 0.000
Thresholds
U11$1 -0.500 -0.6581 0.0000 0.0670 0.0250 0.000 1.000
U21$1 0.000 -0.0149 0.0000 0.0704 0.0002 1.000 0.000
U31$1 0.500 0.4892 0.0000 0.0912 0.0001 1.000 1.000
U12$1 -0.500 -0.6581 0.0000 0.0670 0.0250 0.000 1.000
U22$1 0.000 -0.0149 0.0000 0.0704 0.0002 1.000 0.000
U32$1 0.500 0.4892 0.0000 0.0912 0.0001 1.000 1.000
U13$1 -0.500 -0.6581 0.0000 0.0670 0.0250 0.000 1.000
U23$1 0.000 -0.0149 0.0000 0.0704 0.0002 1.000 0.000
U33$1 0.500 0.4892 0.0000 0.0912 0.0001 1.000 1.000
Variances
I 0.400 0.5268 0.0000 0.1421 0.0161 1.000 1.000
S 0.000 0.0734 0.0000 0.0483 0.0054 0.000 1.000
Residual Variances
F1 0.400 0.4899 0.0000 0.1820 0.0081 1.000 1.000
F2 0.400 0.3995 0.0000 0.1201 0.0000 1.000 1.000
F3 0.400 0.4550 0.0000 0.1968 0.0030 1.000 1.000
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
TAU
U11$1 U21$1 U31$1 U12$1 U22$1
________ ________ ________ ________ ________
10 11 12 10 11
TAU
U32$1 U13$1 U23$1 U33$1
________ ________ ________ ________
12 10 11 12
NU
U11 U21 U31 U12 U22
________ ________ ________ ________ ________
0 0 0 0 0
NU
U32 U13 U23 U33
________ ________ ________ ________
0 0 0 0
LAMBDA
F1 F2 F3 I S
________ ________ ________ ________ ________
U11 0 0 0 0 0
U21 1 0 0 0 0
U31 2 0 0 0 0
U12 0 0 0 0 0
U22 0 1 0 0 0
U32 0 2 0 0 0
U13 0 0 0 0 0
U23 0 0 1 0 0
U33 0 0 2 0 0
THETA
U11 U21 U31 U12 U22
________ ________ ________ ________ ________
U11 0
U21 0 0
U31 0 0 0
U12 0 0 0 0
U22 0 0 0 0 0
U32 0 0 0 0 0
U13 0 0 0 0 0
U23 0 0 0 0 0
U33 0 0 0 0 0
THETA
U32 U13 U23 U33
________ ________ ________ ________
U32 0
U13 0 0
U23 0 0 0
U33 0 0 0 0
ALPHA
F1 F2 F3 I S
________ ________ ________ ________ ________
0 0 0 0 3
BETA
F1 F2 F3 I S
________ ________ ________ ________ ________
F1 0 0 0 0 0
F2 0 0 0 0 0
F3 0 0 0 0 0
I 0 0 0 0 0
S 0 0 0 0 0
PSI
F1 F2 F3 I S
________ ________ ________ ________ ________
F1 4
F2 0 5
F3 0 0 6
I 0 0 0 7
S 0 0 0 8 9
STARTING VALUES
TAU
U11$1 U21$1 U31$1 U12$1 U22$1
________ ________ ________ ________ ________
-0.500 0.000 0.500 -0.500 0.000
TAU
U32$1 U13$1 U23$1 U33$1
________ ________ ________ ________
0.500 -0.500 0.000 0.500
NU
U11 U21 U31 U12 U22
________ ________ ________ ________ ________
0.000 0.000 0.000 0.000 0.000
NU
U32 U13 U23 U33
________ ________ ________ ________
0.000 0.000 0.000 0.000
LAMBDA
F1 F2 F3 I S
________ ________ ________ ________ ________
U11 1.000 0.000 0.000 0.000 0.000
U21 1.000 0.000 0.000 0.000 0.000
U31 1.000 0.000 0.000 0.000 0.000
U12 0.000 1.000 0.000 0.000 0.000
U22 0.000 1.000 0.000 0.000 0.000
U32 0.000 1.000 0.000 0.000 0.000
U13 0.000 0.000 1.000 0.000 0.000
U23 0.000 0.000 1.000 0.000 0.000
U33 0.000 0.000 1.000 0.000 0.000
THETA
U11 U21 U31 U12 U22
________ ________ ________ ________ ________
U11 1.000
U21 0.000 1.000
U31 0.000 0.000 1.000
U12 0.000 0.000 0.000 1.000
U22 0.000 0.000 0.000 0.000 1.000
U32 0.000 0.000 0.000 0.000 0.000
U13 0.000 0.000 0.000 0.000 0.000
U23 0.000 0.000 0.000 0.000 0.000
U33 0.000 0.000 0.000 0.000 0.000
THETA
U32 U13 U23 U33
________ ________ ________ ________
U32 1.000
U13 0.000 1.000
U23 0.000 0.000 1.000
U33 0.000 0.000 0.000 1.000
ALPHA
F1 F2 F3 I S
________ ________ ________ ________ ________
0.000 0.000 0.000 0.000 1.000
BETA
F1 F2 F3 I S
________ ________ ________ ________ ________
F1 0.000 0.000 0.000 1.000 0.000
F2 0.000 0.000 0.000 1.000 1.000
F3 0.000 0.000 0.000 1.000 2.000
I 0.000 0.000 0.000 0.000 0.000
S 0.000 0.000 0.000 0.000 0.000
PSI
F1 F2 F3 I S
________ ________ ________ ________ ________
F1 0.400
F2 0.000 0.400
F3 0.000 0.000 0.400
I 0.000 0.000 0.000 0.400
S 0.000 0.000 0.000 0.000 0.000
POPULATION VALUES
TAU
U11$1 U21$1 U31$1 U12$1 U22$1
________ ________ ________ ________ ________
-0.500 0.000 0.500 -0.500 0.000
TAU
U32$1 U13$1 U23$1 U33$1
________ ________ ________ ________
0.500 -0.500 0.000 0.500
NU
U11 U21 U31 U12 U22
________ ________ ________ ________ ________
0.000 0.000 0.000 0.000 0.000
NU
U32 U13 U23 U33
________ ________ ________ ________
0.000 0.000 0.000 0.000
LAMBDA
F1 F2 F3 I S
________ ________ ________ ________ ________
U11 1.000 0.000 0.000 0.000 0.000
U21 1.000 0.000 0.000 0.000 0.000
U31 1.000 0.000 0.000 0.000 0.000
U12 0.000 1.000 0.000 0.000 0.000
U22 0.000 1.000 0.000 0.000 0.000
U32 0.000 1.000 0.000 0.000 0.000
U13 0.000 0.000 1.000 0.000 0.000
U23 0.000 0.000 1.000 0.000 0.000
U33 0.000 0.000 1.000 0.000 0.000
THETA
U11 U21 U31 U12 U22
________ ________ ________ ________ ________
U11 0.000
U21 0.000 0.000
U31 0.000 0.000 0.000
U12 0.000 0.000 0.000 0.000
U22 0.000 0.000 0.000 0.000 0.000
U32 0.000 0.000 0.000 0.000 0.000
U13 0.000 0.000 0.000 0.000 0.000
U23 0.000 0.000 0.000 0.000 0.000
U33 0.000 0.000 0.000 0.000 0.000
THETA
U32 U13 U23 U33
________ ________ ________ ________
U32 0.000
U13 0.000 0.000
U23 0.000 0.000 0.000
U33 0.000 0.000 0.000 0.000
ALPHA
F1 F2 F3 I S
________ ________ ________ ________ ________
0.000 0.000 0.000 0.000 1.000
BETA
F1 F2 F3 I S
________ ________ ________ ________ ________
F1 0.000 0.000 0.000 1.000 0.000
F2 0.000 0.000 0.000 1.000 1.000
F3 0.000 0.000 0.000 1.000 2.000
I 0.000 0.000 0.000 0.000 0.000
S 0.000 0.000 0.000 0.000 0.000
PSI
F1 F2 F3 I S
________ ________ ________ ________ ________
F1 0.400
F2 0.000 0.400
F3 0.000 0.000 0.400
I 0.000 0.000 0.000 0.400
S 0.000 0.000 0.000 0.000 0.000
PRIORS FOR ALL PARAMETERS PRIOR MEAN PRIOR VARIANCE PRIOR STD. DEV.
Parameter 1~N(0.000,5.000) 0.0000 5.0000 2.2361
Parameter 2~N(0.000,5.000) 0.0000 5.0000 2.2361
Parameter 3~N(0.000,infinity) 0.0000 infinity infinity
Parameter 4~IG(-1.000,0.000) infinity infinity infinity
Parameter 5~IG(-1.000,0.000) infinity infinity infinity
Parameter 6~IG(-1.000,0.000) infinity infinity infinity
Parameter 7~IW(0.000,-3) infinity infinity infinity
Parameter 8~IW(0.000,-3) infinity infinity infinity
Parameter 9~IW(0.000,-3) infinity infinity infinity
Parameter 10~N(0.000,5.000) 0.0000 5.0000 2.2361
Parameter 11~N(0.000,5.000) 0.0000 5.0000 2.2361
Parameter 12~N(0.000,5.000) 0.0000 5.0000 2.2361
TECHNICAL 9 OUTPUT
Error messages for each replication (if any)
SAVEDATA INFORMATION
Order of variables
U11
U21
U31
U12
U22
U32
U13
U23
U33
Save file
ex11.7.dat
Save file format Free
Save file record length 10000
Beginning Time: 22:22:07
Ending Time: 22:22:10
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-2022 Muthen & Muthen
Back to examples