Mplus VERSION 7.4
MUTHEN & MUTHEN
06/10/2016 1:16 PM
INPUT INSTRUCTIONS
title:
NLSY cohort 64
hd82-hd89 (ages 18 - 25)
data:
file is big.dat;
format is
2f5,f2,t14,5f7,t50,f8,t60,6f1.0,t67,2f2.0,t71,8f1.0,t79,f2.0,t82,4f2.0;
variable:
names are
id houseid cohort
weight82 weight83 weight84 weight88 weight89 weight94
hd82 hd83 hd84 hd88 hd89 hd94
dep89 dep94 male black hisp es fh1 fh23 fh123 hsdrp coll
ed89 ed94 cd89 cd94;
useobservations = cohort EQ 64 AND (coll GT 0 AND coll LT 20);
usev are male black hisp es fh123 hsdrp u positive;
missing are .;
categorical = u;
data twopart:
names = hd84;
binary = u;
continuous = positive;
analysis:
estimator = mlr;
link = probit;
mconvergence = 0.00001;
integration = 30;
model:
positive u on male black hisp es fh123 hsdrp;
f by u
positive (lam); f@1;
positive (v);
model constraint:
new(rescorr);
rescorr = lam/(SQRT(lam*lam+v)*SQRT(1+1));
output:
tech1 tech8 standardized;
plot:
type = plot3;
*** WARNING
Data set contains cases with missing on x-variables.
These cases were not included in the analysis.
Number of cases with missing on x-variables: 20
*** WARNING
Data set contains cases with missing on all variables except
x-variables. These cases were not included in the analysis.
Number of cases with missing on all variables except x-variables: 19
2 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS
NLSY cohort 64
hd82-hd89 (ages 18 - 25)
SUMMARY OF ANALYSIS
Number of groups 1
Number of observations 1153
Number of dependent variables 2
Number of independent variables 6
Number of continuous latent variables 1
Observed dependent variables
Continuous
POSITIVE
Binary and ordered categorical (ordinal)
U
Observed independent variables
MALE BLACK HISP ES FH123 HSDRP
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-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 10
Minimum value for logit thresholds -10
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 30
Dimensions of numerical integration 1
Adaptive quadrature ON
Link PROBIT
Cholesky ON
Input data file(s)
big.dat
Input data format
(2F5,F2,T14,5F7,T50,F8,T60,6F1.0,T67,2F2.0,T71,8F1.0,T79,F2.0,T82,4F2.0)
SUMMARY OF DATA
Number of missing data patterns 2
Number of y missing data patterns 2
Number of u missing data patterns 1
COVARIANCE COVERAGE OF DATA
Minimum covariance coverage value 0.100
PROPORTION OF DATA PRESENT FOR Y
Covariance Coverage
POSITIVE MALE BLACK HISP ES
________ ________ ________ ________ ________
POSITIVE 0.400
MALE 0.400 1.000
BLACK 0.400 1.000 1.000
HISP 0.400 1.000 1.000 1.000
ES 0.400 1.000 1.000 1.000 1.000
FH123 0.400 1.000 1.000 1.000 1.000
HSDRP 0.400 1.000 1.000 1.000 1.000
Covariance Coverage
FH123 HSDRP
________ ________
FH123 1.000
HSDRP 1.000 1.000
UNIVARIATE PROPORTIONS AND COUNTS FOR CATEGORICAL VARIABLES
U
Category 1 0.600 692.000
Category 2 0.400 461.000
THE MODEL ESTIMATION TERMINATED NORMALLY
MODEL FIT INFORMATION
Number of Free Parameters 16
Loglikelihood
H0 Value -1088.182
H0 Scaling Correction Factor 0.9583
for MLR
Information Criteria
Akaike (AIC) 2208.364
Bayesian (BIC) 2289.166
Sample-Size Adjusted BIC 2238.344
(n* = (n + 2) / 24)
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
F BY
U 1.000 0.000 999.000 999.000
POSITIVE 0.565 0.328 1.723 0.085
POSITIVE ON
MALE 0.459 0.127 3.606 0.000
BLACK -0.477 0.132 -3.612 0.000
HISP -0.201 0.107 -1.883 0.060
ES 0.187 0.079 2.377 0.017
FH123 0.055 0.097 0.570 0.568
HSDRP 0.138 0.072 1.931 0.053
U ON
MALE 0.992 0.112 8.845 0.000
BLACK -1.009 0.133 -7.567 0.000
HISP -0.528 0.156 -3.395 0.001
ES 0.272 0.154 1.773 0.076
FH123 0.370 0.177 2.096 0.036
HSDRP 0.040 0.141 0.281 0.778
Intercepts
POSITIVE 0.154 0.284 0.542 0.588
Thresholds
U$1 0.644 0.101 6.394 0.000
Variances
F 1.000 0.000 999.000 999.000
Residual Variances
POSITIVE 0.098 0.250 0.393 0.694
New/Additional Parameters
RESCORR 0.618 0.269 2.299 0.022
STANDARDIZED MODEL RESULTS
STDYX Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
F BY
U 0.635 0.010 60.740 0.000
POSITIVE 0.779 0.326 2.387 0.017
POSITIVE ON
MALE 0.316 0.047 6.729 0.000
BLACK -0.292 0.046 -6.302 0.000
HISP -0.102 0.046 -2.215 0.027
ES 0.091 0.039 2.345 0.019
FH123 0.023 0.040 0.592 0.554
HSDRP 0.076 0.038 1.975 0.048
U ON
MALE 0.315 0.033 9.625 0.000
BLACK -0.285 0.035 -8.067 0.000
HISP -0.123 0.036 -3.429 0.001
ES 0.061 0.034 1.777 0.076
FH123 0.072 0.034 2.106 0.035
HSDRP 0.010 0.036 0.281 0.778
Intercepts
POSITIVE 0.212 0.425 0.499 0.618
Thresholds
U$1 0.409 0.063 6.511 0.000
Variances
F 1.000 0.000 999.000 999.000
Residual Variances
POSITIVE 0.186 0.534 0.349 0.727
STDY Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
F BY
U 0.635 0.010 60.740 0.000
POSITIVE 0.779 0.326 2.387 0.017
POSITIVE ON
MALE 0.632 0.093 6.781 0.000
BLACK -0.657 0.104 -6.348 0.000
HISP -0.277 0.125 -2.217 0.027
ES 0.257 0.110 2.348 0.019
FH123 0.076 0.128 0.592 0.554
HSDRP 0.191 0.096 1.977 0.048
U ON
MALE 0.630 0.064 9.783 0.000
BLACK -0.641 0.078 -8.166 0.000
HISP -0.335 0.098 -3.438 0.001
ES 0.173 0.097 1.778 0.075
FH123 0.235 0.112 2.108 0.035
HSDRP 0.025 0.090 0.281 0.778
Intercepts
POSITIVE 0.212 0.425 0.499 0.618
Thresholds
U$1 0.409 0.063 6.511 0.000
Variances
F 1.000 0.000 999.000 999.000
Residual Variances
POSITIVE 0.186 0.534 0.349 0.727
STD Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
F BY
U 1.000 0.000 999.000 999.000
POSITIVE 0.565 0.328 1.723 0.085
POSITIVE ON
MALE 0.459 0.127 3.606 0.000
BLACK -0.477 0.132 -3.612 0.000
HISP -0.201 0.107 -1.883 0.060
ES 0.187 0.079 2.377 0.017
FH123 0.055 0.097 0.570 0.568
HSDRP 0.138 0.072 1.931 0.053
U ON
MALE 0.992 0.112 8.845 0.000
BLACK -1.009 0.133 -7.567 0.000
HISP -0.528 0.156 -3.395 0.001
ES 0.272 0.154 1.773 0.076
FH123 0.370 0.177 2.096 0.036
HSDRP 0.040 0.141 0.281 0.778
Intercepts
POSITIVE 0.154 0.284 0.542 0.588
Thresholds
U$1 0.644 0.101 6.394 0.000
Variances
F 1.000 0.000 999.000 999.000
Residual Variances
POSITIVE 0.098 0.250 0.393 0.694
R-SQUARE
Observed Two-Tailed
Variable Estimate S.E. Est./S.E. P-Value
U 0.597 0.013 44.945 0.000
POSITIVE 0.814 0.534 1.525 0.127
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.186E-03
(ratio of smallest to largest eigenvalue)
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
TAU
U$1
________
1 16
NU
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
1 0 0 0 0 0
NU
ES FH123 HSDRP
________ ________ ________
1 0 0 0
LAMBDA
F U POSITIVE MALE BLACK
________ ________ ________ ________ ________
U 0 0 0 0 0
POSITIVE 0 0 0 0 0
MALE 0 0 0 0 0
BLACK 0 0 0 0 0
HISP 0 0 0 0 0
ES 0 0 0 0 0
FH123 0 0 0 0 0
HSDRP 0 0 0 0 0
LAMBDA
HISP ES FH123 HSDRP
________ ________ ________ ________
U 0 0 0 0
POSITIVE 0 0 0 0
MALE 0 0 0 0
BLACK 0 0 0 0
HISP 0 0 0 0
ES 0 0 0 0
FH123 0 0 0 0
HSDRP 0 0 0 0
THETA
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
U 0
POSITIVE 0 0
MALE 0 0 0
BLACK 0 0 0 0
HISP 0 0 0 0 0
ES 0 0 0 0 0
FH123 0 0 0 0 0
HSDRP 0 0 0 0 0
THETA
ES FH123 HSDRP
________ ________ ________
ES 0
FH123 0 0
HSDRP 0 0 0
ALPHA
F U POSITIVE MALE BLACK
________ ________ ________ ________ ________
1 0 0 1 0 0
ALPHA
HISP ES FH123 HSDRP
________ ________ ________ ________
1 0 0 0 0
BETA
F U POSITIVE MALE BLACK
________ ________ ________ ________ ________
F 0 0 0 0 0
U 0 0 0 2 3
POSITIVE 8 0 0 9 10
MALE 0 0 0 0 0
BLACK 0 0 0 0 0
HISP 0 0 0 0 0
ES 0 0 0 0 0
FH123 0 0 0 0 0
HSDRP 0 0 0 0 0
BETA
HISP ES FH123 HSDRP
________ ________ ________ ________
F 0 0 0 0
U 4 5 6 7
POSITIVE 11 12 13 14
MALE 0 0 0 0
BLACK 0 0 0 0
HISP 0 0 0 0
ES 0 0 0 0
FH123 0 0 0 0
HSDRP 0 0 0 0
PSI
F U POSITIVE MALE BLACK
________ ________ ________ ________ ________
F 0
U 0 0
POSITIVE 0 0 15
MALE 0 0 0 0
BLACK 0 0 0 0 0
HISP 0 0 0 0 0
ES 0 0 0 0 0
FH123 0 0 0 0 0
HSDRP 0 0 0 0 0
PSI
HISP ES FH123 HSDRP
________ ________ ________ ________
HISP 0
ES 0 0
FH123 0 0 0
HSDRP 0 0 0 0
PARAMETER SPECIFICATION FOR THE ADDITIONAL PARAMETERS
NEW/ADDITIONAL PARAMETERS
RESCORR
________
1 17
STARTING VALUES
TAU
U$1
________
1 0.226
NU
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
1 0.000 0.000 0.000 0.000 0.000
NU
ES FH123 HSDRP
________ ________ ________
1 0.000 0.000 0.000
LAMBDA
F U POSITIVE MALE BLACK
________ ________ ________ ________ ________
U 0.000 1.000 0.000 0.000 0.000
POSITIVE 0.000 0.000 1.000 0.000 0.000
MALE 0.000 0.000 0.000 1.000 0.000
BLACK 0.000 0.000 0.000 0.000 1.000
HISP 0.000 0.000 0.000 0.000 0.000
ES 0.000 0.000 0.000 0.000 0.000
FH123 0.000 0.000 0.000 0.000 0.000
HSDRP 0.000 0.000 0.000 0.000 0.000
LAMBDA
HISP ES FH123 HSDRP
________ ________ ________ ________
U 0.000 0.000 0.000 0.000
POSITIVE 0.000 0.000 0.000 0.000
MALE 0.000 0.000 0.000 0.000
BLACK 0.000 0.000 0.000 0.000
HISP 1.000 0.000 0.000 0.000
ES 0.000 1.000 0.000 0.000
FH123 0.000 0.000 1.000 0.000
HSDRP 0.000 0.000 0.000 1.000
THETA
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
U 0.000
POSITIVE 0.000 0.000
MALE 0.000 0.000 0.000
BLACK 0.000 0.000 0.000 0.000
HISP 0.000 0.000 0.000 0.000 0.000
ES 0.000 0.000 0.000 0.000 0.000
FH123 0.000 0.000 0.000 0.000 0.000
HSDRP 0.000 0.000 0.000 0.000 0.000
THETA
ES FH123 HSDRP
________ ________ ________
ES 0.000
FH123 0.000 0.000
HSDRP 0.000 0.000 0.000
ALPHA
F U POSITIVE MALE BLACK
________ ________ ________ ________ ________
1 0.000 0.000 0.773 0.000 0.000
ALPHA
HISP ES FH123 HSDRP
________ ________ ________ ________
1 0.000 0.000 0.000 0.000
BETA
F U POSITIVE MALE BLACK
________ ________ ________ ________ ________
F 0.000 0.000 0.000 0.000 0.000
U 1.000 0.000 0.000 0.000 0.000
POSITIVE 1.000 0.000 0.000 0.000 0.000
MALE 0.000 0.000 0.000 0.000 0.000
BLACK 0.000 0.000 0.000 0.000 0.000
HISP 0.000 0.000 0.000 0.000 0.000
ES 0.000 0.000 0.000 0.000 0.000
FH123 0.000 0.000 0.000 0.000 0.000
HSDRP 0.000 0.000 0.000 0.000 0.000
BETA
HISP ES FH123 HSDRP
________ ________ ________ ________
F 0.000 0.000 0.000 0.000
U 0.000 0.000 0.000 0.000
POSITIVE 0.000 0.000 0.000 0.000
MALE 0.000 0.000 0.000 0.000
BLACK 0.000 0.000 0.000 0.000
HISP 0.000 0.000 0.000 0.000
ES 0.000 0.000 0.000 0.000
FH123 0.000 0.000 0.000 0.000
HSDRP 0.000 0.000 0.000 0.000
PSI
F U POSITIVE MALE BLACK
________ ________ ________ ________ ________
F 1.000
U 0.000 1.000
POSITIVE 0.000 0.000 0.174
MALE 0.000 0.000 0.000 0.125
BLACK 0.000 0.000 0.000 0.000 0.099
HISP 0.000 0.000 0.000 0.000 0.000
ES 0.000 0.000 0.000 0.000 0.000
FH123 0.000 0.000 0.000 0.000 0.000
HSDRP 0.000 0.000 0.000 0.000 0.000
PSI
HISP ES FH123 HSDRP
________ ________ ________ ________
HISP 0.067
ES 0.000 0.062
FH123 0.000 0.000 0.047
HSDRP 0.000 0.000 0.000 0.079
STARTING VALUES FOR THE ADDITIONAL PARAMETERS
NEW/ADDITIONAL PARAMETERS
RESCORR
________
1 0.500
TECHNICAL 8 OUTPUT
E STEP ITER LOGLIKELIHOOD ABS CHANGE REL CHANGE ALGORITHM
1 -0.13635198D+04 0.0000000 0.0000000 EM
2 -0.11955034D+04 168.0164243 0.1232226 EM
3 -0.11423969D+04 53.1065109 0.0444219 EM
4 -0.11168947D+04 25.5022113 0.0223234 EM
5 -0.11035912D+04 13.3035097 0.0119112 EM
6 -0.10964917D+04 7.0994259 0.0064330 EM
7 -0.10926724D+04 3.8193884 0.0034833 EM
8 -0.10906111D+04 2.0612303 0.0018864 EM
9 -0.10894972D+04 1.1139693 0.0010214 EM
10 -0.10888946D+04 0.6025601 0.0005531 EM
11 -0.10885684D+04 0.3261967 0.0002996 EM
12 -0.10883917D+04 0.1767302 0.0001624 EM
13 -0.10882958D+04 0.0958305 0.0000880 EM
14 -0.10882438D+04 0.0520126 0.0000478 EM
15 -0.10882156D+04 0.0282575 0.0000260 EM
16 -0.10882002D+04 0.0153675 0.0000141 EM
17 -0.10881918D+04 0.0083659 0.0000077 EM
18 -0.10881873D+04 0.0045591 0.0000042 EM
19 -0.10881848D+04 0.0024872 0.0000023 EM
20 -0.10881834D+04 0.0013584 0.0000012 EM
21 -0.10881827D+04 0.0007428 0.0000007 EM
22 -0.10881823D+04 0.0004066 0.0000004 EM
23 -0.10881821D+04 0.0002229 0.0000002 EM
24 -0.10881819D+04 0.0001223 0.0000001 EM
25 -0.10881819D+04 0.0000672 0.0000001 EM
26 -0.10881818D+04 0.0000370 0.0000000 EM
27 -0.10881818D+04 0.0000204 0.0000000 EM
28 -0.10881818D+04 0.0000112 0.0000000 EM
29 -0.10881818D+04 0.0000062 0.0000000 EM
30 -0.10881818D+04 0.0000034 0.0000000 EM
31 -0.10881818D+04 0.0000019 0.0000000 EM
32 -0.10881818D+04 0.0000011 0.0000000 EM
33 -0.10881818D+04 0.0000006 0.0000000 EM
34 -0.10881818D+04 0.0000003 0.0000000 EM
SAMPLE STATISTICS FOR ESTIMATED FACTOR SCORES
SAMPLE STATISTICS
Means
F F_SE
________ ________
1 0.000 0.699
Covariances
F F_SE
________ ________
F 0.472
F_SE -0.092 0.039
Correlations
F F_SE
________ ________
F 1.000
F_SE -0.678 1.000
PLOT INFORMATION
The following plots are available:
Histograms (sample values, estimated factor scores, estimated values, residuals)
Scatterplots (sample values, estimated factor scores, estimated values, residuals)
Sample proportions, estimated and conditional estimated probabilities
Item characteristic curves
Information curves
Latent variable distribution plots
DIAGRAM INFORMATION
Use View Diagram under the Diagram menu in the Mplus Editor to view the diagram.
If running Mplus from the Mplus Diagrammer, the diagram opens automatically.
Diagram output
c:\users\gryphon\desktop\chapter7\ex7.5.dgm
Beginning Time: 13:16:10
Ending Time: 13:16:11
Elapsed Time: 00:00:01
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