!==================================================================================== ! SCRIPT maker JJH ! Sex and age as regression variable ! MZ vs DZ, assumptions T1=T2, var DOS=DZ=MZ, Covar DZ=DOS, TT=ST=SS. ! Check ACE. !==================================================================================== ! In this script age and sex effects are estimated for each independent variable. ! Addtive genetics, unique and common environment standard Choleski decomposition. ! Systolic Blood Pressure ! Can be used on any data file as it is supposed not to change per location. The following MX script lines were read for group 1 #DEFINE NVAR 4 #DEFINE NDEF 2 ! NUMBER OF DEFINITION VARIABLES #NGROUPS 3 ! NUMBER OF GROUPS Note: #NGroup set number of groups to 3 G1: CALCULATION GROUP DATA CALCULATION BEGIN MATRICES X FULL NVAR 1 FREE ! ADDTIVE GENETIC Y LOWER NVAR NVAR ! FREE ! COMMON ENVIRONMENT Z DIA NVAR NVAR FREE ! UNIQUE ENVIRONMENT H FULL 1 1 FIX ! HALF-MATRIX G FULL 1 8 FREE ! GENERAL MEANS SAMPLES R FULL NDEF 1 FREE ! DO REGRESSION COEFFICIENTS COVARIATES (DORRET) S FULL NDEF 1 FREE ! HA REGRESSION COEFFICIENTS COVARIATES (HAROLD) T FULL NDEF 1 FREE ! DA REGRESSION COEFFICIENTS COVARIATES (DANIELLE) U FULL NDEF 1 FREE ! NI REGRESSION COEFFICIENTS COVARIATES (NINA) END MATRICES MATRIX H 0.5000; EQ G 1 1 1 G 1 1 5; ! EQUATE MEANS FOR TW1 AND TW2 IN EACH VARIABLE EQ G 1 1 2 G 1 1 6; EQ G 1 1 3 G 1 1 7; EQ G 1 1 4 G 1 1 8; !Starting values ST 9 X 1 1 1 X 1 2 1 X 1 3 1 X 1 4 1 ! st 9 X 1 2 2 X 1 3 3 X 1 4 4 X 1 3 2 X 1 4 2 X 1 4 3; !St 9 Y 1 1 1 Y 1 2 2 Y 1 3 3 Y 1 4 4 Y 1 2 1 Y 1 3 1 Y 1 3 2 Y 1 4 1 Y 1 4 2 Y 1 4 3; ST 9 Z 1 1 1 Z 1 2 2 Z 1 3 3 Z 1 4 4; ST 110 G 1 1 1 G 1 1 2 G 1 1 3 G 1 1 4 G 1 1 5 G 1 1 6 G 1 1 7 G 1 1 8; ST -5.0 R 1 1 1 S 1 1 1 T 1 1 1 U 1 1 1; ST 0.40 R 1 2 1 S 1 2 1 T 1 2 1 U 1 2 1; BEGIN ALGEBRA; A = X*X'; E = Z*Z'; C = Y*Y'; V = A%(A+C+E) | C%(A+C+E) | E%(A+C+E); !TOTAL VARIANCEK K = \STND(A); ! L = \stnd(C); M= \STND(E); END ALGEBRA; END; The following MX script lines were read for group 2 G2: MZM DATA NINPUT_VARS=46 MISSING=-99.0000 NOTE: Missing value *string* set to '-99.0000 ' RECTANGULAR FILE = C11P50.PRN Rectangular continuous data read initiated NOTE: Rectangular file contained 772 records with data that contained a total of 23449 observations LABELS TRAPPREG ID1 ID2 PAIRTP TWZYG DOSEX1 DOAGE1 DOMDBP1 DOMSBP1 DOMED1 HASEX1 HAAGE1 HAMDBP1 HAMSBP1 HAMED1 DASEX1 DAAGE1 DAMDBP1 DAMSBP1 DAMED1 NISEX1 NIAGE1 NIMDBP1 NIMSBP1 NIMED1 DOSEX2 DOAGE2 DOMDBP2 DOMSBP2 DOMED2 HASEX2 HAAGE2 HAMDBP2 HAMSBP2 HAMED2 DASEX2 DAAGE2 DAMDBP2 DAMSBP2 DAMED2 NISEX2 NIAGE2 NIMDBP2 NIMSBP2 NIMED2 PIHAT SELECT IF TWZYG < 4; !MZS NOTE: Select if yields 355 data vectors for analysis NOTE: Vectors contain a total of 10645 observations SELECT IF TWZYG ^= 2; NOTE: Select if yields 283 data vectors for analysis NOTE: Vectors contain a total of 8470 observations SELECT DOSEX1 DOAGE1 HASEX1 HAAGE1 DASEX1 DAAGE1 NISEX1 NIAGE1 DOMSBP1 HAMSBP1 DAMSBP1 NIMSBP1 DOSEX2 DOAGE2 HASEX2 HAAGE2 DASEX2 DAAGE2 NISEX2 NIAGE2 DOMSBP2 HAMSBP2 DAMSBP2 NIMSBP2; DEFINITION NOTE: Selection yields 283 data vectors for analysis NOTE: Vectors contain a total of 5276 observations DOSEX1 DOAGE1 DOSEX2 DOAGE2 HASEX1 HAAGE1 HASEX2 HAAGE2 DASEX1 DAAGE1 DASEX2 DAAGE2 NISEX1 NIAGE1 NISEX2 NIAGE2; MATRICES = GROUP 1 NOTE: Definition yields 283 data vectors for analysis NOTE: Vectors contain a total of 748 observations BEGIN MATRICES B FULL 1 NDEF FIX ! COVARS DO TW1 D FULL 1 NDEF FIX ! COVARS HA TW1 F FULL 1 NDEF FIX ! COVARS DA TW1 K FULL 1 NDEF FIX ! COVARS NI TW1 *** WARNING! *** Matrix was previously defined for this group. Mx will redefine L FULL 1 NDEF FIX ! COVARS DO TW2 M FULL 1 NDEF FIX ! COVARS HA TW2 *** WARNING! *** Matrix was previously defined for this group. Mx will redefine N FULL 1 NDEF FIX ! COVARS DA TW2 O FULL 1 NDEF FIX ! COVARS NI TW1 END MATRICES; SPECIFY B DOSEX1 DOAGE1; ! DEFINITION VARIABLE MATRICES TO ADJUST MEANS TW1. SPECIFY D HASEX1 HAAGE1; SPECIFY F DASEX1 DAAGE1; SPECIFY K NISEX1 NIAGE1; SPECIFY L DOSEX2 DOAGE2; ! DEFINITION VARIABLE MATRICES TO ADJUST MEANS TW2. SPECIFY M HASEX2 HAAGE2; SPECIFY N DASEX2 DAAGE2; SPECIFY O NISEX2 NIAGE2; BEGIN ALGEBRA; J = B*R | D*S | F*T | K*U | L*R | M*S | N*T | O*U; END ALGEBRA; MEANS G+J; COVARIANCES A+C+E | C+A_ C+A | A+C+E; END; The following MX script lines were read for group 3 G3: DZM+SIBS DATA NINPUT_VARS=46 MISSING=-99.0000 NOTE: Missing value *string* set to '-99.0000 ' RECTANGULAR FILE = C11P50.PRN Rectangular continuous data read initiated NOTE: Rectangular file contained 772 records with data that contained a total of 23449 observations LABELS TRAPPREG ID1 ID2 PAIRTP TWZYG DOSEX1 DOAGE1 DOMDBP1 DOMSBP1 DOMED1 HASEX1 HAAGE1 HAMDBP1 HAMSBP1 HAMED1 DASEX1 DAAGE1 DAMDBP1 DAMSBP1 DAMED1 NISEX1 NIAGE1 NIMDBP1 NIMSBP1 NIMED1 DOSEX2 DOAGE2 DOMDBP2 DOMSBP2 DOMED2 HASEX2 HAAGE2 HAMDBP2 HAMSBP2 HAMED2 DASEX2 DAAGE2 DAMDBP2 DAMSBP2 DAMED2 NISEX2 NIAGE2 NIMDBP2 NIMSBP2 NIMED2 PIHAT SELECT IF TWZYG > 1; NOTE: Select if yields 656 data vectors for analysis NOTE: Vectors contain a total of 19901 observations SELECT IF TWZYG ^= 3; NOTE: Select if yields 489 data vectors for analysis NOTE: Vectors contain a total of 14979 observations SELECT DOSEX1 DOAGE1 HASEX1 HAAGE1 DASEX1 DAAGE1 NISEX1 NIAGE1 DOMSBP1 HAMSBP1 DAMSBP1 NIMSBP1 DOSEX2 DOAGE2 HASEX2 HAAGE2 DASEX2 DAAGE2 NISEX2 NIAGE2 DOMSBP2 HAMSBP2 DAMSBP2 NIMSBP2; DEFINITION NOTE: Selection yields 489 data vectors for analysis NOTE: Vectors contain a total of 9231 observations DOSEX1 DOAGE1 DOSEX2 DOAGE2 Note: Global variable previously defined. Updating DOSEX1 Note: Global variable previously defined. Updating DOAGE1 Note: Global variable previously defined. Updating DOSEX2 Note: Global variable previously defined. Updating DOAGE2 HASEX1 HAAGE1 HASEX2 HAAGE2 Note: Global variable previously defined. Updating HASEX1 Note: Global variable previously defined. Updating HAAGE1 Note: Global variable previously defined. Updating HASEX2 Note: Global variable previously defined. Updating HAAGE2 DASEX1 DAAGE1 DASEX2 DAAGE2 Note: Global variable previously defined. Updating DASEX1 Note: Global variable previously defined. Updating DAAGE1 Note: Global variable previously defined. Updating DASEX2 Note: Global variable previously defined. Updating DAAGE2 NISEX1 NIAGE1 NISEX2 NIAGE2; Note: Global variable previously defined. Updating NISEX1 Note: Global variable previously defined. Updating NIAGE1 Note: Global variable previously defined. Updating NISEX2 Note: Global variable previously defined. Updating NIAGE2 MATRICES = GROUP 1 NOTE: Definition yields 489 data vectors for analysis NOTE: Vectors contain a total of 1407 observations BEGIN MATRICES B FULL 1 NDEF FIX ! COVARS DO TW1 D FULL 1 NDEF FIX ! COVARS HA TW1 F FULL 1 NDEF FIX ! COVARS DA TW1 K FULL 1 NDEF FIX ! COVARS NI TW1 *** WARNING! *** Matrix was previously defined for this group. Mx will redefine L FULL 1 NDEF FIX ! COVARS DO TW2 M FULL 1 NDEF FIX ! COVARS HA TW2 *** WARNING! *** Matrix was previously defined for this group. Mx will redefine N FULL 1 NDEF FIX ! COVARS DA TW2 O FULL 1 NDEF FIX ! COVARS NI TW1 END MATRICES; SPECIFY B DOSEX1 DOAGE1; SPECIFY D HASEX1 HAAGE1; SPECIFY F DASEX1 DAAGE1; SPECIFY K NISEX1 NIAGE1; SPECIFY L DOSEX2 DOAGE2; SPECIFY M HASEX2 HAAGE2; SPECIFY N DASEX2 DAAGE2; SPECIFY O NISEX2 NIAGE2; BEGIN ALGEBRA; J = B*R | D*S | F*T | K*U | L*R | M*S | N*T | O*U; END ALGEBRA; MEANS G+J; COVARIANCES A+C+E | H@A+C_ H@A+C | A+C+E; OPTION NDECIMALS=4; OPTION RSIDUALS; END; Summary of VL file data for group 2 NIAGE2 NISEX2 DAAGE2 DASEX2 HAAGE2 HASEX2 DOAGE2 Code -16.0000 -15.0000 -14.0000 -13.0000 -12.0000 -11.0000 -10.0000 Number 283.0000 283.0000 283.0000 283.0000 283.0000 283.0000 283.0000 Mean 8.0536 0.4276 7.9402 0.4311 17.7597 0.7314 18.2968 Variance 218.2635 0.5769 211.3521 0.5845 436.8045 0.7265 444.5409 Minimum 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Maximum 71.4667 2.0000 71.4666 2.0000 69.0000 2.0000 69.0000 DOSEX2 NIAGE1 NISEX1 DAAGE1 DASEX1 HAAGE1 HASEX1 Code -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 -3.0000 Number 283.0000 283.0000 283.0000 283.0000 283.0000 283.0000 283.0000 Mean 0.7420 15.1537 0.5230 14.8392 0.5088 4.0847 0.3781 Variance 0.7214 455.1857 0.6099 450.6234 0.5962 50.8356 0.4896 Minimum 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Maximum 2.0000 62.5000 2.0000 62.5000 2.0000 21.7500 2.0000 DOAGE1 DOSEX1 DOMSBP1 HAMSBP1 DAMSBP1 NIMSBP1 DOMSBP2 Code -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 5.0000 Number 283.0000 283.0000 70.0000 95.0000 135.0000 74.0000 71.0000 Mean 4.0317 0.3710 117.8771 125.2193 126.3259 126.6490 116.9585 Variance 50.4368 0.4807 53.8938 158.4679 184.3793 132.7988 55.4185 Minimum 0.0000 0.0000 105.8500 99.3333 95.6667 104.0000 104.5000 Maximum 21.7500 2.0000 143.6500 165.0000 171.3333 159.8571 137.0000 HAMSBP2 DAMSBP2 NIMSBP2 Code 6.0000 7.0000 8.0000 Number 97.0000 132.0000 74.0000 Mean 124.3797 125.6149 126.2235 Variance 154.5640 262.1247 127.9241 Minimum 97.6667 97.3333 106.5000 Maximum 155.6667 180.0000 162.0000 Summary of VL file data for group 3 NIAGE2 NISEX2 DAAGE2 DASEX2 HAAGE2 HASEX2 DOAGE2 Code -16.0000 -15.0000 -14.0000 -13.0000 -12.0000 -11.0000 -10.0000 Number 489.0000 489.0000 489.0000 489.0000 489.0000 489.0000 489.0000 Mean 20.7177 0.8650 21.1229 0.8671 14.6687 0.5808 15.2106 Variance 435.2160 0.7425 441.0288 0.7247 423.0927 0.6607 434.9679 Minimum 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Maximum 74.9241 2.0000 74.9241 2.0000 65.0000 2.0000 65.0000 DOSEX2 NIAGE1 NISEX1 DAAGE1 DASEX1 HAAGE1 HASEX1 Code -9.0000 -8.0000 -7.0000 -6.0000 -5.0000 -4.0000 -3.0000 Number 489.0000 489.0000 489.0000 489.0000 489.0000 489.0000 489.0000 Mean 0.5419 12.5726 0.4335 14.6033 0.4642 3.4015 0.3006 Variance 0.5714 406.4717 0.5401 445.3344 0.5146 47.4686 0.4147 Minimum 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Maximum 2.0000 59.5000 2.0000 60.5000 2.0000 22.1000 2.0000 DOAGE1 DOSEX1 DOMSBP1 HAMSBP1 DAMSBP1 NIMSBP1 DOMSBP2 Code -2.0000 -1.0000 1.0000 2.0000 3.0000 4.0000 5.0000 Number 489.0000 489.0000 95.0000 162.0000 186.0000 275.0000 97.0000 Mean 3.3500 0.2924 116.6458 126.6903 129.1774 127.7846 117.3933 Variance 47.2802 0.4032 48.2380 235.5634 238.7289 160.1679 68.8937 Minimum 0.0000 0.0000 101.3500 97.3333 97.6667 98.2000 99.3500 Maximum 22.1000 2.0000 136.5000 177.6667 190.6667 165.1667 139.6500 HAMSBP2 DAMSBP2 NIMSBP2 Code 6.0000 7.0000 8.0000 Number 140.0000 182.0000 270.0000 Mean 126.9810 130.8883 127.3998 Variance 146.1917 314.0480 143.5528 Minimum 101.6667 101.3333 103.6000 Maximum 170.0000 182.3333 183.3333 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: CALCULATION GROUP MATRIX A This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 8 1 2 3 4 5 6 7 8 1 9 10 11 12 9 10 11 12 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX R This is a FULL matrix of order 2 by 1 1 1 17 2 18 MATRIX S This is a FULL matrix of order 2 by 1 1 1 19 2 20 MATRIX T This is a FULL matrix of order 2 by 1 1 1 21 2 22 MATRIX U This is a FULL matrix of order 2 by 1 1 1 23 2 24 MATRIX V This is a computed FULL matrix of order 4 by 12 It has no free parameters specified MATRIX X This is a FULL matrix of order 4 by 1 1 1 1 2 2 3 3 4 4 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 It has no free parameters specified MATRIX Z This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 5 2 0 6 3 0 0 7 4 0 0 0 8 GROUP NUMBER: 2 G2: MZM MATRIX A This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 2 1 2 1 -1 -2 MATRIX C This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX D This is a FULL matrix of order 1 by 2 1 2 1 -5 -6 MATRIX E This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 2 1 2 1 -9 -10 MATRIX G This is a FULL matrix of order 1 by 8 1 2 3 4 5 6 7 8 1 9 10 11 12 9 10 11 12 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 8 It has no free parameters specified MATRIX K This is a FULL matrix of order 1 by 2 1 2 1 -13 -14 MATRIX L This is a FULL matrix of order 1 by 2 1 2 1 -3 -4 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 -7 -8 MATRIX N This is a FULL matrix of order 1 by 2 1 2 1 -11 -12 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 -15 -16 MATRIX R This is a FULL matrix of order 2 by 1 1 1 17 2 18 MATRIX S This is a FULL matrix of order 2 by 1 1 1 19 2 20 MATRIX T This is a FULL matrix of order 2 by 1 1 1 21 2 22 MATRIX U This is a FULL matrix of order 2 by 1 1 1 23 2 24 MATRIX V This is a computed FULL matrix of order 4 by 12 It has no free parameters specified MATRIX X This is a FULL matrix of order 4 by 1 1 1 1 2 2 3 3 4 4 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 It has no free parameters specified MATRIX Z This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 5 2 0 6 3 0 0 7 4 0 0 0 8 GROUP NUMBER: 3 G3: DZM+SIBS MATRIX A This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX B This is a FULL matrix of order 1 by 2 1 2 1 -1 -2 MATRIX C This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX D This is a FULL matrix of order 1 by 2 1 2 1 -5 -6 MATRIX E This is a computed FULL matrix of order 4 by 4 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 2 1 2 1 -9 -10 MATRIX G This is a FULL matrix of order 1 by 8 1 2 3 4 5 6 7 8 1 9 10 11 12 9 10 11 12 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 8 It has no free parameters specified MATRIX K This is a FULL matrix of order 1 by 2 1 2 1 -13 -14 MATRIX L This is a FULL matrix of order 1 by 2 1 2 1 -3 -4 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 -7 -8 MATRIX N This is a FULL matrix of order 1 by 2 1 2 1 -11 -12 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 -15 -16 MATRIX R This is a FULL matrix of order 2 by 1 1 1 17 2 18 MATRIX S This is a FULL matrix of order 2 by 1 1 1 19 2 20 MATRIX T This is a FULL matrix of order 2 by 1 1 1 21 2 22 MATRIX U This is a FULL matrix of order 2 by 1 1 1 23 2 24 MATRIX V This is a computed FULL matrix of order 4 by 12 It has no free parameters specified MATRIX X This is a FULL matrix of order 4 by 1 1 1 1 2 2 3 3 4 4 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 It has no free parameters specified MATRIX Z This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 5 2 0 6 3 0 0 7 4 0 0 0 8 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 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 G1: CALCULATION GROUP MATRIX A This is a computed FULL matrix of order 4 by 4 [=X*X'] 1 2 3 4 1 35.1586 45.9446 63.8596 45.9162 2 45.9446 60.0396 83.4505 60.0024 3 63.8596 83.4505 115.9899 83.3989 4 45.9162 60.0024 83.3989 59.9653 MATRIX C This is a computed FULL matrix of order 4 by 4 [=Y*Y'] 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 4 by 4 [=Z*Z'] 1 2 3 4 1 24.0944 0.0000 0.0000 0.0000 2 0.0000 91.8092 0.0000 0.0000 3 0.0000 0.0000 77.1924 0.0000 4 0.0000 0.0000 0.0000 65.0514 MATRIX G This is a FULL matrix of order 1 by 8 1 2 3 4 5 6 7 1 110.6834 124.7202 116.1748 128.0688 110.6834 124.7202 116.1748 8 1 128.0688 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 4 by 4 [=\STND(A)] 1 2 3 4 1 1.0000 1.0000 1.0000 1.0000 2 1.0000 1.0000 1.0000 1.0000 3 1.0000 1.0000 1.0000 1.0000 4 1.0000 1.0000 1.0000 1.0000 MATRIX M This is a computed FULL matrix of order 4 by 4 [=\STND(E)] 1 2 3 4 1 1.0000 0.0000 0.0000 0.0000 2 0.0000 1.0000 0.0000 0.0000 3 0.0000 0.0000 1.0000 0.0000 4 0.0000 0.0000 0.0000 1.0000 MATRIX R This is a FULL matrix of order 2 by 1 1 1 -3.7860 2 0.7127 MATRIX S This is a FULL matrix of order 2 by 1 1 1 -7.2228 2 0.2610 MATRIX T This is a FULL matrix of order 2 by 1 1 1 -6.3103 2 0.5559 MATRIX U This is a FULL matrix of order 2 by 1 1 1 -5.9046 2 0.2267 MATRIX V This is a computed FULL matrix of order 4 by 12 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 7 1 0.5934 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 2 1.0000 0.3954 1.0000 1.0000 0.0000 0.0000 0.0000 3 1.0000 1.0000 0.6004 1.0000 0.0000 0.0000 0.0000 4 1.0000 1.0000 1.0000 0.4797 0.0000 0.0000 0.0000 8 9 10 11 12 1 0.0000 0.4066 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.6046 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.3996 0.0000 4 0.0000 0.0000 0.0000 0.0000 0.5203 MATRIX X This is a FULL matrix of order 4 by 1 1 1 5.9295 2 7.7485 3 10.7699 4 7.7437 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 4 0.0000 0.0000 0.0000 0.0000 MATRIX Z This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 4.9086 2 0.0000 9.5817 3 0.0000 0.0000 8.7859 4 0.0000 0.0000 0.0000 8.0654 GROUP NUMBER: 2 G2: MZM MATRIX A This is a computed FULL matrix of order 4 by 4 [=X*X'] 1 2 3 4 1 35.1586 45.9446 63.8596 45.9162 2 45.9446 60.0396 83.4505 60.0024 3 63.8596 83.4505 115.9899 83.3989 4 45.9162 60.0024 83.3989 59.9653 MATRIX B This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX C This is a computed FULL matrix of order 4 by 4 [=Y*Y'] 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 0.0000 0.0000 MATRIX D This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 4 by 4 [=Z*Z'] 1 2 3 4 1 24.0944 0.0000 0.0000 0.0000 2 0.0000 91.8092 0.0000 0.0000 3 0.0000 0.0000 77.1924 0.0000 4 0.0000 0.0000 0.0000 65.0514 MATRIX F This is a FULL matrix of order 1 by 2 1 2 1 2.0000 34.0000 MATRIX G This is a FULL matrix of order 1 by 8 1 2 3 4 5 6 7 1 110.6834 124.7202 116.1748 128.0688 110.6834 124.7202 116.1748 8 1 128.0688 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 8 [=B*R|D*S|F*T|K*U|L*R|M*S|N*T|O*U] 1 2 3 4 5 6 7 1 0.0000 0.0000 6.2808 0.0000 0.0000 0.0000 6.2808 8 1 0.0000 MATRIX K This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX L This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX N This is a FULL matrix of order 1 by 2 1 2 1 2.0000 34.0000 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX R This is a FULL matrix of order 2 by 1 1 1 -3.7860 2 0.7127 MATRIX S This is a FULL matrix of order 2 by 1 1 1 -7.2228 2 0.2610 MATRIX T This is a FULL matrix of order 2 by 1 1 1 -6.3103 2 0.5559 MATRIX U This is a FULL matrix of order 2 by 1 1 1 -5.9046 2 0.2267 MATRIX V This is a computed FULL matrix of order 4 by 12 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 7 1 0.5934 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 2 1.0000 0.3954 1.0000 1.0000 0.0000 0.0000 0.0000 3 1.0000 1.0000 0.6004 1.0000 0.0000 0.0000 0.0000 4 1.0000 1.0000 1.0000 0.4797 0.0000 0.0000 0.0000 8 9 10 11 12 1 0.0000 0.4066 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.6046 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.3996 0.0000 4 0.0000 0.0000 0.0000 0.0000 0.5203 MATRIX X This is a FULL matrix of order 4 by 1 1 1 5.9295 2 7.7485 3 10.7699 4 7.7437 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 4 0.0000 0.0000 0.0000 0.0000 MATRIX Z This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 4.9086 2 0.0000 9.5817 3 0.0000 0.0000 8.7859 4 0.0000 0.0000 0.0000 8.0654 GROUP NUMBER: 3 G3: DZM+SIBS MATRIX A This is a computed FULL matrix of order 4 by 4 [=X*X'] 1 2 3 4 1 35.1586 45.9446 63.8596 45.9162 2 45.9446 60.0396 83.4505 60.0024 3 63.8596 83.4505 115.9899 83.3989 4 45.9162 60.0024 83.3989 59.9653 MATRIX B This is a FULL matrix of order 1 by 2 1 2 1 1.0000 17.4300 MATRIX C This is a computed FULL matrix of order 4 by 4 [=Y*Y'] 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 0.0000 0.0000 MATRIX D This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX E This is a computed FULL matrix of order 4 by 4 [=Z*Z'] 1 2 3 4 1 24.0944 0.0000 0.0000 0.0000 2 0.0000 91.8092 0.0000 0.0000 3 0.0000 0.0000 77.1924 0.0000 4 0.0000 0.0000 0.0000 65.0514 MATRIX F This is a FULL matrix of order 1 by 2 1 2 1 1.0000 26.0000 MATRIX G This is a FULL matrix of order 1 by 8 1 2 3 4 5 6 7 1 110.6834 124.7202 116.1748 128.0688 110.6834 124.7202 116.1748 8 1 128.0688 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 8 [=B*R|D*S|F*T|K*U|L*R|M*S|N*T|O*U] 1 2 3 4 5 6 7 1 8.6364 0.0000 8.1438 0.0000 8.6364 0.0000 8.1438 8 1 0.0000 MATRIX K This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX L This is a FULL matrix of order 1 by 2 1 2 1 1.0000 17.4300 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX N This is a FULL matrix of order 1 by 2 1 2 1 1.0000 26.0000 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX R This is a FULL matrix of order 2 by 1 1 1 -3.7860 2 0.7127 MATRIX S This is a FULL matrix of order 2 by 1 1 1 -7.2228 2 0.2610 MATRIX T This is a FULL matrix of order 2 by 1 1 1 -6.3103 2 0.5559 MATRIX U This is a FULL matrix of order 2 by 1 1 1 -5.9046 2 0.2267 MATRIX V This is a computed FULL matrix of order 4 by 12 [=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)] 1 2 3 4 5 6 7 1 0.5934 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 2 1.0000 0.3954 1.0000 1.0000 0.0000 0.0000 0.0000 3 1.0000 1.0000 0.6004 1.0000 0.0000 0.0000 0.0000 4 1.0000 1.0000 1.0000 0.4797 0.0000 0.0000 0.0000 8 9 10 11 12 1 0.0000 0.4066 0.0000 0.0000 0.0000 2 0.0000 0.0000 0.6046 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.3996 0.0000 4 0.0000 0.0000 0.0000 0.0000 0.5203 MATRIX X This is a FULL matrix of order 4 by 1 1 1 5.9295 2 7.7485 3 10.7699 4 7.7437 MATRIX Y This is a LOWER TRIANGULAR matrix of order 4 by 4 1 2 3 4 1 0.0000 2 0.0000 0.0000 3 0.0000 0.0000 0.0000 4 0.0000 0.0000 0.0000 0.0000 MATRIX Z This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 4.9086 2 0.0000 9.5817 3 0.0000 0.0000 8.7859 4 0.0000 0.0000 0.0000 8.0654 Vector of OBSERVED means DOMSBP1 HAMSBP1 DAMSBP1 NIMSBP1 DOMSBP2 HAMSBP2 DAMSBP2 Mean 116.6458 126.6903 129.1774 127.7846 117.3933 126.9810 130.8883 NIMSBP2 Mean 127.3998 Vector of EXPECTED means DOMSBP1 HAMSBP1 DAMSBP1 NIMSBP1 DOMSBP2 HAMSBP2 DAMSBP2 Mean 119.3198 124.7202 124.3186 128.0688 119.3198 124.7202 124.3186 NIMSBP2 Mean 128.0688 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DOMSBP1 HAMSBP1 DAMSBP1 NIMSBP1 DOMSBP2 HAMSBP2 DAMSBP2 DOMSBP1 59.2531 HAMSBP1 45.9446 151.8488 DAMSBP1 63.8596 83.4505 193.1823 NIMSBP1 45.9162 60.0024 83.3989 125.0168 DOMSBP2 17.5793 22.9723 31.9298 22.9581 59.2531 HAMSBP2 22.9723 30.0198 41.7252 30.0012 45.9446 151.8488 DAMSBP2 31.9298 41.7252 57.9949 41.6994 63.8596 83.4505 193.1823 NIMSBP2 22.9581 30.0012 41.6994 29.9827 45.9162 60.0024 83.3989 NIMSBP2 DOMSBP1 HAMSBP1 DAMSBP1 NIMSBP1 DOMSBP2 HAMSBP2 DAMSBP2 NIMSBP2 125.0168 Function value of this group: 10734.0855 Where the fit function is -2 * Log-likelihood of raw data *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 20 estimated parameters and 2155 Observed statistics -2 times log-likelihood of data >>> 16270.298 Degrees of freedom >>>>>>>>>>>>>>>> 2135 This problem used 2.6% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 3.47 Execution 0: 0: 0:28.59 TOTAL 0: 0: 0:32.06 Total number of warnings issued: 6 ______________________________________________________________________________