Mplus VERSION 5.2
MUTHEN & MUTHEN
12/02/2008 2:01 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;
INPUT READING TERMINATED NORMALLY
SUMMARY OF ANALYSIS
Number of groups 2
Number of observations
Group MZ 144
Group DZ 178
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 14
Group DZ
Number of missing data patterns 14
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.736
VARA2 0.708 0.743
VARA3 0.333 0.340 0.382
VARA4 0.063 0.069 0.083 0.090
VARA5 0.021 0.028 0.028 0.028 0.028
VARA6 0.007 0.014 0.014 0.014 0.014
VARA7 0.007 0.007 0.007 0.007 0.007
VARA8 0.007 0.007 0.007 0.007 0.007
SEX1 0.736 0.743 0.382 0.090 0.028
SEX2 0.736 0.743 0.382 0.090 0.028
SEX3 0.736 0.743 0.382 0.090 0.028
SEX4 0.736 0.743 0.382 0.090 0.028
SEX5 0.736 0.743 0.382 0.090 0.028
SEX6 0.736 0.743 0.382 0.090 0.028
SEX7 0.736 0.743 0.382 0.090 0.028
SEX8 0.736 0.743 0.382 0.090 0.028
AGE1 0.736 0.743 0.382 0.090 0.028
AGE2 0.736 0.743 0.382 0.090 0.028
AGE3 0.736 0.743 0.382 0.090 0.028
AGE4 0.736 0.743 0.382 0.090 0.028
AGE5 0.736 0.743 0.382 0.090 0.028
AGE6 0.736 0.743 0.382 0.090 0.028
AGE7 0.736 0.743 0.382 0.090 0.028
AGE8 0.736 0.743 0.382 0.090 0.028
Covariance Coverage
VARA6 VARA7 VARA8 SEX1 SEX2
________ ________ ________ ________ ________
VARA6 0.014
VARA7 0.007 0.007
VARA8 0.007 0.007 0.007
SEX1 0.014 0.007 0.007 1.000
SEX2 0.014 0.007 0.007 1.000 1.000
SEX3 0.014 0.007 0.007 1.000 1.000
SEX4 0.014 0.007 0.007 1.000 1.000
SEX5 0.014 0.007 0.007 1.000 1.000
SEX6 0.014 0.007 0.007 1.000 1.000
SEX7 0.014 0.007 0.007 1.000 1.000
SEX8 0.014 0.007 0.007 1.000 1.000
AGE1 0.014 0.007 0.007 1.000 1.000
AGE2 0.014 0.007 0.007 1.000 1.000
AGE3 0.014 0.007 0.007 1.000 1.000
AGE4 0.014 0.007 0.007 1.000 1.000
AGE5 0.014 0.007 0.007 1.000 1.000
AGE6 0.014 0.007 0.007 1.000 1.000
AGE7 0.014 0.007 0.007 1.000 1.000
AGE8 0.014 0.007 0.007 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.820
VARA2 0.736 0.781
VARA3 0.455 0.449 0.506
VARA4 0.090 0.096 0.107 0.107
VARA5 0.011 0.017 0.017 0.017 0.017
VARA6 0.006 0.011 0.011 0.011 0.011
VARA7 0.000 0.000 0.000 0.000 0.000
VARA8 0.000 0.000 0.000 0.000 0.000
SEX1 0.820 0.781 0.506 0.107 0.017
SEX2 0.820 0.781 0.506 0.107 0.017
SEX3 0.820 0.781 0.506 0.107 0.017
SEX4 0.820 0.781 0.506 0.107 0.017
SEX5 0.820 0.781 0.506 0.107 0.017
SEX6 0.820 0.781 0.506 0.107 0.017
SEX7 0.820 0.781 0.506 0.107 0.017
SEX8 0.820 0.781 0.506 0.107 0.017
AGE1 0.820 0.781 0.506 0.107 0.017
AGE2 0.820 0.781 0.506 0.107 0.017
AGE3 0.820 0.781 0.506 0.107 0.017
AGE4 0.820 0.781 0.506 0.107 0.017
AGE5 0.820 0.781 0.506 0.107 0.017
AGE6 0.820 0.781 0.506 0.107 0.017
AGE7 0.820 0.781 0.506 0.107 0.017
AGE8 0.820 0.781 0.506 0.107 0.017
Covariance Coverage
VARA6 VARA7 VARA8 SEX1 SEX2
________ ________ ________ ________ ________
VARA6 0.011
VARA7 0.000 0.000
VARA8 0.000 0.000 0.000
SEX1 0.011 0.000 0.000 1.000
SEX2 0.011 0.000 0.000 1.000 1.000
SEX3 0.011 0.000 0.000 1.000 1.000
SEX4 0.011 0.000 0.000 1.000 1.000
SEX5 0.011 0.000 0.000 1.000 1.000
SEX6 0.011 0.000 0.000 1.000 1.000
SEX7 0.011 0.000 0.000 1.000 1.000
SEX8 0.011 0.000 0.000 1.000 1.000
AGE1 0.011 0.000 0.000 1.000 1.000
AGE2 0.011 0.000 0.000 1.000 1.000
AGE3 0.011 0.000 0.000 1.000 1.000
AGE4 0.011 0.000 0.000 1.000 1.000
AGE5 0.011 0.000 0.000 1.000 1.000
AGE6 0.011 0.000 0.000 1.000 1.000
AGE7 0.011 0.000 0.000 1.000 1.000
AGE8 0.011 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 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 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
TESTS OF MODEL FIT
Degrees of Freedom 291
Loglikelihood
H0 Value -3039.514
Information Criteria
Number of Free Parameters 6
Akaike (AIC) 6091.028
Bayesian (BIC) 6113.675
Sample-Size Adjusted BIC 6094.644
(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.552 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.833 -0.322 0.748
VARA2 -0.590 1.833 -0.322 0.748
VARA3 -0.590 1.833 -0.322 0.748
VARA4 -0.590 1.833 -0.322 0.748
VARA5 -0.590 1.833 -0.322 0.748
VARA6 -0.590 1.833 -0.322 0.748
VARA7 -0.590 1.833 -0.322 0.748
VARA8 -0.590 1.833 -0.322 0.748
Residual Variances
VARA1 168.158 9.992 16.828 0.000
VARA2 168.158 9.992 16.828 0.000
VARA3 168.158 9.992 16.828 0.000
VARA4 168.158 9.992 16.828 0.000
VARA5 168.158 9.992 16.828 0.000
VARA6 168.158 9.992 16.828 0.000
VARA7 168.158 9.992 16.828 0.000
VARA8 168.158 9.992 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.833 -0.322 0.748
VARA2 -0.590 1.833 -0.322 0.748
VARA3 -0.590 1.833 -0.322 0.748
VARA4 -0.590 1.833 -0.322 0.748
VARA5 -0.590 1.833 -0.322 0.748
VARA6 -0.590 1.833 -0.322 0.748
VARA7 -0.590 1.833 -0.322 0.748
VARA8 -0.590 1.833 -0.322 0.748
Residual Variances
VARA1 168.158 9.992 16.828 0.000
VARA2 168.158 9.992 16.828 0.000
VARA3 168.158 9.992 16.828 0.000
VARA4 168.158 9.992 16.828 0.000
VARA5 168.158 9.992 16.828 0.000
VARA6 168.158 9.992 16.828 0.000
VARA7 168.158 9.992 16.828 0.000
VARA8 168.158 9.992 16.828 0.000
New/Additional Parameters
A 49.328 37.138 1.328 0.184
D 67.224 37.482 1.793 0.073
E 51.606 7.152 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.408E-04
(ratio of smallest to largest eigenvalue)
Beginning Time: 14:01:35
Ending Time: 14:01:35
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-2008 Muthen & Muthen
Back to the list of genetics examples