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

INPUT INSTRUCTIONS

  ! SCRIPT NAME        : ctSATut6 (cb)
  ! GOAL                : To calculate tetrachoric correlations
  ! DATA                : ordinal
  ! INPUT                : contingency tables
  ! UNI/BI/MULTI        : uni
  ! DATA-GROUPS        : MZM DZM MZF DZF DOSMF DOSFM
  ! MEANS MODEL        : -
  ! VARIANCE COVARIANCE MODEL(S)        :
  ! 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 ct6.dat;

  variable: names are g y1 y2 weight;
            categorical=y1 y2;
            grouping=g(1=MZM 2=DZM 3=MZF 4=DZF 5=DOSMF 6=DOSFM);  ! specify the groups
            freqweight=weight;

  analysis: conv=0.000001;

  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                                                   243
   Group DZM                                                   137
   Group MZF                                                   620
   Group DZF                                                   317
   Group DOSMF                                                 173
   Group DOSFM                                                 140
Number of patterns
   Group MZM                                                     4
   Group DZM                                                     4
   Group MZF                                                     4
   Group DZF                                                     4
   Group DOSMF                                                   4
   Group DOSFM                                                   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.100D-05
Maximum number of steepest descent iterations                   20
Parameterization                                             DELTA

Input data file(s)
  ct6.dat

Input data format  FREE


UNIVARIATE PROPORTIONS AND COUNTS FOR CATEGORICAL VARIABLES

  Group MZM
    Y1
      Category 1    0.951      231.000
      Category 2    0.049       12.000
    Y2
      Category 1    0.934      227.000
      Category 2    0.066       16.000

  Group DZM
    Y1
      Category 1    0.942      129.000
      Category 2    0.058        8.000
    Y2
      Category 1    0.927      127.000
      Category 2    0.073       10.000

  Group MZF
    Y1
      Category 1    0.821      509.000
      Category 2    0.179      111.000
    Y2
      Category 1    0.832      516.000
      Category 2    0.168      104.000

  Group DZF
    Y1
      Category 1    0.864      274.000
      Category 2    0.136       43.000
    Y2
      Category 1    0.817      259.000
      Category 2    0.183       58.000

  Group DOSMF
    Y1
      Category 1    0.908      157.000
      Category 2    0.092       16.000
    Y2
      Category 1    0.809      140.000
      Category 2    0.191       33.000

  Group DOSFM
    Y1
      Category 1    0.879      123.000
      Category 2    0.121       17.000
    Y2
      Category 1    0.914      128.000
      Category 2    0.086       12.000



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                            173.209
          Degrees of Freedom                     6
          P-Value                           0.0000

CFI/TLI

          CFI                                1.000
          TLI                                1.000

Wald Test of Parameter Constraints

          Value                              7.572
          Degrees of Freedom                     6
          P-Value                           0.2712

Number of Free Parameters                       18

RMSEA (Root Mean Square Error Of Approximation)

          Estimate                           0.000

WRMR (Weighted Root Mean Square Residual)

          Value                              0.000



MODEL RESULTS

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

Group MZM

 Y1       WITH
    Y2                 0.651      0.138      4.712      0.000

 Thresholds
    Y1$1               1.651      0.136     12.129      0.000
    Y2$1               1.507      0.124     12.135      0.000

Group DZM

 Y1       WITH
    Y2                 0.158      0.291      0.543      0.587

 Thresholds
    Y1$1               1.568      0.172      9.130      0.000
    Y2$1               1.454      0.160      9.070      0.000

Group MZF

 Y1       WITH
    Y2                 0.646      0.055     11.682      0.000

 Thresholds
    Y1$1               0.919      0.059     15.610      0.000
    Y2$1               0.963      0.060     16.103      0.000

Group DZF

 Y1       WITH
    Y2                 0.333      0.112      2.970      0.003

 Thresholds
    Y1$1               1.100      0.088     12.460      0.000
    Y2$1               0.904      0.082     11.037      0.000

Group DOSMF

 Y1       WITH
    Y2                 0.319      0.167      1.909      0.056

 Thresholds
    Y1$1               1.326      0.133      9.973      0.000
    Y2$1               0.875      0.110      7.969      0.000

Group DOSFM

 Y1       WITH
    Y2                 0.283      0.211      1.338      0.181

 Thresholds
    Y1$1               1.168      0.137      8.534      0.000
    Y2$1               1.368      0.151      9.052      0.000


QUALITY OF NUMERICAL RESULTS

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


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

Back to the list of genetics examples