Mplus VERSION 6
MUTHEN & MUTHEN
04/25/2010 10:57 PM
INPUT INSTRUCTIONS
! SCRIPT NAME : rawVC1a (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.ACE 2.AE 3.CE 4.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 ACE
new(a c e x y z);
a=x*x;
c=y*y;
e=z*z;
v=a+c+e;
c1=a+c;
c2=0.5*a+c;
! Uncomment for Model AE
! c=0;
! Uncomment for Model CE
! a=0;
! Uncomment for Model E
! a=0; c=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
THE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES MAY NOT BE
TRUSTWORTHY FOR SOME PARAMETERS DUE TO A NON-POSITIVE DEFINITE
FIRST-ORDER DERIVATIVE PRODUCT MATRIX. THIS MAY BE DUE TO THE STARTING
VALUES BUT MAY ALSO BE AN INDICATION OF MODEL NONIDENTIFICATION. THE
CONDITION NUMBER IS 0.420D-13. PROBLEM INVOLVING PARAMETER 11.
TESTS OF MODEL FIT
Chi-Square Test of Model Fit
Value 30.813
Degrees of Freedom 20
P-Value 0.0577
Chi-Square Contributions From Each Group
MZ 11.353
DZ 19.460
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.937
TLI 0.943
Loglikelihood
H0 Value -4327.091
H1 Value -4311.684
Information Criteria
Number of Free Parameters 6
Akaike (AIC) 8666.182
Bayesian (BIC) 8688.484
Sample-Size Adjusted BIC 8669.455
(n* = (n + 2) / 24)
RMSEA (Root Mean Square Error Of Approximation)
Estimate 0.060
90 Percent C.I. 0.000 0.099
Probability RMSEA <= .05 0.321
SRMR (Standardized Root Mean Square Residual)
Value 0.120
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
Group MZ
BMI1 ON
SEX1 0.837 0.340 2.463 0.014
AGE1 0.097 0.014 6.773 0.000
BMI2 ON
SEX2 0.837 0.340 2.463 0.014
AGE2 0.097 0.014 6.773 0.000
BMI1 WITH
BMI2 9.655 0.991 9.739 0.000
Intercepts
BMI1 20.220 0.590 34.254 0.000
BMI2 20.220 0.590 34.254 0.000
Residual Variances
BMI1 13.197 0.882 14.966 0.000
BMI2 13.197 0.882 14.966 0.000
Group DZ
BMI1 ON
SEX1 0.837 0.340 2.463 0.014
AGE1 0.097 0.014 6.773 0.000
BMI2 ON
SEX2 0.837 0.340 2.463 0.014
AGE2 0.097 0.014 6.773 0.000
BMI1 WITH
BMI2 4.827 0.496 9.739 0.000
Intercepts
BMI1 20.220 0.590 34.254 0.000
BMI2 20.220 0.590 34.254 0.000
Residual Variances
BMI1 13.197 0.882 14.966 0.000
BMI2 13.197 0.882 14.966 0.000
New/Additional Parameters
A 9.655 0.991 9.739 0.000
C 0.000 0.000 0.000 1.000
E 3.542 0.456 7.764 0.000
X 3.107 0.160 19.477 0.000
Y 0.000 0.755 0.000 1.000
Z 1.882 0.121 15.529 0.000
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.205E-04
(ratio of smallest to largest eigenvalue)
TECHNICAL 1 OUTPUT
PARAMETER SPECIFICATION FOR MZ
NU
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
1 0 0 0 0 0
NU
AGE2
________
1 0
LAMBDA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0 0 0 0 0
BMI2 0 0 0 0 0
SEX1 0 0 0 0 0
AGE1 0 0 0 0 0
SEX2 0 0 0 0 0
AGE2 0 0 0 0 0
LAMBDA
AGE2
________
BMI1 0
BMI2 0
SEX1 0
AGE1 0
SEX2 0
AGE2 0
THETA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0
BMI2 0 0
SEX1 0 0 0
AGE1 0 0 0 0
SEX2 0 0 0 0 0
AGE2 0 0 0 0 0
THETA
AGE2
________
AGE2 0
ALPHA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
1 1 1 0 0 0
ALPHA
AGE2
________
1 0
BETA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0 0 2 3 0
BMI2 0 0 0 0 2
SEX1 0 0 0 0 0
AGE1 0 0 0 0 0
SEX2 0 0 0 0 0
AGE2 0 0 0 0 0
BETA
AGE2
________
BMI1 0
BMI2 3
SEX1 0
AGE1 0
SEX2 0
AGE2 0
PSI
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 4
BMI2 5 4
SEX1 0 0 0
AGE1 0 0 0 0
SEX2 0 0 0 0 0
AGE2 0 0 0 0 0
PSI
AGE2
________
AGE2 0
PARAMETER SPECIFICATION FOR DZ
NU
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
1 0 0 0 0 0
NU
AGE2
________
1 0
LAMBDA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0 0 0 0 0
BMI2 0 0 0 0 0
SEX1 0 0 0 0 0
AGE1 0 0 0 0 0
SEX2 0 0 0 0 0
AGE2 0 0 0 0 0
LAMBDA
AGE2
________
BMI1 0
BMI2 0
SEX1 0
AGE1 0
SEX2 0
AGE2 0
THETA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0
BMI2 0 0
SEX1 0 0 0
AGE1 0 0 0 0
SEX2 0 0 0 0 0
AGE2 0 0 0 0 0
THETA
AGE2
________
AGE2 0
ALPHA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
1 1 1 0 0 0
ALPHA
AGE2
________
1 0
BETA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0 0 2 3 0
BMI2 0 0 0 0 2
SEX1 0 0 0 0 0
AGE1 0 0 0 0 0
SEX2 0 0 0 0 0
AGE2 0 0 0 0 0
BETA
AGE2
________
BMI1 0
BMI2 3
SEX1 0
AGE1 0
SEX2 0
AGE2 0
PSI
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 4
BMI2 6 4
SEX1 0 0 0
AGE1 0 0 0 0
SEX2 0 0 0 0 0
AGE2 0 0 0 0 0
PSI
AGE2
________
AGE2 0
PARAMETER SPECIFICATION FOR THE ADDITIONAL PARAMETERS
New/Additional Parameters
A C E X Y
________ ________ ________ ________ ________
1 7 8 9 10 11
New/Additional Parameters
Z
________
1 12
STARTING VALUES FOR MZ
NU
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
1 0.000 0.000 0.000 0.000 0.000
NU
AGE2
________
1 0.000
LAMBDA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 1.000 0.000 0.000 0.000 0.000
BMI2 0.000 1.000 0.000 0.000 0.000
SEX1 0.000 0.000 1.000 0.000 0.000
AGE1 0.000 0.000 0.000 1.000 0.000
SEX2 0.000 0.000 0.000 0.000 1.000
AGE2 0.000 0.000 0.000 0.000 0.000
LAMBDA
AGE2
________
BMI1 0.000
BMI2 0.000
SEX1 0.000
AGE1 0.000
SEX2 0.000
AGE2 1.000
THETA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.000
BMI2 0.000 0.000
SEX1 0.000 0.000 0.000
AGE1 0.000 0.000 0.000 0.000
SEX2 0.000 0.000 0.000 0.000 0.000
AGE2 0.000 0.000 0.000 0.000 0.000
THETA
AGE2
________
AGE2 0.000
ALPHA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
1 0.000 0.000 0.348 36.909 0.355
ALPHA
AGE2
________
1 36.568
BETA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.000 0.000 0.000 0.000 0.000
BMI2 0.000 0.000 0.000 0.000 0.000
SEX1 0.000 0.000 0.000 0.000 0.000
AGE1 0.000 0.000 0.000 0.000 0.000
SEX2 0.000 0.000 0.000 0.000 0.000
AGE2 0.000 0.000 0.000 0.000 0.000
BETA
AGE2
________
BMI1 0.000
BMI2 0.000
SEX1 0.000
AGE1 0.000
SEX2 0.000
AGE2 0.000
PSI
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 6.683
BMI2 0.000 6.326
SEX1 0.000 0.000 0.444
AGE1 0.000 0.000 1.958 208.238
SEX2 0.000 0.000 0.210 -0.983 0.446
AGE2 0.000 0.000 -0.993 154.076 2.931
PSI
AGE2
________
AGE2 220.236
STARTING VALUES FOR DZ
NU
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
1 0.000 0.000 0.000 0.000 0.000
NU
AGE2
________
1 0.000
LAMBDA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 1.000 0.000 0.000 0.000 0.000
BMI2 0.000 1.000 0.000 0.000 0.000
SEX1 0.000 0.000 1.000 0.000 0.000
AGE1 0.000 0.000 0.000 1.000 0.000
SEX2 0.000 0.000 0.000 0.000 1.000
AGE2 0.000 0.000 0.000 0.000 0.000
LAMBDA
AGE2
________
BMI1 0.000
BMI2 0.000
SEX1 0.000
AGE1 0.000
SEX2 0.000
AGE2 1.000
THETA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.000
BMI2 0.000 0.000
SEX1 0.000 0.000 0.000
AGE1 0.000 0.000 0.000 0.000
SEX2 0.000 0.000 0.000 0.000 0.000
AGE2 0.000 0.000 0.000 0.000 0.000
THETA
AGE2
________
AGE2 0.000
ALPHA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
1 0.000 0.000 0.325 35.344 0.181
ALPHA
AGE2
________
1 33.743
BETA
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.000 0.000 0.000 0.000 0.000
BMI2 0.000 0.000 0.000 0.000 0.000
SEX1 0.000 0.000 0.000 0.000 0.000
AGE1 0.000 0.000 0.000 0.000 0.000
SEX2 0.000 0.000 0.000 0.000 0.000
AGE2 0.000 0.000 0.000 0.000 0.000
BETA
AGE2
________
BMI1 0.000
BMI2 0.000
SEX1 0.000
AGE1 0.000
SEX2 0.000
AGE2 0.000
PSI
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 6.683
BMI2 0.000 6.326
SEX1 0.000 0.000 0.509
AGE1 0.000 0.000 4.240 229.798
SEX2 0.000 0.000 -0.011 -1.991 0.618
AGE2 0.000 0.000 -0.830 134.283 6.002
PSI
AGE2
________
AGE2 261.295
STARTING VALUES FOR THE ADDITIONAL PARAMETERS
New/Additional Parameters
A C E X Y
________ ________ ________ ________ ________
1 0.500 0.500 0.500 0.500 0.500
New/Additional Parameters
Z
________
1 0.500
Beginning Time: 22:57:41
Ending Time: 22:57:41
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-2010 Muthen & Muthen
Back to the list of genetics examples