Mplus VERSION 7
MUTHEN & MUTHEN
09/22/2012  10:01 PM

INPUT INSTRUCTIONS

  title:	this is an example of a 2-group twin model
      for a categorical outcome using probit WLSMV
  	
  montecarlo:
  	names = u1 u2;
  	ngroups = 2;
  	nobs = 1000 1000;
      generate = u1-u2(1);
  	categorical = u1 u2;
  	nreps = 1;
  	save = ex5.22.dat;

  model population:
      [u1$1-u2$1*0];
      u1-u2@1; ! u* variances are fixed at 1 by default in
               ! the analysis
      u1 WITH u2*0.73; ! Value for MZ twins
  	
  model population-g2:
  	u1 WITH u2*0.41;  ! Value for DZ twins

  model:
      [u1$1-u2$1*0](1);
      u1 WITH u2*0.73(covmz);
  		
  model g2:
      u1 WITH u2*0.41(covdz);

  model constraint:
      NEW(a*.8 c*.3 h*.64);
      ! e is not a free parameter but obtained as 1 - (a**2 + c**2)
      covmz = a**2 + c**2;
      covdz = 0.5*a**2 + c**2;
      h = a**2/1;

  output:
  	tech9;




INPUT READING TERMINATED NORMALLY



this is an example of a 2-group twin model
for a categorical outcome using probit WLSMV

SUMMARY OF ANALYSIS

Number of groups                                                 2
Number of observations
   Group G1                                                   1000
   Group G2                                                   1000

Number of replications
    Requested                                                    1
    Completed                                                    1
Value of seed                                                    0

Number of dependent variables                                    2
Number of independent variables                                  0
Number of continuous latent variables                            0

Observed dependent variables

  Binary and ordered categorical (ordinal)
   U1          U2


Estimator                                                    WLSMV
Maximum number of iterations                                  1000
Convergence criterion                                    0.500D-04
Maximum number of steepest descent iterations                   20
Parameterization                                             DELTA



SAMPLE STATISTICS FOR THE FIRST REPLICATION


     SAMPLE STATISTICS FOR G1


           Means/Intercepts/Thresholds
              U1$1          U2$1
              ________      ________
      1         0.005        -0.015


           Covariances/Correlations/Residual Correlations
              U1            U2
              ________      ________
 U1
 U2             0.771


     SAMPLE STATISTICS FOR G2


           Means/Intercepts/Thresholds
              U1$1          U2$1
              ________      ________
      1         0.010         0.028


           Covariances/Correlations/Residual Correlations
              U1            U2
              ________      ________
 U1
 U2             0.411


MODEL FIT INFORMATION

Number of Free Parameters                        3

Chi-Square Test of Model Fit

        Degrees of freedom                       3

        Mean                                 0.905
        Std Dev                              0.000
        Number of successful computations        1

             Proportions                   Percentiles
        Expected    Observed         Expected       Observed
           0.990       1.000            0.115          0.905
           0.980       1.000            0.185          0.905
           0.950       1.000            0.352          0.905
           0.900       1.000            0.584          0.905
           0.800       0.000            1.005          0.905
           0.700       0.000            1.424          0.905
           0.500       0.000            2.366          0.905
           0.300       0.000            3.665          0.905
           0.200       0.000            4.642          0.905
           0.100       0.000            6.251          0.905
           0.050       0.000            7.815          0.905
           0.020       0.000            9.837          0.905
           0.010       0.000           11.345          0.905

RMSEA (Root Mean Square Error Of Approximation)

        Mean                                 0.000
        Std Dev                              0.000
        Number of successful computations        1

             Proportions                   Percentiles
        Expected    Observed         Expected       Observed
           0.990       0.000            0.000          0.000
           0.980       0.000            0.000          0.000
           0.950       0.000            0.000          0.000
           0.900       0.000            0.000          0.000
           0.800       0.000            0.000          0.000
           0.700       0.000            0.000          0.000
           0.500       0.000            0.000          0.000
           0.300       0.000            0.000          0.000
           0.200       0.000            0.000          0.000
           0.100       0.000            0.000          0.000
           0.050       0.000            0.000          0.000
           0.020       0.000            0.000          0.000
           0.010       0.000            0.000          0.000

