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

INPUT INSTRUCTIONS

  ! SCRIPT NAME        : rawVC1a  (dp)
  ! GOAL                : To evaluate best model for variance components
  ! DATA                : continuous
  ! INPUT                : raw data
  ! UNI/BI/MULTI        : uni
  ! DATA-GROUPS        : MZ DZ
  ! MEANS MODEL        : grand mean, age effect, sex effect
  ! VARIANCE COVARIANCE MODEL(S)        ## 1.ACE  2.AE  3.CE  4.E

  data: file is example.dat;

  variable:

    names are country famid zygos sex1 age1 height1 weight1 bmi1
                                  sex2 age2 height2 weight2 bmi2;

    usevar are bmi1 sex1 age1 bmi2 sex2 age2 g;

    grouping=g(1=MZ 2=DZ);  ! specify the two groups MZ and DZ

    missing=all(-1); ! specify missing data symbol

  define: if (zygos==1 .or. zygos==3) then g=1 else g=2; ! defines the two groups

  model :
    bmi1 on sex1 (b1)
            age1 (b2);
    bmi2 on sex2 (b1)
            age2 (b2);
    [bmi1 bmi2] (m);
    bmi1 bmi2 (v);
    bmi1 with bmi2 (c1);

  model DZ:
    bmi1 with bmi2 (c2);

  model constraint:

  ! Uncomment for Model ACE
    new(a c e x y z);
    a=x*x;
    c=y*y;
    e=z*z;
    v=a+c+e;
    c1=a+c;
    c2=0.5*a+c;

  ! Uncomment for Model AE
  ! c=0;

  ! Uncomment for Model CE
  ! a=0;

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



*** WARNING
  Data set contains cases with missing on all variables except
  x-variables.  These cases were not included in the analysis.
  Number of cases with missing on all variables except x-variables:  3
   1 WARNING(S) FOUND IN THE INPUT INSTRUCTIONS




SUMMARY OF ANALYSIS

Number of groups                                                 2
Number of observations
   Group MZ                                                    138
   Group DZ                                                    166

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

Observed dependent variables

  Continuous
   BMI1        BMI2

Observed independent variables
   SEX1        AGE1        SEX2        AGE2

Variables with special functions

  Grouping variable     G

Estimator                                                       ML
Information matrix                                        OBSERVED
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

Input data file(s)
  example.dat

Input data format  FREE


SUMMARY OF DATA

   Group MZ
     Number of missing data patterns             3

   Group DZ
     Number of missing data patterns             3


COVARIANCE COVERAGE OF DATA

Minimum covariance coverage value   0.100


     PROPORTION OF DATA PRESENT FOR MZ


           Covariance Coverage
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           0.935
 BMI2           0.891         0.957
 SEX1           0.935         0.957         1.000
 AGE1           0.935         0.957         1.000         1.000
 SEX2           0.935         0.957         1.000         1.000         1.000
 AGE2           0.935         0.957         1.000         1.000         1.000


           Covariance Coverage
              AGE2
              ________
 AGE2           1.000


     PROPORTION OF DATA PRESENT FOR DZ


           Covariance Coverage
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           0.940
 BMI2           0.843         0.904
 SEX1           0.940         0.904         1.000
 AGE1           0.940         0.904         1.000         1.000
 SEX2           0.940         0.904         1.000         1.000         1.000
 AGE2           0.940         0.904         1.000         1.000         1.000


           Covariance Coverage
              AGE2
              ________
 AGE2           1.000



THE MODEL ESTIMATION TERMINATED NORMALLY

     THE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES MAY NOT BE
     TRUSTWORTHY FOR SOME PARAMETERS DUE TO A NON-POSITIVE DEFINITE
     FIRST-ORDER DERIVATIVE PRODUCT MATRIX.  THIS MAY BE DUE TO THE STARTING
     VALUES BUT MAY ALSO BE AN INDICATION OF MODEL NONIDENTIFICATION.  THE
     CONDITION NUMBER IS       0.420D-13.  PROBLEM INVOLVING PARAMETER 11.




TESTS OF MODEL FIT

Chi-Square Test of Model Fit

          Value                             30.813
          Degrees of Freedom                    20
          P-Value                           0.0577

Chi-Square Contributions From Each Group

          MZ                                11.353
          DZ                                19.460

Chi-Square Test of Model Fit for the Baseline Model

          Value                            188.497
          Degrees of Freedom                    18
          P-Value                           0.0000

CFI/TLI

          CFI                                0.937
          TLI                                0.943

Loglikelihood

          H0 Value                       -4327.091
          H1 Value                       -4311.684

Information Criteria

          Number of Free Parameters              6
          Akaike (AIC)                    8666.182
          Bayesian (BIC)                  8688.484
          Sample-Size Adjusted BIC        8669.455
            (n* = (n + 2) / 24)

RMSEA (Root Mean Square Error Of Approximation)

          Estimate                           0.060
          90 Percent C.I.                    0.000  0.099
          Probability RMSEA <= .05           0.321

