Mplus VERSION 6
MUTHEN & MUTHEN
04/25/2010 10:57 PM
INPUT INSTRUCTIONS
! SCRIPT NAME : rawVC1b (dp)
! GOAL : To evaluate best model for variance components
! DATA : continuous
! INPUT : raw data
! UNI/BI/MULTI : uni
! DATA-GROUPS : MZ DZ
! MEANS MODEL : grand mean, age effect, sex effect
! VARIANCE COVARIANCE MODEL(S) : 1.ADE 2.AE 3.E
data: file is example.dat;
variable:
names are country famid zygos sex1 age1 height1 weight1 bmi1
sex2 age2 height2 weight2 bmi2;
usevar are bmi1 sex1 age1 bmi2 sex2 age2 g;
grouping=g(1=MZ 2=DZ); ! specify the two groups MZ and DZ
missing=all(-1); ! specify missing data symbol
define: if (zygos==1 .or. zygos==3) then g=1 else g=2; ! defines the two groups
model :
bmi1 on sex1 (b1)
age1 (b2);
bmi2 on sex2 (b1)
age2 (b2);
[bmi1 bmi2] (m);
bmi1 bmi2 (v);
bmi1 with bmi2 (c1);
model DZ:
bmi1 with bmi2 (c2);
model constraint:
! Uncomment for Model ADE
new(a d e x w z);
a=x*x;
d=w*w;
e=z*z;
v=a+d+e;
c1=a+d;
c2=0.5*a+0.25*d;
! Uncomment for Model AE
! d=0;
! Uncomment for Model E
! a=0; d=0;
*** 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: 3
1 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS
SUMMARY OF ANALYSIS
Number of groups 2
Number of observations
Group MZ 138
Group DZ 166
Number of dependent variables 2
Number of independent variables 4
Number of continuous latent variables 0
Observed dependent variables
Continuous
BMI1 BMI2
Observed independent variables
SEX1 AGE1 SEX2 AGE2
Variables with special functions
Grouping variable G
Estimator ML
Information matrix OBSERVED
Maximum number of iterations 1000
Convergence criterion 0.500D-04
Maximum number of steepest descent iterations 20
Maximum number of iterations for H1 2000
Convergence criterion for H1 0.100D-03
Input data file(s)
example.dat
Input data format FREE
SUMMARY OF DATA
Group MZ
Number of missing data patterns 3
Group DZ
Number of missing data patterns 3
COVARIANCE COVERAGE OF DATA
Minimum covariance coverage value 0.100
PROPORTION OF DATA PRESENT FOR MZ
Covariance Coverage
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.935
BMI2 0.891 0.957
SEX1 0.935 0.957 1.000
AGE1 0.935 0.957 1.000 1.000
SEX2 0.935 0.957 1.000 1.000 1.000
AGE2 0.935 0.957 1.000 1.000 1.000
Covariance Coverage
AGE2
________
AGE2 1.000
PROPORTION OF DATA PRESENT FOR DZ
Covariance Coverage
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.940
BMI2 0.843 0.904
SEX1 0.940 0.904 1.000
AGE1 0.940 0.904 1.000 1.000
SEX2 0.940 0.904 1.000 1.000 1.000
AGE2 0.940 0.904 1.000 1.000 1.000
Covariance Coverage
AGE2
________
AGE2 1.000
THE MODEL ESTIMATION TERMINATED NORMALLY
TESTS OF MODEL FIT
Chi-Square Test of Model Fit
Value 28.544
Degrees of Freedom 20
P-Value 0.0971
Chi-Square Contributions From Each Group
MZ 11.178
DZ 17.366
Chi-Square Test of Model Fit for the Baseline Model
Value 188.497
Degrees of Freedom 18
P-Value 0.0000
CFI/TLI
CFI 0.950
TLI 0.955
Loglikelihood
H0 Value -4325.956
H1 Value -4311.684
Information Criteria
Number of Free Parameters 6
Akaike (AIC) 8663.913
Bayesian (BIC) 8686.215
Sample-Size Adjusted BIC 8667.186
(n* = (n + 2) / 24)
RMSEA (Root Mean Square Error Of Approximation)
Estimate 0.053
90 Percent C.I. 0.000 0.094
Probability RMSEA <= .05 0.418
SRMR (Standardized Root Mean Square Residual)
Value 0.119
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
Group MZ
BMI1 ON
SEX1 0.866 0.343 2.529 0.011
AGE1 0.097 0.014 6.980 0.000
BMI2 ON
SEX2 0.866 0.343 2.529 0.011
AGE2 0.097 0.014 6.980 0.000
BMI1 WITH
BMI2 9.665 0.929 10.399 0.000
Intercepts
BMI1 20.203 0.577 35.002 0.000
BMI2 20.203 0.577 35.002 0.000
Residual Variances
BMI1 13.019 0.857 15.194 0.000
BMI2 13.019 0.857 15.194 0.000
Group DZ
BMI1 ON
SEX1 0.866 0.343 2.529 0.011
AGE1 0.097 0.014 6.980 0.000
BMI2 ON
SEX2 0.866 0.343 2.529 0.011
AGE2 0.097 0.014 6.980 0.000
BMI1 WITH
BMI2 3.490 0.988 3.532 0.000
Intercepts
BMI1 20.203 0.577 35.002 0.000
BMI2 20.203 0.577 35.002 0.000
Residual Variances
BMI1 13.019 0.857 15.194 0.000
BMI2 13.019 0.857 15.194 0.000
New/Additional Parameters
A 4.296 3.725 1.153 0.249
D 5.370 3.723 1.442 0.149
E 3.354 0.430 7.807 0.000
X 2.073 0.899 2.306 0.021
W -2.317 0.803 -2.885 0.004
Z 1.831 0.117 15.613 0.000
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.160E-04
(ratio of smallest to largest eigenvalue)
Beginning Time: 22:57:41
Ending Time: 22:57:42
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-2010 Muthen & Muthen
Back to the list of genetics examples