Mplus VERSION 5.2
MUTHEN & MUTHEN
12/02/2008 3:05 PM
INPUT INSTRUCTIONS
! SCRIPT NAME : rawVC3b (dp)
! GOAL : To evaluate whether different genes are important for males and fe
! DATA : continuous
! INPUT : raw data
! UNI/BI/MULTI : uni
! DATA-GROUPS : MZM DZM MZF DZF DOSMF DOSFM
! MEANS MODEL : grand mean, age effect, sex effect
! VARIANCE COVARIANCE MODEL(S) : 1.ADE males + ADE females + DOS A/D-correlation is
! : 2.ADE males = ADE females + DOS A/D-correlation is not
! : 3.ADE males = ADE females + DOS A-correlation = DZ A-co
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;
grouping=zygos(1=MZM 2=DZM 3=MZF 4=DZF 5=DOSMF 6=DOSFM); ! specify the groups
missing=all(-1); ! specify missing data symbol
model :
bmi1 on sex1 (b1)
age1 (b2);
bmi2 on sex2 (b1)
age2 (b2);
[bmi1 bmi2] (m);
model MZM :
bmi1 bmi2 (mv);
bmi1 with bmi2 (mc1);
model DZM :
bmi1 bmi2 (mv);
bmi1 with bmi2 (mc2);
model MZF :
bmi1 bmi2 (fv);
bmi1 with bmi2 (fc1);
model DZF :
bmi1 bmi2 (fv);
bmi1 with bmi2 (fc2);
model DOSMF :
bmi1 (mv); bmi2 (fv);
bmi1 with bmi2 (mfc);
model DOSFM :
bmi1 (fv); bmi2 (mv);
bmi1 with bmi2 (mfc);
model constraint:
new(ma md me mx mw mz);
ma=mx*mx;
md=mw*mw;
me=mz*mz;
mv=ma+md+me;
mc1=ma+md;
mc2=0.5*ma+0.25*md;
new(fa fd fe fx fw fz);
fa=fx*fx;
fd=fw*fw;
fe=fz*fz;
fv=fa+fd+fe;
fc1=fa+fd;
fc2=0.5*fa+0.25*fd;
new(f*0.1 j*0.1); j>0; j<0.25;
mfc=f*mx*fx+j*mw*fw; ! f and j can not be identified simultaniously
! in this example they will get fixed at their boundaries
f=0.5; ! we fix this parameter to identify the model
! Uncomment for Model ADE with all components equal
! ma=fa;
! md=fd;
! me=fe;
! Uncomment for Model ACE with all components equal and DOS correlations = DZ correlations
! ma=fa;
! md=fd;
! me=fe;
! f=0.5;
! j=0.25
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! GOAL : To evaluate whether different genes are important for males and females, ADE m
1 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS
SUMMARY OF ANALYSIS
Number of groups 6
Number of observations
Group MZM 61
Group DZM 40
Group MZF 77
Group DZF 68
Group DOSMF 36
Group DOSFM 25
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 ZYGOS
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 MZM
Number of missing data patterns 3
Group DZM
Number of missing data patterns 4
Group MZF
Number of missing data patterns 3
Group DZF
Number of missing data patterns 4
Group DOSMF
Number of missing data patterns 3
Group DOSFM
Number of missing data patterns 4
COVARIANCE COVERAGE OF DATA
Minimum covariance coverage value 0.100
PROPORTION OF DATA PRESENT FOR MZM
Covariance Coverage
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.902
BMI2 0.852 0.951
SEX1 0.902 0.951 1.000
AGE1 0.902 0.951 1.000 1.000
SEX2 0.902 0.951 1.000 1.000 1.000
AGE2 0.902 0.951 1.000 1.000 1.000
Covariance Coverage
AGE2
________
AGE2 1.000
PROPORTION OF DATA PRESENT FOR DZM
Covariance Coverage
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.875
BMI2 0.725 0.825
SEX1 0.875 0.825 1.000
AGE1 0.875 0.825 1.000 1.000
SEX2 0.875 0.825 1.000 1.000 1.000
AGE2 0.875 0.825 1.000 1.000 1.000
Covariance Coverage
AGE2
________
AGE2 1.000
PROPORTION OF DATA PRESENT FOR MZF
Covariance Coverage
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.961
BMI2 0.922 0.961
SEX1 0.961 0.961 1.000
AGE1 0.961 0.961 1.000 1.000
SEX2 0.961 0.961 1.000 1.000 1.000
AGE2 0.961 0.961 1.000 1.000 1.000
Covariance Coverage
AGE2
________
AGE2 1.000
PROPORTION OF DATA PRESENT FOR DZF
Covariance Coverage
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.956
BMI2 0.926 0.956
SEX1 0.956 0.956 1.000
AGE1 0.956 0.956 1.000 1.000
SEX2 0.956 0.956 1.000 1.000 1.000
AGE2 0.956 0.956 1.000 1.000 1.000
Covariance Coverage
AGE2
________
AGE2 1.000
PROPORTION OF DATA PRESENT FOR DOSMF
Covariance Coverage
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.917
BMI2 0.806 0.889
SEX1 0.917 0.889 1.000
AGE1 0.917 0.889 1.000 1.000
SEX2 0.917 0.889 1.000 1.000 1.000
AGE2 0.917 0.889 1.000 1.000 1.000
Covariance Coverage
AGE2
________
AGE2 1.000
PROPORTION OF DATA PRESENT FOR DOSFM
Covariance Coverage
BMI1 BMI2 SEX1 AGE1 SEX2
________ ________ ________ ________ ________
BMI1 0.920
BMI2 0.760 0.800
SEX1 0.920 0.800 1.000
AGE1 0.920 0.800 1.000 1.000
SEX2 0.920 0.800 1.000 1.000 1.000
AGE2 0.920 0.800 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.444D-11.
TESTS OF MODEL FIT
Chi-Square Test of Model Fit
Value 43.805
Degrees of Freedom 68
P-Value 0.9901
Chi-Square Contributions From Each Group
MZM 14.536
DZM 5.440
MZF 4.341
DZF 9.485
DOSMF 6.286
DOSFM 3.717
Chi-Square Test of Model Fit for the Baseline Model
Value 204.730
Degrees of Freedom 54
P-Value 0.0000
CFI/TLI
CFI 1.000
TLI 1.127
Loglikelihood
H0 Value -3973.212
H1 Value -3951.310
Information Criteria
Number of Free Parameters 10
Akaike (AIC) 7966.425
Bayesian (BIC) 8003.693
Sample-Size Adjusted BIC 7971.978
(n* = (n + 2) / 24)
RMSEA (Root Mean Square Error Of Approximation)
Estimate 0.000
90 Percent C.I. 0.000 0.000
SRMR (Standardized Root Mean Square Residual)
Value 0.252
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
Group MZM
BMI1 ON
SEX1 0.886 0.344 2.576 0.010
AGE1 0.097 0.014 6.935 0.000
BMI2 ON
SEX2 0.886 0.344 2.576 0.010
AGE2 0.097 0.014 6.935 0.000
BMI2 WITH
BMI1 10.182 1.303 7.814 0.000
Intercepts
BMI1 20.204 0.582 34.696 0.000
BMI2 20.204 0.582 34.696 0.000
Residual Variances
BMI1 12.091 1.256 9.625 0.000
BMI2 12.091 1.256 9.625 0.000
Group DZM
BMI1 ON
SEX1 0.886 0.344 2.576 0.010
AGE1 0.097 0.014 6.935 0.000
BMI2 ON
SEX2 0.886 0.344 2.576 0.010
AGE2 0.097 0.014 6.935 0.000
BMI2 WITH
BMI1 3.847 1.177 3.269 0.001
Intercepts
BMI1 20.204 0.582 34.696 0.000
BMI2 20.204 0.582 34.696 0.000
Residual Variances
BMI1 12.091 1.256 9.625 0.000
BMI2 12.091 1.256 9.625 0.000
Group MZF
BMI1 ON
SEX1 0.886 0.344 2.576 0.010
AGE1 0.097 0.014 6.935 0.000
BMI2 ON
SEX2 0.886 0.344 2.576 0.010
AGE2 0.097 0.014 6.935 0.000
BMI2 WITH
BMI1 9.335 1.317 7.087 0.000
Intercepts
BMI1 20.204 0.582 34.696 0.000
BMI2 20.204 0.582 34.696 0.000
Residual Variances
BMI1 13.769 1.165 11.821 0.000
BMI2 13.769 1.165 11.821 0.000
Group DZF
BMI1 ON
SEX1 0.886 0.344 2.576 0.010
AGE1 0.097 0.014 6.935 0.000
BMI2 ON
SEX2 0.886 0.344 2.576 0.010
AGE2 0.097 0.014 6.935 0.000
BMI2 WITH
BMI1 4.193 1.253 3.346 0.001
Intercepts
BMI1 20.204 0.582 34.696 0.000
BMI2 20.204 0.582 34.696 0.000
Residual Variances
BMI1 13.769 1.165 11.821 0.000
BMI2 13.769 1.165 11.821 0.000
Group DOSMF
BMI1 ON
SEX1 0.886 0.344 2.576 0.010
AGE1 0.097 0.014 6.935 0.000
BMI2 ON
SEX2 0.886 0.344 2.576 0.010
AGE2 0.097 0.014 6.935 0.000
BMI2 WITH
BMI1 3.110 1.364 2.281 0.023
Intercepts
BMI1 20.204 0.582 34.696 0.000
BMI2 20.204 0.582 34.696 0.000
Residual Variances
BMI1 12.091 1.256 9.625 0.000
BMI2 13.769 1.165 11.821 0.000
Group DOSFM
BMI1 ON
SEX1 0.886 0.344 2.576 0.010
AGE1 0.097 0.014 6.935 0.000
BMI2 ON
SEX2 0.886 0.344 2.576 0.010
AGE2 0.097 0.014 6.935 0.000
BMI2 WITH
BMI1 3.110 1.364 2.281 0.023
Intercepts
BMI1 20.204 0.582 34.696 0.000
BMI2 20.204 0.582 34.696 0.000
Residual Variances
BMI1 13.769 1.165 11.821 0.000
BMI2 12.091 1.256 9.625 0.000
New/Additional Parameters
MA 5.204 4.207 1.237 0.216
MD 4.978 4.079 1.220 0.222
ME 1.908 0.381 5.011 0.000
MX 2.281 0.922 2.474 0.013
MW 2.231 0.914 2.441 0.015
MZ 1.381 0.138 10.023 0.000
FA 7.438 4.631 1.606 0.108
FD 1.898 4.609 0.412 0.681
FE 4.434 0.742 5.976 0.000
FX 2.727 0.849 3.212 0.001
FW -1.378 1.673 -0.823 0.410
FZ 2.106 0.176 11.951 0.000
F 0.500 0.000 0.000 1.000
J 0.000 0.000 1752.142 0.000
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.473E-08
(ratio of smallest to largest eigenvalue)
Beginning Time: 15:05:25
Ending Time: 15:05:27
Elapsed Time: 00:00:02
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-2008 Muthen & Muthen
Back to the list of genetics examples