** 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=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 1 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 1 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.0001 5.5000 1.5000 1.0000 0.5000 Variance 0.0000 0.0000 2.2500 0.2500 0.0000 0.2500 Minimum 2.4089 0.0001 4.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.0001 7.0000 2.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.0001 13.5000 1.5000 1.0000 0.5000 Variance 0.0000 0.0000 90.2500 0.2500 0.0000 0.2500 Minimum 2.4089 0.0001 4.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.0001 23.0000 2.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5559 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9420 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 7.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5559 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9420 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.0001 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.0001 2 0.0001 3 0.0002 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9420 1.0000 Function value of this group: 29.6691 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 4.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5559 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9420 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.0001 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.0001 2 0.0001 3 0.0002 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5559 1.0000 Function value of this group: 60.0816 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 2 estimated parameters and 8 Observed statistics -2 times log-likelihood of data >>> 89.751 Degrees of freedom >>>>>>>>>>>>>>>> 6 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.9420 0.7858 0.9900 0 0 1 0 L 1 2 1 95.0 0.5559 0.2978 0.7573 0 0 1 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 0.10 TOTAL 0: 0: 0: 0.16 Total number of warnings issued: 0 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 2 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 2 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.1110 13.6667 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 77.5556 0.6667 0.0000 0.2222 Minimum 2.4089 0.1110 6.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.1110 26.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.1110 20.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 366.2222 0.6667 0.0000 0.2222 Minimum 2.4089 0.1110 3.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.1110 47.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.6846 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9813 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 6.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0017 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.6846 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9813 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.1110 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.1110 2 0.1110 3 0.2097 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9813 1.0000 Function value of this group: 88.7069 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 3.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0018 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.6846 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9813 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.1110 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.1110 2 0.1110 3 0.2097 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.6846 1.0000 Function value of this group: 144.8680 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 2 estimated parameters and 12 Observed statistics -2 times log-likelihood of data >>> 233.575 Degrees of freedom >>>>>>>>>>>>>>>> 10 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.9813 0.9471 0.9900 0 0 1 0 L 1 2 1 95.0 0.6846 0.5472 0.7957 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.39 TOTAL 0: 0: 0: 0.45 Total number of warnings issued: 0 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 3 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 3 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.4240 58.0000 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 88.6667 0.6667 0.0000 0.2222 Minimum 2.4089 0.4240 49.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.4240 71.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.4240 172.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.000042554.8889 0.6667 0.0000 0.2222 Minimum 2.4089 0.4240 21.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.4240 464.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5159 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9706 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 71.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0057 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5159 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9706 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4240 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.4240 2 0.4240 3 0.6682 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9706 1.0000 Function value of this group: 334.9821 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 32.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0066 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5159 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9706 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4240 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.4240 2 0.4240 3 0.6682 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5159 1.0000 Function value of this group: 1004.1585 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 2 estimated parameters and 12 Observed statistics -2 times log-likelihood of data >>> 1339.141 Degrees of freedom >>>>>>>>>>>>>>>> 10 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.9706 0.9535 0.9824 0 0 0 0 L 1 2 1 95.0 0.5159 0.4518 0.5761 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.77 TOTAL 0: 0: 0: 0.83 Total number of warnings issued: 0 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 4 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 4 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4066 0.9999 34.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 12.2500 1.0000 0.0000 0.2500 Minimum 2.4066 0.9999 31.0000 1.0000 1.0000 0.0000 Maximum 2.4066 0.9999 38.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4066 0.9999 48.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 1156.0000 1.0000 0.0000 0.2500 Minimum 2.4066 0.9999 14.0000 1.0000 1.0000 0.0000 Maximum 2.4066 0.9999 82.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7007 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9363 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 31.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0111 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7007 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9363 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4066 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4066 2.4066 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9363 1.0000 Function value of this group: 147.6221 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 14.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0141 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7007 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9363 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4066 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4066 2.4066 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.7007 1.0000 Function value of this group: 193.4354 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 2 estimated parameters and 8 Observed statistics -2 times log-likelihood of data >>> 341.058 Degrees of freedom >>>>>>>>>>>>>>>> 6 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.9363 0.8860 0.9671 0 0 0 0 L 1 2 1 95.0 0.7007 0.5695 0.8030 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.07 Execution 0: 0: 0: 0.10 TOTAL 0: 0: 0: 0.17 Total number of warnings issued: 2 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 5 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 5 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4135 0.9999 11.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 6.2500 1.0000 0.0000 0.2500 Minimum 2.4135 0.9999 9.0000 1.0000 1.0000 0.0000 Maximum 2.4135 0.9999 14.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4135 0.9999 21.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 90.2500 1.0000 0.0000 0.2500 Minimum 2.4135 0.9999 12.0000 1.0000 1.0000 0.0000 Maximum 2.4135 0.9999 31.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8500 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9748 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 14.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0098 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8500 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9748 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4135 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4135 2.4135 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9748 1.0000 Function value of this group: 43.2657 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 12.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0124 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8500 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9748 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4135 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4135 2.4135 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8500 1.0000 Function value of this group: 93.8934 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 2 estimated parameters and 8 Observed statistics -2 times log-likelihood of data >>> 137.159 Degrees of freedom >>>>>>>>>>>>>>>> 6 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.9748 0.9203 0.9900 0 0 1 0 L 1 2 1 95.0 0.8500 0.7246 0.9274 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.32 TOTAL 0: 0: 0: 0.38 Total number of warnings issued: 2 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 6 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 6 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 1.8250 0.9999 10.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 12.2500 1.0000 0.0000 0.2500 Minimum 1.8250 0.9999 7.0000 1.0000 1.0000 0.0000 Maximum 1.8250 0.9999 14.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 1.8250 0.9999 21.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 225.0000 1.0000 0.0000 0.2500 Minimum 1.8250 0.9999 6.0000 1.0000 1.0000 0.0000 Maximum 1.8250 0.9999 36.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5653 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8290 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 7.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0510 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5653 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8290 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 1.8250 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 1.8250 1.8250 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8290 1.0000 Function value of this group: 46.1417 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 6.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0595 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5653 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8290 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 1.8250 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 1.8250 1.8250 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5653 1.0000 Function value of this group: 84.3564 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 2 estimated parameters and 8 Observed statistics -2 times log-likelihood of data >>> 130.498 Degrees of freedom >>>>>>>>>>>>>>>> 6 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.8290 0.5999 0.9437 0 0 0 0 L 1 2 1 95.0 0.5653 0.2954 0.7664 1 0 0 0 This problem used 0.3% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.06 Execution 0: 0: 0: 0.08 TOTAL 0: 0: 0: 0.14 Total number of warnings issued: 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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 7 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! population sample, so no weight required ! Weight \part(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 7 NOTE: Select if yields 6 data vectors for analysis NOTE: Vectors contain a total of 54 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! population sample, so no weight required (could estimate threshold) ! Weight \part(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.0465 0.9999 1977.6667 1.6667 0.6667 0.3333 Variance 0.0000 0.0000********** 0.8889 0.2222 0.2222 Minimum 2.0465 0.9999 60.0000 1.0000 0.0000 0.0000 Maximum 2.0465 0.9999 5739.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.0465 0.9999 2518.0000 1.6667 0.6667 0.3333 Variance 0.0000 0.0000********** 0.8889 0.2222 0.2222 Minimum 2.0465 0.9999 18.0000 1.0000 0.0000 0.0000 Maximum 2.0465 0.9999 7277.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.4107 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8298 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 60.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0313 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.4107 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8298 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0465 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.0465 2.0465 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8298 1.0000 Function value of this group: 2134.9071 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 18.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0381 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.4107 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8298 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0465 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.0465 2.0465 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.4107 1.0000 Function value of this group: 2867.8783 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 2 estimated parameters and 12 Observed statistics -2 times log-likelihood of data >>> 5002.785 Degrees of freedom >>>>>>>>>>>>>>>> 10 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.8298 0.7737 0.8750 0 0 0 0 L 1 2 1 95.0 0.4107 0.2891 0.5205 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.14 TOTAL 0: 0: 0: 0.20 Total number of warnings issued: 6 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 8 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 8 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.9999 15.5000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 0.2500 1.0000 0.0000 0.2500 Minimum 2.4089 0.9999 15.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.9999 16.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.9999 14.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 169.0000 1.0000 0.0000 0.2500 Minimum 2.4089 0.9999 1.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.9999 27.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.3728 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9472 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 15.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0108 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.3728 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9472 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9472 1.0000 Function value of this group: 65.1236 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0154 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.3728 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9472 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.3728 1.0000 Function value of this group: 46.0582 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 2 estimated parameters and 8 Observed statistics -2 times log-likelihood of data >>> 111.182 Degrees of freedom >>>>>>>>>>>>>>>> 6 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.9472 0.8738 0.9816 0 0 0 0 L 1 2 1 95.0 0.3728 -0.0835 0.7052 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.09 TOTAL 0: 0: 0: 0.15 Total number of warnings issued: 6 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 9 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 9 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is not 1 so weight is needed WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.0001 7.0000 1.5000 1.0000 0.5000 Variance 0.0000 0.0000 16.0000 0.2500 0.0000 0.2500 Minimum 2.4089 0.0001 3.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.0001 11.0000 2.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.4089 0.0001 6.0000 1.5000 1.0000 0.5000 Variance 0.0000 0.0000 16.0000 0.2500 0.0000 0.2500 Minimum 2.4089 0.0001 2.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.0001 10.0000 2.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5866 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9805 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 11.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5866 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9805 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.0001 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.0001 2 0.0001 3 0.0002 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9805 1.0000 Function value of this group: 33.9553 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 2.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.5866 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9805 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.0001 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.0001 2 0.0001 3 0.0002 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.5866 1.0000 Function value of this group: 27.4488 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 2 estimated parameters and 8 Observed statistics -2 times log-likelihood of data >>> 61.404 Degrees of freedom >>>>>>>>>>>>>>>> 6 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.9805 0.9013 0.9900 0 0 0 0 L 1 2 1 95.0 0.5866 0.2036 0.8468 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.21 TOTAL 0: 0: 0: 0.27 Total number of warnings issued: 6 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 10 NOTE: Select if yields 8 data vectors for analysis NOTE: Vectors contain a total of 72 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 4 data vectors for analysis NOTE: Vectors contain a total of 24 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 4 data vectors for analysis NOTE: Vectors contain a total of 8 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! population sample, so no weight required ! Weight \part(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 10 NOTE: Select if yields 8 data vectors for analysis NOTE: Vectors contain a total of 72 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 4 data vectors for analysis NOTE: Vectors contain a total of 36 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 4 data vectors for analysis NOTE: Vectors contain a total of 24 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 4 data vectors for analysis NOTE: Vectors contain a total of 8 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! population sample, so no weight required (could estimate threshold) ! Weight \part(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 Mean 2.0990 0.9999 623.7500 1.5000 0.5000 0.5000 Variance 0.0000 0.0000********** 0.7500 0.2500 0.2500 Minimum 2.0990 0.9999 20.0000 1.0000 0.0000 0.0000 Maximum 2.0990 0.9999 2428.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 4.0000 4.0000 4.0000 4.0000 4.0000 4.0000 Mean 2.0990 0.9999 1344.5000 1.5000 0.5000 0.5000 Variance 0.0000 0.0000********** 0.7500 0.2500 0.2500 Minimum 2.0990 0.9999 9.0000 1.0000 0.0000 0.0000 Maximum 2.0990 0.9999 5192.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.3400 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8362 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 20.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0275 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.3400 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8362 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0990 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.0990 2.0990 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8362 1.0000 Function value of this group: 763.7952 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 9.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0342 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.3400 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8362 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.9999 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.9999 2 0.9999 3 1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0990 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.0990 2.0990 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.3400 1.0000 Function value of this group: 1934.4931 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 2 estimated parameters and 16 Observed statistics -2 times log-likelihood of data >>> 2698.288 Degrees of freedom >>>>>>>>>>>>>>>> 14 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.8362 0.7368 0.9047 0 0 0 0 L 1 2 1 95.0 0.3400 0.1777 0.4849 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.15 TOTAL 0: 0: 0: 0.21 Total number of warnings issued: 6 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 11 NOTE: Select if yields 5 data vectors for analysis NOTE: Vectors contain a total of 45 observations SELECT IF ZYGOSITY = 1; NOTE: Select if yields 3 data vectors for analysis NOTE: Vectors contain a total of 27 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 3 data vectors for analysis NOTE: Vectors contain a total of 18 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 3 data vectors for analysis NOTE: Vectors contain a total of 6 observations BEGIN MATRICES= GROUP 1; F FULL 1 1 ! FREQUENCY G FULL 1 4 ! TO SELECT ASCERTAINMENT PROB P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 11 NOTE: Select if yields 5 data vectors for analysis NOTE: Vectors contain a total of 45 observations SELECT IF ZYGOSITY = 2 NOTE: Select if yields 2 data vectors for analysis NOTE: Vectors contain a total of 18 observations SELECT TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ; DEFINITION ASC FREQ PIHAT THRESHOLD ; NOTE: Selection yields 2 data vectors for analysis NOTE: Vectors contain a total of 12 observations Note: Global variable previously defined. Updating ASC Note: Global variable previously defined. Updating FREQ Note: Global variable previously defined. Updating PIHAT Note: Global variable previously defined. Updating THRESHOLD NOTE: Definition yields 2 data vectors for analysis NOTE: Vectors contain a total of 4 observations BEGIN MATRICES= GROUP 1 ; F FULL 1 1 ! FREQUENCY G FULL 1 4 =G2 ! TO SELECT ASCERTAINMENT PROBABILITY P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.2603 0.4670 16.6667 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 80.2222 0.6667 0.0000 0.2222 Minimum 2.2603 0.4670 8.0000 1.0000 1.0000 0.0000 Maximum 2.2603 0.4670 29.0000 3.0000 1.0000 1.0000 Summary of VL file data for group 3 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Mean 2.2603 0.4670 29.0000 2.0000 1.0000 0.5000 Variance 0.0000 0.0000 784.0000 1.0000 0.0000 0.2500 Minimum 2.2603 0.4670 1.0000 1.0000 1.0000 0.0000 Maximum 2.2603 0.4670 57.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.1007 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8655 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 13.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0099 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.1007 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8655 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.4670 2 0.4670 3 0.7159 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.2603 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.2603 2.2603 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.8655 1.0000 Function value of this group: 115.0674 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0111 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.1007 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8655 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX S This is a computed FULL matrix of order 3 by 1 [=(P)_(P)_P+P-P*P] 1 1 0.4670 2 0.4670 3 0.7159 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.2603 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.2603 2.2603 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.1007 1.0000 Function value of this group: 89.1223 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 2 estimated parameters and 10 Observed statistics -2 times log-likelihood of data >>> 204.190 Degrees of freedom >>>>>>>>>>>>>>>> 8 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.8655 0.7670 0.9308 0 0 0 0 L 1 2 1 95.0 0.1007 -0.2554 0.4006 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.10 TOTAL 0: 0: 0: 0.16 Total number of warnings issued: 8 ______________________________________________________________________________ !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 ! ! Ascertainment correction, estimating pi ! The following MX script lines were read for group 1 G1: MODEL PARAMETERS DATA CALC NGROUPS=3 BEGIN MATRICES; H FULL 1 1 Q FULL 1 1 U UNIT 1 1 M STAN 2 2 FREE L STAN 2 2 FREE W LOWER 1 1 O ZERO 1 1 END MATRICES; MATRIX H .5 MATRIX Q .25 START .6 ALL BOUND -.99 .99 M 2 1 L 2 1 BEGIN ALGEBRA; ! A= X*X'; ! C= Y*Y'; ! E= Z*Z'; ! D= W*W'; ! L= A+C+D+E | H@A+C+Q@D _ ! H@A+C+Q@D | A+C+D+E ; ! DZ cov matrix ! M= A+C+D+E | A+C+D _ ! A+C+D | A+C+D+E ; ! MZ cov matrix END ALGEBRA: END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 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 P FULL 1 1 T FULL 1 1 END MATRICES; BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES J= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U); END ALGEBRA; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 COVARIANCES M; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%J ; OPTIONS RSIDUAL END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS DATA NINPUT=9 LABELS STUDYID STUDYYR ZYGOSITY TWIN1 TWIN2 ASC FREQ PIHAT THRESHOLD ORDINAL FILE=SCZ_MX.DAT Ordinal data read initiated NOTE: Rectangular file contained 61 records with data that contained a total of 549 observations SELECT IF STUDYID = 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 P FULL 1 1 T FULL 1 1 END MATRICES; SPECIFY F FREQ SPECIFY P PIHAT SPECIFY T THRESHOLD MATRIX G 1 1 1 1 SPECIFY G ASC 0 ASC 0 BEGIN ALGEBRA; S= (P) _ (P) _ P+P-P*P ; ! VECTOR OF ASCERTAINMENT PROBABILITIES K= (P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U); END ALGEBRA; COVARIANCES L; THRESHOLDS T|T ; FREQUENCY F ; ! pihat is 1 but -- not sampled WEIGHT \PART(S,G)%K ; ! Drop @0.9248 M 1 2 1 ! Drop @0.5216 L 1 2 1 INTERVALS M 1 2 1 L 1 2 1 OPTIONS RSIDUAL NDECIMALS=4 END Summary of VL file data for group 2 THRESHOLD PIHAT FREQ ASC TWIN1 TWIN2 Code -4.0000 -3.0000 -2.0000 -1.0000 1.0000 2.0000 Number 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 Mean 2.4089 0.2940 4.3333 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 1.5556 0.6667 0.0000 0.2222 Minimum 2.4089 0.2940 3.0000 1.0000 1.0000 0.0000 Maximum 2.4089 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.4089 0.2940 10.0000 2.0000 1.0000 0.6667 Variance 0.0000 0.0000 180.6667 0.6667 0.0000 0.2222 Minimum 2.4089 0.2940 0.0000 1.0000 1.0000 0.0000 Maximum 2.4089 0.2940 29.0000 3.0000 1.0000 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2: MZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 GROUP NUMBER: 3 G3: DZ twin pairs MATRIX F This is a FULL matrix of order 1 by 1 1 1 -2 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 0 -1 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 2 0 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 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 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Model parameters MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.2544 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9204 1.0000 MATRIX O This is a NULL matrix of order 1 by 1 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 0.2500 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 GROUP NUMBER: 2 G2: MZ twin pairs 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(U_O_T_T_U)-P*P*\MNOR(M_O|O_T|T_T|T_U|U)] 1 1 0.0043 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.2544 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9204 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.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.9204 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 F This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX G This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX K This is a computed FULL matrix of order 1 by 1 [=(P+P)*\MNOR(U_O_T_T_U)-P*P*\MNOR(L_O|O_T|T_T|T_U|U)] 1 1 0.0047 MATRIX L This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.2544 1.0000 MATRIX M This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.9204 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.4089 MATRIX U This is a UNIT matrix of order 1 by 1 MATRIX W This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 Matrix of EXPECTED thresholds TWIN1 TWIN2 Threshold 1 2.4089 2.4089 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX TWIN1 TWIN2 TWIN1 1.0000 TWIN2 0.2544 1.0000 Function value of this group: 50.0395 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 2 estimated parameters and 12 Observed statistics -2 times log-likelihood of data >>> 80.575 Degrees of freedom >>>>>>>>>>>>>>>> 10 2 Confidence intervals requested in group 3 Matrix Element Int. Estimate Lower Upper Lfail Ufail M 1 2 1 95.0 0.9204 0.7576 0.9846 0 0 0 0 L 1 2 1 95.0 0.2544 -0.1481 0.5745 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.11 TOTAL 0: 0: 0: 0.17 Total number of warnings issued: 8 ______________________________________________________________________________ !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 Note: repeat section completed ______________________________________________________________________________