!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 Y DIAG NVAR NVAR FREE !COMMON ENVIRONMENT Z DIAG NVAR NVAR FREE !SPECIFIC ENVIRONMENTAL INNOVATIONS N LOWER NVAR NVAR FREE !GENETIC TRANSMISSION PATHS O LOWER NVAR NVAR FREE !COMMON ENVIRONMENTAL 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 END MATRICES; MATRIX H .5 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'); C= (I-O)~ & (Y*Y'); E= (I-P)~ & (Z*Z') + R*R'; Q= L*L'; !VARIANCE DUE TO QTL V= A+Q+C+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 Rectangular continuous data read initiated NOTE: Rectangular file contained 519 records with data that contained a total of 5624 observations 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+C+Q+E | A+C+Q _ A+C+Q | A+C+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 Rectangular continuous data read initiated NOTE: Rectangular file contained 519 records with data that contained a total of 5624 observations 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+C; !IBD 0 COVARIANCE (=NON-QTL COVARIANCE) S= U+H@Q; !IBD 1 COVARIANCE T= U+Q; !IBD 2 COVARIANCE D= 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 D; 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 MULTIPLE OPTION ISSAT 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 C 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 C 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 18 by 6 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 starting optimization; number of parameters = 22 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.142 0.134 0.054 2 0.134 0.179 0.072 3 0.054 0.072 0.029 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.000 0.000 0.000 2 0.000 0.007 0.007 3 0.000 0.007 0.008 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.002 2 0.010 0.063 0.012 3 0.002 0.012 0.059 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.350 2 0.276 3 0.420 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.506 2 2.879 2.690 2.506 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.945 0.000 3 0.000 0.402 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.639 0.000 3 0.000 1.106 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.222 0.000 3 0.000 0.185 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.097 0.147 2 0.097 0.076 0.116 3 0.147 0.116 0.176 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.000 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.310 0.241 0.203 2 0.241 0.325 0.207 3 0.203 0.207 0.272 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.377 2 0.000 0.229 3 0.000 0.000 0.000 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.000 2 0.000 0.081 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.212 2 0.000 0.246 3 0.000 0.000 0.238 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.142 0.134 0.054 2 0.134 0.179 0.072 3 0.054 0.072 0.029 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.000 0.000 0.000 2 0.000 0.007 0.007 3 0.000 0.007 0.008 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.002 2 0.010 0.063 0.012 3 0.002 0.012 0.059 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.350 2 0.276 3 0.420 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.506 2 2.879 2.690 2.506 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.945 0.000 3 0.000 0.402 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.639 0.000 3 0.000 1.106 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.222 0.000 3 0.000 0.185 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.097 0.147 2 0.097 0.076 0.116 3 0.147 0.116 0.176 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.000 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.310 0.241 0.203 2 0.241 0.325 0.207 3 0.203 0.207 0.272 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.377 2 0.000 0.229 3 0.000 0.000 0.000 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.000 2 0.000 0.081 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.212 2 0.000 0.246 3 0.000 0.000 0.238 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.142 0.134 0.054 2 0.134 0.179 0.072 3 0.054 0.072 0.029 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.000 0.000 0.000 2 0.000 0.007 0.007 3 0.000 0.007 0.008 MATRIX D 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.310 0.241 0.203 0.071 0.067 0.027 2 0.241 0.325 0.207 0.067 0.096 0.043 3 0.203 0.207 0.272 0.027 0.043 0.022 4 0.071 0.067 0.027 0.310 0.241 0.203 5 0.067 0.096 0.043 0.241 0.325 0.207 6 0.027 0.043 0.022 0.203 0.207 0.272 7 0.310 0.241 0.203 0.132 0.116 0.100 8 0.241 0.325 0.207 0.116 0.134 0.101 9 0.203 0.207 0.272 0.100 0.101 0.111 10 0.132 0.116 0.100 0.310 0.241 0.203 11 0.116 0.134 0.101 0.241 0.325 0.207 12 0.100 0.101 0.111 0.203 0.207 0.272 13 0.310 0.241 0.203 0.193 0.164 0.174 14 0.241 0.325 0.207 0.164 0.173 0.159 15 0.203 0.207 0.272 0.174 0.159 0.199 16 0.193 0.164 0.174 0.310 0.241 0.203 17 0.164 0.173 0.159 0.241 0.325 0.207 18 0.174 0.159 0.199 0.203 0.207 0.272 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.002 2 0.010 0.063 0.012 3 0.002 0.012 0.059 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.350 2 0.276 3 0.420 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.879 2.690 2.506 2 2.879 2.690 2.506 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.945 0.000 3 0.000 0.402 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.639 0.000 3 0.000 1.106 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.222 0.000 3 0.000 0.185 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.122 0.097 0.147 2 0.097 0.076 0.116 3 0.147 0.116 0.176 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.000 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.132 0.116 0.100 2 0.116 0.134 0.101 3 0.100 0.101 0.111 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.193 0.164 0.174 2 0.164 0.173 0.159 3 0.174 0.159 0.199 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+C] 1 2 3 1 0.071 0.067 0.027 2 0.067 0.096 0.043 3 0.027 0.043 0.022 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.310 0.241 0.203 2 0.241 0.325 0.207 3 0.203 0.207 0.272 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.377 2 0.000 0.229 3 0.000 0.000 0.000 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.000 2 0.000 0.081 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.212 2 0.000 0.246 3 0.000 0.000 0.238 Your model has 22 estimated parameters and 1991 Observed statistics -2 times log-likelihood of data >>> 1979.626 Degrees of freedom >>>>>>>>>>>>>>>> 1969 Model number: 1 Model number: 2 Model number: 3 Model number: 1 Model number: 2 Model number: 3 This problem used 14.7% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 2.60 Execution 0: 0: 0:21.66 TOTAL 0: 0: 0:24.26 Total number of warnings issued: 0 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: OPTION TH=-1 EQ 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 C 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 C 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 18 by 6 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 starting optimization; number of parameters = 20 *** 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.117 0.092 0.043 2 0.092 0.124 0.058 3 0.043 0.058 0.050 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.005 0.000 0.009 2 0.000 0.000 0.000 3 0.009 0.000 0.019 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.004 2 0.010 0.064 0.009 3 0.004 0.009 0.059 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.378 2 0.378 3 0.378 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.510 2 2.880 2.691 2.510 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.785 0.000 3 0.000 0.468 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.007 0.000 3 0.000 303.978 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.325 0.000 3 0.000 0.379 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.143 0.143 0.143 2 0.143 0.143 0.143 3 0.143 0.143 0.143 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.203 2 0.000 0.203 3 0.000 0.000 0.203 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.309 0.244 0.199 2 0.244 0.330 0.209 3 0.199 0.209 0.270 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.342 2 0.000 0.229 3 0.000 0.000 0.150 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.068 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.065 2 0.000 0.000 3 0.000 0.000 0.121 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.117 0.092 0.043 2 0.092 0.124 0.058 3 0.043 0.058 0.050 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.005 0.000 0.009 2 0.000 0.000 0.000 3 0.009 0.000 0.019 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.004 2 0.010 0.064 0.009 3 0.004 0.009 0.059 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.378 2 0.378 3 0.378 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.510 2 2.880 2.691 2.510 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.785 0.000 3 0.000 0.468 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.007 0.000 3 0.000 303.978 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.325 0.000 3 0.000 0.379 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.143 0.143 0.143 2 0.143 0.143 0.143 3 0.143 0.143 0.143 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.203 2 0.000 0.203 3 0.000 0.000 0.203 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.309 0.244 0.199 2 0.244 0.330 0.209 3 0.199 0.209 0.270 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.342 2 0.000 0.229 3 0.000 0.000 0.150 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.068 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.065 2 0.000 0.000 3 0.000 0.000 0.121 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.117 0.092 0.043 2 0.092 0.124 0.058 3 0.043 0.058 0.050 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.005 0.000 0.009 2 0.000 0.000 0.000 3 0.009 0.000 0.019 MATRIX D 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.309 0.244 0.199 0.063 0.046 0.031 2 0.244 0.330 0.209 0.046 0.062 0.029 3 0.199 0.209 0.270 0.031 0.029 0.044 4 0.063 0.046 0.031 0.309 0.244 0.199 5 0.046 0.062 0.029 0.244 0.330 0.209 6 0.031 0.029 0.044 0.199 0.209 0.270 7 0.309 0.244 0.199 0.134 0.117 0.102 8 0.244 0.330 0.209 0.117 0.133 0.100 9 0.199 0.209 0.270 0.102 0.100 0.115 10 0.134 0.117 0.102 0.309 0.244 0.199 11 0.117 0.133 0.100 0.244 0.330 0.209 12 0.102 0.100 0.115 0.199 0.209 0.270 13 0.309 0.244 0.199 0.206 0.188 0.173 14 0.244 0.330 0.209 0.188 0.205 0.172 15 0.199 0.209 0.270 0.173 0.172 0.186 16 0.206 0.188 0.173 0.309 0.244 0.199 17 0.188 0.205 0.172 0.244 0.330 0.209 18 0.173 0.172 0.186 0.199 0.209 0.270 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.004 2 0.010 0.064 0.009 3 0.004 0.009 0.059 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.378 2 0.378 3 0.378 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.510 2 2.880 2.691 2.510 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.785 0.000 3 0.000 0.468 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.007 0.000 3 0.000 303.978 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.325 0.000 3 0.000 0.379 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.143 0.143 0.143 2 0.143 0.143 0.143 3 0.143 0.143 0.143 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.203 2 0.000 0.203 3 0.000 0.000 0.203 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.134 0.117 0.102 2 0.117 0.133 0.100 3 0.102 0.100 0.115 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.206 0.188 0.173 2 0.188 0.205 0.172 3 0.173 0.172 0.186 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+C] 1 2 3 1 0.063 0.046 0.031 2 0.046 0.062 0.029 3 0.031 0.029 0.044 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.309 0.244 0.199 2 0.244 0.330 0.209 3 0.199 0.209 0.270 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.342 2 0.000 0.229 3 0.000 0.000 0.150 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.068 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.065 2 0.000 0.000 3 0.000 0.000 0.121 *** 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 >>> 1984.982 Degrees of freedom >>>>>>>>>>>>>>>> 1971 Saturated model fit* >>>>>>>>>>> 1979.626 Saturated model df* >>>>>>>>>>> 1969 Difference Chi-squared >>>>>>>> 5.356 Difference d.f. >>>>>>>>>>>>>>> 2 Probability >>>>>>>>>>>>>>>>>>>> 0.069 Akaike's Information Criterion > 1.356 * Saturated model statistic computed earlier in this job Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 20 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.117 0.092 0.043 2 0.092 0.124 0.058 3 0.043 0.058 0.050 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.005 0.000 0.009 2 0.000 0.000 0.000 3 0.009 0.000 0.019 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.004 2 0.010 0.064 0.009 3 0.004 0.009 0.059 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.378 2 0.378 3 0.378 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.510 2 2.880 2.691 2.510 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.785 0.000 3 0.000 0.468 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.007 0.000 3 0.000 303.978 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.325 0.000 3 0.000 0.379 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.143 0.143 0.143 2 0.143 0.143 0.143 3 0.143 0.143 0.143 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.203 2 0.000 0.203 3 0.000 0.000 0.203 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.309 0.244 0.199 2 0.244 0.330 0.209 3 0.199 0.209 0.270 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.342 2 0.000 0.229 3 0.000 0.000 0.150 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.068 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.065 2 0.000 0.000 3 0.000 0.000 0.121 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.117 0.092 0.043 2 0.092 0.124 0.058 3 0.043 0.058 0.050 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.005 0.000 0.009 2 0.000 0.000 0.000 3 0.009 0.000 0.019 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.004 2 0.010 0.064 0.009 3 0.004 0.009 0.059 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.378 2 0.378 3 0.378 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.510 2 2.880 2.691 2.510 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.785 0.000 3 0.000 0.468 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.007 0.000 3 0.000 303.978 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.325 0.000 3 0.000 0.379 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.143 0.143 0.143 2 0.143 0.143 0.143 3 0.143 0.143 0.143 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.203 2 0.000 0.203 3 0.000 0.000 0.203 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.309 0.244 0.199 2 0.244 0.330 0.209 3 0.199 0.209 0.270 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.342 2 0.000 0.229 3 0.000 0.000 0.150 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.068 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.065 2 0.000 0.000 3 0.000 0.000 0.121 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.117 0.092 0.043 2 0.092 0.124 0.058 3 0.043 0.058 0.050 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.005 0.000 0.009 2 0.000 0.000 0.000 3 0.009 0.000 0.019 MATRIX D 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.309 0.244 0.199 0.063 0.046 0.031 2 0.244 0.330 0.209 0.046 0.062 0.029 3 0.199 0.209 0.270 0.031 0.029 0.044 4 0.063 0.046 0.031 0.309 0.244 0.199 5 0.046 0.062 0.029 0.244 0.330 0.209 6 0.031 0.029 0.044 0.199 0.209 0.270 7 0.309 0.244 0.199 0.134 0.117 0.102 8 0.244 0.330 0.209 0.117 0.133 0.100 9 0.199 0.209 0.270 0.102 0.100 0.115 10 0.134 0.117 0.102 0.309 0.244 0.199 11 0.117 0.133 0.100 0.244 0.330 0.209 12 0.102 0.100 0.115 0.199 0.209 0.270 13 0.309 0.244 0.199 0.206 0.188 0.173 14 0.244 0.330 0.209 0.188 0.205 0.172 15 0.199 0.209 0.270 0.173 0.172 0.186 16 0.206 0.188 0.173 0.309 0.244 0.199 17 0.188 0.205 0.172 0.244 0.330 0.209 18 0.173 0.172 0.186 0.199 0.209 0.270 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.010 0.004 2 0.010 0.064 0.009 3 0.004 0.009 0.059 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.378 2 0.378 3 0.378 MATRIX M This is a FULL matrix of order 2 by 3 1 2 3 1 2.880 2.691 2.510 2 2.880 2.691 2.510 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.785 0.000 3 0.000 0.468 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.007 0.000 3 0.000 303.978 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.325 0.000 3 0.000 0.379 0.000 MATRIX Q This is a computed FULL matrix of order 3 by 3 [=L*L'] 1 2 3 1 0.143 0.143 0.143 2 0.143 0.143 0.143 3 0.143 0.143 0.143 MATRIX R This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.203 2 0.000 0.203 3 0.000 0.000 0.203 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.134 0.117 0.102 2 0.117 0.133 0.100 3 0.102 0.100 0.115 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.206 0.188 0.173 2 0.188 0.205 0.172 3 0.173 0.172 0.186 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+C] 1 2 3 1 0.063 0.046 0.031 2 0.046 0.062 0.029 3 0.031 0.029 0.044 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.309 0.244 0.199 2 0.244 0.330 0.209 3 0.199 0.209 0.270 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.342 2 0.000 0.229 3 0.000 0.000 0.150 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.068 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.065 2 0.000 0.000 3 0.000 0.000 0.121 Your model has 20 estimated parameters and 1991 Observed statistics -2 times log-likelihood of data >>> 1984.982 Degrees of freedom >>>>>>>>>>>>>>>> 1971 Saturated model fit* >>>>>>>>>>> 1979.626 Saturated model df* >>>>>>>>>>> 1969 Difference Chi-squared >>>>>>>> 5.356 Difference d.f. >>>>>>>>>>>>>>> 2 Probability >>>>>>>>>>>>>>>>>>>> 0.069 Akaike's Information Criterion > 1.356 * Saturated model statistic computed earlier in this job Model number: 1 Model number: 2 Model number: 3 Model number: 1 Model number: 2 Model number: 3 This problem used 14.7% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.00 Execution 0: 0: 1:21.62 TOTAL 0: 0: 1:21.62 Total number of warnings issued: 2 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: DROP @.0 L 1 1 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 C 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 C 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 C This is a computed FULL matrix of order 3 by 3 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 18 by 6 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 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 DIAGONAL matrix of order 3 by 3 1 2 3 1 4 2 0 5 3 0 0 6 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 starting optimization; number of parameters = 19 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.257 0.234 0.182 2 0.234 0.260 0.203 3 0.182 0.203 0.188 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.009 0.000 0.016 2 0.000 0.000 0.000 3 0.016 0.000 0.027 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.009 0.003 2 0.009 0.065 0.007 3 0.003 0.007 0.059 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.509 2 2.879 2.689 2.509 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.911 0.000 3 0.000 0.779 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.006 0.000 3 0.000 303.977 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.625 0.000 3 0.000 0.288 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.204 2 0.000 0.204 3 0.000 0.000 0.204 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.311 0.243 0.200 2 0.243 0.325 0.209 3 0.200 0.209 0.274 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.507 2 0.000 0.217 3 0.000 0.000 0.174 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.096 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.058 2 0.000 0.000 3 0.000 0.000 0.124 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.257 0.234 0.182 2 0.234 0.260 0.203 3 0.182 0.203 0.188 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.009 0.000 0.016 2 0.000 0.000 0.000 3 0.016 0.000 0.027 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.009 0.003 2 0.009 0.065 0.007 3 0.003 0.007 0.059 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.509 2 2.879 2.689 2.509 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.911 0.000 3 0.000 0.779 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.006 0.000 3 0.000 303.977 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.625 0.000 3 0.000 0.288 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.204 2 0.000 0.204 3 0.000 0.000 0.204 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.311 0.243 0.200 2 0.243 0.325 0.209 3 0.200 0.209 0.274 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.507 2 0.000 0.217 3 0.000 0.000 0.174 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.096 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.058 2 0.000 0.000 3 0.000 0.000 0.124 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.257 0.234 0.182 2 0.234 0.260 0.203 3 0.182 0.203 0.188 MATRIX C This is a computed FULL matrix of order 3 by 3 [=(I-O)~&(Y*Y')] 1 2 3 1 0.009 0.000 0.016 2 0.000 0.000 0.000 3 0.016 0.000 0.027 MATRIX D 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.311 0.243 0.200 0.138 0.117 0.107 2 0.243 0.325 0.209 0.117 0.130 0.101 3 0.200 0.209 0.274 0.107 0.101 0.121 4 0.138 0.117 0.107 0.311 0.243 0.200 5 0.117 0.130 0.101 0.243 0.325 0.209 6 0.107 0.101 0.121 0.200 0.209 0.274 7 0.311 0.243 0.200 0.138 0.117 0.107 8 0.243 0.325 0.209 0.117 0.130 0.101 9 0.200 0.209 0.274 0.107 0.101 0.121 10 0.138 0.117 0.107 0.311 0.243 0.200 11 0.117 0.130 0.101 0.243 0.325 0.209 12 0.107 0.101 0.121 0.200 0.209 0.274 13 0.311 0.243 0.200 0.138 0.117 0.107 14 0.243 0.325 0.209 0.117 0.130 0.101 15 0.200 0.209 0.274 0.107 0.101 0.121 16 0.138 0.117 0.107 0.311 0.243 0.200 17 0.117 0.130 0.101 0.243 0.325 0.209 18 0.107 0.101 0.121 0.200 0.209 0.274 MATRIX E This is a computed FULL matrix of order 3 by 3 [=(I-P)~&(Z*Z')+R*R'] 1 2 3 1 0.045 0.009 0.003 2 0.009 0.065 0.007 3 0.003 0.007 0.059 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.509 2 2.879 2.689 2.509 MATRIX N This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.911 0.000 3 0.000 0.779 0.000 MATRIX O This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 0.006 0.000 3 0.000 303.977 0.000 MATRIX P This is a LOWER TRIANGULAR matrix of order 3 by 3 1 2 3 1 0.000 2 2.625 0.000 3 0.000 0.288 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.204 2 0.000 0.204 3 0.000 0.000 0.204 MATRIX S This is a computed FULL matrix of order 3 by 3 [=U+H@Q] 1 2 3 1 0.138 0.117 0.107 2 0.117 0.130 0.101 3 0.107 0.101 0.121 MATRIX T This is a computed FULL matrix of order 3 by 3 [=U+Q] 1 2 3 1 0.138 0.117 0.107 2 0.117 0.130 0.101 3 0.107 0.101 0.121 MATRIX U This is a computed FULL matrix of order 3 by 3 [=H@A+C] 1 2 3 1 0.138 0.117 0.107 2 0.117 0.130 0.101 3 0.107 0.101 0.121 MATRIX V This is a computed FULL matrix of order 3 by 3 [=A+Q+C+E] 1 2 3 1 0.311 0.243 0.200 2 0.243 0.325 0.209 3 0.200 0.209 0.274 MATRIX X This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.507 2 0.000 0.217 3 0.000 0.000 0.174 MATRIX Y This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 0.096 2 0.000 0.000 3 0.000 0.000 0.000 MATRIX Z This is a DIAGONAL matrix of order 3 by 3 1 2 3 1 -0.058 2 0.000 0.000 3 0.000 0.000 0.124 Your model has 19 estimated parameters and 1991 Observed statistics -2 times log-likelihood of data >>> 1998.690 Degrees of freedom >>>>>>>>>>>>>>>> 1972 Saturated model fit* >>>>>>>>>>> 1979.626 Saturated model df* >>>>>>>>>>> 1969 Difference Chi-squared >>>>>>>> 19.064 Difference d.f. >>>>>>>>>>>>>>> 3 Probability >>>>>>>>>>>>>>>>>>>> 0.000 Akaike's Information Criterion > 13.064 * Saturated model statistic computed earlier in this job Model number: 1 Model number: 2 Model number: 3 Model number: 1 Model number: 2 Model number: 3 This problem used 14.7% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.00 Execution 0: 0: 0:10.83 TOTAL 0: 0: 0:10.83 Total number of warnings issued: 2 ______________________________________________________________________________