Mplus VERSION 5.2
MUTHEN & MUTHEN
12/02/2008 2:01 PM
INPUT INSTRUCTIONS
! SCRIPT NAME : ctVCut4c (cb)
! GOAL : To evaluate best model for variance components A C E
! DATA : ordinal
! INPUT : contingency tables
! UNI/BI/MULTI : uni
! DATA-GROUPS : MZM DZM MZF DZF
! MEANS MODEL : -
! VARIANCE COVARIANCE MODEL(S) :
! 1. different prevalences for males and females, ACE for males, KLM for females
! 2. same prevalences for males and females, ACE for males, KLM for females
! 3. same prevalences for males and females, ACE= KLM for males and females
! 4. same prevalences for males and females, AE= KM for males and females
! 5. same prevalences for males and females, CE= LM for males and females
! 6. same prevalences for males and females, E= M for males and females
! 7. different prevalences for males and females, ACE= KLM for males and females
! 8. different prevalences for males and females, AE= KM for males and females
! 9. different prevalences for males and females, CE= LM for males and females
! 10. different prevalences for males and females, E= M for males and females
data: file is ct4.dat;
variable: names are g y1 y2 weight;
categorical=y1 y2;
grouping=g(1=MZM 2=DZM 3=MZF 4=DZF); ! specify the groups
freqweight=weight;
model:
[y1$1] (mt);
[y2$1] (mt);
y1 with y2 (mzmc);
model dzm:
[y1$1] (mt);
[y2$1] (mt);
y1 with y2 (dzmc);
model mzf:
[y1$1] (ft);
[y2$1] (ft);
y1 with y2 (mzfc);
model dzf:
[y1$1] (ft);
[y2$1] (ft);
y1 with y2 (dzfc);
model constraint:
new(a c e x y z);
a=x*x;
c=y*y;
e=1-x*x-y*y;
z=sqrt(1-x*x-y*y);
mzmc=x*x+y*y;
dzmc=0.5*x*x+y*y;
new(k l m s t u);
k=s*s;
l=t*t;
m=1-s*s-t*t;
u=sqrt(1-s*s-t*t);
mzfc=s*s+t*t;
dzfc=0.5*s*s+t*t;
! Uncomment for same prevalences for males and females
! mt=ft;
! Uncomment for Model ACE=KLM
! a=k;
! c=l;
! Uncomment for Model AE=KM
! c=0;
! Uncomment for Model CE=LM
! a=0;
! Uncomment for Model E=M
! a=0;
! c=0;
INPUT READING TERMINATED NORMALLY
SUMMARY OF ANALYSIS
Number of groups 4
Number of observations
Group MZM 243
Group DZM 137
Group MZF 620
Group DZF 317
Number of patterns
Group MZM 4
Group DZM 4
Group MZF 4
Group DZF 4
Number of dependent variables 2
Number of independent variables 0
Number of continuous latent variables 0
Observed dependent variables
Binary and ordered categorical (ordinal)
Y1 Y2
Variables with special functions
Grouping variable G
Weight variable WEIGHT
Estimator WLSMV
Maximum number of iterations 1000
Convergence criterion 0.500D-04
Maximum number of steepest descent iterations 20
Parameterization DELTA
Input data file(s)
ct4.dat
Input data format FREE
SUMMARY OF CATEGORICAL DATA PROPORTIONS
Group MZM
Y1
Category 1 0.951
Category 2 0.049
Y2
Category 1 0.934
Category 2 0.066
Group DZM
Y1
Category 1 0.942
Category 2 0.058
Y2
Category 1 0.927
Category 2 0.073
Group MZF
Y1
Category 1 0.821
Category 2 0.179
Y2
Category 1 0.832
Category 2 0.168
Group DZF
Y1
Category 1 0.864
Category 2 0.136
Y2
Category 1 0.817
Category 2 0.183
THE MODEL ESTIMATION TERMINATED NORMALLY
TESTS OF MODEL FIT
Chi-Square Test of Model Fit
Value 5.294*
Degrees of Freedom 6**
P-Value 0.5067
Chi-Square Contributions From Each Group
MZM 0.767
DZM 0.716
MZF 0.498
DZF 3.313
* The chi-square value for MLM, MLMV, MLR, ULSMV, WLSM and WLSMV cannot be used
for chi-square difference tests. MLM, MLR and WLSM chi-square difference
testing is described in the Mplus Technical Appendices at www.statmodel.com.
See chi-square difference testing in the index of the Mplus User's Guide.
** The degrees of freedom for MLMV, ULSMV and WLSMV are estimated according to
a formula given in the Mplus Technical Appendices at www.statmodel.com.
See degrees of freedom in the index of the Mplus User's Guide.
Chi-Square Test of Model Fit for the Baseline Model
Value 167.776
Degrees of Freedom 4
P-Value 0.0000
CFI/TLI
CFI 1.000
TLI 1.003
Number of Free Parameters 6
RMSEA (Root Mean Square Error Of Approximation)
Estimate 0.000
WRMR (Weighted Root Mean Square Residual)
Value 1.267
MODEL RESULTS
Two-Tailed
Estimate S.E. Est./S.E. P-Value
Group MZM
Y1 WITH
Y2 0.633 0.138 4.582 0.000
Thresholds
Y1$1 1.547 0.080 19.409 0.000
Y2$1 1.547 0.080 19.409 0.000
Group DZM
Y1 WITH
Y2 0.316 0.291 1.088 0.277
Thresholds
Y1$1 1.547 0.080 19.409 0.000
Y2$1 1.547 0.080 19.409 0.000
Group MZF
Y1 WITH
Y2 0.646 0.055 11.680 0.000
Thresholds
Y1$1 0.958 0.040 24.225 0.000
Y2$1 0.958 0.040 24.225 0.000
Group DZF
Y1 WITH
Y2 0.333 0.112 2.976 0.003
Thresholds
Y1$1 0.958 0.040 24.225 0.000
Y2$1 0.958 0.040 24.225 0.000
New/Additional Parameters
A 0.633 0.644 0.982 0.326
C 0.000 0.598 0.000 1.000
E 0.367 0.138 2.661 0.008
X 0.795 0.405 1.965 0.049
Y -0.003 86.876 0.000 1.000
Z 0.606 0.114 5.322 0.000
K 0.625 0.250 2.502 0.012
L 0.021 0.231 0.091 0.928
M 0.354 0.055 6.404 0.000
S 0.791 0.158 5.003 0.000
T -0.145 0.798 -0.181 0.856
U 0.595 0.046 12.809 0.000
QUALITY OF NUMERICAL RESULTS
Condition Number for the Information Matrix 0.101E-06
(ratio of smallest to largest eigenvalue)
Beginning Time: 14:01:23
Ending Time: 14:01:24
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-2008 Muthen & Muthen
Back to the list of genetics examples