WRMR (Weighted Root Mean Square Residual)

        Mean                                 0.442
        Std Dev                              0.000
        Number of successful computations        1

             Proportions                   Percentiles
        Expected    Observed         Expected       Observed
           0.990       0.000            0.442          0.442
           0.980       0.000            0.442          0.442
           0.950       0.000            0.442          0.442
           0.900       0.000            0.442          0.442
           0.800       0.000            0.442          0.442
           0.700       0.000            0.442          0.442
           0.500       0.000            0.442          0.442
           0.300       0.000            0.442          0.442
           0.200       0.000            0.442          0.442
           0.100       0.000            0.442          0.442
           0.050       0.000            0.442          0.442
           0.020       0.000            0.442          0.442
           0.010       0.000            0.442          0.442


MODEL RESULTS

                           ESTIMATES              S. E.     M. S. E.  95%  % Sig
              Population   Average   Std. Dev.   Average             Cover Coeff
Group G1

 U1       WITH
  U2               0.730     0.7707     0.0000     0.0262     0.0017 1.000 1.000

 Thresholds
  U1$1             0.000     0.0069     0.0000     0.0236     0.0000 1.000 0.000
  U2$1             0.000     0.0069     0.0000     0.0236     0.0000 1.000 0.000

Group G2

 U1       WITH
  U2               0.410     0.4114     0.0000     0.0436     0.0000 1.000 1.000

 Thresholds
  U1$1             0.000     0.0069     0.0000     0.0236     0.0000 1.000 0.000
  U2$1             0.000     0.0069     0.0000     0.0236     0.0000 1.000 0.000

 New/Additional Parameters
  A                0.800     0.8477     0.0000     0.0600     0.0023 1.000 1.000
  C                0.300     0.2284     0.0000     0.1993     0.0051 1.000 0.000
  H                0.640     0.7185     0.0000     0.1018     0.0062 1.000 1.000


QUALITY OF NUMERICAL RESULTS

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


TECHNICAL OUTPUT


     PARAMETER SPECIFICATION FOR G1


           TAU
              U1$1          U2$1
              ________      ________
      1           1             1


           NU
              U1            U2
              ________      ________
      1           0             0


           THETA
              U1            U2
              ________      ________
 U1                 0
 U2                 2             0


           DELTA
              U1            U2
              ________      ________
      1           0             0


     PARAMETER SPECIFICATION FOR G2


           TAU
              U1$1          U2$1
              ________      ________
      1           1             1


           NU
              U1            U2
              ________      ________
      1           0             0


           THETA
              U1            U2
              ________      ________
 U1                 0
 U2                 3             0


           DELTA
              U1            U2
              ________      ________
      1           0             0


     PARAMETER SPECIFICATION FOR THE ADDITIONAL PARAMETERS


           New/Additional Parameters
              A             C             H
              ________      ________      ________
 1                  4             5             6


     STARTING VALUES FOR G1


           TAU
              U1$1          U2$1
              ________      ________
      1         0.000         0.000


           NU
              U1            U2
              ________      ________
      1         0.000         0.000


           THETA
              U1            U2
              ________      ________
 U1             1.000
 U2             0.730         1.000


           DELTA
              U1            U2
              ________      ________
      1         1.000         1.000


     STARTING VALUES FOR G2


           TAU
              U1$1          U2$1
              ________      ________
      1         0.000         0.000


           NU
              U1            U2
              ________      ________
      1         0.000         0.000


           THETA
              U1            U2
              ________      ________
 U1             1.000
 U2             0.410         1.000


           DELTA
              U1            U2
              ________      ________
      1         1.000         1.000


     STARTING VALUES FOR THE ADDITIONAL PARAMETERS


           New/Additional Parameters
              A             C             H
              ________      ________      ________
 1              0.800         0.300         0.640


     POPULATION VALUES FOR G1


           TAU
              U1$1          U2$1
              ________      ________
 1              0.000         0.000


           NU
              U1            U2
              ________      ________
 1              0.000         0.000


           THETA
              U1            U2
              ________      ________
 U1             1.000
 U2             0.730         1.000


     POPULATION VALUES FOR G2


           TAU
              U1$1          U2$1
              ________      ________
 1              0.000         0.000


           NU
              U1            U2
              ________      ________
 1              0.000         0.000


           THETA
              U1            U2
              ________      ________
 U1             1.000
 U2             0.410         1.000


TECHNICAL 9 OUTPUT

  Error messages for each replication (if any)



SAVEDATA INFORMATION

  Order of variables

    U1
    U2
    GROUP

  Save file
    ex5.22.dat

  Save file format           Free
  Save file record length    10000


     Beginning Time:  22:01:55
        Ending Time:  22:01:55
       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-2012 Muthen & Muthen

Back to examples