Mplus VERSION 6
MUTHEN & MUTHEN
04/25/2010 10:57 PM
INPUT INSTRUCTIONS
! SCRIPT NAME : rawSAT2 (dp)
! GOAL : To calculate saturated statistics & test assumptions, including se
! DATA : continuous
! INPUT : raw data
! UNI/BI/MULTI : uni
! DATA-GROUPS : MZM DZM MZF DZF DOSMF DOSFM
! MEANS MODEL : grand mean (mztwin1, mztwin2, dztwin1, dztwin2, dostwin1, dostwin2)
! VARIANCE COVARIANCE MODEL(S) :
! saturated model:
! variances: (mzmt1, mzmt2, dzmt1, dzmt2, mzft1, mzft2, dzft1, dzft2, dosm1, dosf2, dosf1,
! covariances: mzm, dzm, mzf. dzf, dosmf, dosfm
! 1.Test whether gm_mz1 = gm_mz2 AND gm_dz1 = gm_dz2 AND gm_dos1 = gm_dos2 (birth order ef
! 2.Test whether grand means MZ = grand means DZ (zygosity effect on means)
! 3.Test significance of sex effect
! 4.Test significance of age effect
! 5.Test whether variance twin1=variance twin2, within zygosity (birth order effect on var
! 6.Test whether variances are equal across zygosity (zygosity effect on variance )
! 7.Test whether variances are equal across sexes (sex effect on variance )
! 8.Test whether covariances are equal across sexes
data: file is example.dat;
variable:
names = country famid zygos sex1 age1 height1 weight1 bmi1
sex2 age2 height2 weight2 bmi2;
usevar = 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);
model MZM : [bmi1] (m11); [bmi2] (m12);
bmi1 (v11); bmi2 (v12);
bmi1 with bmi2 (c1);
model DZM : [bmi1] (m21); [bmi2] (m22);
bmi1 (v21); bmi2 (v22);
bmi1 with bmi2 (c2);
model MZF : [bmi1] (m31); [bmi2] (m32);
bmi1 (v31); bmi2 (v32);
bmi1 with bmi2 (c3);
model DZF : [bmi1] (m41); [bmi2] (m42);
bmi1 (v41); bmi2 (v42);
bmi1 with bmi2 (c4);
model DOSMF : [bmi1] (m51); [bmi2] (m52);
bmi1 (v51); bmi2 (v52);
bmi1 with bmi2 (c5);
model DOSFM : [bmi1] (m61); [bmi2] (m62);
bmi1 (v61); bmi2 (v62);
bmi1 with bmi2 (c6);
model constraint:
m31=m11;
m32=m12;
m41=m21;
m42=m22;
m61=m51;
m62=m52;
! Uncomment for test 1
! m11=m12;
! m21=m22;
! m51=m52;
! Uncomment for test 2
! m11=m12;
! m21=m22;
! m51=m52;
! m21=m11;
! m51=m11;
! Uncomment for test 3
! b1=0;
! Uncomment for test 4
! b2=0;
! Uncomment for test 5
! v11=v12;
! v21=v22;
! v31=v32;
! v41=v42;
! v51=v62;
! v61=v52;
! Uncomment for test 6
! v21=v11;
! v22=v12;
! v41=v31;
! v42=v32;
! v51=v11;
! v52=v32;
! v61=v31;
! v62=v12;
! Uncomment for test 7
! v31=v11;
! v32=v12;
! v41=v21;
! v42=v22;
! v61=v51;
! v62=v52;
! Uncomment for test 8
! c3=c1;
! c4=c2;
! c5=c2;
! c6=c2;
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! GOAL : To calculate saturated statistics & test assumptions, including sex effects on va
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! MEANS MODEL : grand mean (mztwin1, mztwin2, dztwin1, dztwin2, dostwin1, dostwin2), age ef
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! 1.Test whether gm_mz1 = gm_mz2 AND gm_dz1 = gm_dz2 AND gm_dos1 = gm_dos2 (birth order eff
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! 5.Test whether variance twin1=variance twin2, within zygosity (birth order effect on vari
*** 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
5 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS
SUMMARY OF ANALYSIS
Number of groups 6
Number of observations
Group MZM 61
Group DZM 39
Group MZF 77
Group DZF 67
Group DOSMF 36
Group DOSFM 24
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 3
Group MZF
Number of missing data patterns 3
Group DZF
Number of missing data patterns 3
Group DOSMF
Number of missing data patterns 3
Group DOSFM
Number of missing data patterns 3
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.897
BMI2 0.744 0.846
SEX1 0.897 0.846 1.000
AGE1 0.897 0.846 1.000 1.000
SEX2 0.897 0.846 1.000 1.000 1.000
AGE2 0.897 0.846 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.970
BMI2 0.940 0.970
SEX1 0.970 0.970 1.000
AGE1 0.970 0.970 1.000 1.000
SEX2 0.970 0.970 1.000 1.000 1.000
AGE2 0.970 0.970 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.958
BMI2 0.792 0.833
SEX1 0.958 0.833 1.000
AGE1 0.958 0.833 1.000 1.000
SEX2 0.958 0.833 1.000 1.000 1.000
AGE2 0.958 0.833 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 24.045
Degrees of Freedom 52
P-Value 0.9997
Chi-Square Contributions From Each Group
MZM 4.943
DZM 5.363
MZF 5.755
DZF 2.052
DOSMF 2.472
DOSFM 3.460
Chi-Square Test of Model Fit for the Baseline Model
Value 251.891
Degrees of Freedom 54
P-Value 0.0000
CFI/TLI
CFI 1.000
TLI 1.147
Loglikelihood
H0 Value -3757.818
H1 Value -3745.795
Information Criteria
Number of Free Parameters 26
Akaike (AIC) 7567.635
Bayesian (BIC) 7664.278
Sample-Size Adjusted BIC 7581.819
(n* = (n + 2) / 24)
RMSEA (Root Mean Square Error Of Approximation)
Estimate 0.000
90 Percent C.I. 0.000 0.000
Probability RMSEA <= .05 1.000
SRMR (Standardized Root Mean Square Residual)
Value 0.189
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
Group MZM
BMI1 ON
SEX1 0.862 0.343 2.511 0.012
AGE1 0.097 0.014 6.999 0.000
BMI2 ON
SEX2 0.862 0.343 2.511 0.012
AGE2 0.097 0.014 6.999 0.000
BMI2 WITH
BMI1 9.007 1.822 4.944 0.000
Intercepts
BMI1 20.100 0.639 31.466 0.000
BMI2 19.608 0.632 31.031 0.000
Residual Variances
BMI1 11.366 2.181 5.213 0.000
BMI2 9.822 1.787 5.498 0.000
Group DZM
BMI1 ON
SEX1 0.862 0.343 2.511 0.012
AGE1 0.097 0.014 6.999 0.000
BMI2 ON
SEX2 0.862 0.343 2.511 0.012
AGE2 0.097 0.014 6.999 0.000
BMI2 WITH
BMI1 2.661 2.197 1.211 0.226
Intercepts
BMI1 20.520 0.679 30.212 0.000
BMI2 20.142 0.637 31.606 0.000
Residual Variances
BMI1 12.333 2.944 4.189 0.000
BMI2 12.437 3.105 4.005 0.000
Group MZF
BMI1 ON
SEX1 0.862 0.343 2.511 0.012
AGE1 0.097 0.014 6.999 0.000
BMI2 ON
SEX2 0.862 0.343 2.511 0.012
AGE2 0.097 0.014 6.999 0.000
BMI2 WITH
BMI1 8.564 1.800 4.757 0.000
Intercepts
BMI1 20.100 0.639 31.466 0.000
BMI2 19.608 0.632 31.031 0.000
Residual Variances
BMI1 13.224 2.161 6.120 0.000
BMI2 12.877 2.123 6.066 0.000
Group DZF
BMI1 ON
SEX1 0.862 0.343 2.511 0.012
AGE1 0.097 0.014 6.999 0.000
BMI2 ON
SEX2 0.862 0.343 2.511 0.012
AGE2 0.097 0.014 6.999 0.000
BMI2 WITH
BMI1 4.212 1.801 2.339 0.019
Intercepts
BMI1 20.520 0.679 30.212 0.000
BMI2 20.142 0.637 31.606 0.000
Residual Variances
BMI1 18.748 3.289 5.701 0.000
BMI2 10.483 1.841 5.695 0.000
Group DOSMF
BMI1 ON
SEX1 0.862 0.343 2.511 0.012
AGE1 0.097 0.014 6.999 0.000
BMI2 ON
SEX2 0.862 0.343 2.511 0.012
AGE2 0.097 0.014 6.999 0.000
BMI2 WITH
BMI1 5.359 2.480 2.161 0.031
Intercepts
BMI1 20.550 0.761 27.016 0.000
BMI2 20.596 0.694 29.694 0.000
Residual Variances
BMI1 15.427 3.771 4.091 0.000
BMI2 11.803 2.938 4.017 0.000
Group DOSFM
BMI1 ON
SEX1 0.862 0.343 2.511 0.012
AGE1 0.097 0.014 6.999 0.000
BMI2 ON
SEX2 0.862 0.343 2.511 0.012
AGE2 0.097 0.014 6.999 0.000
BMI2 WITH
BMI1 1.522 2.541 0.599 0.549
Intercepts
BMI1 20.550 0.761 27.016 0.000
BMI2 20.596 0.694 29.694 0.000
Residual Variances
BMI1 13.859 4.091 3.388 0.001
BMI2 10.095 3.225 3.130 0.002
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.300E-05
(ratio of smallest to largest eigenvalue)
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