Mplus VERSION 6
MUTHEN & MUTHEN
04/25/2010 10:57 PM
INPUT INSTRUCTIONS
! SCRIPT NAME : rawVC5b (dp)
! GOAL : To evaluate best model for variance components
! DATA : continuous
! INPUT : raw data
! UNI/BI/MULTI : uni
! DATA-GROUPS : MZ+sibs DZ+sibs
! MEANS MODEL : grand mean, age effect, sex effect
! VARIANCE COVARIANCE MODEL(S) : 1.ADE 2.AE 3.E
data: file is example2.dat; variance=nocheck;
variable:
names are
famnr zyg
sex1 age1 vara1 varb1
sex2 age2 vara2 varb2
sex3 age3 vara3 varb3
sex4 age4 vara4 varb4
sex5 age5 vara5 varb5
sex6 age6 vara6 varb6
sex7 age7 vara7 varb7
sex8 age8 vara8 varb8;
usevar are
sex1 sex2 sex3 sex4 sex5 sex6 sex7 sex8
age1 age2 age3 age4 age5 age6 age7 age8
vara1 vara2 vara3 vara4 vara5 vara6 vara7 vara8 g;
grouping=g(1=MZ 2=DZ); ! specify the two groups MZ and DZ
missing=all(-200); ! specify missing data symbol
define: if (zyg==1 .or. zyg==2) then g=1 else g=2; ! defines the two groups
model :
vara1 on sex1 (b1)
age1 (b2);
vara2 on sex2 (b1)
age2 (b2);
vara3 on sex3 (b1)
age3 (b2);
vara4 on sex4 (b1)
age4 (b2);
vara5 on sex5 (b1)
age5 (b2);
vara6 on sex6 (b1)
age6 (b2);
vara7 on sex7 (b1)
age7 (b2);
vara8 on sex8 (b1)
age8 (b2);
[vara1-vara8] (m);
vara1-vara8 (v);
vara1 with vara2 (c1);
vara1 with vara3-vara8 (c2);
vara2 with vara3-vara8 (c2);
vara3 with vara4-vara8 (c2);
vara4 with vara5-vara8 (c2);
vara5 with vara6-vara8 (c2);
vara6 with vara7-vara8 (c2);
vara7 with vara8 (c2);
model DZ:
vara1 with vara2 (c2);
model constraint:
! Uncomment for Model ACE
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 DE
! a=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: 51
1 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS
SUMMARY OF ANALYSIS
Number of groups 2
Number of observations
Group MZ 115
Group DZ 156
Number of dependent variables 8
Number of independent variables 16
Number of continuous latent variables 0
Observed dependent variables
Continuous
VARA1 VARA2 VARA3 VARA4 VARA5 VARA6
VARA7 VARA8
Observed independent variables
SEX1 SEX2 SEX3 SEX4 SEX5 SEX6
SEX7 SEX8 AGE1 AGE2 AGE3 AGE4
AGE5 AGE6 AGE7 AGE8
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)
example2.dat
Input data format FREE
SUMMARY OF DATA
Group MZ
Number of missing data patterns 13
Group DZ
Number of missing data patterns 13
COVARIANCE COVERAGE OF DATA
Minimum covariance coverage value 0.100
PROPORTION OF DATA PRESENT FOR MZ
Covariance Coverage
VARA1 VARA2 VARA3 VARA4 VARA5
________ ________ ________ ________ ________
VARA1 0.922
VARA2 0.887 0.930
VARA3 0.417 0.426 0.478
VARA4 0.078 0.087 0.104 0.113
VARA5 0.026 0.035 0.035 0.035 0.035
VARA6 0.009 0.017 0.017 0.017 0.017
VARA7 0.009 0.009 0.009 0.009 0.009
VARA8 0.009 0.009 0.009 0.009 0.009
SEX1 0.922 0.930 0.478 0.113 0.035
SEX2 0.922 0.930 0.478 0.113 0.035
SEX3 0.922 0.930 0.478 0.113 0.035
SEX4 0.922 0.930 0.478 0.113 0.035
SEX5 0.922 0.930 0.478 0.113 0.035
SEX6 0.922 0.930 0.478 0.113 0.035
SEX7 0.922 0.930 0.478 0.113 0.035
SEX8 0.922 0.930 0.478 0.113 0.035
AGE1 0.922 0.930 0.478 0.113 0.035
AGE2 0.922 0.930 0.478 0.113 0.035
AGE3 0.922 0.930 0.478 0.113 0.035
AGE4 0.922 0.930 0.478 0.113 0.035
AGE5 0.922 0.930 0.478 0.113 0.035
AGE6 0.922 0.930 0.478 0.113 0.035
AGE7 0.922 0.930 0.478 0.113 0.035
AGE8 0.922 0.930 0.478 0.113 0.035
Covariance Coverage
VARA6 VARA7 VARA8 SEX1 SEX2
________ ________ ________ ________ ________
VARA6 0.017
VARA7 0.009 0.009
VARA8 0.009 0.009 0.009
SEX1 0.017 0.009 0.009 1.000
SEX2 0.017 0.009 0.009 1.000 1.000
SEX3 0.017 0.009 0.009 1.000 1.000
SEX4 0.017 0.009 0.009 1.000 1.000
SEX5 0.017 0.009 0.009 1.000 1.000
SEX6 0.017 0.009 0.009 1.000 1.000
SEX7 0.017 0.009 0.009 1.000 1.000
SEX8 0.017 0.009 0.009 1.000 1.000
AGE1 0.017 0.009 0.009 1.000 1.000
AGE2 0.017 0.009 0.009 1.000 1.000
AGE3 0.017 0.009 0.009 1.000 1.000
AGE4 0.017 0.009 0.009 1.000 1.000
AGE5 0.017 0.009 0.009 1.000 1.000
AGE6 0.017 0.009 0.009 1.000 1.000
AGE7 0.017 0.009 0.009 1.000 1.000
AGE8 0.017 0.009 0.009 1.000 1.000
Covariance Coverage
SEX3 SEX4 SEX5 SEX6 SEX7
________ ________ ________ ________ ________
SEX3 1.000
SEX4 1.000 1.000
SEX5 1.000 1.000 1.000
SEX6 1.000 1.000 1.000 1.000
SEX7 1.000 1.000 1.000 1.000 1.000
SEX8 1.000 1.000 1.000 1.000 1.000
AGE1 1.000 1.000 1.000 1.000 1.000
AGE2 1.000 1.000 1.000 1.000 1.000
AGE3 1.000 1.000 1.000 1.000 1.000
AGE4 1.000 1.000 1.000 1.000 1.000
AGE5 1.000 1.000 1.000 1.000 1.000
AGE6 1.000 1.000 1.000 1.000 1.000
AGE7 1.000 1.000 1.000 1.000 1.000
AGE8 1.000 1.000 1.000 1.000 1.000
Covariance Coverage
SEX8 AGE1 AGE2 AGE3 AGE4
________ ________ ________ ________ ________
SEX8 1.000
AGE1 1.000 1.000
AGE2 1.000 1.000 1.000
AGE3 1.000 1.000 1.000 1.000
AGE4 1.000 1.000 1.000 1.000 1.000
AGE5 1.000 1.000 1.000 1.000 1.000
AGE6 1.000 1.000 1.000 1.000 1.000
AGE7 1.000 1.000 1.000 1.000 1.000
AGE8 1.000 1.000 1.000 1.000 1.000
Covariance Coverage
AGE5 AGE6 AGE7 AGE8
________ ________ ________ ________
AGE5 1.000
AGE6 1.000 1.000
AGE7 1.000 1.000 1.000
AGE8 1.000 1.000 1.000 1.000
THE COVARIANCE COVERAGE FALLS BELOW THE SPECIFIED LIMIT.
THE MISSING DATA EM ALGORITHM WILL NOT BE INITIATED.
CHECK YOUR DATA OR LOWER THE COVARIANCE COVERAGE LIMIT.
PROPORTION OF DATA PRESENT FOR DZ
Covariance Coverage
VARA1 VARA2 VARA3 VARA4 VARA5
________ ________ ________ ________ ________
VARA1 0.936
VARA2 0.840 0.891
VARA3 0.519 0.513 0.577
VARA4 0.103 0.109 0.122 0.122
VARA5 0.013 0.019 0.019 0.019 0.019
VARA6 0.006 0.013 0.013 0.013 0.013
VARA7 0.000 0.000 0.000 0.000 0.000
VARA8 0.000 0.000 0.000 0.000 0.000
SEX1 0.936 0.891 0.577 0.122 0.019
SEX2 0.936 0.891 0.577 0.122 0.019
SEX3 0.936 0.891 0.577 0.122 0.019
SEX4 0.936 0.891 0.577 0.122 0.019
SEX5 0.936 0.891 0.577 0.122 0.019
SEX6 0.936 0.891 0.577 0.122 0.019
SEX7 0.936 0.891 0.577 0.122 0.019
SEX8 0.936 0.891 0.577 0.122 0.019
AGE1 0.936 0.891 0.577 0.122 0.019
AGE2 0.936 0.891 0.577 0.122 0.019
AGE3 0.936 0.891 0.577 0.122 0.019
AGE4 0.936 0.891 0.577 0.122 0.019
AGE5 0.936 0.891 0.577 0.122 0.019
AGE6 0.936 0.891 0.577 0.122 0.019
AGE7 0.936 0.891 0.577 0.122 0.019
AGE8 0.936 0.891 0.577 0.122 0.019
Covariance Coverage
VARA6 VARA7 VARA8 SEX1 SEX2
________ ________ ________ ________ ________
VARA6 0.013
VARA7 0.000 0.000
VARA8 0.000 0.000 0.000
SEX1 0.013 0.000 0.000 1.000
SEX2 0.013 0.000 0.000 1.000 1.000
SEX3 0.013 0.000 0.000 1.000 1.000
SEX4 0.013 0.000 0.000 1.000 1.000
SEX5 0.013 0.000 0.000 1.000 1.000
SEX6 0.013 0.000 0.000 1.000 1.000
SEX7 0.013 0.000 0.000 1.000 1.000
SEX8 0.013 0.000 0.000 1.000 1.000
AGE1 0.013 0.000 0.000 1.000 1.000
AGE2 0.013 0.000 0.000 1.000 1.000
AGE3 0.013 0.000 0.000 1.000 1.000
AGE4 0.013 0.000 0.000 1.000 1.000
AGE5 0.013 0.000 0.000 1.000 1.000
AGE6 0.013 0.000 0.000 1.000 1.000
AGE7 0.013 0.000 0.000 1.000 1.000
AGE8 0.013 0.000 0.000 1.000 1.000
Covariance Coverage
SEX3 SEX4 SEX5 SEX6 SEX7
________ ________ ________ ________ ________
SEX3 1.000
SEX4 1.000 1.000
SEX5 1.000 1.000 1.000
SEX6 1.000 1.000 1.000 1.000
SEX7 1.000 1.000 1.000 1.000 1.000
SEX8 1.000 1.000 1.000 1.000 1.000
AGE1 1.000 1.000 1.000 1.000 1.000
AGE2 1.000 1.000 1.000 1.000 1.000
AGE3 1.000 1.000 1.000 1.000 1.000
AGE4 1.000 1.000 1.000 1.000 1.000
AGE5 1.000 1.000 1.000 1.000 1.000
AGE6 1.000 1.000 1.000 1.000 1.000
AGE7 1.000 1.000 1.000 1.000 1.000
AGE8 1.000 1.000 1.000 1.000 1.000
Covariance Coverage
SEX8 AGE1 AGE2 AGE3 AGE4
________ ________ ________ ________ ________
SEX8 1.000
AGE1 1.000 1.000
AGE2 1.000 1.000 1.000
AGE3 1.000 1.000 1.000 1.000
AGE4 1.000 1.000 1.000 1.000 1.000
AGE5 1.000 1.000 1.000 1.000 1.000
AGE6 1.000 1.000 1.000 1.000 1.000
AGE7 1.000 1.000 1.000 1.000 1.000
AGE8 1.000 1.000 1.000 1.000 1.000
Covariance Coverage
AGE5 AGE6 AGE7 AGE8
________ ________ ________ ________
AGE5 1.000
AGE6 1.000 1.000
AGE7 1.000 1.000 1.000
AGE8 1.000 1.000 1.000 1.000
THE MINIMUM COVARIANCE COVERAGE WAS NOT FULFILLED FOR ALL GROUPS.
WARNING: THE SAMPLE COVARIANCE OF THE INDEPENDENT VARIABLES IN GROUP MZ
IS SINGULAR.
WARNING: THE SAMPLE CORRELATION OF SEX8 AND SEX7
IN GROUP MZ IS 1.000.
WARNING: THE SAMPLE CORRELATION OF AGE7 AND SEX7
IN GROUP MZ IS 1.000.
WARNING: THE SAMPLE CORRELATION OF AGE8 AND SEX8
IN GROUP MZ IS 1.000.
WARNING: THE SAMPLE COVARIANCE OF THE INDEPENDENT VARIABLES IN GROUP DZ
IS SINGULAR.
WARNING: THE SAMPLE VARIANCE OF SEX7 IN GROUP DZ IS 0.000.
WARNING: THE SAMPLE VARIANCE OF SEX8 IN GROUP DZ IS 0.000.
WARNING: THE SAMPLE VARIANCE OF AGE7 IN GROUP DZ IS 0.000.
WARNING: THE SAMPLE VARIANCE OF AGE8 IN GROUP DZ IS 0.000.
WARNING: THE SAMPLE COVARIANCE OF THE INDEPENDENT VARIABLES IN GROUP DZ
IS SINGULAR. PROBLEM INVOLVING VARIABLE AGE6.
THE MODEL ESTIMATION TERMINATED NORMALLY
THE CHI-SQUARE COULD NOT BE COMPUTED.
PROBLEMS OCCURRED DURING THE H1 ESTIMATION.
TESTS OF MODEL FIT
Degrees of Freedom 291
Loglikelihood
H0 Value -3198.779
Information Criteria
Number of Free Parameters 6
Akaike (AIC) 6409.557
Bayesian (BIC) 6431.170
Sample-Size Adjusted BIC 6412.146
(n* = (n + 2) / 24)
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
Group MZ
VARA1 ON
SEX1 0.170 1.063 0.160 0.873
AGE1 0.006 0.045 0.127 0.899
VARA2 ON
SEX2 0.170 1.063 0.160 0.873
AGE2 0.006 0.045 0.127 0.899
VARA3 ON
SEX3 0.170 1.063 0.160 0.873
AGE3 0.006 0.045 0.127 0.899
VARA4 ON
SEX4 0.170 1.063 0.160 0.873
AGE4 0.006 0.045 0.127 0.899
VARA5 ON
SEX5 0.170 1.063 0.160 0.873
AGE5 0.006 0.045 0.127 0.899
VARA6 ON
SEX6 0.170 1.063 0.160 0.873
AGE6 0.006 0.045 0.127 0.899
VARA7 ON
SEX7 0.170 1.063 0.160 0.873
AGE7 0.006 0.045 0.127 0.899
VARA8 ON
SEX8 0.170 1.063 0.160 0.873
AGE8 0.006 0.045 0.127 0.899
VARA1 WITH
VARA2 116.553 11.748 9.921 0.000
VARA3 41.470 10.092 4.109 0.000
VARA4 41.470 10.092 4.109 0.000
VARA5 41.470 10.092 4.109 0.000
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA2 WITH
VARA3 41.470 10.092 4.109 0.000
VARA4 41.470 10.092 4.109 0.000
VARA5 41.470 10.092 4.109 0.000
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA3 WITH
VARA4 41.470 10.092 4.109 0.000
VARA5 41.470 10.092 4.109 0.000
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA4 WITH
VARA5 41.470 10.092 4.109 0.000
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA5 WITH
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA6 WITH
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA7 WITH
VARA8 41.470 10.092 4.109 0.000
Intercepts
VARA1 -0.590 1.829 -0.322 0.747
VARA2 -0.590 1.829 -0.322 0.747
VARA3 -0.590 1.829 -0.322 0.747
VARA4 -0.590 1.829 -0.322 0.747
VARA5 -0.590 1.829 -0.322 0.747
VARA6 -0.590 1.829 -0.322 0.747
VARA7 -0.590 1.829 -0.322 0.747
VARA8 -0.590 1.829 -0.322 0.747
Residual Variances
VARA1 168.158 9.993 16.828 0.000
VARA2 168.158 9.993 16.828 0.000
VARA3 168.158 9.993 16.828 0.000
VARA4 168.158 9.993 16.828 0.000
VARA5 168.158 9.993 16.828 0.000
VARA6 168.158 9.993 16.828 0.000
VARA7 168.158 9.993 16.828 0.000
VARA8 168.158 9.993 16.828 0.000
Group DZ
VARA1 ON
SEX1 0.170 1.063 0.160 0.873
AGE1 0.006 0.045 0.127 0.899
VARA2 ON
SEX2 0.170 1.063 0.160 0.873
AGE2 0.006 0.045 0.127 0.899
VARA3 ON
SEX3 0.170 1.063 0.160 0.873
AGE3 0.006 0.045 0.127 0.899
VARA4 ON
SEX4 0.170 1.063 0.160 0.873
AGE4 0.006 0.045 0.127 0.899
VARA5 ON
SEX5 0.170 1.063 0.160 0.873
AGE5 0.006 0.045 0.127 0.899
VARA6 ON
SEX6 0.170 1.063 0.160 0.873
AGE6 0.006 0.045 0.127 0.899
VARA7 ON
SEX7 0.170 1.063 0.160 0.873
AGE7 0.006 0.045 0.127 0.899
VARA8 ON
SEX8 0.170 1.063 0.160 0.873
AGE8 0.006 0.045 0.127 0.899
VARA1 WITH
VARA2 41.470 10.092 4.109 0.000
VARA3 41.470 10.092 4.109 0.000
VARA4 41.470 10.092 4.109 0.000
VARA5 41.470 10.092 4.109 0.000
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA2 WITH
VARA3 41.470 10.092 4.109 0.000
VARA4 41.470 10.092 4.109 0.000
VARA5 41.470 10.092 4.109 0.000
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA3 WITH
VARA4 41.470 10.092 4.109 0.000
VARA5 41.470 10.092 4.109 0.000
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA4 WITH
VARA5 41.470 10.092 4.109 0.000
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA5 WITH
VARA6 41.470 10.092 4.109 0.000
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA6 WITH
VARA7 41.470 10.092 4.109 0.000
VARA8 41.470 10.092 4.109 0.000
VARA7 WITH
VARA8 41.470 10.092 4.109 0.000
Intercepts
VARA1 -0.590 1.829 -0.322 0.747
VARA2 -0.590 1.829 -0.322 0.747
VARA3 -0.590 1.829 -0.322 0.747
VARA4 -0.590 1.829 -0.322 0.747
VARA5 -0.590 1.829 -0.322 0.747
VARA6 -0.590 1.829 -0.322 0.747
VARA7 -0.590 1.829 -0.322 0.747
VARA8 -0.590 1.829 -0.322 0.747
Residual Variances
VARA1 168.158 9.993 16.828 0.000
VARA2 168.158 9.993 16.828 0.000
VARA3 168.158 9.993 16.828 0.000
VARA4 168.158 9.993 16.828 0.000
VARA5 168.158 9.993 16.828 0.000
VARA6 168.158 9.993 16.828 0.000
VARA7 168.158 9.993 16.828 0.000
VARA8 168.158 9.993 16.828 0.000
New/Additional Parameters
A 49.326 37.136 1.328 0.184
D 67.227 37.481 1.794 0.073
E 51.605 7.151 7.216 0.000
X 7.023 2.644 2.656 0.008
W 8.199 2.286 3.587 0.000
Z 7.184 0.498 14.432 0.000
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.410E-04
(ratio of smallest to largest eigenvalue)
Beginning Time: 22:57:49
Ending Time: 22:57:50
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