SRMR (Standardized Root Mean Square Residual)

          Value                              0.120



MODEL RESULTS

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

Group MZ

 BMI1     ON
    SEX1               0.837      0.340      2.463      0.014
    AGE1               0.097      0.014      6.773      0.000

 BMI2     ON
    SEX2               0.837      0.340      2.463      0.014
    AGE2               0.097      0.014      6.773      0.000

 BMI1     WITH
    BMI2               9.655      0.991      9.739      0.000

 Intercepts
    BMI1              20.220      0.590     34.254      0.000
    BMI2              20.220      0.590     34.254      0.000

 Residual Variances
    BMI1              13.197      0.882     14.966      0.000
    BMI2              13.197      0.882     14.966      0.000

Group DZ

 BMI1     ON
    SEX1               0.837      0.340      2.463      0.014
    AGE1               0.097      0.014      6.773      0.000

 BMI2     ON
    SEX2               0.837      0.340      2.463      0.014
    AGE2               0.097      0.014      6.773      0.000

 BMI1     WITH
    BMI2               4.827      0.496      9.739      0.000

 Intercepts
    BMI1              20.220      0.590     34.254      0.000
    BMI2              20.220      0.590     34.254      0.000

 Residual Variances
    BMI1              13.197      0.882     14.966      0.000
    BMI2              13.197      0.882     14.966      0.000

 New/Additional Parameters
    A                  9.655      0.991      9.739      0.000
    C                  0.000      0.000      0.000      1.000
    E                  3.542      0.456      7.764      0.000
    X                  3.107      0.160     19.477      0.000
    Y                  0.000      0.755      0.000      1.000
    Z                  1.882      0.121     15.529      0.000


QUALITY OF NUMERICAL RESULTS

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


