Mplus VERSION 7.4
MUTHEN & MUTHEN
06/02/2016 12:11 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 hd84 male black hisp es fh123 hsdrp;
censored = hd84 (bi);
missing are .;
analysis:
estimator = mlr;
model:
hd84 on male (beta1)
black hisp es fh123 hsdrp;
[hd84] (beta0);
hd84 (v);
hd84#1 on male (gamma1)
black hisp es fh123 hsdrp;
[hd84#1] (gamma0);
model constraint:
! expectation for covs=1,0,0,0,0,0
new(pi z tobitp p0 pgt0 fi tobit predict);
pi = 1/(1+EXP(-gamma0-gamma1)); ! logit to probability
z = (beta0+beta1)/SQRT(v);
tobitp = PHI(z);
p0 = pi + (1-pi)*(1-tobitp);
pgt0 = (1-pi)*tobitp;
fi = 0.3989* EXP(-z*z/2);
tobit = beta0 + beta1 + SQRT(v)*fi/tobitp;
predict = (1-p0)*tobitp*tobit;
output:
tech1 tech8 standardized sampstat;
*** 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 1
Number of independent variables 6
Number of continuous latent variables 0
Observed dependent variables
Censored
HD84
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 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
Cholesky OFF
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
COVARIANCE COVERAGE OF DATA
Minimum covariance coverage value 0.100
SUMMARY OF CENSORED LIMITS
HD84 0.000
SAMPLE STATISTICS
SAMPLE STATISTICS
Means
MALE BLACK HISP ES FH123
________ ________ ________ ________ ________
1 0.513 0.271 0.160 0.146 0.106
Means
HSDRP
________
1 0.195
Covariances
MALE BLACK HISP ES FH123
________ ________ ________ ________ ________
MALE 0.250
BLACK -0.004 0.197
HISP 0.004 -0.043 0.135
ES 0.015 -0.018 0.005 0.124
FH123 -0.011 -0.009 0.003 0.004 0.095
HSDRP 0.008 -0.004 0.018 0.005 -0.002
Covariances
HSDRP
________
HSDRP 0.157
Correlations
MALE BLACK HISP ES FH123
________ ________ ________ ________ ________
MALE 1.000
BLACK -0.019 1.000
HISP 0.020 -0.266 1.000
ES 0.083 -0.113 0.040 1.000
FH123 -0.071 -0.064 0.026 0.034 1.000
HSDRP 0.038 -0.024 0.125 0.032 -0.020
Correlations
HSDRP
________
HSDRP 1.000
UNIVARIATE SAMPLE STATISTICS
UNIVARIATE HIGHER-ORDER MOMENT DESCRIPTIVE STATISTICS
Variable/ Mean/ Skewness/ Minimum/ % with Percentiles
Sample Size Variance Kurtosis Maximum Min/Max 20%/60% 40%/80% Median
MALE 0.513 -0.050 0.000 48.74% 0.000 0.000 1.000
1153.000 0.250 -1.997 1.000 51.26% 1.000 1.000
BLACK 0.271 1.033 0.000 72.94% 0.000 0.000 0.000
1153.000 0.197 -0.934 1.000 27.06% 0.000 1.000
HISP 0.160 1.850 0.000 83.95% 0.000 0.000 0.000
1153.000 0.135 1.424 1.000 16.05% 0.000 0.000
ES 0.146 2.008 0.000 85.43% 0.000 0.000 0.000
1153.000 0.124 2.034 1.000 14.57% 0.000 0.000
FH123 0.106 2.563 0.000 89.42% 0.000 0.000 0.000
1153.000 0.095 4.569 1.000 10.58% 0.000 0.000
HSDRP 0.195 1.538 0.000 80.49% 0.000 0.000 0.000
1153.000 0.157 0.367 1.000 19.51% 0.000 0.000
THE MODEL ESTIMATION TERMINATED NORMALLY
MODEL FIT INFORMATION
Number of Free Parameters 15
Loglikelihood
H0 Value -1499.409
H0 Scaling Correction Factor 0.9951
for MLR
Information Criteria
Akaike (AIC) 3028.818
Bayesian (BIC) 3104.570
Sample-Size Adjusted BIC 3056.925
(n* = (n + 2) / 24)
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
HD84 ON
MALE 0.957 0.236 4.057 0.000
BLACK -1.150 0.282 -4.073 0.000
HISP -0.405 0.320 -1.264 0.206
ES 0.585 0.276 2.120 0.034
FH123 -0.031 0.329 -0.095 0.924
HSDRP 0.390 0.263 1.487 0.137
HD84#1 ON
MALE -1.025 0.166 -6.157 0.000
BLACK 0.962 0.208 4.621 0.000
HISP 0.570 0.215 2.651 0.008
ES -0.204 0.198 -1.032 0.302
FH123 -0.512 0.273 -1.876 0.061
HSDRP 0.040 0.188 0.213 0.831
Intercepts
HD84#1 0.412 0.145 2.848 0.004
HD84 1.567 0.189 8.290 0.000
Residual Variances
HD84 3.168 0.286 11.089 0.000
New/Additional Parameters
PI 0.352 0.034 10.490 0.000
Z 1.418 0.126 11.258 0.000
TOBITP 0.922 0.018 50.143 0.000
P0 0.402 0.028 14.564 0.000
PGT0 0.598 0.028 21.650 0.000
FI 0.146 0.026 5.598 0.000
TOBIT 2.806 0.132 21.182 0.000
PREDICT 1.547 0.115 13.446 0.000
STANDARDIZED MODEL RESULTS
STDYX Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
HD84 ON
MALE 0.247 0.057 4.304 0.000
BLACK -0.264 0.059 -4.482 0.000
HISP -0.077 0.060 -1.273 0.203
ES 0.107 0.050 2.122 0.034
FH123 -0.005 0.052 -0.095 0.924
HSDRP 0.080 0.053 1.494 0.135
HD84#1 ON
MALE -0.264 0.040 -6.633 0.000
BLACK 0.220 0.045 4.850 0.000
HISP 0.108 0.040 2.679 0.007
ES -0.037 0.036 -1.032 0.302
FH123 -0.081 0.043 -1.890 0.059
HSDRP 0.008 0.038 0.213 0.831
Intercepts
HD84#1 0.212 0.074 2.885 0.004
HD84 0.808 0.119 6.803 0.000
Residual Variances
HD84 0.843 0.041 20.573 0.000
STDY Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
HD84 ON
MALE 0.494 0.114 4.319 0.000
BLACK -0.593 0.132 -4.499 0.000
HISP -0.209 0.164 -1.274 0.203
ES 0.302 0.142 2.124 0.034
FH123 -0.016 0.170 -0.095 0.924
HSDRP 0.201 0.135 1.495 0.135
HD84#1 ON
MALE -0.528 0.079 -6.688 0.000
BLACK 0.495 0.102 4.873 0.000
HISP 0.293 0.109 2.684 0.007
ES -0.105 0.102 -1.032 0.302
FH123 -0.263 0.139 -1.891 0.059
HSDRP 0.021 0.097 0.213 0.831
Intercepts
HD84#1 0.212 0.074 2.885 0.004
HD84 0.808 0.119 6.803 0.000
Residual Variances
HD84 0.843 0.041 20.573 0.000
STD Standardization
Two-Tailed
Estimate S.E. Est./S.E. P-Value
HD84 ON
MALE 0.957 0.236 4.057 0.000
BLACK -1.150 0.282 -4.073 0.000
HISP -0.405 0.320 -1.264 0.206
ES 0.585 0.276 2.120 0.034
FH123 -0.031 0.329 -0.095 0.924
HSDRP 0.390 0.263 1.487 0.137
HD84#1 ON
MALE -1.025 0.166 -6.157 0.000
BLACK 0.962 0.208 4.621 0.000
HISP 0.570 0.215 2.651 0.008
ES -0.204 0.198 -1.032 0.302
FH123 -0.512 0.273 -1.876 0.061
HSDRP 0.040 0.188 0.213 0.831
Intercepts
HD84#1 0.412 0.145 2.848 0.004
HD84 1.567 0.189 8.290 0.000
Residual Variances
HD84 3.168 0.286 11.089 0.000
R-SQUARE
Observed Two-Tailed
Variable Estimate S.E. Est./S.E. P-Value
HD84 0.128 0.030 4.281 0.000
HD84 0.157 0.041 3.822 0.000
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.290E-02
(ratio of smallest to largest eigenvalue)
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION
NU
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
1 0 0 0 0 0
NU
ES FH123 HSDRP
________ ________ ________
1 0 0 0
LAMBDA
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
HD84#1 0 0 0 0 0
HD84 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
________ ________ ________
HD84#1 0 0 0
HD84 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
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
HD84#1 0
HD84 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
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
1 1 2 0 0 0
ALPHA
ES FH123 HSDRP
________ ________ ________
1 0 0 0
BETA
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
HD84#1 0 0 3 4 5
HD84 0 0 9 10 11
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
________ ________ ________
HD84#1 6 7 8
HD84 12 13 14
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
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
HD84#1 0
HD84 0 15
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
PI Z TOBITP P0 PGT0
________ ________ ________ ________ ________
1 16 17 18 19 20
NEW/ADDITIONAL PARAMETERS
FI TOBIT PREDICT
________ ________ ________
1 21 22 23
STARTING VALUES
NU
HD84#1 HD84 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
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
HD84#1 1.000 0.000 0.000 0.000 0.000
HD84 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
________ ________ ________
HD84#1 0.000 0.000 0.000
HD84 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
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
HD84#1 0.000
HD84 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
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
1 0.424 1.029 0.000 0.000 0.000
ALPHA
ES FH123 HSDRP
________ ________ ________
1 0.000 0.000 0.000
BETA
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
HD84#1 0.000 0.000 0.000 0.000 0.000
HD84 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
________ ________ ________
HD84#1 0.000 0.000 0.000
HD84 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
HD84#1 HD84 MALE BLACK HISP
________ ________ ________ ________ ________
HD84#1 0.000
HD84 0.000 1.272
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
PI Z TOBITP P0 PGT0
________ ________ ________ ________ ________
1 0.500 0.500 0.500 0.500 0.500
NEW/ADDITIONAL PARAMETERS
FI TOBIT PREDICT
________ ________ ________
1 0.500 0.500 0.500
TECHNICAL 8 OUTPUT
E STEP ITER LOGLIKELIHOOD ABS CHANGE REL CHANGE ALGORITHM
1 -0.20423926D+04 0.0000000 0.0000000 EM
2 -0.14994474D+04 542.9452444 0.2658378 EM
3 -0.14994092D+04 0.0382116 0.0000255 EM
4 -0.14994092D+04 0.0000228 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:11:37
Ending Time: 12:11:38
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