Mplus VERSION 7.3
MUTHEN & MUTHEN
09/22/2014   5:17 PM

INPUT INSTRUCTIONS

  title:	this is an example of a 2-group twin model
  	
  montecarlo:
  	names = y1 y2;
  	ngroups = 2;
  	nobs = 1000 1000;
  	nreps = 1;
  	save = ex5.21.dat;

  model population:
      [y1-y2*0];
      y1-y2*1;
      y1 WITH y2*0.73; ! Value for MZ twins
  	
  model population-g2:
  	y1 WITH y2*0.41;  ! Value for DZ twins

  model:
      [y1-y2*0] (1);
      y1-y2*1 (var);
      y1 WITH y2*0.73(covmz);
  		
  model g2:
      y1 WITH y2*0.41(covdz);

  model constraint:
      NEW(a*.8 c*.3 e*.52 h*.64);
      var = a**2 + c**2 + e**2;
      covmz = a**2 + c**2;
      covdz = 0.5*a**2 + c**2;
      h = a**2/(a**2 + c**2 + e**2);

  output:
  	tech9;




INPUT READING TERMINATED NORMALLY



this is an example of a 2-group twin model

SUMMARY OF ANALYSIS

Number of groups                                                 2
Number of observations
   Group G1                                                   1000
   Group G2                                                   1000
   Total sample size                                          2000

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

  Continuous
   Y1          Y2


Estimator                                                       ML
Information matrix                                        OBSERVED
Maximum number of iterations                                  1000
Convergence criterion                                    0.500D-04
Maximum number of steepest descent iterations                   20



SAMPLE STATISTICS FOR THE FIRST REPLICATION


     SAMPLE STATISTICS FOR G1


           Means
              Y1            Y2
              ________      ________
      1         0.013         0.035


           Covariances
              Y1            Y2
              ________      ________
 Y1             0.992
 Y2             0.759         1.055


           Correlations
              Y1            Y2
              ________      ________
 Y1             1.000
 Y2             0.741         1.000


     SAMPLE STATISTICS FOR G2


           Means
              Y1            Y2
              ________      ________
      1        -0.017         0.006


           Covariances
              Y1            Y2
              ________      ________
 Y1             0.984
 Y2             0.429         1.042


           Correlations
              Y1            Y2
              ________      ________
 Y1             1.000
 Y2             0.423         1.000


MODEL FIT INFORMATION

Number of Free Parameters                        4

Loglikelihood

    H0 Value

        Mean                             -5216.004
        Std Dev                              0.000
        Number of successful computations        1

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

    H1 Value

        Mean                             -5213.423
        Std Dev                              0.000
        Number of successful computations        1

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

Information Criteria

    Akaike (AIC)

        Mean                             10440.008
        Std Dev                              0.000
        Number of successful computations        1

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

    Bayesian (BIC)

        Mean                             10462.411
        Std Dev                              0.000
        Number of successful computations        1

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

    Sample-Size Adjusted BIC (n* = (n + 2) / 24)

        Mean                             10449.703
        Std Dev                              0.000
        Number of successful computations        1

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

Chi-Square Test of Model Fit

        Degrees of freedom                       6

        Mean                                 5.161
        Std Dev                              0.000
        Number of successful computations        1

             Proportions                   Percentiles
        Expected    Observed         Expected       Observed
           0.990       1.000            0.872          5.161
           0.980       1.000            1.134          5.161
           0.950       1.000            1.635          5.161
           0.900       1.000            2.204          5.161
           0.800       1.000            3.070          5.161
           0.700       1.000            3.828          5.161
           0.500       0.000            5.348          5.161
           0.300       0.000            7.231          5.161
           0.200       0.000            8.558          5.161
           0.100       0.000           10.645          5.161
           0.050       0.000           12.592          5.161
           0.020       0.000           15.033          5.161
           0.010       0.000           16.812          5.161

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

SRMR (Standardized Root Mean Square Residual)

        Mean                                 0.103
        Std Dev                              0.000
        Number of successful computations        1

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


MODEL RESULTS

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

 Y1       WITH
  Y2                  0.730     0.7530     0.0000     0.0278     0.0005 1.000 1.000

 Means
  Y1                  0.000     0.0079     0.0000     0.0200     0.0001 1.000 0.000
  Y2                  0.000     0.0079     0.0000     0.0200     0.0001 1.000 0.000

 Variances
  Y1                  1.000     1.0180     0.0000     0.0263     0.0003 1.000 1.000
  Y2                  1.000     1.0180     0.0000     0.0263     0.0003 1.000 1.000

Group G2

 Y1       WITH
  Y2                  0.410     0.4320     0.0000     0.0309     0.0005 1.000 1.000

 Means
  Y1                  0.000     0.0079     0.0000     0.0200     0.0001 1.000 0.000
  Y2                  0.000     0.0079     0.0000     0.0200     0.0001 1.000 0.000

 Variances
  Y1                  1.000     1.0180     0.0000     0.0263     0.0003 1.000 1.000
  Y2                  1.000     1.0180     0.0000     0.0263     0.0003 1.000 1.000

New/Additional Parameters
  A                   0.800     0.8012     0.0000     0.0346     0.0000 1.000 1.000
  C                   0.300     0.3332     0.0000     0.0776     0.0011 1.000 1.000
  E                   0.520     0.5148     0.0000     0.0114     0.0000 1.000 1.000
  H                   0.640     0.6306     0.0000     0.0533     0.0001 1.000 1.000


QUALITY OF NUMERICAL RESULTS

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


TECHNICAL OUTPUT


     PARAMETER SPECIFICATION FOR G1


           NU
              Y1            Y2
              ________      ________
      1           1             1


           THETA
              Y1            Y2
              ________      ________
 Y1                 2
 Y2                 3             2


     PARAMETER SPECIFICATION FOR G2


           NU
              Y1            Y2
              ________      ________
      1           1             1


           THETA
              Y1            Y2
              ________      ________
 Y1                 2
 Y2                 4             2


     PARAMETER SPECIFICATION FOR THE ADDITIONAL PARAMETERS


           New/Additional Parameters
              A             C             E             H
              ________      ________      ________      ________
 1                  5             6             7             8


     STARTING VALUES FOR G1


           NU
              Y1            Y2
              ________      ________
      1         0.000         0.000


           THETA
              Y1            Y2
              ________      ________
 Y1             1.000
 Y2             0.730         1.000


     STARTING VALUES FOR G2


           NU
              Y1            Y2
              ________      ________
      1         0.000         0.000


           THETA
              Y1            Y2
              ________      ________
 Y1             1.000
 Y2             0.410         1.000


     STARTING VALUES FOR THE ADDITIONAL PARAMETERS


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


     POPULATION VALUES FOR G1


           NU
              Y1            Y2
              ________      ________
 1              0.000         0.000


           THETA
              Y1            Y2
              ________      ________
 Y1             1.000
 Y2             0.730         1.000


     POPULATION VALUES FOR G2


           NU
              Y1            Y2
              ________      ________
 1              0.000         0.000


           THETA
              Y1            Y2
              ________      ________
 Y1             1.000
 Y2             0.410         1.000


TECHNICAL 9 OUTPUT

  Error messages for each replication (if any)



SAVEDATA INFORMATION

  Order of variables

    Y1
    Y2
    GROUP

  Save file
    ex5.21.dat

  Save file format           Free
  Save file record length    10000


     Beginning Time:  17:17:51
        Ending Time:  17:17:51
       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-2014 Muthen & Muthen

Back to examples