** Mx startup successful ** **MX-Linux version 1.52b** #repeat 11 Note: beginning of #Repeat section ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 1 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 1 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.0001 5.5000 1.5000 1.0000 0.5000 Variance 0.0000 0.0000 2.2500 0.2500 0.0000 0.2500 Minimum 2.4089 0.0001 4.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.0001 7.0000 2.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.0001 13.5000 1.5000 1.0000 0.5000 Variance 0.0000 0.0000 90.2500 0.2500 0.0000 0.2500 Minimum 2.4089 0.0001 4.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.0001 23.0000 2.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7722 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1698 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0580 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5559 2 0.5559 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9420 2 0.9420 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8788 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4120 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2408 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7722 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1698 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0580 MATRIX F This is a FULL matrix of order 1 by 1 1 1 7.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5559 2 0.5559 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9420 2 0.9420 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.0001 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.0001 2 0.0001 3 0.0002 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8788 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4120 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2408 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9420 1.0000 Function value of this group: 29.6691 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7722 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1698 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0580 MATRIX F This is a FULL matrix of order 1 by 1 1 1 4.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5559 2 0.5559 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9420 2 0.9420 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.0001 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.0001 2 0.0001 3 0.0002 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8788 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4120 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2408 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5559 1.0000 Function value of this group: 60.0816 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7722 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1698 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0580 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5559 2 0.5559 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9420 2 0.9420 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8788 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4120 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2408 Your model has 3 estimated parameters and 9 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 89.751 Degrees of freedom >>>>>>>>>>>>>>>> 6 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.7722 0.2790 0.8245 0 1 0 0 C 1 1 1 95.0 0.1698 0.0000 0.5924 6 5 0 0 E 1 1 1 95.0 0.0580 0.0100 0.2142 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 1.72 TOTAL 0: 0: 0: 1.78 Total number of warnings issued: 0 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 2 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 2 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.1110 13.6667 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 77.5556 0.6667 0.0000 0.2222 Minimum 2.4089 0.1110 6.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.1110 26.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.1110 20.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 366.2222 0.6667 0.0000 0.2222 Minimum 2.4089 0.1110 3.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.1110 47.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5933 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3880 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0187 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.6846 2 0.6846 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9813 2 0.9813 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7703 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6229 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1368 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5933 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3880 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0187 MATRIX F This is a FULL matrix of order 1 by 1 1 1 6.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0017 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.6846 2 0.6846 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9813 2 0.9813 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.1110 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.1110 2 0.1110 3 0.2097 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7703 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6229 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1368 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9813 1.0000 Function value of this group: 88.7069 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5933 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3880 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0187 MATRIX F This is a FULL matrix of order 1 by 1 1 1 3.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0018 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.6846 2 0.6846 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9813 2 0.9813 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.1110 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.1110 2 0.1110 3 0.2097 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7703 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6229 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1368 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.6846 1.0000 Function value of this group: 144.8680 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5933 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3880 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0187 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.6846 2 0.6846 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9813 2 0.9813 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7703 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6229 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1368 Your model has 3 estimated parameters and 13 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 233.575 Degrees of freedom >>>>>>>>>>>>>>>> 10 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.5933 0.3646 0.8715 0 0 1 0 C 1 1 1 95.0 0.3880 0.1119 0.6117 0 0 0 0 E 1 1 1 95.0 0.0187 0.0100 0.0529 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 1.74 TOTAL 0: 0: 0: 1.80 Total number of warnings issued: 0 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 3 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 3 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.4240 58.0000 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 88.6667 0.6667 0.0000 0.2222 Minimum 2.4089 0.4240 49.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.4240 71.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.4240 172.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.000042554.8889 0.6667 0.0000 0.2222 Minimum 2.4089 0.4240 21.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.4240 464.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 *** 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: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.9094 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0612 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0294 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5159 2 0.5159 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9706 2 0.9706 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9536 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2474 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1714 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.9094 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0612 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0294 MATRIX F This is a FULL matrix of order 1 by 1 1 1 71.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0057 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5159 2 0.5159 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9706 2 0.9706 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4240 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.4240 2 0.4240 3 0.6682 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9536 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2474 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1714 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9706 1.0000 Function value of this group: 334.9821 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.9094 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0612 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0294 MATRIX F This is a FULL matrix of order 1 by 1 1 1 32.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0066 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5159 2 0.5159 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9706 2 0.9706 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4240 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.4240 2 0.4240 3 0.6682 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9536 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2474 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1714 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5159 1.0000 Function value of this group: 1004.1585 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.9094 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0612 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0294 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5159 2 0.5159 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9706 2 0.9706 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9536 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2474 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1714 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 3 estimated parameters and 13 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 1339.141 Degrees of freedom >>>>>>>>>>>>>>>> 10 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.9094 0.7852 0.9816 0 0 0 1 C 1 1 1 95.0 0.0612 0.0000 0.1825 0 0 0 0 E 1 1 1 95.0 0.0294 0.0176 0.0465 0 0 6 1 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 9.64 TOTAL 0: 0: 0: 9.70 Total number of warnings issued: 2 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 4 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 4 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4066 0.9999 34.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 12.2500 1.0000 0.0000 0.2500 Minimum 2.4066 0.9999 31.0000 1.0000 1.0000 0.0000 Maximum 2.4066 0.9999 38.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4066 0.9999 48.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 1156.0000 1.0000 0.0000 0.2500 Minimum 2.4066 0.9999 14.0000 1.0000 1.0000 0.0000 Maximum 2.4066 0.9999 82.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 *** 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: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.4711 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4652 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0637 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.7007 2 0.7007 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9363 2 0.9363 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6864 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6820 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2524 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.4711 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4652 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0637 MATRIX F This is a FULL matrix of order 1 by 1 1 1 31.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0111 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.7007 2 0.7007 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9363 2 0.9363 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4066 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6864 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6820 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2524 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4066 2.4066 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9363 1.0000 Function value of this group: 147.6221 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.4711 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4652 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0637 MATRIX F This is a FULL matrix of order 1 by 1 1 1 14.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0141 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.7007 2 0.7007 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9363 2 0.9363 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4066 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6864 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6820 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2524 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4066 2.4066 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.7007 1.0000 Function value of this group: 193.4354 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.4711 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4652 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0637 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.7007 2 0.7007 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9363 2 0.9363 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6864 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6820 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2524 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 3 estimated parameters and 9 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 341.058 Degrees of freedom >>>>>>>>>>>>>>>> 6 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.4711 0.2456 0.7434 0 0 0 0 C 1 1 1 95.0 0.4652 0.1998 0.6748 0 0 0 0 E 1 1 1 95.0 0.0637 0.0329 0.1140 0 1 0 1 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 2.81 TOTAL 0: 0: 0: 2.87 Total number of warnings issued: 4 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 5 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 5 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4135 0.9999 11.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 6.2500 1.0000 0.0000 0.2500 Minimum 2.4135 0.9999 9.0000 1.0000 1.0000 0.0000 Maximum 2.4135 0.9999 14.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4135 0.9999 21.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 90.2500 1.0000 0.0000 0.2500 Minimum 2.4135 0.9999 12.0000 1.0000 1.0000 0.0000 Maximum 2.4135 0.9999 31.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 *** 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: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.2496 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7252 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0252 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.8500 2 0.8500 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9748 2 0.9748 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4996 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8516 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1587 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.2496 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7252 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0252 MATRIX F This is a FULL matrix of order 1 by 1 1 1 14.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0098 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.8500 2 0.8500 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9748 2 0.9748 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4135 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4996 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8516 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1587 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4135 2.4135 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9748 1.0000 Function value of this group: 43.2657 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.2496 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7252 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0252 MATRIX F This is a FULL matrix of order 1 by 1 1 1 12.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0124 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.8500 2 0.8500 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9748 2 0.9748 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4135 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4996 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8516 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1587 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4135 2.4135 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8500 1.0000 Function value of this group: 93.8934 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.2496 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7252 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0252 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.8500 2 0.8500 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9748 2 0.9748 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4996 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8516 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1587 *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 3 estimated parameters and 9 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 137.159 Degrees of freedom >>>>>>>>>>>>>>>> 6 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.2496 0.0646 0.5066 1 0 0 0 C 1 1 1 95.0 0.7252 0.4724 0.8864 0 0 0 0 E 1 1 1 95.0 0.0252 0.0100 0.0797 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 1.40 TOTAL 0: 0: 0: 1.46 Total number of warnings issued: 6 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 6 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 6 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 1.8250 0.9999 10.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 12.2500 1.0000 0.0000 0.2500 Minimum 1.8250 0.9999 7.0000 1.0000 1.0000 0.0000 Maximum 1.8250 0.9999 14.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 1.8250 0.9999 21.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 225.0000 1.0000 0.0000 0.2500 Minimum 1.8250 0.9999 6.0000 1.0000 1.0000 0.0000 Maximum 1.8250 0.9999 36.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5274 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3016 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1710 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5653 2 0.5653 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8290 2 0.8290 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7262 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5492 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4136 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5274 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3016 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1710 MATRIX F This is a FULL matrix of order 1 by 1 1 1 7.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0510 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5653 2 0.5653 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8290 2 0.8290 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 1.8250 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7262 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5492 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4136 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 1.8250 1.8250 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8290 1.0000 Function value of this group: 46.1417 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5274 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3016 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1710 MATRIX F This is a FULL matrix of order 1 by 1 1 1 6.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0595 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5653 2 0.5653 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8290 2 0.8290 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 1.8250 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7262 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5492 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4136 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 1.8250 1.8250 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5653 1.0000 Function value of this group: 84.3564 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5274 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3016 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1710 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5653 2 0.5653 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8290 2 0.8290 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7262 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5492 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4136 Your model has 3 estimated parameters and 9 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 130.498 Degrees of freedom >>>>>>>>>>>>>>>> 6 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.5274 0.0000 0.7073 1 0 0 0 C 1 1 1 95.0 0.3016 0.0000 0.7533 0 0 0 0 E 1 1 1 95.0 0.1710 0.0563 0.4001 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.05 Execution 0: 0: 0: 3.39 TOTAL 0: 0: 0: 3.44 Total number of warnings issued: 6 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 7 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! population sample, so no weight required ! Weight \part(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 7 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! population sample, so no weight required (could estimate threshold) ! Weight \part(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.0465 0.9999 1977.6667 1.6667 0.6667 0.3333 Variance 0.0000 0.0000********** 0.8889 0.2222 0.2222 Minimum 2.0465 0.9999 60.0000 1.0000 0.0000 0.0000 Maximum 2.0465 0.9999 5739.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.0465 0.9999 2518.0000 1.6667 0.6667 0.3333 Variance 0.0000 0.0000********** 0.8889 0.2222 0.2222 Minimum 2.0465 0.9999 18.0000 1.0000 0.0000 0.0000 Maximum 2.0465 0.9999 7277.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8294 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1706 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4147 2 0.4147 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8294 2 0.8294 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9107 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4130 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8294 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1706 MATRIX F This is a FULL matrix of order 1 by 1 1 1 60.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0313 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4147 2 0.4147 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8294 2 0.8294 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0465 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9107 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4130 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.0465 2.0465 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8294 1.0000 Function value of this group: 2134.9073 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8294 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1706 MATRIX F This is a FULL matrix of order 1 by 1 1 1 18.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0381 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4147 2 0.4147 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8294 2 0.8294 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0465 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9107 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4130 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.0465 2.0465 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4147 1.0000 Function value of this group: 2867.8828 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8294 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1706 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4147 2 0.4147 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8294 2 0.8294 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9107 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4130 Your model has 3 estimated parameters and 13 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 5002.790 Degrees of freedom >>>>>>>>>>>>>>>> 10 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.8294 0.5927 0.8740 1 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.0000 0 0 0 0 E 1 1 1 95.0 0.1706 0.1260 0.2263 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 0.80 TOTAL 0: 0: 0: 0.86 Total number of warnings issued: 6 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 8 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 8 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.9999 15.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 0.2500 1.0000 0.0000 0.2500 Minimum 2.4089 0.9999 15.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.9999 16.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.9999 14.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 169.0000 1.0000 0.0000 0.2500 Minimum 2.4089 0.9999 1.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.9999 27.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.9463 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0537 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4731 2 0.4731 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9463 2 0.9463 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9728 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2318 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.9463 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0537 MATRIX F This is a FULL matrix of order 1 by 1 1 1 15.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0108 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4731 2 0.4731 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9463 2 0.9463 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9728 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2318 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9463 1.0000 Function value of this group: 65.1248 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.9463 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0537 MATRIX F This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0151 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4731 2 0.4731 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9463 2 0.9463 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9728 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2318 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4731 1.0000 Function value of this group: 46.3111 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.9463 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0537 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4731 2 0.4731 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9463 2 0.9463 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9728 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2318 Your model has 3 estimated parameters and 9 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 111.436 Degrees of freedom >>>>>>>>>>>>>>>> 6 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.9463 0.8717 0.9814 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.0000 0 0 0 0 E 1 1 1 95.0 0.0537 0.0186 0.1283 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 1.55 TOTAL 0: 0: 0: 1.61 Total number of warnings issued: 6 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 9 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 9 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.0001 7.0000 1.5000 1.0000 0.5000 Variance 0.0000 0.0000 16.0000 0.2500 0.0000 0.2500 Minimum 2.4089 0.0001 3.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.0001 11.0000 2.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.0001 6.0000 1.5000 1.0000 0.5000 Variance 0.0000 0.0000 16.0000 0.2500 0.0000 0.2500 Minimum 2.4089 0.0001 2.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.0001 10.0000 2.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7878 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1927 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0195 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5866 2 0.5866 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9805 2 0.9805 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8876 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4389 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1397 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7878 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1927 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0195 MATRIX F This is a FULL matrix of order 1 by 1 1 1 11.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5866 2 0.5866 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9805 2 0.9805 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.0001 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.0001 2 0.0001 3 0.0002 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8876 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4389 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1397 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9805 1.0000 Function value of this group: 33.9553 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7878 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1927 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0195 MATRIX F This is a FULL matrix of order 1 by 1 1 1 2.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5866 2 0.5866 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9805 2 0.9805 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.0001 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.0001 2 0.0001 3 0.0002 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8876 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4389 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1397 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5866 1.0000 Function value of this group: 27.4488 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.7878 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1927 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.0195 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.5866 2 0.5866 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.9805 2 0.9805 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8876 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4389 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1397 Your model has 3 estimated parameters and 9 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 61.404 Degrees of freedom >>>>>>>>>>>>>>>> 6 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.7878 0.2538 0.9102 0 0 0 0 C 1 1 1 95.0 0.1927 0.0000 0.7164 6 5 0 0 E 1 1 1 95.0 0.0195 0.0100 0.0987 1 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 1.23 TOTAL 0: 0: 0: 1.29 Total number of warnings issued: 6 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 10 NOTE: Select if yields 8 data vectors for analysis NOTE: Vectors contain a total of 72 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 4 data vectors for analysis NOTE: Vectors contain a total of 24 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 4 data vectors for analysis NOTE: Vectors contain a total of 8 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! population sample, so no weight required ! Weight \part(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 10 NOTE: Select if yields 8 data vectors for analysis NOTE: Vectors contain a total of 72 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 4 data vectors for analysis NOTE: Vectors contain a total of 24 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 4 data vectors for analysis NOTE: Vectors contain a total of 8 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! population sample, so no weight required (could estimate threshold) ! Weight \part(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 Mean 2.0990 0.9999 623.7500 1.5000 0.5000 0.5000 Variance 0.0000 0.0000********** 0.7500 0.2500 0.2500 Minimum 2.0990 0.9999 20.0000 1.0000 0.0000 0.0000 Maximum 2.0990 0.9999 2428.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 Mean 2.0990 0.9999 1344.5000 1.5000 0.5000 0.5000 Variance 0.0000 0.0000********** 0.7500 0.2500 0.2500 Minimum 2.0990 0.9999 9.0000 1.0000 0.0000 0.0000 Maximum 2.0990 0.9999 5192.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 *** 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 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8239 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1761 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4119 2 0.4119 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8239 2 0.8239 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9077 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4197 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8239 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1761 MATRIX F This is a FULL matrix of order 1 by 1 1 1 20.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0278 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4119 2 0.4119 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8239 2 0.8239 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0990 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9077 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4197 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.0990 2.0990 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8239 1.0000 Function value of this group: 763.8755 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8239 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1761 MATRIX F This is a FULL matrix of order 1 by 1 1 1 9.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0337 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4119 2 0.4119 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8239 2 0.8239 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0990 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9077 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4197 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.0990 2.0990 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4119 1.0000 Function value of this group: 1935.3774 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8239 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1761 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4119 2 0.4119 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8239 2 0.8239 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9077 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4197 *** WARNING! *** Minimization may not be successful. See above CODE RED - Hessian/precision problem Your model has 3 estimated parameters and 17 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 2699.253 Degrees of freedom >>>>>>>>>>>>>>>> 14 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.8239 0.6047 0.8960 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.1816 0 0 0 0 E 1 1 1 95.0 0.1761 0.1040 0.2760 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 0.99 TOTAL 0: 0: 0: 1.05 Total number of warnings issued: 8 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 X LOWER 1 1 FREE Y LOWER 1 1 FREE Z LOWER 1 1 FREE O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND 0 1 X 1 1 1 Y 1 1 1 BOUND .1 1 Z 1 1 1 BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; L= A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; ! DZ COV MATRIX M= A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; ! MZ COV MATRIX END ALGEBRA: INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 11 NOTE: Select if yields 5 data vectors for analysis NOTE: Vectors contain a total of 45 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 11 NOTE: Select if yields 5 data vectors for analysis NOTE: Vectors contain a total of 45 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END The following MX script lines were read for group 4 GROUP 4: CONSTRAIN VARIANCE TO 1 CONSTRAINT NI=1 BEGIN MATRICES = GROUP 1 ; I UNIT 1 1 END MATRICES; CONSTRAINT I = A+C+E+D ; OPTION END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.2603 0.4670 16.6667 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 80.2222 0.6667 0.0000 0.2222 Minimum 2.2603 0.4670 8.0000 1.0000 1.0000 0.0000 Maximum 2.2603 0.4670 29.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.2603 0.4670 29.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 784.0000 1.0000 0.0000 0.2500 Minimum 2.2603 0.4670 1.0000 1.0000 1.0000 0.0000 Maximum 2.2603 0.4670 57.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 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 1 by 1 It has no free parameters specified MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 -3 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a computed FULL matrix of order 3 by 1 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 1 1 1 -4 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 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 a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 3 Mx starting optimization; number of parameters = 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8497 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1503 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4249 2 0.4249 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8497 2 0.8497 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9218 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3877 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8497 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1503 MATRIX F This is a FULL matrix of order 1 by 1 1 1 13.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0099 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4249 2 0.4249 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8497 2 0.8497 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.4670 2 0.4670 3 0.7159 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.2603 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9218 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3877 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.2603 2.2603 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8497 1.0000 Function value of this group: 115.2021 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8497 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1503 MATRIX F This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 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 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0108 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4249 2 0.4249 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8497 2 0.8497 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.4670 2 0.4670 3 0.7159 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.2603 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9218 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3877 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.2603 2.2603 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4249 1.0000 Function value of this group: 93.7344 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Group 4: constrain variance to 1 MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.8497 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX D This is a computed FULL matrix of order 1 by 1 [=W*W'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.1503 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX L This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|H@A+C+Q@D_H@A+C+Q@D|A+C+D+E] 1 2 1 1.0000 0.4249 2 0.4249 1.0000 MATRIX M This is a computed FULL matrix of order 2 by 2 [=A+C+D+E|A+C+D_A+C+D|A+C+D+E] 1 2 1 1.0000 0.8497 2 0.8497 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9218 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3877 Your model has 3 estimated parameters and 11 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 208.936 Degrees of freedom >>>>>>>>>>>>>>>> 8 3 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.8497 0.6356 0.9222 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.1902 0 0 0 0 E 1 1 1 95.0 0.1503 0.0778 0.2566 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 2.14 TOTAL 0: 0: 0: 2.20 Total number of warnings issued: 8 ______________________________________________________________________________ Note: repeat section completed ______________________________________________________________________________