** Mx startup successful ** **MX-SGI version 1.50** !Genetic Simplex Model !Analysis of platelet count !Platelet count rescaled by dividing count by 100 to aid optimisation The following MX script lines were read for group 1 #DEFINE NVAR 3 #DEFINE NSIB 2 ! NUMBER OF SIBLINGS #NGROUPS 3 Note: #NGroup set number of groups to 3 TITLE GENETIC STRUCTURE - SIMPLEX CALCULATION BEGIN MATRICES X DIAG NVAR NVAR FREE !GENETIC INNOVATIONS W DIAG NVAR NVAR FREE !DOMINANCE Z DIAG NVAR NVAR FREE !SPECIFIC ENVIRONMENTAL INNOVATIONS N LOWER NVAR NVAR FREE !GENETIC TRANSMISSION PATHS O LOWER NVAR NVAR FREE !DOMINANCE TRANSMISSION PATHS P LOWER NVAR NVAR FREE !SPECIFIC ENVIRONMENTAL PATHS R DIAG NVAR NVAR FREE !MEASUREMENT ERROR L FULL NVAR 1 FREE !QTL EFFECT I IDEN NVAR NVAR M FULL NSIB NVAR FREE H FULL 1 1 !.5 B FULL 1 1 !.25 END MATRICES; MATRIX H .5 MATRIX B .25 SPECIFY N 0 13 0 0 14 0 SPECIFY O 0 15 0 0 16 0 SPECIFY P 0 17 0 0 18 0 SPECIFY R 21 21 21 !EQUATE MEASUREMENT ERROR VARIANCES FOR IDENTIFICATION SPECIFY M 40 41 42 40 41 42 BEGIN ALGEBRA; A= (I-N)~ & (X*X'); D= (I-O)~ & (W*W'); E= (I-P)~ & (Z*Z') + R*R'; Q= L*L'; !VARIANCE DUE TO QTL V= A+Q+D+E; !TOTAL VARIANCE END ALGEBRA; END The following MX script lines were read for group 2 G2: MONOZYGOTIC TWINS DATA NINPUT=13 RECTANGULAR FILE=MARKER6.DAT NOTE: Rectangular file contained 519 records with data LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 SELECT IF ZYGOSITY <3; NOTE: Select if yields 287 data vectors for analysis NOTE: Vectors contain a total of 2962 observations SELECT PLT11 PLT12 PLT13 PLT21 PLT22 PLT23; BEGIN MATRICES = GROUP 1; NOTE: Selection yields 277 data vectors for analysis NOTE: Vectors contain a total of 953 observations MEANS M; COVARIANCE A+D+Q+E | A+D+Q _ A+D+Q | A+D+Q+E; END The following MX script lines were read for group 3 G3: DIZYGOTIC TWINS DATA NINPUT=13 NMODEL=3 RECTANGULAR FILE=MARKER6.DAT NOTE: Rectangular file contained 519 records with data LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 SELECT IF ZYGOSITY >2; NOTE: Select if yields 232 data vectors for analysis NOTE: Vectors contain a total of 2662 observations SELECT PIBD0 PIBD1 PIBD2 PLT11 PLT12 PLT13 PLT21 PLT22 PLT23; DEFINITION PIBD0 PIBD1 PIBD2; NOTE: Selection yields 232 data vectors for analysis NOTE: Vectors contain a total of 1734 observations BEGIN MATRICES = GROUP 1; NOTE: Definition yields 230 data vectors for analysis NOTE: Vectors contain a total of 1038 observations K FULL 3 1 !IBD PROBABILITIES (FROM GENEHUNTER) END MATRICES; SPECIFY K PIBD0 PIBD1 PIBD2 BEGIN ALGEBRA; U= H@A+B@D; !IBD 0 COVARIANCE (=NON-QTL COVARIANCE) S= U+H@Q; !IBD 1 COVARIANCE T= U+Q; !IBD 2 COVARIANCE Y= V|U_ U|V_ !IBD 0 MATRIX V|S_ S|V_ !IBD 1 MATRIX V|T_ T|V; !IBD 2 MATRIX END ALGEBRA; MEANS M_M_M; COVARIANCE Y; WEIGHTS K; START 0.5 ALL START 2.8 M 1 1 1 M 1 1 2 M 1 1 3 OPTION NDECIMALS=3 OPTION ITERATIONS=10000 OPTION MULTIPLE OPTION ISSAT OPTION TH=-1 END Summary of VL file data for group 2 PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 Code 1.000 2.000 3.000 4.000 5.000 6.000 Number 165.000 145.000 171.000 163.000 140.000 169.000 Mean 2.849 2.696 2.541 2.862 2.711 2.568 Variance 0.293 0.356 0.272 0.302 0.307 0.268 Minimum 1.800 1.150 1.240 1.800 1.220 1.560 Maximum 4.480 4.360 4.050 4.580 4.280 4.070 Summary of VL file data for group 3 PIBD2 PIBD1 PIBD0 PLT11 PLT12 PLT13 PLT21 Code -3.000 -2.000 -1.000 1.000 2.000 3.000 4.000 Number 230.000 230.000 230.000 193.000 205.000 123.000 190.000 Mean 0.241 0.479 0.279 2.866 2.671 2.477 2.796 Variance 0.134 0.138 0.132 0.294 0.297 0.232 0.306 Minimum 0.000 0.000 0.000 1.580 1.310 1.520 1.480 Maximum 1.000 1.000 1.000 4.410 4.350 4.210 4.590 PLT22 PLT23 Code 5.000 6.000 Number 205.000 122.000 Mean 2.634 2.491 Variance 0.357 0.314 Minimum 1.380 1.400 Maximum 4.780 4.110 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Title Genetic structure - Simplex MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 31 2 32 3 33 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 GROUP NUMBER: 2 G2: Monozygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 31 2 32 3 33 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 GROUP NUMBER: 3 G3: Dizygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX K This is a FULL matrix of order 3 by 1 1 1 -1 2 -2 3 -3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 31 2 32 3 33 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX S This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Y This is a computed FULL matrix of order 18 by 6 It has no free parameters specified MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title Genetic structure - Simplex MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.111 0.130 0.053 2 0.130 0.182 0.075 3 0.053 0.075 0.031 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.031 0.001 -0.003 2 0.001 0.000 0.000 3 -0.003 0.000 0.000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.003 2 0.011 0.062 0.011 3 0.003 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.349 2 0.283 3 0.427 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.505 2 2.879 2.690 2.505 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.173 0.000 3 0.000 0.410 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.024 0.000 3 0.000 -4.406 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.410 0.000 3 0.000 0.249 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.099 0.149 2 0.099 0.080 0.121 3 0.149 0.121 0.183 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.132 2 0.000 0.132 3 0.000 0.000 0.132 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.240 0.202 2 0.240 0.325 0.207 3 0.202 0.207 0.272 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.177 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.333 2 0.000 0.173 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.163 2 0.000 -0.201 3 0.000 0.000 -0.195 GROUP NUMBER: 2 G2: Monozygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.111 0.130 0.053 2 0.130 0.182 0.075 3 0.053 0.075 0.031 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.031 0.001 -0.003 2 0.001 0.000 0.000 3 -0.003 0.000 0.000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.003 2 0.011 0.062 0.011 3 0.003 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.349 2 0.283 3 0.427 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.505 2 2.879 2.690 2.505 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.173 0.000 3 0.000 0.410 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.024 0.000 3 0.000 -4.406 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.410 0.000 3 0.000 0.249 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.099 0.149 2 0.099 0.080 0.121 3 0.149 0.121 0.183 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.132 2 0.000 0.132 3 0.000 0.000 0.132 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.240 0.202 2 0.240 0.325 0.207 3 0.202 0.207 0.272 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.177 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.333 2 0.000 0.173 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.163 2 0.000 -0.201 3 0.000 0.000 -0.195 GROUP NUMBER: 3 G3: Dizygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.111 0.130 0.053 2 0.130 0.182 0.075 3 0.053 0.075 0.031 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.031 0.001 -0.003 2 0.001 0.000 0.000 3 -0.003 0.000 0.000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.003 2 0.011 0.062 0.011 3 0.003 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX K This is a FULL matrix of order 3 by 1 1 1 0.001 2 0.184 3 0.814 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.349 2 0.283 3 0.427 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.505 2 2.879 2.690 2.505 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.173 0.000 3 0.000 0.410 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.024 0.000 3 0.000 -4.406 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.410 0.000 3 0.000 0.249 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.099 0.149 2 0.099 0.080 0.121 3 0.149 0.121 0.183 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.132 2 0.000 0.132 3 0.000 0.000 0.132 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.124 0.115 0.100 2 0.115 0.131 0.098 3 0.100 0.098 0.107 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.185 0.164 0.175 2 0.164 0.171 0.158 3 0.175 0.158 0.198 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+B@D] 1 2 3 1 0.063 0.065 0.026 2 0.065 0.091 0.037 3 0.026 0.037 0.015 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.240 0.202 2 0.240 0.325 0.207 3 0.202 0.207 0.272 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.177 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.333 2 0.000 0.173 3 0.000 0.000 0.000 MATRIX Y This is a computed FULL matrix of order 18 by 6 [=V|U_U|V_V|S_S|V_V|T_T|V] 1 2 3 4 5 6 1 0.308 0.240 0.202 0.063 0.065 0.026 2 0.240 0.325 0.207 0.065 0.091 0.037 3 0.202 0.207 0.272 0.026 0.037 0.015 4 0.063 0.065 0.026 0.308 0.240 0.202 5 0.065 0.091 0.037 0.240 0.325 0.207 6 0.026 0.037 0.015 0.202 0.207 0.272 7 0.308 0.240 0.202 0.124 0.115 0.100 8 0.240 0.325 0.207 0.115 0.131 0.098 9 0.202 0.207 0.272 0.100 0.098 0.107 10 0.124 0.115 0.100 0.308 0.240 0.202 11 0.115 0.131 0.098 0.240 0.325 0.207 12 0.100 0.098 0.107 0.202 0.207 0.272 13 0.308 0.240 0.202 0.185 0.164 0.175 14 0.240 0.325 0.207 0.164 0.171 0.158 15 0.202 0.207 0.272 0.175 0.158 0.198 16 0.185 0.164 0.175 0.308 0.240 0.202 17 0.164 0.171 0.158 0.240 0.325 0.207 18 0.175 0.158 0.198 0.202 0.207 0.272 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.163 2 0.000 -0.201 3 0.000 0.000 -0.195 Your model has 22 estimated parameters and 1991 Observed statistics -2 times log-likelihood of data >>> 1979.020 Degrees of freedom >>>>>>>>>>>>>>>> 1969 Now I'm trying to improve on the current solution for you... *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title Genetic structure - Simplex MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.111 0.130 0.053 2 0.130 0.182 0.075 3 0.053 0.075 0.031 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.031 0.001 -0.003 2 0.001 0.000 0.000 3 -0.003 0.000 0.000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.003 2 0.011 0.062 0.011 3 0.003 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.349 2 0.283 3 0.427 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.505 2 2.879 2.690 2.505 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.173 0.000 3 0.000 0.410 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.024 0.000 3 0.000 -4.406 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.410 0.000 3 0.000 0.249 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.099 0.149 2 0.099 0.080 0.121 3 0.149 0.121 0.183 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.132 2 0.000 0.132 3 0.000 0.000 0.132 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.240 0.202 2 0.240 0.325 0.207 3 0.202 0.207 0.272 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.177 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.333 2 0.000 0.173 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.163 2 0.000 -0.201 3 0.000 0.000 -0.195 GROUP NUMBER: 2 G2: Monozygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.111 0.130 0.053 2 0.130 0.182 0.075 3 0.053 0.075 0.031 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.031 0.001 -0.003 2 0.001 0.000 0.000 3 -0.003 0.000 0.000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.003 2 0.011 0.062 0.011 3 0.003 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.349 2 0.283 3 0.427 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.505 2 2.879 2.690 2.505 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.173 0.000 3 0.000 0.410 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.024 0.000 3 0.000 -4.406 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.410 0.000 3 0.000 0.249 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.099 0.149 2 0.099 0.080 0.121 3 0.149 0.121 0.183 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.132 2 0.000 0.132 3 0.000 0.000 0.132 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.240 0.202 2 0.240 0.325 0.207 3 0.202 0.207 0.272 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.177 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.333 2 0.000 0.173 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.163 2 0.000 -0.201 3 0.000 0.000 -0.195 GROUP NUMBER: 3 G3: Dizygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.111 0.130 0.053 2 0.130 0.182 0.075 3 0.053 0.075 0.031 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.031 0.001 -0.003 2 0.001 0.000 0.000 3 -0.003 0.000 0.000 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.003 2 0.011 0.062 0.011 3 0.003 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX K This is a FULL matrix of order 3 by 1 1 1 0.001 2 0.184 3 0.814 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.349 2 0.283 3 0.427 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.505 2 2.879 2.690 2.505 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.173 0.000 3 0.000 0.410 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.024 0.000 3 0.000 -4.406 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.410 0.000 3 0.000 0.249 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.099 0.149 2 0.099 0.080 0.121 3 0.149 0.121 0.183 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.132 2 0.000 0.132 3 0.000 0.000 0.132 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.124 0.115 0.100 2 0.115 0.131 0.098 3 0.100 0.098 0.107 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.185 0.164 0.175 2 0.164 0.171 0.158 3 0.175 0.158 0.198 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+B@D] 1 2 3 1 0.063 0.065 0.026 2 0.065 0.091 0.037 3 0.026 0.037 0.015 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.240 0.202 2 0.240 0.325 0.207 3 0.202 0.207 0.272 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.177 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.333 2 0.000 0.173 3 0.000 0.000 0.000 MATRIX Y This is a computed FULL matrix of order 18 by 6 [=V|U_U|V_V|S_S|V_V|T_T|V] 1 2 3 4 5 6 1 0.308 0.240 0.202 0.063 0.065 0.026 2 0.240 0.325 0.207 0.065 0.091 0.037 3 0.202 0.207 0.272 0.026 0.037 0.015 4 0.063 0.065 0.026 0.308 0.240 0.202 5 0.065 0.091 0.037 0.240 0.325 0.207 6 0.026 0.037 0.015 0.202 0.207 0.272 7 0.308 0.240 0.202 0.124 0.115 0.100 8 0.240 0.325 0.207 0.115 0.131 0.098 9 0.202 0.207 0.272 0.100 0.098 0.107 10 0.124 0.115 0.100 0.308 0.240 0.202 11 0.115 0.131 0.098 0.240 0.325 0.207 12 0.100 0.098 0.107 0.202 0.207 0.272 13 0.308 0.240 0.202 0.185 0.164 0.175 14 0.240 0.325 0.207 0.164 0.171 0.158 15 0.202 0.207 0.272 0.175 0.158 0.198 16 0.185 0.164 0.175 0.308 0.240 0.202 17 0.164 0.171 0.158 0.240 0.325 0.207 18 0.175 0.158 0.198 0.202 0.207 0.272 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.163 2 0.000 -0.201 3 0.000 0.000 -0.195 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 22 estimated parameters and 1991 Observed statistics -2 times log-likelihood of data >>> 1979.020 Degrees of freedom >>>>>>>>>>>>>>>> 1969 This problem used 1.5% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.63 Execution 0: 0: 1:14.27 TOTAL 0: 0: 1:14.90 Total number of warnings issued: 2 ______________________________________________________________________________ ! Equate loadings on QTL Multiple fit option in effect. The following MX script lines have been read: OPTION TH=-1 EQUATE L 1 1 1 L 1 2 1 L 1 3 1 END Summary of VL file data for group 2 PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 Code 1.000 2.000 3.000 4.000 5.000 6.000 Number 165.000 145.000 171.000 163.000 140.000 169.000 Mean 2.849 2.696 2.541 2.862 2.711 2.568 Variance 0.293 0.356 0.272 0.302 0.307 0.268 Minimum 1.800 1.150 1.240 1.800 1.220 1.560 Maximum 4.480 4.360 4.050 4.580 4.280 4.070 Summary of VL file data for group 3 PIBD2 PIBD1 PIBD0 PLT11 PLT12 PLT13 PLT21 Code -3.000 -2.000 -1.000 1.000 2.000 3.000 4.000 Number 230.000 230.000 230.000 193.000 205.000 123.000 190.000 Mean 0.241 0.479 0.279 2.866 2.671 2.477 2.796 Variance 0.134 0.138 0.132 0.294 0.297 0.232 0.306 Minimum 0.000 0.000 0.000 1.580 1.310 1.520 1.480 Maximum 1.000 1.000 1.000 4.410 4.350 4.210 4.590 PLT22 PLT23 Code 5.000 6.000 Number 205.000 122.000 Mean 2.634 2.491 Variance 0.357 0.314 Minimum 1.380 1.400 Maximum 4.780 4.110 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Title Genetic structure - Simplex MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 31 2 31 3 31 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 GROUP NUMBER: 2 G2: Monozygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 31 2 31 3 31 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 GROUP NUMBER: 3 G3: Dizygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX K This is a FULL matrix of order 3 by 1 1 1 -1 2 -2 3 -3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 31 2 31 3 31 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX S This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Y This is a computed FULL matrix of order 18 by 6 It has no free parameters specified MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 6 Looks like I got stuck here. Check the following: 1. The model is correctly specified 2. Starting values are good 3. You are not already at the solution The error can arise if the Hessian is ill-conditioned You can try resetting it to an identity matrix and fit from the solution by putting TH=-n on the OU line where n is the number of refits that you want to do If all else fails try putting NAG=30 on the OU line and examine the file NAGDUMP.OUT and the NAG manual MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title Genetic structure - Simplex MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.082 0.091 0.041 2 0.091 0.125 0.056 3 0.041 0.056 0.065 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.042 0.000 0.017 2 0.000 0.000 0.000 3 0.017 0.000 0.007 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.373 2 0.373 3 0.373 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.508 2 2.880 2.691 2.508 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.110 0.000 3 0.000 0.448 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.002 0.000 3 0.000 -191.655 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.349 0.000 3 0.000 0.216 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.139 0.139 0.139 2 0.139 0.139 0.139 3 0.139 0.139 0.139 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.109 2 0.000 -0.109 3 0.000 0.000 -0.109 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.241 0.200 2 0.241 0.327 0.206 3 0.200 0.206 0.269 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.205 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.287 2 0.000 0.156 3 0.000 0.000 0.200 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.179 2 0.000 -0.216 3 0.000 0.000 -0.209 GROUP NUMBER: 2 G2: Monozygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.082 0.091 0.041 2 0.091 0.125 0.056 3 0.041 0.056 0.065 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.042 0.000 0.017 2 0.000 0.000 0.000 3 0.017 0.000 0.007 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.373 2 0.373 3 0.373 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.508 2 2.880 2.691 2.508 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.110 0.000 3 0.000 0.448 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.002 0.000 3 0.000 -191.655 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.349 0.000 3 0.000 0.216 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.139 0.139 0.139 2 0.139 0.139 0.139 3 0.139 0.139 0.139 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.109 2 0.000 -0.109 3 0.000 0.000 -0.109 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.241 0.200 2 0.241 0.327 0.206 3 0.200 0.206 0.269 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.205 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.287 2 0.000 0.156 3 0.000 0.000 0.200 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.179 2 0.000 -0.216 3 0.000 0.000 -0.209 GROUP NUMBER: 3 G3: Dizygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.082 0.091 0.041 2 0.091 0.125 0.056 3 0.041 0.056 0.065 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.042 0.000 0.017 2 0.000 0.000 0.000 3 0.017 0.000 0.007 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX K This is a FULL matrix of order 3 by 1 1 1 0.001 2 0.184 3 0.814 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.373 2 0.373 3 0.373 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.508 2 2.880 2.691 2.508 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.110 0.000 3 0.000 0.448 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.002 0.000 3 0.000 -191.655 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.349 0.000 3 0.000 0.216 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.139 0.139 0.139 2 0.139 0.139 0.139 3 0.139 0.139 0.139 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.109 2 0.000 -0.109 3 0.000 0.000 -0.109 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.121 0.115 0.094 2 0.115 0.132 0.098 3 0.094 0.098 0.104 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.191 0.185 0.164 2 0.185 0.202 0.167 3 0.164 0.167 0.173 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+B@D] 1 2 3 1 0.052 0.046 0.025 2 0.046 0.063 0.028 3 0.025 0.028 0.034 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.241 0.200 2 0.241 0.327 0.206 3 0.200 0.206 0.269 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.205 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.287 2 0.000 0.156 3 0.000 0.000 0.200 MATRIX Y This is a computed FULL matrix of order 18 by 6 [=V|U_U|V_V|S_S|V_V|T_T|V] 1 2 3 4 5 6 1 0.308 0.241 0.200 0.052 0.046 0.025 2 0.241 0.327 0.206 0.046 0.063 0.028 3 0.200 0.206 0.269 0.025 0.028 0.034 4 0.052 0.046 0.025 0.308 0.241 0.200 5 0.046 0.063 0.028 0.241 0.327 0.206 6 0.025 0.028 0.034 0.200 0.206 0.269 7 0.308 0.241 0.200 0.121 0.115 0.094 8 0.241 0.327 0.206 0.115 0.132 0.098 9 0.200 0.206 0.269 0.094 0.098 0.104 10 0.121 0.115 0.094 0.308 0.241 0.200 11 0.115 0.132 0.098 0.241 0.327 0.206 12 0.094 0.098 0.104 0.200 0.206 0.269 13 0.308 0.241 0.200 0.191 0.185 0.164 14 0.241 0.327 0.206 0.185 0.202 0.167 15 0.200 0.206 0.269 0.164 0.167 0.173 16 0.191 0.185 0.164 0.308 0.241 0.200 17 0.185 0.202 0.167 0.241 0.327 0.206 18 0.164 0.167 0.173 0.200 0.206 0.269 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.179 2 0.000 -0.216 3 0.000 0.000 -0.209 *** WARNING! *** Minimization may not be successful. See above CODE RED - Hessian/precision problem Your model has 20 estimated parameters and 1991 Observed statistics -2 times log-likelihood of data >>> 1982.347 Degrees of freedom >>>>>>>>>>>>>>>> 1971 Saturated model fit* >>>>>>>>>>> 1979.020 Saturated model df* >>>>>>>>>>> 1969 Difference Chi-squared >>>>>>>> 3.327 Difference d.f. >>>>>>>>>>>>>>> 2 Probability >>>>>>>>>>>>>>>>>>>> 0.189 Akaike's Information Criterion > -0.673 * Saturated model statistic computed earlier in this job Now I'm trying to improve on the current solution for you... MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title Genetic structure - Simplex MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.082 0.091 0.041 2 0.091 0.125 0.056 3 0.041 0.056 0.065 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.042 0.000 0.017 2 0.000 0.000 0.000 3 0.017 0.000 0.007 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.373 2 0.373 3 0.373 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.508 2 2.880 2.691 2.508 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.110 0.000 3 0.000 0.448 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.002 0.000 3 0.000 -191.655 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.349 0.000 3 0.000 0.216 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.139 0.139 0.139 2 0.139 0.139 0.139 3 0.139 0.139 0.139 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.109 2 0.000 -0.109 3 0.000 0.000 -0.109 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.241 0.200 2 0.241 0.327 0.206 3 0.200 0.206 0.269 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.205 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.287 2 0.000 0.156 3 0.000 0.000 0.200 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.179 2 0.000 -0.216 3 0.000 0.000 -0.209 GROUP NUMBER: 2 G2: Monozygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.082 0.091 0.041 2 0.091 0.125 0.056 3 0.041 0.056 0.065 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.042 0.000 0.017 2 0.000 0.000 0.000 3 0.017 0.000 0.007 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.373 2 0.373 3 0.373 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.508 2 2.880 2.691 2.508 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.110 0.000 3 0.000 0.448 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.002 0.000 3 0.000 -191.655 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.349 0.000 3 0.000 0.216 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.139 0.139 0.139 2 0.139 0.139 0.139 3 0.139 0.139 0.139 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.109 2 0.000 -0.109 3 0.000 0.000 -0.109 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.241 0.200 2 0.241 0.327 0.206 3 0.200 0.206 0.269 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.205 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.287 2 0.000 0.156 3 0.000 0.000 0.200 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.179 2 0.000 -0.216 3 0.000 0.000 -0.209 GROUP NUMBER: 3 G3: Dizygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.082 0.091 0.041 2 0.091 0.125 0.056 3 0.041 0.056 0.065 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.042 0.000 0.017 2 0.000 0.000 0.000 3 0.017 0.000 0.007 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX K This is a FULL matrix of order 3 by 1 1 1 0.001 2 0.184 3 0.814 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.373 2 0.373 3 0.373 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.508 2 2.880 2.691 2.508 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.110 0.000 3 0.000 0.448 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.002 0.000 3 0.000 -191.655 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.349 0.000 3 0.000 0.216 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.139 0.139 0.139 2 0.139 0.139 0.139 3 0.139 0.139 0.139 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.109 2 0.000 -0.109 3 0.000 0.000 -0.109 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.121 0.115 0.094 2 0.115 0.132 0.098 3 0.094 0.098 0.104 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.191 0.185 0.164 2 0.185 0.202 0.167 3 0.164 0.167 0.173 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+B@D] 1 2 3 1 0.052 0.046 0.025 2 0.046 0.063 0.028 3 0.025 0.028 0.034 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.308 0.241 0.200 2 0.241 0.327 0.206 3 0.200 0.206 0.269 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.205 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.287 2 0.000 0.156 3 0.000 0.000 0.200 MATRIX Y This is a computed FULL matrix of order 18 by 6 [=V|U_U|V_V|S_S|V_V|T_T|V] 1 2 3 4 5 6 1 0.308 0.241 0.200 0.052 0.046 0.025 2 0.241 0.327 0.206 0.046 0.063 0.028 3 0.200 0.206 0.269 0.025 0.028 0.034 4 0.052 0.046 0.025 0.308 0.241 0.200 5 0.046 0.063 0.028 0.241 0.327 0.206 6 0.025 0.028 0.034 0.200 0.206 0.269 7 0.308 0.241 0.200 0.121 0.115 0.094 8 0.241 0.327 0.206 0.115 0.132 0.098 9 0.200 0.206 0.269 0.094 0.098 0.104 10 0.121 0.115 0.094 0.308 0.241 0.200 11 0.115 0.132 0.098 0.241 0.327 0.206 12 0.094 0.098 0.104 0.200 0.206 0.269 13 0.308 0.241 0.200 0.191 0.185 0.164 14 0.241 0.327 0.206 0.185 0.202 0.167 15 0.200 0.206 0.269 0.164 0.167 0.173 16 0.191 0.185 0.164 0.308 0.241 0.200 17 0.185 0.202 0.167 0.241 0.327 0.206 18 0.164 0.167 0.173 0.200 0.206 0.269 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.179 2 0.000 -0.216 3 0.000 0.000 -0.209 Your model has 20 estimated parameters and 1991 Observed statistics -2 times log-likelihood of data >>> 1982.347 Degrees of freedom >>>>>>>>>>>>>>>> 1971 Saturated model fit* >>>>>>>>>>> 1979.020 Saturated model df* >>>>>>>>>>> 1969 Difference Chi-squared >>>>>>>> 3.327 Difference d.f. >>>>>>>>>>>>>>> 2 Probability >>>>>>>>>>>>>>>>>>>> 0.189 Akaike's Information Criterion > -0.673 * Saturated model statistic computed earlier in this job This problem used 1.5% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.00 Execution 0: 0: 2: 7.57 TOTAL 0: 0: 2: 7.58 Total number of warnings issued: 4 ______________________________________________________________________________ ! Test significance of QTL effect Multiple fit option in effect. The following MX script lines have been read: DROP L 1 1 1 TO L 1 3 1 END Summary of VL file data for group 2 PLT11 PLT12 PLT13 PLT21 PLT22 PLT23 Code 1.000 2.000 3.000 4.000 5.000 6.000 Number 165.000 145.000 171.000 163.000 140.000 169.000 Mean 2.849 2.696 2.541 2.862 2.711 2.568 Variance 0.293 0.356 0.272 0.302 0.307 0.268 Minimum 1.800 1.150 1.240 1.800 1.220 1.560 Maximum 4.480 4.360 4.050 4.580 4.280 4.070 Summary of VL file data for group 3 PIBD2 PIBD1 PIBD0 PLT11 PLT12 PLT13 PLT21 Code -3.000 -2.000 -1.000 1.000 2.000 3.000 4.000 Number 230.000 230.000 230.000 193.000 205.000 123.000 190.000 Mean 0.241 0.479 0.279 2.866 2.671 2.477 2.796 Variance 0.134 0.138 0.132 0.294 0.297 0.232 0.306 Minimum 0.000 0.000 0.000 1.580 1.310 1.520 1.480 Maximum 1.000 1.000 1.000 4.410 4.350 4.210 4.590 PLT22 PLT23 Code 5.000 6.000 Number 205.000 122.000 Mean 2.634 2.491 Variance 0.357 0.314 Minimum 1.380 1.400 Maximum 4.780 4.110 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Title Genetic structure - Simplex MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 GROUP NUMBER: 2 G2: Monozygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 GROUP NUMBER: 3 G3: Dizygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX K This is a FULL matrix of order 3 by 1 1 1 -1 2 -2 3 -3 MATRIX L This is a FULL matrix of order 3 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 40 41 42 2 40 41 42 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 13 0 3 0 14 0 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 15 0 3 0 16 0 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0 2 17 0 3 0 18 0 MATRIX Q This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 21 2 0 21 3 0 0 21 MATRIX S This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX T This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX V This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 1 2 0 2 3 0 0 3 MATRIX Y This is a computed FULL matrix of order 18 by 6 It has no free parameters specified MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 7 2 0 8 3 0 0 9 MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title Genetic structure - Simplex MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.215 0.229 0.171 2 0.229 0.262 0.195 3 0.171 0.195 0.198 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.048 0.000 0.027 2 0.000 0.000 0.000 3 0.027 0.000 0.015 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.000 2 0.000 3 0.000 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.689 2.507 2 2.879 2.689 2.507 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.064 0.000 3 0.000 0.746 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.003 0.000 3 0.000 -191.654 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.332 0.000 3 0.000 0.204 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.000 0.000 0.000 2 0.000 0.000 0.000 3 0.000 0.000 0.000 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.101 2 0.000 -0.101 3 0.000 0.000 -0.101 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.307 0.240 0.200 2 0.240 0.324 0.206 3 0.200 0.206 0.271 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.220 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.464 2 0.000 0.135 3 0.000 0.000 0.228 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.183 2 0.000 -0.220 3 0.000 0.000 -0.213 GROUP NUMBER: 2 G2: Monozygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.215 0.229 0.171 2 0.229 0.262 0.195 3 0.171 0.195 0.198 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.048 0.000 0.027 2 0.000 0.000 0.000 3 0.027 0.000 0.015 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.000 2 0.000 3 0.000 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.689 2.507 2 2.879 2.689 2.507 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.064 0.000 3 0.000 0.746 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.003 0.000 3 0.000 -191.654 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.332 0.000 3 0.000 0.204 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.000 0.000 0.000 2 0.000 0.000 0.000 3 0.000 0.000 0.000 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.101 2 0.000 -0.101 3 0.000 0.000 -0.101 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.307 0.240 0.200 2 0.240 0.324 0.206 3 0.200 0.206 0.271 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.220 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.464 2 0.000 0.135 3 0.000 0.000 0.228 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.183 2 0.000 -0.220 3 0.000 0.000 -0.213 GROUP NUMBER: 3 G3: Dizygotic twins MATRIX A This is a computed FULL matrix of order 3 by 3 [=(I-N)~&(X*X')] 1 2 3 1 0.215 0.229 0.171 2 0.229 0.262 0.195 3 0.171 0.195 0.198 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.250 MATRIX D This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(W*W')] 1 2 3 1 0.048 0.000 0.027 2 0.000 0.000 0.000 3 0.027 0.000 0.015 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.044 0.011 0.002 2 0.011 0.062 0.011 3 0.002 0.011 0.058 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.500 MATRIX I This is an IDENTITY matrix of order 3 by 3 MATRIX K This is a FULL matrix of order 3 by 1 1 1 0.001 2 0.184 3 0.814 MATRIX L This is a FULL matrix of order 3 by 1 1 1 0.000 2 0.000 3 0.000 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.689 2.507 2 2.879 2.689 2.507 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 1.064 0.000 3 0.000 0.746 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 -0.003 0.000 3 0.000 -191.654 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.332 0.000 3 0.000 0.204 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.000 0.000 0.000 2 0.000 0.000 0.000 3 0.000 0.000 0.000 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.101 2 0.000 -0.101 3 0.000 0.000 -0.101 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.120 0.114 0.092 2 0.114 0.131 0.098 3 0.092 0.098 0.103 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.120 0.114 0.092 2 0.114 0.131 0.098 3 0.092 0.098 0.103 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+B@D] 1 2 3 1 0.120 0.114 0.092 2 0.114 0.131 0.098 3 0.092 0.098 0.103 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+D+E] 1 2 3 1 0.307 0.240 0.200 2 0.240 0.324 0.206 3 0.200 0.206 0.271 MATRIX W This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.220 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.464 2 0.000 0.135 3 0.000 0.000 0.228 MATRIX Y This is a computed FULL matrix of order 18 by 6 [=V|U_U|V_V|S_S|V_V|T_T|V] 1 2 3 4 5 6 1 0.307 0.240 0.200 0.120 0.114 0.092 2 0.240 0.324 0.206 0.114 0.131 0.098 3 0.200 0.206 0.271 0.092 0.098 0.103 4 0.120 0.114 0.092 0.307 0.240 0.200 5 0.114 0.131 0.098 0.240 0.324 0.206 6 0.092 0.098 0.103 0.200 0.206 0.271 7 0.307 0.240 0.200 0.120 0.114 0.092 8 0.240 0.324 0.206 0.114 0.131 0.098 9 0.200 0.206 0.271 0.092 0.098 0.103 10 0.120 0.114 0.092 0.307 0.240 0.200 11 0.114 0.131 0.098 0.240 0.324 0.206 12 0.092 0.098 0.103 0.200 0.206 0.271 13 0.307 0.240 0.200 0.120 0.114 0.092 14 0.240 0.324 0.206 0.114 0.131 0.098 15 0.200 0.206 0.271 0.092 0.098 0.103 16 0.120 0.114 0.092 0.307 0.240 0.200 17 0.114 0.131 0.098 0.240 0.324 0.206 18 0.092 0.098 0.103 0.200 0.206 0.271 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.183 2 0.000 -0.220 3 0.000 0.000 -0.213 Your model has 19 estimated parameters and 1991 Observed statistics -2 times log-likelihood of data >>> 1994.522 Degrees of freedom >>>>>>>>>>>>>>>> 1972 Saturated model fit* >>>>>>>>>>> 1979.020 Saturated model df* >>>>>>>>>>> 1969 Difference Chi-squared >>>>>>>> 15.503 Difference d.f. >>>>>>>>>>>>>>> 3 Probability >>>>>>>>>>>>>>>>>>>> 0.001 Akaike's Information Criterion > 9.503 * Saturated model statistic computed earlier in this job This problem used 1.5% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.00 Execution 0: 0: 0:17.27 TOTAL 0: 0: 0:17.27 Total number of warnings issued: 4 ______________________________________________________________________________ ______________________________________________________________________________