** 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_MX005.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_MX005.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.5758 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.5758 0.0001 4.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.0001 4.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.7088 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2390 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.0522 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.5934 2 0.5934 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.9478 2 0.9478 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.5758 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.8419 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4889 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2284 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.7088 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2390 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.0522 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.5934 2 0.5934 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.9478 2 0.9478 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.5758 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.8419 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4889 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2284 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9478 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.7088 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2390 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.0522 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.5934 2 0.5934 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.9478 2 0.9478 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.5758 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.8419 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4889 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2284 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5934 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.7088 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2390 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.0522 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.5934 2 0.5934 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.9478 2 0.9478 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.5758 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.8419 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4889 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2284 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.7088 0.2535 0.8580 0 0 0 0 C 1 1 1 95.0 0.2390 0.0000 0.6294 0 0 0 0 E 1 1 1 95.0 0.0522 0.0064 0.1941 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.66 TOTAL 0: 0: 0: 0.72 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_MX005.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_MX005.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.5758 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.5758 0.1110 6.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.1110 3.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5387 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4436 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.0177 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.0011 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.0011 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.7130 2 0.7130 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.9823 2 0.9823 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.5758 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.7339 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6661 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1330 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.5387 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4436 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.0177 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.0011 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.0011 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.7130 2 0.7130 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.9823 2 0.9823 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.5758 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.7339 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6661 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1330 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9823 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.5387 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4436 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.0177 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.0011 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.0011 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.7130 2 0.7130 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.9823 2 0.9823 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.5758 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.7339 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6661 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1330 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.7130 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.5387 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.4436 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.0177 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.0011 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.0011 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.7130 2 0.7130 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.9823 2 0.9823 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.5758 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.7339 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6661 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1330 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.5387 0.3283 0.7972 0 0 0 0 C 1 1 1 95.0 0.4436 0.1872 0.6492 0 0 0 0 E 1 1 1 95.0 0.0177 0.0045 0.0499 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.05 Execution 0: 0: 0: 0.83 TOTAL 0: 0: 0: 0.88 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_MX005.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_MX005.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.5758 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.5758 0.4240 43.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 0.4240 172.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.000042554.8889 0.6667 0.0000 0.2222 Minimum 2.5758 0.4240 21.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.8509 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1305 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.0185 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.0035 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.0041 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.5560 2 0.5560 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.9815 2 0.9815 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.5758 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.9225 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3613 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1361 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.8509 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1305 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.0185 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.0035 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.0041 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.5560 2 0.5560 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.9815 2 0.9815 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.5758 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.9225 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3613 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1361 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9815 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.8509 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1305 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.0185 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.0035 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.0041 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.5560 2 0.5560 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.9815 2 0.9815 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.5758 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.9225 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3613 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1361 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5560 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.8509 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.1305 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.0185 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.0035 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.0041 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.5560 2 0.5560 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.9815 2 0.9815 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.5758 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.9225 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3613 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1361 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.8509 0.7362 0.9733 0 0 0 0 C 1 1 1 95.0 0.1305 0.0093 0.2436 0 0 0 0 E 1 1 1 95.0 0.0185 0.0103 0.0312 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: 0.79 TOTAL 0: 0: 0: 0.85 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_MX005.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_MX005.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.5758 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.5758 0.9999 31.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.9999 14.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.4291 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5136 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.0573 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.0069 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.0087 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.7282 2 0.7282 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.9427 2 0.9427 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.5758 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.6550 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7167 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2393 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.4291 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5136 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.0573 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.0069 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.0087 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.7282 2 0.7282 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.9427 2 0.9427 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.5758 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.6550 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7167 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2393 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9427 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.4291 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5136 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.0573 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.0069 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.0087 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.7282 2 0.7282 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.9427 2 0.9427 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.5758 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.6550 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7167 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2393 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.7282 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.4291 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5136 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.0573 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.0069 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.0087 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.7282 2 0.7282 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.9427 2 0.9427 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.5758 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.6550 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7167 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2393 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.4291 0.2230 0.6808 0 0 1 0 C 1 1 1 95.0 0.5136 0.2680 0.7056 0 0 0 0 E 1 1 1 95.0 0.0573 0.0295 0.1026 6 1 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.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_MX005.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_MX005.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.5758 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.5758 0.9999 9.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.9999 12.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.2264 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7509 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.0227 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.0062 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.0078 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.8641 2 0.8641 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.9773 2 0.9773 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.5758 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.4758 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8665 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1507 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.2264 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7509 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.0227 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.0062 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.0078 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.8641 2 0.8641 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.9773 2 0.9773 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.5758 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.4758 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8665 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1507 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9773 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.2264 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7509 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.0227 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.0062 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.0078 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.8641 2 0.8641 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.9773 2 0.9773 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.5758 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.4758 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8665 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1507 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8641 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.2264 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.7509 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.0227 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.0062 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.0078 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.8641 2 0.8641 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.9773 2 0.9773 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.5758 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.4758 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.8665 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1507 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.2264 0.0585 0.4618 0 0 0 0 C 1 1 1 95.0 0.7509 0.5192 0.8973 0 0 0 0 E 1 1 1 95.0 0.0227 0.0053 0.0720 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.74 TOTAL 0: 0: 0: 0.80 Total number of warnings issued: 0 ______________________________________________________________________________ ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=4 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 W LOWER 1 1 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_MX005.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_MX005.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.5758 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.5758 0.9999 7.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.9999 6.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.3455 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5509 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.1037 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.0075 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.0087 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.7236 2 0.7236 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.8963 2 0.8963 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.5758 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.5878 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7422 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3220 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.3455 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5509 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.1037 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.0075 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.0087 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.7236 2 0.7236 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.8963 2 0.8963 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.5758 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.5878 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7422 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3220 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8963 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.3455 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5509 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.1037 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.0075 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.0087 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.7236 2 0.7236 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.8963 2 0.8963 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.5758 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.5878 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7422 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3220 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.7236 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.3455 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5509 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.1037 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.0075 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.0087 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.7236 2 0.7236 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.8963 2 0.8963 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.5758 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.5878 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7422 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.3220 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.3455 0.0000 0.7747 0 0 0 0 C 1 1 1 95.0 0.5509 0.1471 0.8483 0 0 0 0 E 1 1 1 95.0 0.1037 0.0334 0.2525 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.58 TOTAL 0: 0: 0: 0.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_MX005.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_MX005.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.5758 0.9999 1977.6667 1.6667 0.6667 0.3333 Variance 0.0000 0.0000********** 0.8889 0.2222 0.2222 Minimum 2.5758 0.9999 60.0000 1.0000 0.0000 0.0000 Maximum 2.5758 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.5758 0.9999 2518.0000 1.6667 0.6667 0.3333 Variance 0.0000 0.0000********** 0.8889 0.2222 0.2222 Minimum 2.5758 0.9999 18.0000 1.0000 0.0000 0.0000 Maximum 2.5758 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.7364 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0677 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.1959 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.0096 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.4359 2 0.4359 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.8041 2 0.8041 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.5758 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.8582 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2602 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4426 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.7364 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0677 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.1959 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.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.0096 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.4359 2 0.4359 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.8041 2 0.8041 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.5758 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.8582 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2602 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4426 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8041 1.0000 Function value of this group: 2392.2205 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.7364 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0677 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.1959 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.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.0096 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.4359 2 0.4359 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.8041 2 0.8041 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.5758 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.8582 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2602 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4426 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4359 1.0000 Function value of this group: 3209.6591 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.7364 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0677 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.1959 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.0096 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.4359 2 0.4359 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.8041 2 0.8041 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.5758 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.8582 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2602 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4426 Your model has 3 estimated parameters and 13 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 5601.880 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.7364 0.5267 0.8503 0 1 0 0 C 1 1 1 95.0 0.0677 0.0000 0.2554 6 1 0 1 E 1 1 1 95.0 0.1959 0.1494 0.2505 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: 4.11 TOTAL 0: 0: 0: 4.17 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_MX005.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_MX005.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.5758 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.5758 0.9999 15.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.9999 1.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.9520 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.0480 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.0067 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.0096 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.4760 2 0.4760 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.9520 2 0.9520 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.5758 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.9757 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.2190 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.9520 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.0480 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.0067 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.0096 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.4760 2 0.4760 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.9520 2 0.9520 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.5758 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.9757 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.2190 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9520 1.0000 Function value of this group: 65.1239 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.9520 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.0480 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.0067 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.0096 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.4760 2 0.4760 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.9520 2 0.9520 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.5758 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.9757 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.2190 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4760 1.0000 Function value of this group: 46.1399 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.9520 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.0480 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.0067 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.0096 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.4760 2 0.4760 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.9520 2 0.9520 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.5758 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.9757 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.2190 Your model has 3 estimated parameters and 9 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 111.264 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.9520 0.4245 0.9834 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.5196 8 0 0 0 E 1 1 1 95.0 0.0480 0.0166 0.1147 0 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.05 Execution 0: 0: 0: 1.89 TOTAL 0: 0: 0: 1.94 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_MX005.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_MX005.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.5758 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.5758 0.0001 3.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.0001 2.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.7208 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2616 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.0175 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.6220 2 0.6220 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.9825 2 0.9825 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.5758 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.8490 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5115 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1324 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.7208 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2616 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.0175 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.6220 2 0.6220 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.9825 2 0.9825 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.5758 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.8490 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5115 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1324 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9825 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.7208 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2616 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.0175 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.6220 2 0.6220 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.9825 2 0.9825 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.5758 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.8490 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5115 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1324 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.6220 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.7208 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2616 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.0175 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.6220 2 0.6220 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.9825 2 0.9825 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.5758 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.8490 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5115 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.1324 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.7208 0.2299 0.9240 0 0 0 0 C 1 1 1 95.0 0.2616 0.0000 0.7435 0 0 1 0 E 1 1 1 95.0 0.0175 0.0013 0.0889 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.62 TOTAL 0: 0: 0: 0.68 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_MX005.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_MX005.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.5758 0.9999 623.7500 1.5000 0.5000 0.5000 Variance 0.0000 0.0000********** 0.7500 0.2500 0.2500 Minimum 2.5758 0.9999 20.0000 1.0000 0.0000 0.0000 Maximum 2.5758 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.5758 0.9999 1344.5000 1.5000 0.5000 0.5000 Variance 0.0000 0.0000********** 0.7500 0.2500 0.2500 Minimum 2.5758 0.9999 9.0000 1.0000 0.0000 0.0000 Maximum 2.5758 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.7964 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.2036 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.0097 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.3982 2 0.3982 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.7964 2 0.7964 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.5758 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.8924 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.4512 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.7964 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.2036 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.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.0097 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.3982 2 0.3982 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.7964 2 0.7964 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.5758 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.8924 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.4512 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.7964 1.0000 Function value of this group: 832.6716 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.7964 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.2036 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.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.0097 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.3982 2 0.3982 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.7964 2 0.7964 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.5758 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.8924 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.4512 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.3982 1.0000 Function value of this group: 2144.7892 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.7964 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.2036 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.0097 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.3982 2 0.3982 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.7964 2 0.7964 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.5758 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.8924 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.4512 Your model has 3 estimated parameters and 17 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 2977.461 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.7964 0.5344 0.8705 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.2150 8 0 0 1 E 1 1 1 95.0 0.2036 0.1295 0.2972 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.70 TOTAL 0: 0: 0: 2.76 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_MX005.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_MX005.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.5758 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.5758 0.4670 8.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.4670 1.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.8818 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.1182 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.0042 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.0046 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.4409 2 0.4409 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.8818 2 0.8818 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.5758 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.9391 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.3437 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.8818 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.1182 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.0042 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.0046 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.4409 2 0.4409 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.8818 2 0.8818 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.5758 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.9391 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.3437 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8818 1.0000 Function value of this group: 115.1166 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.8818 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.1182 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.0042 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.0046 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.4409 2 0.4409 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.8818 2 0.8818 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.5758 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.9391 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.3437 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4409 1.0000 Function value of this group: 91.6269 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.8818 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.1182 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.0042 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.0046 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.4409 2 0.4409 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.8818 2 0.8818 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.5758 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.9391 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.3437 Your model has 3 estimated parameters and 11 Observed statistics Observed statistics include 1 constraints. -2 times log-likelihood of data >>> 206.744 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.8818 0.6399 0.9393 0 0 0 0 C 1 1 1 95.0 0.0000 0.0000 0.2257 0 0 0 0 E 1 1 1 95.0 0.1182 0.0607 0.2045 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_MX005.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_MX005.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.5758 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.5758 0.2940 3.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.5758 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.5758 0.2940 2.0000 1.0000 1.0000 0.0000 Maximum 2.5758 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.4075 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5208 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.0717 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.0027 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.0028 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.7245 2 0.7245 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.9283 2 0.9283 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.5758 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.6383 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7217 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2678 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.4075 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5208 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.0717 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.0027 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.0028 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.7245 2 0.7245 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.9283 2 0.9283 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.5758 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.6383 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7217 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2678 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9283 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.4075 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5208 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.0717 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.0027 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.0028 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.7245 2 0.7245 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.9283 2 0.9283 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.5758 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.6383 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7217 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2678 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.5758 2.5758 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.7245 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.4075 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.5208 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.0717 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.0027 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.0028 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.7245 2 0.7245 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.9283 2 0.9283 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.5758 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.6383 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7217 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.2678 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.4075 0.0000 0.9411 0 0 6 2 C 1 1 1 95.0 0.5208 0.0047 0.7940 0 1 0 0 E 1 1 1 95.0 0.0717 0.0138 0.1248 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.39 TOTAL 0: 0: 0: 1.45 Total number of warnings issued: 0 ______________________________________________________________________________ Note: repeat section completed ______________________________________________________________________________