Mplus VERSION 7.4
MUTHEN & MUTHEN
06/02/2016 12:17 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;
model:
positive on male (pmale)
black hisp es fh123 hsdrp;
u on male (male)
black hisp es fh123 hsdrp;
[u$1] (thresh);
[positive] (int);
positive (v);
model constraint:
! expectation for covs=1,0,0,0,0,0
new(p0 fi expcond predict);
p0 = Phi(thresh - male);
fi = exp(v/2);
expcond = exp(int+pmale);
predict = (1-p0)*fi*expcond;
output:
tech1 tech8 standardized;
*** 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 0
Observed dependent variables
Continuous
POSITIVE
Binary and ordered categorical (ordinal)
U
Observed independent variables
MALE BLACK HISP ES FH123 HSDRP
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 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 15
Dimensions of numerical integration 0
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 15
Loglikelihood
H0 Value -1088.400
H0 Scaling Correction Factor 0.9678
for MLR
Information Criteria
Akaike (AIC) 2206.801
Bayesian (BIC) 2282.552
Sample-Size Adjusted BIC 2234.908
(n* = (n + 2) / 24)
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
POSITIVE ON
MALE 0.271 0.057 4.747 0.000
BLACK -0.285 0.068 -4.190 0.000
HISP -0.102 0.080 -1.271 0.204
ES 0.141 0.073 1.931 0.053
FH123 -0.013 0.082 -0.156 0.876
HSDRP 0.128 0.065 1.961 0.050
U ON
MALE 0.704 0.079 8.938 0.000
BLACK -0.713 0.095 -7.528 0.000
HISP -0.378 0.109 -3.456 0.001
ES 0.190 0.108 1.763 0.078
FH123 0.264 0.126 2.102 0.036
HSDRP 0.034 0.099 0.349 0.727
Intercepts
POSITIVE 0.601 0.050 12.105 0.000
Thresholds
U$1 0.458 0.071 6.460 0.000
Residual Variances
POSITIVE 0.315 0.015 20.791 0.000
New/Additional Parameters
P0 0.403 0.027 14.969 0.000
FI 1.170 0.009 132.153 0.000
EXPCOND 2.392 0.105 22.691 0.000
PREDICT 1.671 0.105 15.862 0.000
STANDARDIZED MODEL RESULTS
STDYX Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
POSITIVE ON
MALE 0.227 0.047 4.868 0.000
BLACK -0.212 0.049 -4.326 0.000
HISP -0.063 0.049 -1.274 0.203
ES 0.083 0.043 1.933 0.053
FH123 -0.007 0.042 -0.156 0.876
HSDRP 0.085 0.043 1.966 0.049
U ON
MALE 0.316 0.032 9.743 0.000
BLACK -0.284 0.035 -8.030 0.000
HISP -0.125 0.036 -3.493 0.000
ES 0.060 0.034 1.768 0.077
FH123 0.073 0.034 2.112 0.035
HSDRP 0.012 0.035 0.349 0.727
Intercepts
POSITIVE 1.007 0.089 11.253 0.000
Thresholds
U$1 0.411 0.062 6.581 0.000
Residual Variances
POSITIVE 0.883 0.030 29.485 0.000
STDY Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
POSITIVE ON
MALE 0.454 0.093 4.890 0.000
BLACK -0.477 0.110 -4.342 0.000
HISP -0.171 0.134 -1.274 0.203
ES 0.236 0.122 1.934 0.053
FH123 -0.021 0.137 -0.156 0.876
HSDRP 0.214 0.109 1.967 0.049
U ON
MALE 0.632 0.064 9.907 0.000
BLACK -0.640 0.079 -8.128 0.000
HISP -0.339 0.097 -3.502 0.000
ES 0.171 0.097 1.769 0.077
FH123 0.237 0.112 2.114 0.035
HSDRP 0.031 0.089 0.349 0.727
Intercepts
POSITIVE 1.007 0.089 11.253 0.000
Thresholds
U$1 0.411 0.062 6.581 0.000
Residual Variances
POSITIVE 0.883 0.030 29.485 0.000
STD Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
POSITIVE ON
MALE 0.271 0.057 4.747 0.000
BLACK -0.285 0.068 -4.190 0.000
HISP -0.102 0.080 -1.271 0.204
ES 0.141 0.073 1.931 0.053
FH123 -0.013 0.082 -0.156 0.876
HSDRP 0.128 0.065 1.961 0.050
U ON
MALE 0.704 0.079 8.938 0.000
BLACK -0.713 0.095 -7.528 0.000
HISP -0.378 0.109 -3.456 0.001
ES 0.190 0.108 1.763 0.078
FH123 0.264 0.126 2.102 0.036
HSDRP 0.034 0.099 0.349 0.727
Intercepts
POSITIVE 0.601 0.050 12.105 0.000
Thresholds
U$1 0.458 0.071 6.460 0.000
Residual Variances
POSITIVE 0.315 0.015 20.791 0.000
R-SQUARE
Observed Two-Tailed
Variable Estimate S.E. Est./S.E. P-Value
U 0.194 0.027 7.242 0.000
POSITIVE 0.117 0.030 3.906 0.000
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.341E-03
(ratio of smallest to largest eigenvalue)
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
TAU
U$1
________
1 15
NU
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
1 0 0 0 0 0
NU
ES FH123 HSDRP
________ ________ ________
1 0 0 0
LAMBDA
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
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
ES FH123 HSDRP
________ ________ ________
U 0 0 0
POSITIVE 0 0 0
MALE 0 0 0
BLACK 0 0 0
HISP 0 0 0
ES 0 0 0
FH123 0 0 0
HSDRP 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
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
1 0 1 0 0 0
ALPHA
ES FH123 HSDRP
________ ________ ________
1 0 0 0
BETA
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
U 0 0 2 3 4
POSITIVE 0 0 8 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
ES FH123 HSDRP
________ ________ ________
U 5 6 7
POSITIVE 11 12 13
MALE 0 0 0
BLACK 0 0 0
HISP 0 0 0
ES 0 0 0
FH123 0 0 0
HSDRP 0 0 0
PSI
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
U 0
POSITIVE 0 14
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
PSI
ES FH123 HSDRP
________ ________ ________
ES 0
FH123 0 0
HSDRP 0 0 0
PARAMETER SPECIFICATION FOR THE ADDITIONAL PARAMETERS
NEW/ADDITIONAL PARAMETERS
P0 FI EXPCOND PREDICT
________ ________ ________ ________
1 16 17 18 19
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
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
U 1.000 0.000 0.000 0.000 0.000
POSITIVE 0.000 1.000 0.000 0.000 0.000
MALE 0.000 0.000 1.000 0.000 0.000
BLACK 0.000 0.000 0.000 1.000 0.000
HISP 0.000 0.000 0.000 0.000 1.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
ES FH123 HSDRP
________ ________ ________
U 0.000 0.000 0.000
POSITIVE 0.000 0.000 0.000
MALE 0.000 0.000 0.000
BLACK 0.000 0.000 0.000
HISP 0.000 0.000 0.000
ES 1.000 0.000 0.000
FH123 0.000 1.000 0.000
HSDRP 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
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
1 0.000 0.773 0.000 0.000 0.000
ALPHA
ES FH123 HSDRP
________ ________ ________
1 0.000 0.000 0.000
BETA
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
U 0.000 0.000 0.000 0.000 0.000
POSITIVE 0.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
ES FH123 HSDRP
________ ________ ________
U 0.000 0.000 0.000
POSITIVE 0.000 0.000 0.000
MALE 0.000 0.000 0.000
BLACK 0.000 0.000 0.000
HISP 0.000 0.000 0.000
ES 0.000 0.000 0.000
FH123 0.000 0.000 0.000
HSDRP 0.000 0.000 0.000
PSI
U POSITIVE MALE BLACK HISP
________ ________ ________ ________ ________
U 1.000
POSITIVE 0.000 0.174
MALE 0.000 0.000 0.125
BLACK 0.000 0.000 0.000 0.099
HISP 0.000 0.000 0.000 0.000 0.067
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
ES FH123 HSDRP
________ ________ ________
ES 0.062
FH123 0.000 0.047
HSDRP 0.000 0.000 0.079
STARTING VALUES FOR THE ADDITIONAL PARAMETERS
NEW/ADDITIONAL PARAMETERS
P0 FI EXPCOND PREDICT
________ ________ ________ ________
1 0.500 0.500 0.500 0.500
TECHNICAL 8 OUTPUT
E STEP ITER LOGLIKELIHOOD ABS CHANGE REL CHANGE ALGORITHM
1 -0.12572710D+04 0.0000000 0.0000000 EM
2 -0.11241862D+04 133.0847888 0.1058521 FS
3 -0.10936506D+04 30.5355972 0.0271624 FS
4 -0.10886079D+04 5.0427285 0.0046109 FS
5 -0.10884166D+04 0.1912835 0.0001757 FS
6 -0.10884131D+04 0.0034302 0.0000032 FS
7 -0.10884003D+04 0.0128736 0.0000118 EM
8 -0.10884003D+04 0.0000000 0.0000000 EM
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\bengt 2013\documents\bengt\mplus runs\a book - topic 2 mplus runs\nlsy hd tobit cens-in
Beginning Time: 12:17:40
Ending Time: 12:17:40
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-2015 Muthen & Muthen
Back to examples