** Mx startup successful **
  
  **MX-PC 1.54** Job started on 03/03/03 at 12:28:52
 
 
 
 
 The following MX script lines were read for group  1
 
 ANALYSIS OF AUSTRALIAN BMI DATA-YOUNG FEMALE MZ TWINS PAIRS-MX DIY VERSION
 !
 DA  NG=1 NI=2 NO=0
 BEGIN MATRICES;
 A LO 1 1 FR ! ADDITIVE GENETIC VARIANCE
 C LO 1 1 FR ! SHARED ENVIRONMENTAL VARIANCE
 E LO 1 1 FR ! NON-SHARED ENVIRONMENTAL VARIANCE
 M FU 2 2  ! THIS WILL BE OBSERVED MZ COVARIANCE MATRIX
 D FU 2 2  ! THIS WILL BE OBSERVED DZ COVARIANCE MATRIX
 G FU 1 1  ! COEFFICIENT OF 0.5 FOR DZ PAIRS
 N FU 1 1  ! SAMPLE SIZE FOR MZ PAIRS (FEMALE IN THIS ILLUSTRATION)
 K FU 1 1  ! SAMPLE SIZE FOR DZ PAIRS (FEMALE IN THIS ILLUSTRATION)
 P FU 1 1  ! ORDER OF MATRICES (I.E. NUMBER OF VARIABLES =2 IN THIS CASE)
 END MATRICES;
 MAT G 0.5
 MAT P 2
 MAT M
 0.7247 0.5891
 0.5891 0.7915
 MAT N 532
 MAT D
 0.7786 0.2461
 0.2461 0.8365
 MAT K 326
 MAT A 0.25
 MAT C 0.25
 MAT E 0.25
 BEGIN ALGEBRA;
 T=N+K; ! TOTAL SAMPLE SIZE
 U=(A+C+E | A+C _ A+C | A+C+E); ! EXPECTED MZ COVARIANCE MATRIX
 V=(A+C+E | G*A+C _ G*A+C | A+C+E); ! EXPECTED DZ COVARIANCE MATRIX
 H=N*(\LN (\DET(U))-\LN (\DET(M)) + \TR((U~ *M))-P); ! FIT-FUNCTION FOR MZ GROUP
 J=K*(\LN (\DET(V))-\LN (\DET(D)) + \TR((V~ *D))-P); ! FIT-FUNCTION FOR DZ GROUP
 F=H+J;
 END ALGEBRA;
 BO 0.01 1.0 E(1,1)
 BO 0.0  1.0 C(1,1) A(1,1)
 CO F ;
 OPTION USER DF=6
 END
  
  
  PARAMETER SPECIFICATIONS
  
  GROUP NUMBER: 1
  
Analysis of Australian BMI data-young female MZ twins pairs-MX DIY version                                                      
  
  MATRIX A
 This is a LOWER TRIANGULAR matrix of order    1 by    1
    1
 1  1
  
  MATRIX C
 This is a LOWER TRIANGULAR matrix of order    1 by    1
    1
 1  2
  
  MATRIX D
 This is a FULL matrix of order    2 by    2
  It has no free parameters specified
  
  MATRIX E
 This is a LOWER TRIANGULAR matrix of order    1 by    1
    1
 1  3
  
  MATRIX F
 This is a computed FULL matrix of order    1 by    1
  It has no free parameters specified
  
  MATRIX G
 This is a FULL matrix of order    1 by    1
  It has no free parameters specified
  
  MATRIX H
 This is a computed FULL matrix of order    1 by    1
  It has no free parameters specified
  
  MATRIX J
 This is a computed FULL matrix of order    1 by    1
  It has no free parameters specified
  
  MATRIX K
 This is a FULL matrix of order    1 by    1
  It has no free parameters specified
  
  MATRIX M
 This is a FULL matrix of order    2 by    2
  It has no free parameters specified
  
  MATRIX N
 This is a FULL matrix of order    1 by    1
  It has no free parameters specified
  
  MATRIX P
 This is a FULL matrix of order    1 by    1
  It has no free parameters specified
  
  MATRIX T
 This is a computed FULL matrix of order    1 by    1
  It has no free parameters specified
  
  MATRIX U
 This is a computed FULL matrix of order    2 by    2
  It has no free parameters specified
  
  MATRIX V
 This is a computed FULL matrix of order    2 by    2
  It has no free parameters specified
  
  Mx starting optimization; number of parameters =  3
  
  
  
  MX PARAMETER ESTIMATES
  
  GROUP NUMBER: 1
  
Analysis of Australian BMI data-young female MZ twins pairs-MX DIY version                                                      
  
  MATRIX A
 This is a LOWER TRIANGULAR matrix of order    1 by    1
          1
 1   0.6196
  
  MATRIX C
 This is a LOWER TRIANGULAR matrix of order    1 by    1
          1
 1   0.0000
  
  MATRIX D
 This is a FULL matrix of order    2 by    2
          1       2
 1   0.7786  0.2461
 2   0.2461  0.8365
  
  MATRIX E
 This is a LOWER TRIANGULAR matrix of order    1 by    1
          1
 1   0.1734
  
  MATRIX F
 This is a computed FULL matrix of order    1 by    1
  [=H+J]
             1
 1      8.0344
  
  MATRIX G
 This is a FULL matrix of order    1 by    1
          1
 1   0.5000
  
  MATRIX H
 This is a computed FULL matrix of order    1 by    1
  [=N*(\LN(\DET(U))-\LN(\DET(M))+\TR((U~*M))-P)]
             1
 1      3.3785
  
  MATRIX J
 This is a computed FULL matrix of order    1 by    1
  [=K*(\LN(\DET(V))-\LN(\DET(D))+\TR((V~*D))-P)]
             1
 1      4.6559
  
  MATRIX K
 This is a FULL matrix of order    1 by    1
             1
 1    326.0000
  
  MATRIX M
 This is a FULL matrix of order    2 by    2
          1       2
 1   0.7247  0.5891
 2   0.5891  0.7915
  
  MATRIX N
 This is a FULL matrix of order    1 by    1
             1
 1    532.0000
  
  MATRIX P
 This is a FULL matrix of order    1 by    1
             1
 1      2.0000
  
  MATRIX T
 This is a computed FULL matrix of order    1 by    1
  [=N+K]
             1
 1    858.0000
  
  MATRIX U
 This is a computed FULL matrix of order    2 by    2
  [=(A+C+E|A+C_A+C|A+C+E)]
          1       2
 1   0.7930  0.6196
 2   0.6196  0.7930
  
  MATRIX V
 This is a computed FULL matrix of order    2 by    2
  [=(A+C+E|G*A+C_G*A+C|A+C+E)]
          1       2
 1   0.7930  0.3098
 2   0.3098  0.7930
  
 Your model has    3 estimated parameters and      6 Observed statistics
 Note that Observed statistics includes your DF= adjustment of   6
  
 User defined function value =     8.034
 'Degrees of freedom' >>>>>>>>>>>>>>>>         3
  
  
  
  
 This problem used  0.3% of my workspace
  
 Task                     Time elapsed (DD:HH:MM:SS)
 Reading script & data      0: 0: 0: 0.20
 Execution                  0: 0: 0:-0.03
 TOTAL                      0: 0: 0: 0.17
  
 Total number of warnings issued: 0
 ______________________________________________________________________________