```Mplus VERSION 6
MUTHEN & MUTHEN
04/25/2010  10:57 PM

INPUT INSTRUCTIONS

! SCRIPT NAME        : ordSATut6  (cvb)
! GOAL                : calculation of tetrachoric correlations, check assumptions of twin
! DATA                : ordinal
! INPUT                : raw data
! UNI/BI/MULTI        : uni
! DATA-GROUPS        : MZM, DZM, MZF, DZF, DOSMF, DOSFM
! MEANS MODEL        : n.a.
! VARIANCE COVARIANCE MODEL(S)        : n.a.
!
!
! Evaluated models:
! 1. no model
! 2. as 1, plus thresholds same for twin 1 and twin 2 (but different for DOS twins)
! 3. as 2, plus thresholds same for DZ-ss and DZ-os of the same sex
! 4. as 3, plus thresholds same for MZ and DZ
! 5. as 4, plus thresholds same for males and females but correlations different for 6 gro
! 6. as 5, plus correlations of DOS-MF = DOS-FM
! 7. as 6, plus correlations MZM = MZF and DZM = DZF = DOS-MF = DOS-FM (no sex differences
! 8. as 4, plus correlations of DOS-MF = DOS-FM, but thresholds differ for males and femal
! 9. as 8, plus correlations MZM = MZF and DZM = DZF = DOS-MF = DOS-FM but thresholds diff

data: file is ordraw1.dat;

variable: names are id y1 y2 zygot age;
categorical=y1 y2;
usevar are y1 y2;
grouping=zygot(1=MZM 2=DZM 3=MZF 4=DZF 5=DOSMF 6=DOSFM);  ! specify the groups
missing=all(-9); ! specify missing data symbol

model:
[y1\$1] (mzmt1);
[y2\$1] (mzmt2);
y1 with y2 (mzmc);

model dzm:
[y1\$1] (dzmt1);
[y2\$1] (dzmt2);
y1 with y2 (dzmc);

model mzf:
[y1\$1] (mzft1);
[y2\$1] (mzft2);
y1 with y2 (mzfc);

model dzf:
[y1\$1] (dzft1);
[y2\$1] (dzft2);
y1 with y2 (dzfc);

model dosmf:
[y1\$1] (dosmft1);
[y2\$1] (dosmft2);
y1 with y2 (dosmfc);

model dosfm:
[y1\$1] (dosfmt1);
[y2\$1] (dosfmt2);
y1 with y2 (dosfmc);

! model test will produce Wald's test for these constraints
model test:

! Uncomment to test equal thresholds for for twin 1 and twin 2
mzmt1 = mzmt2;
dzmt1 = dzmt2;
mzft1 = mzft2;
dzft1 = dzft2;
dosfmt1 = dosmft2;
dosfmt2 = dosmft1;

! Uncomment to test equal thresholds for same for DZ-ss and DZ-os of the same sex
! dzmt1 = dosmft1;
! dzft1 = dosfmt1;

! Uncomment to test equal thresholds for MZ and DZ
! mzmt1 = dzmt1;
! mzft1 = dzft1;

! Uncomment to test equal thresholds across gender
! mzmt1 = mzft1;

! Uncomment to test correlations of DOS-MF = DOS-FM
! dosfmc = dosmfc;

! Uncomment to test no sex differences on correlation
! correlations MZM = MZF and DZM = DZF = DOS-MF = DOS-FM (no sex differences)
! mzmc = mzfc;
! dzmc = dzfc;
! dzmc = dosfmc;

*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! 5. as 4, plus thresholds same for males and females but correlations different for 6 grou
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! 7. as 6, plus correlations MZM = MZF and DZM = DZF = DOS-MF = DOS-FM (no sex differences)
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! 8. as 4, plus correlations of DOS-MF = DOS-FM, but thresholds differ for males and female
*** WARNING
Input line exceeded 90 characters. Some input may be truncated.
! 9. as 8, plus correlations MZM = MZF and DZM = DZF = DOS-MF = DOS-FM but thresholds diffe
4 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS

SUMMARY OF ANALYSIS

Number of groups                                                 6
Number of observations
Group MZM                                                   399
Group DZM                                                   273
Group MZF                                                   891
Group DZF                                                   577
Group DOSMF                                                 381
Group DOSFM                                                 330

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     ZYGOT

Estimator                                                    WLSMV
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
Parameterization                                             DELTA

Input data file(s)
ordraw1.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
Y1            Y2
________      ________
Y1             0.827
Y2             0.609         0.782

PROPORTION OF DATA PRESENT FOR DZM

Covariance Coverage
Y1            Y2
________      ________
Y1             0.744
Y2             0.498         0.755

PROPORTION OF DATA PRESENT FOR MZF

Covariance Coverage
Y1            Y2
________      ________
Y1             0.856
Y2             0.694         0.837

PROPORTION OF DATA PRESENT FOR DZF

Covariance Coverage
Y1            Y2
________      ________
Y1             0.792
Y2             0.549         0.757

PROPORTION OF DATA PRESENT FOR DOSMF

Covariance Coverage
Y1            Y2
________      ________
Y1             0.633
Y2             0.454         0.822

PROPORTION OF DATA PRESENT FOR DOSFM

Covariance Coverage
Y1            Y2
________      ________
Y1             0.845
Y2             0.424         0.579

UNIVARIATE PROPORTIONS AND COUNTS FOR CATEGORICAL VARIABLES

Group MZM
Y1
Category 1    0.952      314.000
Category 2    0.048       16.000
Y2
Category 1    0.933      291.000
Category 2    0.067       21.000

Group DZM
Y1
Category 1    0.946      192.000
Category 2    0.054       11.000
Y2
Category 1    0.937      193.000
Category 2    0.063       13.000

Group MZF
Y1
Category 1    0.814      621.000
Category 2    0.186      142.000
Y2
Category 1    0.822      613.000
Category 2    0.178      133.000

Group DZF
Y1
Category 1    0.838      383.000
Category 2    0.162       74.000
Y2
Category 1    0.812      355.000
Category 2    0.188       82.000

Group DOSMF
Y1
Category 1    0.909      219.000
Category 2    0.091       22.000
Y2
Category 1    0.812      254.000
Category 2    0.188       59.000

Group DOSFM
Y1
Category 1    0.857      239.000
Category 2    0.143       40.000
Y2
Category 1    0.916      175.000
Category 2    0.084       16.000

WARNING:  THE BIVARIATE TABLE OF Y2 AND Y1 HAS AN EMPTY CELL.

THE MODEL ESTIMATION TERMINATED NORMALLY

TESTS OF MODEL FIT

Chi-Square Test of Model Fit

Value                              0.000*
Degrees of Freedom                     0
P-Value                           0.0000

Chi-Square Contributions From Each Group

MZM                                0.000
DZM                                0.000
MZF                                0.000
DZF                                0.000
DOSMF                              0.000
DOSFM                              0.000

*   The chi-square value for MLM, MLMV, MLR, ULSMV, WLSM and WLSMV cannot be used
for chi-square difference testing in the regular way.  MLM, MLR and WLSM
chi-square difference testing is described on the Mplus website.  MLMV, WLSMV,
and ULSMV difference testing is done using the DIFFTEST option.

Chi-Square Test of Model Fit for the Baseline Model

Value                            176.182
Degrees of Freedom                     6
P-Value                           0.0000

CFI/TLI

CFI                                1.000
TLI                                1.000

Wald Test of Parameter Constraints

Value                              5.084
Degrees of Freedom                     6
P-Value                           0.5331

Number of Free Parameters                       18

RMSEA (Root Mean Square Error Of Approximation)

Estimate                           0.000

WRMR (Weighted Root Mean Square Residual)

Value                              0.019

MODEL RESULTS

Two-Tailed
Estimate       S.E.  Est./S.E.    P-Value

Group MZM

Y1       WITH
Y2                 0.650      0.137      4.733      0.000

Thresholds
Y1\$1               1.660      0.117     14.126      0.000
Y2\$1               1.496      0.109     13.740      0.000

Group DZM

Y1       WITH
Y2                 0.007      0.581      0.012      0.991

Thresholds
Y1\$1               1.606      0.145     11.109      0.000
Y2\$1               1.529      0.137     11.185      0.000

Group MZF

Y1       WITH
Y2                 0.649      0.055     11.802      0.000

Thresholds
Y1\$1               0.892      0.053     16.968      0.000
Y2\$1               0.922      0.054     17.159      0.000

Group DZF

Y1       WITH
Y2                 0.341      0.113      3.018      0.003

Thresholds
Y1\$1               0.987      0.070     14.039      0.000
Y2\$1               0.887      0.069     12.783      0.000

Group DOSMF

Y1       WITH
Y2                 0.318      0.166      1.907      0.056

Thresholds
Y1\$1               1.333      0.113     11.790      0.000
Y2\$1               0.883      0.082     10.792      0.000

Group DOSFM

Y1       WITH
Y2                 0.279      0.212      1.317      0.188

Thresholds
Y1\$1               1.065      0.093     11.485      0.000
Y2\$1               1.380      0.130     10.597      0.000

QUALITY OF NUMERICAL RESULTS

Condition Number for the Information Matrix              0.820E-02
(ratio of smallest to largest eigenvalue)

Beginning Time:  22:57:39
Ending Time:  22:57:39
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
```