The following MX script lines were read for group 1 #DEFINE NVAR 2 G1: CALCULATION DATA CALC NGROUPS=3 BEGIN MATRICES; X LOWER NVAR NVAR FREE ! GENETIC STRUCTURE Y LOWER NVAR NVAR FREE ! COMMON ENVIRONMENTAL STRUCTURE Z LOWER NVAR NVAR FREE ! SPECIFIC ENVIRONMENTAL STRUCTURE G FULL 1 NVAR FREE ! MEANS H FULL 1 1 END MATRICES; MATRIX H .5 ST 1 ALL START 400 G 1 1 START 70 G 1 2 ST 18 X 1 1 Y 1 1 Z 1 1 ST 4 X 2 2 Y 2 2 Z 2 2 BEGIN ALGEBRA; A= X*X' ; C= Y*Y' ; E= Z*Z' ; P = A%(A+C+E)| C%(A+C+E) | E%(A+C+E) ; !STANDARDIZED ESTIMATES K=\STND(A) ; !GENETIC CORRELATION L=\STND(C); !COMMON ENVIRONMENTAL CORRELATION M=\STND(E); !SPECIFIC ENVIRONMENTAL CORRELATION END ALGEBRA; END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS AND SIB DA NINPUT_VARS=14 MISSING=-200.00 NOTE: Missing value *string* set to '-200.00' RECTANGULAR FILE=MRI4.DAT Rectangular continuous data read initiated NOTE: Rectangular file contained 322 records with data that contained a total of 1835 observations LABELS FAMNR ZYG WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 WM3 BBGMT3 BBWMT3 WM4 BBGMT4 BBWMT4 SELECT IF ZYG < 3 NOTE: Select if yields 144 data vectors for analysis NOTE: Vectors contain a total of 809 observations SELECT BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 BBGMT4 WM4 ; MATRICES = G 1 NOTE: Selection yields 144 data vectors for analysis NOTE: Vectors contain a total of 401 observations MEANS G|G|G|G; !MODEL FOR MEANS COVARIANCES ! MODEL FOR COVARIANCES A+C+E | A+C | H@A+C |H@A+C_ A+C | A+C+E | H@A+C |H@A+C_ H@A+C | H@A+C | A+C+E |H@A+C _ H@A+C | H@A+C | H@A+C | A+C+E; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS AND SIB DA NINPUT_VARS=14 MISSING=-200.00 NOTE: Missing value *string* set to '-200.00' RECTANGULAR FILE=MRI4.DAT Rectangular continuous data read initiated NOTE: Rectangular file contained 322 records with data that contained a total of 1835 observations LABELS FAMNR ZYG WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 WM3 BBGMT3 BBWMT3 WM4 BBGMT4 BBWMT4 SELECT IF ZYG > 2 NOTE: Select if yields 178 data vectors for analysis NOTE: Vectors contain a total of 1026 observations SELECT BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 BBGMT4 WM4 ; MATRICES = G 1 NOTE: Selection yields 178 data vectors for analysis NOTE: Vectors contain a total of 532 observations MEANS G|G|G|G; !MODEL FOR MEANS COVARIANCES ! MODEL FOR COVARIANCES A+C+E |H@A+C | H@A+C |H@A+C_ H@A+C | A+C+E | H@A+C |H@A+C_ H@A+C | H@A+C | A+C+E |H@A+C _ H@A+C | H@A+C | H@A+C | A+C+E; OPTIONS RSIDUAL END Summary of VL file data for group 2 BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 54.0000 106.0000 54.0000 107.0000 12.0000 55.0000 Mean 705.8689 62.3434 697.0530 61.3925 708.9617 61.8987 Variance 3234.4423 157.6016 3779.2227 201.7276 2227.0761 196.3816 Minimum 570.6100 26.1900 523.6600 27.7400 608.9900 30.1200 Maximum 860.5200 90.4800 836.6400 91.7700 805.5700 97.2700 WM4 Code 8.0000 Number 13.0000 Mean 56.6969 Variance 191.3088 Minimum 30.9500 Maximum 75.3600 Summary of VL file data for group 3 BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 Number 57.0000 146.0000 58.0000 139.0000 19.0000 90.0000 Mean 699.6267 62.3736 682.9521 61.8701 706.4447 66.1716 Variance 2234.7938 181.6186 2733.4405 127.3058 2679.9778 135.6165 Minimum 562.2500 24.9400 568.7000 28.5700 606.5200 32.5700 Maximum 799.9600 94.3700 775.9000 86.6600 798.1000 95.3000 BBGMT4 WM4 Code 7.0000 8.0000 Number 4.0000 19.0000 Mean 694.2000 61.3105 Variance 3169.6562 161.9674 Minimum 630.2600 32.5000 Maximum 779.0800 80.1200 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: calculation MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 10 11 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 2 by 6 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 2 G2: MZ twin pairs and sib MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 10 11 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 2 by 6 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 GROUP NUMBER: 3 G3: DZ twin pairs and sib MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 10 11 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 2 by 6 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 4 2 5 6 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 7 2 8 9 Mx starting optimization; number of parameters = 11 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: calculation MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 2030.0534 209.1538 2 209.1538 114.1243 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 309.7520 5.2243 2 5.2243 0.0881 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 529.3391 -28.7396 2 -28.7396 56.4224 MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 696.5500 62.4723 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.4345 2 0.4345 1.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 1.0000 2 1.0000 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 -0.1663 2 -0.1663 1.0000 MATRIX P This is a computed FULL matrix of order 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 1 7.0755E-01 1.1267E+00 1.0796E-01 2.8142E-02 1.8449E-01 2 1.1267E+00 6.6882E-01 2.8142E-02 5.1638E-04 -1.5481E-01 6 1 -1.5481E-01 2 3.3066E-01 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 -45.0561 2 -4.6421 9.6216 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.7600E+01 2 2.9684E-01 -2.7363E-06 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 23.0074 2 -1.2491 7.4069 GROUP NUMBER: 2 G2: MZ twin pairs and sib MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 2030.0534 209.1538 2 209.1538 114.1243 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 309.7520 5.2243 2 5.2243 0.0881 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 529.3391 -28.7396 2 -28.7396 56.4224 MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 696.5500 62.4723 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.4345 2 0.4345 1.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 1.0000 2 1.0000 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 -0.1663 2 -0.1663 1.0000 MATRIX P This is a computed FULL matrix of order 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 1 7.0755E-01 1.1267E+00 1.0796E-01 2.8142E-02 1.8449E-01 2 1.1267E+00 6.6882E-01 2.8142E-02 5.1638E-04 -1.5481E-01 6 1 -1.5481E-01 2 3.3066E-01 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 -45.0561 2 -4.6421 9.6216 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.7600E+01 2 2.9684E-01 -2.7363E-06 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 23.0074 2 -1.2491 7.4069 Vector of OBSERVED means BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 Mean 705.8689 62.3434 697.0530 61.3925 708.9617 61.8987 BBGMT4 WM4 Mean 0.0000 56.6969 Vector of EXPECTED means BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 Mean 696.5500 62.4723 696.5500 62.4723 696.5500 62.4723 BBGMT4 WM4 Mean 696.5500 62.4723 EXPECTED COVARIANCE MATRIX BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 BBGMT1 2869.1445 WM1 185.6385 170.6348 BBGMT2 2339.8054 214.3780 2869.1445 WM2 214.3780 114.2124 185.6385 170.6348 BBGMT3 1324.7787 109.8012 1324.7787 109.8012 2869.1445 WM3 109.8012 57.1503 109.8012 57.1503 185.6385 170.6348 BBGMT4 1324.7787 109.8012 1324.7787 109.8012 1324.7787 109.8012 WM4 109.8012 57.1503 109.8012 57.1503 109.8012 57.1503 BBGMT4 WM4 BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 BBGMT4 2869.1445 WM4 185.6385 170.6348 Function value of this group: 3422.2464 Where the fit function is -2 * Log-likelihood of raw data GROUP NUMBER: 3 G3: DZ twin pairs and sib MATRIX A This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 2030.0534 209.1538 2 209.1538 114.1243 MATRIX C This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 309.7520 5.2243 2 5.2243 0.0881 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Z*Z'] 1 2 1 529.3391 -28.7396 2 -28.7396 56.4224 MATRIX G This is a FULL matrix of order 1 by 2 1 2 1 696.5500 62.4723 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.4345 2 0.4345 1.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 1.0000 2 1.0000 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 -0.1663 2 -0.1663 1.0000 MATRIX P This is a computed FULL matrix of order 2 by 6 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 1 7.0755E-01 1.1267E+00 1.0796E-01 2.8142E-02 1.8449E-01 2 1.1267E+00 6.6882E-01 2.8142E-02 5.1638E-04 -1.5481E-01 6 1 -1.5481E-01 2 3.3066E-01 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 -45.0561 2 -4.6421 9.6216 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.7600E+01 2 2.9684E-01 -2.7363E-06 MATRIX Z This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 23.0074 2 -1.2491 7.4069 Vector of OBSERVED means BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 Mean 699.6267 62.3736 682.9521 61.8701 706.4447 66.1716 BBGMT4 WM4 Mean 694.2000 61.3105 Vector of EXPECTED means BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 Mean 696.5500 62.4723 696.5500 62.4723 696.5500 62.4723 BBGMT4 WM4 Mean 696.5500 62.4723 EXPECTED COVARIANCE MATRIX BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 BBGMT1 2869.1445 WM1 185.6385 170.6348 BBGMT2 1324.7787 109.8012 2869.1445 WM2 109.8012 57.1503 185.6385 170.6348 BBGMT3 1324.7787 109.8012 1324.7787 109.8012 2869.1445 WM3 109.8012 57.1503 109.8012 57.1503 185.6385 170.6348 BBGMT4 1324.7787 109.8012 1324.7787 109.8012 1324.7787 109.8012 WM4 109.8012 57.1503 109.8012 57.1503 109.8012 57.1503 BBGMT4 WM4 BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 BBGMT4 2869.1445 WM4 185.6385 170.6348 Function value of this group: 4550.9224 Where the fit function is -2 * Log-likelihood of raw data Your model has 11 estimated parameters and 933 Observed statistics -2 times log-likelihood of data >>> 7973.169 Degrees of freedom >>>>>>>>>>>>>>>> 922 This problem used 5.9% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.87 Execution 0: 0: 0: 8.53 TOTAL 0: 0: 0: 9.40 Total number of warnings issued: 0 ______________________________________________________________________________