Mplus VERSION 6
MUTHEN & MUTHEN
04/25/2010 10:57 PM
INPUT INSTRUCTIONS
! SCRIPT NAME : rawVCQ1 (dp)
! GOAL : To evaluate significance of linkage using pi-hat method
! 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. AEQ
! 2. AE
data: file is example3.dat;
variable:
names are
famnr zygos
fata1 fata2 mota1 mota2 !fatherallele1 fatherallele2 motherallele1 mothereallele
pheno1 age1 sex1 tw1a1 tw1a2 !fenotypetwin1 agetwin1 sextwin1 twin1allele1 twin1
pheno2 age2 sex2 tw2a1 tw2a2 !fenotypetwin2 agetwin2 sextwin2 twin2allele1 twin2
z0_0 z1_0 z2_0 z0_1 z1_1 z2_1 z0_2 z1_2 z2_2 z0_3 z1_3 z2_3 z0_4 z1_4 z2_4
z0_5 z1_5 z2_5 z0_6 z1_6 z2_6 z0_7 z1_7 z2_7 z0_8 z1_8 z2_8 z0_9 z1_9 z2_9
z0_10 z1_10 z2_10 z0_11 z1_11 z2_11 z0_12 z1_12 z2_12 z0_13 z1_13 z2_13
z0_14 z1_14 z2_14 z0_15 z1_15 z2_15 z0_16 z1_16 z2_16 z0_17 z1_17 z2_17
z0_18 z1_18 z2_18 z0_19 z1_19 z2_19 z0_20 z1_20 z2_20 z0_21 z1_21 z2_21
z0_22 z1_22 z2_22 z0_23 z1_23 z2_23 z0_24 z1_24 z2_24 z0_25 z1_25 z2_25
z0_26 z1_26 z2_26 z0_27 z1_27 z2_27 z0_28 z1_28 z2_28 z0_29 z1_29 z2_29
z0_30 z1_30 z2_30 z0_31 z1_31 z2_31 z0_32 z1_32 z2_32 z0_33 z1_33 z2_33
z0_34 z1_34 z2_34 z0_35 z1_35 z2_35;
usevar are pheno1 pheno2 age1 age2 sex1 sex2;
constraint=pihat h;
missing =all(-99.00);
define: if (zygos>=3) then pihat=0.5*z1_20+z2_20 else pihat=1;
if (zygos>=3) then h=0.5 else h=1;
model:
pheno1 on age1 (b1);
pheno2 on age2 (b1);
pheno1 on sex1 (b2);
pheno2 on sex2 (b2);
[pheno1 pheno2] (m);
pheno1 pheno2 (v);
pheno1 with pheno2 (c);
model constraint:
! Uncomment for AEQ model
new(a e q x z d);
a=x*x;
e=z*z;
q=d*d;
v=x*x+z*z+d*d;
c=h*x*x+pihat*d*d;
! Uncomment for AE model
! new(a e x z);
! a=x*x;
! e=z*z;
! v=x*x+z*z;
! c=h*x*x;
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
fata1 fata2 mota1 mota2 !fatherallele1 fatherallele2 motherallele1 mothereallele2, allel
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
pheno1 age1 sex1 tw1a1 tw1a2 !fenotypetwin1 agetwin1 sextwin1 twin1allele1 twin1allele2,
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
pheno2 age2 sex2 tw2a1 tw2a2 !fenotypetwin2 agetwin2 sextwin2 twin2allele1 twin2allele2,
3 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS
SUMMARY OF ANALYSIS
Number of groups 1
Number of observations 152
Number of dependent variables 2
Number of independent variables 4
Number of continuous latent variables 0
Observed dependent variables
Continuous
PHENO1 PHENO2
Observed independent variables
AGE1 AGE2 SEX1 SEX2
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)
example3.dat
Input data format FREE
SUMMARY OF DATA
Number of missing data patterns 2
COVARIANCE COVERAGE OF DATA
Minimum covariance coverage value 0.100
PROPORTION OF DATA PRESENT
Covariance Coverage
PHENO1 PHENO2 AGE1 AGE2 SEX1
________ ________ ________ ________ ________
PHENO1 0.967
PHENO2 0.967 0.967
AGE1 0.967 0.967 1.000
AGE2 0.967 0.967 1.000 1.000
SEX1 0.967 0.967 1.000 1.000 1.000
SEX2 0.967 0.967 1.000 1.000 1.000
Covariance Coverage
SEX2
________
SEX2 1.000
THE MODEL ESTIMATION TERMINATED NORMALLY
TESTS OF MODEL FIT
Loglikelihood
H0 Value -941.506
Information Criteria
Number of Free Parameters 6
Akaike (AIC) 1895.012
Bayesian (BIC) 1913.155
Sample-Size Adjusted BIC 1894.165
(n* = (n + 2) / 24)
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
PHENO1 ON
AGE1 -0.234 0.081 -2.878 0.004
SEX1 1.475 0.295 5.003 0.000
PHENO2 ON
AGE2 -0.234 0.081 -2.878 0.004
SEX2 1.475 0.295 5.003 0.000
PHENO1 WITH
PHENO2 999.000 0.000 999.000 999.000
Intercepts
PHENO1 9.956 1.390 7.162 0.000
PHENO2 9.956 1.390 7.162 0.000
Residual Variances
PHENO1 5.166 0.491 10.527 0.000
PHENO2 5.166 0.491 10.527 0.000
New/Additional Parameters
A 4.217 1.496 2.818 0.005
E 0.712 0.129 5.501 0.000
Q 0.237 1.424 0.166 0.868
X 2.054 0.364 5.637 0.000
Z 0.844 0.077 11.003 0.000
D -0.487 1.463 -0.333 0.739
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.398E-04
(ratio of smallest to largest eigenvalue)
Beginning Time: 22:57:52
Ending Time: 22:57:53
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