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

INPUT INSTRUCTIONS

  ! SCRIPT NAME        : ctVCut2c  (cb)
  ! GOAL                : To evaluate best model for variance components A C E
  ! DATA                : ordinal
  ! INPUT                : contingency tables
  ! UNI/BI/MULTI        : uni
  ! DATA-GROUPS        : MZ DZ
  ! MEANS MODEL        : -
  ! VARIANCE COVARIANCE MODEL(S)        :
  ! 1. ACE
  ! 2. AE
  ! 3. CE
  ! 4. E

  data: file is ct.dat;

  variable: names are g y1 y2 weight;
            categorical=y1 y2;
            grouping=g(1=MZ 2=DZ);  ! specify the two groups MZ and DZ
            freqweight=weight;

  model:
       [y1$1 y2$1] (t);
       y1 with y2 (mzc);

  model dz:
       y1 with y2 (dzc);

  model constraint:

  ! Uncomment for Model ACE
    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);
    mzc=x*x+y*y;
    dzc=0.5*x*x+y*y;

  ! Uncomment for Model AE
  ! c=0;

  ! Uncomment for Model CE
  ! a=0;

  ! Uncomment for Model E
  ! a=0; c=0;



INPUT READING TERMINATED NORMALLY




SUMMARY OF ANALYSIS

Number of groups                                                 2
Number of observations
   Group MZ                                                    702
   Group DZ                                                    726
Number of patterns
   Group MZ                                                      4
   Group DZ                                                      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)
  ct.dat

Input data format  FREE


UNIVARIATE PROPORTIONS AND COUNTS FOR CATEGORICAL VARIABLES

  Group MZ
    Y1
      Category 1    0.744      522.000
      Category 2    0.256      180.000
    Y2
      Category 1    0.761      534.000
      Category 2    0.239      168.000

  Group DZ
    Y1
      Category 1    0.753      547.000
      Category 2    0.247      179.000
    Y2
      Category 1    0.773      561.000
      Category 2    0.227      165.000



THE MODEL ESTIMATION TERMINATED NORMALLY



TESTS OF MODEL FIT

Chi-Square Test of Model Fit

          Value                              2.185*
          Degrees of Freedom                     3
          P-Value                           0.5349

Chi-Square Contributions From Each Group

          MZ                                 0.885
          DZ                                 1.300

*   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                             56.109
          Degrees of Freedom                     2
          P-Value                           0.0000

CFI/TLI

          CFI                                1.000
          TLI                                1.010

Number of Free Parameters                        3

RMSEA (Root Mean Square Error Of Approximation)

          Estimate                           0.000

WRMR (Weighted Root Mean Square Residual)

          Value                              0.826



MODEL RESULTS

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

Group MZ

 Y1       WITH
    Y2                 0.408      0.060      6.843      0.000

 Thresholds
    Y1$1               0.699      0.028     25.145      0.000
    Y2$1               0.699      0.028     25.145      0.000

Group DZ

 Y1       WITH
    Y2                 0.204      0.068      3.001      0.003

 Thresholds
    Y1$1               0.699      0.028     25.145      0.000
    Y2$1               0.699      0.028     25.145      0.000

 New/Additional Parameters
    A                  0.408      0.181      2.256      0.024
    C                  0.000      0.148      0.000      1.000
    E                  0.592      0.060      9.945      0.000
    X                  0.638      0.141      4.512      0.000
    Y                  0.002     42.196      0.000      1.000
    Z                  0.770      0.039     19.890      0.000


QUALITY OF NUMERICAL RESULTS

     Condition Number for the Information Matrix              0.370E-06
       (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