** Mx startup successful ** **MX-Linux version 1.52b** #repeat 12 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 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.3263 0.0001 4.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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.3263 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.3263 0.0001 4.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.8063 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1325 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.0612 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_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.5356 2 0.5356 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.9388 2 0.9388 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.3263 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.8979 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3640 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2474 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.8063 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1325 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.0612 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_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.5356 2 0.5356 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.9388 2 0.9388 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.3263 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.8979 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3640 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2474 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9388 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.8063 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1325 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.0612 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_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.5356 2 0.5356 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.9388 2 0.9388 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.3263 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.8979 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3640 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2474 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5356 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.8063 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1325 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.0612 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_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.5356 2 0.5356 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.9388 2 0.9388 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.3263 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.8979 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3640 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2474 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 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.8063 0.2918 0.8072 0 2 0 0 C 1 1 1 95.0 0.1325 0.0000 0.5723 0 0 0 0 E 1 1 1 95.0 0.0612 0.0075 0.2253 1 0 1 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.07 Execution 0: 0: 0: 0.75 TOTAL 0: 0: 0: 0.82 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 0.1110 13.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 69.5556 0.6667 0.0000 0.2222 Minimum 2.3263 0.1110 6.0000 1.0000 1.0000 0.0000 Maximum 2.3263 0.1110 25.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.3263 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.3263 0.1110 3.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.6200 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3592 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.0208 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0021 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0022 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.6692 2 0.6692 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.9792 2 0.9792 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.3263 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.7874 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5993 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1442 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.6200 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3592 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.0208 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0021 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0022 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.6692 2 0.6692 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.9792 2 0.9792 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.3263 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.7874 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5993 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1442 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9792 1.0000 Function value of this group: 86.9321 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.6200 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3592 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.0208 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0021 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0022 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.6692 2 0.6692 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.9792 2 0.9792 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.3263 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.7874 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5993 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1442 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.6692 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.6200 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.3592 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.0208 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0021 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0022 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.6692 2 0.6692 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.9792 2 0.9792 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.3263 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.7874 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5993 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1442 Your model has 3 estimated parameters and 13 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 231.800 Degrees of freedom >>>>>>>>>>>>>>>> 10 3 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.6200 0.3804 0.9090 0 0 0 0 C 1 1 1 95.0 0.3592 0.0727 0.5929 0 0 0 0 E 1 1 1 95.0 0.0208 0.0053 0.0585 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.71 TOTAL 0: 0: 0: 0.77 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 0.4240 54.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 144.8889 0.6667 0.0000 0.2222 Minimum 2.3263 0.4240 43.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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.3263 0.4240 172.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.000042554.8889 0.6667 0.0000 0.2222 Minimum 2.3263 0.4240 21.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.9677 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0105 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.0218 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0071 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0083 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.4944 2 0.4944 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.9782 2 0.9782 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.3263 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.9837 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1025 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1476 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.9677 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0105 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.0218 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0071 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0083 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.4944 2 0.4944 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.9782 2 0.9782 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.3263 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.9837 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1025 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1476 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9782 1.0000 Function value of this group: 292.2912 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.9677 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0105 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.0218 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0071 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0083 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.4944 2 0.4944 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.9782 2 0.9782 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.3263 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.9837 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1025 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1476 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4944 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.9677 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0105 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.0218 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0071 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0083 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.4944 2 0.4944 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.9782 2 0.9782 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.3263 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.9837 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1025 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1476 Your model has 3 estimated parameters and 13 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 1296.450 Degrees of freedom >>>>>>>>>>>>>>>> 10 3 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.9677 0.8407 0.9879 0 0 0 0 C 1 1 1 95.0 0.0105 0.0000 0.1356 0 0 0 0 E 1 1 1 95.0 0.0218 0.0121 0.0366 6 1 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: 1.09 TOTAL 0: 0: 0: 1.15 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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 non-random 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but non-random so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 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.3263 0.9999 31.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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.3263 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.3263 0.9999 14.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.4930 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4399 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.0671 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0138 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0175 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.6864 2 0.6864 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.9329 2 0.9329 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.3263 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.7021 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6632 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2591 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.4930 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4399 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.0671 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0138 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0175 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.6864 2 0.6864 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.9329 2 0.9329 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.3263 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.7021 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6632 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2591 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9329 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.4930 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4399 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.0671 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0138 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0175 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.6864 2 0.6864 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.9329 2 0.9329 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.3263 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.7021 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6632 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2591 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.6864 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.4930 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4399 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.0671 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0138 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0175 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.6864 2 0.6864 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.9329 2 0.9329 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.3263 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.7021 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6632 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2591 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 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.4930 0.2574 0.7757 0 0 0 0 C 1 1 1 95.0 0.4399 0.1646 0.6586 0 0 0 0 E 1 1 1 95.0 0.0671 0.0347 0.1199 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: 0.86 TOTAL 0: 0: 0: 0.92 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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 non-random 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but non-random so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 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.3263 0.9999 9.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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.3263 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.3263 0.9999 12.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.2634 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7099 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.0267 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0124 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0156 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.8416 2 0.8416 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.9733 2 0.9733 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.3263 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.5133 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8425 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1634 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.2634 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7099 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.0267 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0124 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0156 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.8416 2 0.8416 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.9733 2 0.9733 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.3263 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.5133 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8425 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1634 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9733 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.2634 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7099 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.0267 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0124 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0156 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.8416 2 0.8416 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.9733 2 0.9733 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.3263 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.5133 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8425 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1634 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8416 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.2634 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7099 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.0267 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0124 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0156 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.8416 2 0.8416 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.9733 2 0.9733 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.3263 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.5133 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8425 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1634 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 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.2634 0.0682 0.5331 0 0 1 0 C 1 1 1 95.0 0.7099 0.4447 0.8798 0 0 0 0 E 1 1 1 95.0 0.0267 0.0062 0.0843 1 0 1 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.69 TOTAL 0: 0: 0: 0.75 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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 non-random 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but non-random so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 0.9999 10.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 12.2500 1.0000 0.0000 0.2500 Minimum 2.3263 0.9999 7.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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.3263 0.9999 21.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 225.0000 1.0000 0.0000 0.2500 Minimum 2.3263 0.9999 6.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.3954 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4835 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.1211 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0150 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0175 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.6812 2 0.6812 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.8789 2 0.8789 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.3263 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.6288 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6953 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3480 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.3954 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4835 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.1211 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0150 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0175 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.6812 2 0.6812 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.8789 2 0.8789 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.3263 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.6288 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6953 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3480 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8789 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.3954 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4835 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.1211 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0150 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0175 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.6812 2 0.6812 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.8789 2 0.8789 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.3263 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.6288 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6953 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3480 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.6812 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.3954 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4835 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.1211 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0150 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0175 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.6812 2 0.6812 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.8789 2 0.8789 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.3263 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.6288 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6953 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3480 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 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.3954 0.0000 0.8744 0 0 0 2 C 1 1 1 95.0 0.4835 0.0352 0.8234 0 1 0 0 E 1 1 1 95.0 0.1211 0.0393 0.2919 0 0 1 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.17 TOTAL 0: 0: 0: 1.23 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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 ; ! random sample so no 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! random sample so no weight is needed ! Weight \part(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 0.9999 1977.6667 1.6667 0.6667 0.3333 Variance 0.0000 0.0000********** 0.8889 0.2222 0.2222 Minimum 2.3263 0.9999 60.0000 1.0000 0.0000 0.0000 Maximum 2.3263 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.3263 0.9999 2518.0000 1.6667 0.6667 0.3333 Variance 0.0000 0.0000********** 0.8889 0.2222 0.2222 Minimum 2.3263 0.9999 18.0000 1.0000 0.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.8021 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.1979 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0162 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0191 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.4011 2 0.4011 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.8021 2 0.8021 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.3263 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.8956 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.4448 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.8021 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.1979 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0162 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0191 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.4011 2 0.4011 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.8021 2 0.8021 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.3263 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.8956 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.4448 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8021 1.0000 Function value of this group: 2218.4961 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.8021 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.1979 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0162 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0191 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.4011 2 0.4011 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.8021 2 0.8021 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.3263 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.8956 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.4448 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4011 1.0000 Function value of this group: 2967.4342 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.8021 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.1979 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0162 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0191 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.4011 2 0.4011 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.8021 2 0.8021 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.3263 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.8956 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.4448 Your model has 3 estimated parameters and 13 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 5185.930 Degrees of freedom >>>>>>>>>>>>>>>> 10 3 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.8021 0.5790 0.8498 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.2006 0 0 6 1 E 1 1 1 95.0 0.1979 0.1502 0.2558 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: 2.58 TOTAL 0: 0: 0: 2.64 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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 non-random 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but non-random so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 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.3263 0.9999 15.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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.3263 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.3263 0.9999 1.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.9430 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.0570 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0135 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0188 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.4715 2 0.4715 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.9430 2 0.9430 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.3263 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.9711 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.2388 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.9430 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.0570 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0135 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0188 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.4715 2 0.4715 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.9430 2 0.9430 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.3263 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.9711 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.2388 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9430 1.0000 Function value of this group: 65.1256 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.9430 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.0570 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0135 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0188 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.4715 2 0.4715 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.9430 2 0.9430 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.3263 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.9711 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.2388 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4715 1.0000 Function value of this group: 46.4340 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.9430 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.0570 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0135 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0188 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.4715 2 0.4715 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.9430 2 0.9430 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.3263 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.9711 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.2388 Your model has 3 estimated parameters and 9 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 111.560 Degrees of freedom >>>>>>>>>>>>>>>> 6 3 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.9430 0.4697 0.9802 6 2 0 0 C 1 1 1 95.0 0.0000 0.0000 0.4639 0 0 0 1 E 1 1 1 95.0 0.0570 0.0198 0.1358 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.13 TOTAL 0: 0: 0: 1.19 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 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.3263 0.0001 3.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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.3263 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.3263 0.0001 2.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.8240 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1554 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.0206 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_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.5674 2 0.5674 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.9794 2 0.9794 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.3263 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.3942 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1436 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.8240 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1554 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.0206 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_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.5674 2 0.5674 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.9794 2 0.9794 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.3263 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.3942 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1436 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9794 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.8240 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1554 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.0206 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_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.5674 2 0.5674 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.9794 2 0.9794 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.3263 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.3942 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1436 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5674 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.8240 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1554 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.0206 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_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.5674 2 0.5674 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.9794 2 0.9794 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.3263 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.3942 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1436 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 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.8240 0.2669 0.9984 0 0 0 0 C 1 1 1 95.0 0.1554 0.0000 0.7015 0 0 0 0 E 1 1 1 95.0 0.0206 0.0206 0.1040 8 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.73 TOTAL 0: 0: 0: 0.79 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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 ; ! random sample so no 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! random sample so no weight is needed ! Weight \part(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 0.9999 623.7500 1.5000 0.5000 0.5000 Variance 0.0000 0.0000********** 0.7500 0.2500 0.2500 Minimum 2.3263 0.9999 20.0000 1.0000 0.0000 0.0000 Maximum 2.3263 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.3263 0.9999 1344.5000 1.5000 0.5000 0.5000 Variance 0.0000 0.0000********** 0.7500 0.2500 0.2500 Minimum 2.3263 0.9999 9.0000 1.0000 0.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.7924 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.2076 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0163 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0191 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.3962 2 0.3962 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.7924 2 0.7924 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.3263 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.8901 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.4557 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.7924 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.2076 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0163 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0191 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.3962 2 0.3962 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.7924 2 0.7924 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.3263 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.8901 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.4557 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.7924 1.0000 Function value of this group: 780.1262 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.7924 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.2076 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0163 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0191 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.3962 2 0.3962 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.7924 2 0.7924 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.3263 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.8901 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.4557 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.3962 1.0000 Function value of this group: 1989.7380 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.7924 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.2076 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0163 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0191 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.3962 2 0.3962 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.7924 2 0.7924 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.3263 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.8901 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.4557 Your model has 3 estimated parameters and 17 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 2769.864 Degrees of freedom >>>>>>>>>>>>>>>> 14 3 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.7924 0.5717 0.8714 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.1765 1 0 0 0 E 1 1 1 95.0 0.2076 0.1286 0.3099 0 0 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.90 TOTAL 0: 0: 0: 2.96 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.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 END MATRICES; 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 non-random 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_MX01.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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but non-random so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 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.3263 0.4670 8.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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.3263 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.3263 0.4670 1.0000 1.0000 1.0000 0.0000 Maximum 2.3263 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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.8574 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.1426 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0083 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0091 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.4287 2 0.4287 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.8574 2 0.8574 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.3263 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.9260 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.3776 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.8574 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.1426 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0083 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0091 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.4287 2 0.4287 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.8574 2 0.8574 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.3263 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.9260 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.3776 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8574 1.0000 Function value of this group: 115.1778 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.8574 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.1426 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0083 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0091 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.4287 2 0.4287 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.8574 2 0.8574 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.3263 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.9260 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.3776 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4287 1.0000 Function value of this group: 93.2325 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.8574 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.1426 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0083 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0091 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.4287 2 0.4287 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.8574 2 0.8574 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.3263 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.9260 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.3776 Your model has 3 estimated parameters and 11 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 208.410 Degrees of freedom >>>>>>>>>>>>>>>> 8 3 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.8574 0.6388 0.9263 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.1967 0 0 0 0 E 1 1 1 95.0 0.1426 0.0737 0.2441 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.28 TOTAL 0: 0: 0: 1.34 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 O ZERO 1 1 P FULL 1 1 ! ASCERTAINMENT PROBABILITY P T FULL 1 1 ! THRESHOLD X LOWER 1 1 FREE ! A Y LOWER 1 1 FREE ! C Z LOWER 1 1 FREE ! E END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BEGIN ALGEBRA; A= X*X'; C= Y*Y'; E= Z*Z'; D= W*W'; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES 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 J= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); K= (P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA: 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_MX01.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 12 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 END MATRICES; 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 non-random 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_MX01.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 12 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 END MATRICES; SPECIFY F FREQ MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but non-random so weight is needed WEIGHT \PART(S,G)%K ; INTERVALS A 1 1 1 C 1 1 1 E 1 1 1 OPTION JIGGLE OPTION FUNC=1.E-12 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 Multiple issat nag=10 db=1 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.3263 0.2940 4.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 1.5556 0.6667 0.0000 0.2222 Minimum 2.3263 0.2940 3.0000 1.0000 1.0000 0.0000 Maximum 2.3263 0.2940 6.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.3263 0.2940 6.0000 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 32.0000 0.6667 0.0000 0.2222 Minimum 2.3263 0.2940 2.0000 1.0000 1.0000 0.0000 Maximum 2.3263 0.2940 14.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 J This is a computed 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: 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 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: 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 J This is a computed 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 J This is a computed 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 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.4675 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4485 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.0840 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0054 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_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.6823 2 0.6823 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.9160 2 0.9160 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.2940 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.2940 2 0.2940 3 0.5016 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.3263 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.6837 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6697 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2898 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.4675 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4485 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.0840 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0054 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_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.6823 2 0.6823 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.9160 2 0.9160 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.2940 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.2940 2 0.2940 3 0.5016 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.3263 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.6837 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6697 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2898 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9160 1.0000 Function value of this group: 30.5358 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.4675 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4485 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.0840 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 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(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0054 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_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.6823 2 0.6823 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.9160 2 0.9160 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.2940 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.2940 2 0.2940 3 0.5016 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.3263 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.6837 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6697 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2898 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.3263 2.3263 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.6823 1.0000 Function value of this group: 40.6142 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.4675 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4485 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.0840 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 J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0054 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(A+C+D+E_O_T_T_U)-P*P*\MNOR(L_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.6823 2 0.6823 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.9160 2 0.9160 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.2940 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.2940 2 0.2940 3 0.5016 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.3263 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.6837 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6697 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2898 Your model has 3 estimated parameters and 13 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 71.150 Degrees of freedom >>>>>>>>>>>>>>>> 10 3 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail A 1 1 1 95.0 0.4675 0.0000 0.8591 0 0 0 0 C 1 1 1 95.0 0.4485 0.0000 0.8551 0 0 0 1 E 1 1 1 95.0 0.0840 0.0162 0.2547 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: 3.18 TOTAL 0: 0: 0: 3.24 Total number of warnings issued: 0 ______________________________________________________________________________ Note: repeat section completed ______________________________________________________________________________