TECHNICAL 1 OUTPUT


     PARAMETER SPECIFICATION FOR MZ


           NU
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
      1           0             0             0             0             0


           NU
              AGE2
              ________
      1           0


           LAMBDA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1               0             0             0             0             0
 BMI2               0             0             0             0             0
 SEX1               0             0             0             0             0
 AGE1               0             0             0             0             0
 SEX2               0             0             0             0             0
 AGE2               0             0             0             0             0


           LAMBDA
              AGE2
              ________
 BMI1               0
 BMI2               0
 SEX1               0
 AGE1               0
 SEX2               0
 AGE2               0


           THETA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1               0
 BMI2               0             0
 SEX1               0             0             0
 AGE1               0             0             0             0
 SEX2               0             0             0             0             0
 AGE2               0             0             0             0             0


           THETA
              AGE2
              ________
 AGE2               0


           ALPHA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
      1           1             1             0             0             0


           ALPHA
              AGE2
              ________
      1           0


           BETA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1               0             0             2             3             0
 BMI2               0             0             0             0             2
 SEX1               0             0             0             0             0
 AGE1               0             0             0             0             0
 SEX2               0             0             0             0             0
 AGE2               0             0             0             0             0


           BETA
              AGE2
              ________
 BMI1               0
 BMI2               3
 SEX1               0
 AGE1               0
 SEX2               0
 AGE2               0


           PSI
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1               4
 BMI2               5             4
 SEX1               0             0             0
 AGE1               0             0             0             0
 SEX2               0             0             0             0             0
 AGE2               0             0             0             0             0


           PSI
              AGE2
              ________
 AGE2               0


     PARAMETER SPECIFICATION FOR DZ


           NU
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
      1           0             0             0             0             0


           NU
              AGE2
              ________
      1           0


           LAMBDA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1               0             0             0             0             0
 BMI2               0             0             0             0             0
 SEX1               0             0             0             0             0
 AGE1               0             0             0             0             0
 SEX2               0             0             0             0             0
 AGE2               0             0             0             0             0


           LAMBDA
              AGE2
              ________
 BMI1               0
 BMI2               0
 SEX1               0
 AGE1               0
 SEX2               0
 AGE2               0


           THETA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1               0
 BMI2               0             0
 SEX1               0             0             0
 AGE1               0             0             0             0
 SEX2               0             0             0             0             0
 AGE2               0             0             0             0             0


           THETA
              AGE2
              ________
 AGE2               0


           ALPHA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
      1           1             1             0             0             0


           ALPHA
              AGE2
              ________
      1           0


           BETA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1               0             0             2             3             0
 BMI2               0             0             0             0             2
 SEX1               0             0             0             0             0
 AGE1               0             0             0             0             0
 SEX2               0             0             0             0             0
 AGE2               0             0             0             0             0


           BETA
              AGE2
              ________
 BMI1               0
 BMI2               3
 SEX1               0
 AGE1               0
 SEX2               0
 AGE2               0


           PSI
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1               4
 BMI2               6             4
 SEX1               0             0             0
 AGE1               0             0             0             0
 SEX2               0             0             0             0             0
 AGE2               0             0             0             0             0


           PSI
              AGE2
              ________
 AGE2               0


     PARAMETER SPECIFICATION FOR THE ADDITIONAL PARAMETERS


           New/Additional Parameters
              A             C             E             X             Y
              ________      ________      ________      ________      ________
 1                  7             8             9            10            11


           New/Additional Parameters
              Z
              ________
 1                 12


     STARTING VALUES FOR MZ


           NU
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
      1         0.000         0.000         0.000         0.000         0.000


           NU
              AGE2
              ________
      1         0.000


           LAMBDA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           1.000         0.000         0.000         0.000         0.000
 BMI2           0.000         1.000         0.000         0.000         0.000
 SEX1           0.000         0.000         1.000         0.000         0.000
 AGE1           0.000         0.000         0.000         1.000         0.000
 SEX2           0.000         0.000         0.000         0.000         1.000
 AGE2           0.000         0.000         0.000         0.000         0.000


           LAMBDA
              AGE2
              ________
 BMI1           0.000
 BMI2           0.000
 SEX1           0.000
 AGE1           0.000
 SEX2           0.000
 AGE2           1.000


           THETA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           0.000
 BMI2           0.000         0.000
 SEX1           0.000         0.000         0.000
 AGE1           0.000         0.000         0.000         0.000
 SEX2           0.000         0.000         0.000         0.000         0.000
 AGE2           0.000         0.000         0.000         0.000         0.000


           THETA
              AGE2
              ________
 AGE2           0.000


           ALPHA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
      1         0.000         0.000         0.348        36.909         0.355


           ALPHA
              AGE2
              ________
      1        36.568


           BETA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           0.000         0.000         0.000         0.000         0.000
 BMI2           0.000         0.000         0.000         0.000         0.000
 SEX1           0.000         0.000         0.000         0.000         0.000
 AGE1           0.000         0.000         0.000         0.000         0.000
 SEX2           0.000         0.000         0.000         0.000         0.000
 AGE2           0.000         0.000         0.000         0.000         0.000


           BETA
              AGE2
              ________
 BMI1           0.000
 BMI2           0.000
 SEX1           0.000
 AGE1           0.000
 SEX2           0.000
 AGE2           0.000


           PSI
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           6.683
 BMI2           0.000         6.326
 SEX1           0.000         0.000         0.444
 AGE1           0.000         0.000         1.958       208.238
 SEX2           0.000         0.000         0.210        -0.983         0.446
 AGE2           0.000         0.000        -0.993       154.076         2.931


           PSI
              AGE2
              ________
 AGE2         220.236


     STARTING VALUES FOR DZ


           NU
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
      1         0.000         0.000         0.000         0.000         0.000


           NU
              AGE2
              ________
      1         0.000


           LAMBDA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           1.000         0.000         0.000         0.000         0.000
 BMI2           0.000         1.000         0.000         0.000         0.000
 SEX1           0.000         0.000         1.000         0.000         0.000
 AGE1           0.000         0.000         0.000         1.000         0.000
 SEX2           0.000         0.000         0.000         0.000         1.000
 AGE2           0.000         0.000         0.000         0.000         0.000


           LAMBDA
              AGE2
              ________
 BMI1           0.000
 BMI2           0.000
 SEX1           0.000
 AGE1           0.000
 SEX2           0.000
 AGE2           1.000


           THETA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           0.000
 BMI2           0.000         0.000
 SEX1           0.000         0.000         0.000
 AGE1           0.000         0.000         0.000         0.000
 SEX2           0.000         0.000         0.000         0.000         0.000
 AGE2           0.000         0.000         0.000         0.000         0.000


           THETA
              AGE2
              ________
 AGE2           0.000


           ALPHA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
      1         0.000         0.000         0.325        35.344         0.181


           ALPHA
              AGE2
              ________
      1        33.743


           BETA
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           0.000         0.000         0.000         0.000         0.000
 BMI2           0.000         0.000         0.000         0.000         0.000
 SEX1           0.000         0.000         0.000         0.000         0.000
 AGE1           0.000         0.000         0.000         0.000         0.000
 SEX2           0.000         0.000         0.000         0.000         0.000
 AGE2           0.000         0.000         0.000         0.000         0.000


           BETA
              AGE2
              ________
 BMI1           0.000
 BMI2           0.000
 SEX1           0.000
 AGE1           0.000
 SEX2           0.000
 AGE2           0.000


           PSI
              BMI1          BMI2          SEX1          AGE1          SEX2
              ________      ________      ________      ________      ________
 BMI1           6.683
 BMI2           0.000         6.326
 SEX1           0.000         0.000         0.509
 AGE1           0.000         0.000         4.240       229.798
 SEX2           0.000         0.000        -0.011        -1.991         0.618
 AGE2           0.000         0.000        -0.830       134.283         6.002


           PSI
              AGE2
              ________
 AGE2         261.295


     STARTING VALUES FOR THE ADDITIONAL PARAMETERS


           New/Additional Parameters
              A             C             E             X             Y
              ________      ________      ________      ________      ________
 1              0.500         0.500         0.500         0.500         0.500


           New/Additional Parameters
              Z
              ________
 1              0.500


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