The following MX script lines were read for group 1 G1 CALCULATION GROUP CONTROL CORRELATION MATRIX A FACTOR DA CALC NG=6 MATRICES A LO 1 1 FIXED C LO 1 1 FREE H FU 1 1 I ID 1 1 BEGIN ALGEBRA; X = A*A' ; Y = C*C' ; END ALGEBRA; COMPUTE I|H@X+Y_ H@X+Y|I; MATRIX A 0 MATRIX H .5 OPTION RS END The following MX script lines were read for group 2 G2 CALCULATION GROUP CLINICAL CORRELATION MATRIX A FACTOR DA CALC MATRICES A LO 1 1 =A1 C LO 1 1 =C1 H FU 1 1 I ID 1 1 BEGIN ALGEBRA; X = A*A' ; Y = C*C' ; END ALGEBRA ; COMPUTE I|H@X+Y_ H@X+Y|I ; MATRIX H .5 OPTION RS END The following MX script lines were read for group 3 G3 CALCULATION GROUP CONTROL CORRELATION MATRIX B FACTOR DA CALC MATRICES A LO 1 1 FIXED C LO 1 1 FREE H FU 1 1 I ID 1 1 BEGIN ALGEBRA; X = A*A' ; Y = C*C' ; END ALGEBRA; COMPUTE I|H@X+Y_ H@X+Y|I; MATRIX H .5 MATRIX A 0 OPTION RS END The following MX script lines were read for group 4 G4 CALCULATION GROUP CLINICAL CORRELATION MATRIX B FACTOR DA CALC MATRICES A LO 1 1 =A3 C LO 1 1 =C3 H FU 1 1 I ID 1 1 BEGIN ALGEBRA; X = A*A' ; Y = C*C' ; END ALGEBRA ; COMPUTE I|H@X+Y_ H@X+Y|I ; MATRIX H .5 OPTION RS END The following MX script lines were read for group 5 EPIPHENOMENA MODEL - FIT MODEL TO CONTROL A/D DATA DATA NI=6 ORDINAL FI=CONTROL.DAT Ordinal data read initiated NOTE: Rectangular file contained 422 records with data that contained a total of 2532 observations LABELS AGE1 AGE2 SEX1 SEX2 TYPE DUMMY DEFINITION_VARIABLES AGE1 AGE2 SEX1 SEX2 TYPE / NOTE: Definition yields 422 data vectors for analysis NOTE: Vectors contain a total of 422 observations MATRICES A FULL 2 2 =%E1 ! CORR BETWEEN A FACTORS B FULL 2 2 =%E3 ! CORR BETWEEN B FACTORS P FULL 1 1 FREE ! PROBABILITY OF BEING COMORBID GIVEN A R FULL 1 1 FREE ! PROBABILITY OF BEING COMORBID GIVEN B O FULL 1 2 ! MATRIX FOR AGE DEFINITION VARIABLE M FULL 1 2 ! MATRIX FOR SEX DEFINITION VARIABLE T FULL 1 2 ! ESTIMATED THRESHOLD FOR A Q FULL 1 2 ! AGE DIFFERENCE IN THRESHOLD FOR A X FULL 1 2 ! SEX DIFFERENCE IN THRESHOLD FOR A U FULL 1 2 ! ESTIMATED THRESHOLD FOR B S FULL 1 2 ! AGE DIFFERENCE IN THRESHOLD FOR B Y FULL 1 2 ! SEX DIFFERENCE IN THRESHOLD FOR B V FULL 1 1 ! VARIANCE OF DUMMY VARIABLE L FULL 4 1 ! FOR EXTRACTING RIGHT PART OF K W FULL 1 1 ! 0 + + THESE ARE TO CONTROL INTEGRAL TYPE Z FULL 1 1 ! 1 - - BEGIN ALGEBRA; C = \MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z)) ; ! LOWER/LOWER A Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 D = \MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W)) ; ! LOWER/UPPER A Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 E = \MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z)) ; ! UPPER/LOWER A Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 F = \MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W)) ; ! UPPER/UPPER A Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 G = \MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z)) ; ! LOWER/LOWER B Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 H = \MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W)) ; ! LOWER/UPPER B Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 I = \MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z)) ; ! UPPER/LOWER B Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 J = \MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W)) ; ! UPPER/UPPER B Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 K = C.G_ ! 1 (EXPECTED PROBABILITY OF 0000) C.I.(Z-R)_ ! 2 (EXPECTED PROBABILITY OF 0100) C.H.(Z-R)_ ! 3 (EXPECTED PROBABILITY OF 0001) E.G.(Z-P)_ ! 4 (EXPECTED PROBABILITY OF 1000) D.G.(Z-P)_ ! 5 (EXPECTED PROBABILITY OF 0010) E.I + E.G.P + C.I.R_ ! 6 (EXPECTED PROBABILITY OF 1100) D.H + D.G.P + C.H.R_ ! 7 (EXPECTED PROBABILITY OF 0011) C.J.(Z-R).(Z-R)_ ! 8 (EXPECTED PROBABILITY OF 0101) E.H.(Z-P).(Z-R)_ ! 9 (EXPECTED PROBABILITY OF 1001) D.I.(Z-P).(Z-R)_ ! 10 (EXPECTED PROBABILITY OF 0110) E.J.(Z-R) + E.H.P.(Z-R) + C.J.R.(Z-R)_ ! 11 (EXPECTED PROBABILITY OF 1101) D.J.(Z-R) + D.I.P.(Z-R) + C.J.R.(Z-R)_ ! 12 (EXPECTED PROBABILITY OF 0111) F.G.(Z-P).(Z-P)_ ! 13 (EXPECTED PROBABILITY OF 1010) F.I.(Z-P) + F.G.P.(Z-P) + D.I.R.(Z-P)_ ! 14 (EXPECTED PROBABILITY OF 1110) F.H.(Z-P) + F.G.P.(Z-P) + E.H.R.(Z-P)_ ! 15 (EXPECTED PROBABILITY OF 1011) F.J + F.H.P + F.I.P + F.G.P.P + D.J.R + D.I.P.R + E.J.R + E.H.P.R + C.J.R.R; ! 16 (EXPECTED PROBABILITY OF 1111) END ALGEBRA; MATRIX L 1 1 1 1 MATRIX V .159154943 MATRIX P .5 MATRIX R .5 MATRIX W 0 MATRIX Z 1 SPECIFY T 100 100 SPECIFY Q 101 101 SPECIFY X 102 102 SPECIFY U 103 103 SPECIFY S 104 104 SPECIFY Y 105 105 MATRIX T 1.5 1.5 MATRIX U 1.5 1.5 SPECIFY O AGE1 AGE2; SPECIFY M SEX1 SEX2; SPECIFY L -5 0 -5 0; MEANS W; COVARIANCE V; WEIGHT (\PART(K,L)) / BOUND 0 1 3 4 OPTION RS END The following MX script lines were read for group 6 EPIPHENOMENA MODEL - FIT MODEL TO CLINICAL A/D DATA DATA NI=6 ORDINAL FI=CLINICAL.DAT Ordinal data read initiated NOTE: Rectangular file contained 253 records with data that contained a total of 1518 observations LABELS AGE1 AGE2 SEX1 SEX2 TYPE DUMMY DEFINITION_VARIABLES AGE1 AGE2 SEX1 SEX2 TYPE / Note: Global variable previously defined. Updating AGE1 Note: Global variable previously defined. Updating AGE2 Note: Global variable previously defined. Updating SEX1 Note: Global variable previously defined. Updating SEX2 Note: Global variable previously defined. Updating TYPE NOTE: Definition yields 253 data vectors for analysis NOTE: Vectors contain a total of 253 observations MATRICES A FULL 2 2 =%E2 ! CORR BETWEEN A FACTORS B FULL 2 2 =%E4 ! CORR BETWEEN B FACTORS P FULL 1 1 =P5 ! PROBABILITY OF BEING COMORBID GIVEN A R FULL 1 1 =R5 ! PROBABILITY OF BEING COMORBID GIVEN B O FULL 1 2 ! MATRIX FOR AGE DEFINITION VARIABLE M FULL 1 2 ! MATRIX FOR SEX DEFINITION VARIABLE T FULL 1 2 =T5! ESTIMATED THRESHOLD FOR A Q FULL 1 2 =Q5! AGE DIFFERENCE IN THRESHOLD FOR A X FULL 1 2 =X5! SEX DIFFERENCE IN THRESHOLD FOR A U FULL 1 2 =U5! ESTIMATED THRESHOLD FOR B S FULL 1 2 =S5! AGE DIFFERENCE IN THRESHOLD FOR B Y FULL 1 2 =Y5! SEX DIFFERENCE IN THRESHOLD FOR B V FULL 1 1 ! VARIANCE OF DUMMY VARIABLE L FULL 4 1 ! FOR EXTRACTING RIGHT PART OF K W FULL 1 1 ! 0 + + THESE ARE TO CONTROL INTEGRAL TYPE Z FULL 1 1 ! 1 - - ! ALSO ASCERTAINMENT PARAMETER FOR ILLICIT DRUG DEPENDENCE (FIXED TO 1) N FULL 1 1 FREE ! ASCERTAINMENT PARAMETER FOR ALCOHOL DEPENDENCE BEGIN ALGEBRA; C = \MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z)) ; ! LOWER/LOWER A Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 D = \MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W)) ; ! LOWER/UPPER A Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 E = \MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z)) ; ! UPPER/LOWER A Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 F = \MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W)) ; ! UPPER/UPPER A Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 G = \MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z)) ; ! LOWER/LOWER B Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 H = \MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W)) ; ! LOWER/UPPER B Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 I = \MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z)) ; ! UPPER/LOWER B Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 J = \MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W)) ; ! UPPER/UPPER B Problem with correlation between variables 2 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 K = (Z).(C.I.(Z-R))_ ! 2 (EXPECTED PROBABILITY OF 0100) (N).(E.G.(Z-P))_ ! 4 (EXPECTED PROBABILITY OF 1000) (Z+N).(E.I + E.G.P + C.I.R)_ ! 6 (EXPECTED PROBABILITY OF 1100) (Z).(C.J.(Z-R).(Z-R))_ ! 8 (EXPECTED PROBABILITY OF 0101) (N).(E.H.(Z-P).(Z-R))_ ! 9 (EXPECTED PROBABILITY OF 1001) (Z).(D.I.(Z-P).(Z-R))_ ! 10 (EXPECTED PROBABILITY OF 0110) (Z+N).(E.J.(Z-R) + E.H.P.(Z-R) + C.J.R.(Z-R))_ ! 11 (EXPECTED PROBABILITY OF 1101) (Z).(D.J.(Z-R) + D.I.P.(Z-R) + C.J.R.(Z-R))_ ! 12 (EXPECTED PROBABILITY OF 0111) (N).(F.G.(Z-P).(Z-P))_ ! 13 (EXPECTED PROBABILITY OF 1010) (Z+N).(F.I.(Z-P) + F.G.P.(Z-P) + D.I.R.(Z-P))_ ! 14 (EXPECTED PROBABILITY OF 1110) (N).(F.H.(Z-P) + F.G.P.(Z-P) + E.H.R.(Z-P))_ ! 15 (EXPECTED PROBABILITY OF 1011) (Z+N).(F.J + F.H.P + F.I.P + F.G.P.P + D.J.R + D.I.P.R + E.J.R + E.H.P.R + C.J.R.R); ! 16 (EXPECTED PROBABILITY OF 1111) END ALGEBRA; MATRIX L 1 1 1 1 MATRIX V .159154943 MATRIX P .5 MATRIX R .5 MATRIX W 0 MATRIX Z 1 MATRIX N .5 SPECIFY O AGE1 AGE2; SPECIFY M SEX1 SEX2; SPECIFY L -5 0 -5 0; MEANS W; COVARIANCE V; WEIGHT (\PART((K%( \SUM(K)_\SUM(K)_\SUM(K)_\SUM(K)_\SUM(K)_\SUM(K)_\SUM(K)_\SUM(K)_\SUM(K)_\SUM(K)_\SUM(K)_\SUM(K))),L)) / OPTION FUNC=1.E-7 ND=7 !DIFF=.001 OPTION RS MU CI 90 OPTION TH=-3 END Summary of VL file data for group 5 TYPE SEX2 SEX1 AGE2 AGE1 Code -5.0000000 -4.0000000 -3.0000000 -2.0000000 -1.0000000 Number 422.0000000 422.0000000 422.0000000 422.0000000 422.0000000 Mean 1.6303318 0.5497630 0.9194313 16.6899526 16.4779858 Variance 3.8017340 0.2475236 0.0740774 11.3756730 2.2599095 Minimum 1.0000000 0.0000000 0.0000000 10.9600000 11.9600000 Maximum 16.0000000 1.0000000 1.0000000 24.9500000 21.0500000 DUMMY Code 1.0000000 Number 422.0000000 Mean 1.0000000 Variance 0.0000000 Minimum 1.0000000 Maximum 1.0000000 Summary of VL file data for group 6 TYPE SEX2 SEX1 AGE2 AGE1 Code -5.0000000 -4.0000000 -3.0000000 -2.0000000 -1.0000000 Number 253.0000000 253.0000000 253.0000000 253.0000000 253.0000000 Mean 3.1778656 0.5375494 0.9288538 17.0966403 16.4343874 Variance 7.7825931 0.2485900 0.0660845 10.7494634 1.4493661 Minimum 1.0000000 0.0000000 0.0000000 11.1500000 13.0900000 Maximum 12.0000000 1.0000000 1.0000000 24.5500000 19.7100000 DUMMY Code 1.0000000 Number 253.0000000 Mean 1.0000000 Variance 0.0000000 Minimum 1.0000000 Maximum 1.0000000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1 CALCULATION GROUP Control Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a computed FULL matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a computed FULL matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 3 G3 CALCULATION GROUP Control Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a computed FULL matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 4 G4 CALCULATION GROUP Clinical Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Y This is a computed FULL matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 5 Epiphenomena model - fit model to Control A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 It has no free parameters specified MATRIX B This is a constrained FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX G This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX H This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 16 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 4 by 1 1 1 -5 2 0 3 -5 4 0 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 -3 -4 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 -1 -2 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 2 1 2 1 101 101 MATRIX R This is a FULL matrix of order 1 by 1 1 1 4 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 104 104 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 100 100 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 103 103 MATRIX V This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 102 102 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 105 105 MATRIX Z This is a FULL matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 6 Epiphenomena model - fit model to Clinical A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 It has no free parameters specified MATRIX B This is a constrained FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX D This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX F This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX G This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX H This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX K This is a computed FULL matrix of order 12 by 1 It has no free parameters specified MATRIX L This is a FULL matrix of order 4 by 1 1 1 -5 2 0 3 -5 4 0 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 -3 -4 MATRIX N This is a FULL matrix of order 1 by 1 1 1 106 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 -1 -2 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 2 1 2 1 101 101 MATRIX R This is a FULL matrix of order 1 by 1 1 1 4 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 104 104 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 100 100 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 103 103 MATRIX V This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX W This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 102 102 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 105 105 MATRIX Z This is a FULL matrix of order 1 by 1 It has no free parameters specified Mx starting optimization; number of parameters = 11 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 6 Looks like I got stuck here. Check the following: 1. The model is correctly specified 2. Starting values are good 3. You are not already at the solution The error can arise if the Hessian is ill-conditioned You can try resetting it to an identity matrix and fit from the solution by putting TH=-n on the OU line where n is the number of refits that you want to do If all else fails try putting NAG=30 on the OU line and examine the file NAGDUMP.OUT and the NAG manual MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5798509 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3362271 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3362271 2 0.3362271 1.0000000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5798509 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3362271 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3362271 2 0.3362271 1.0000000 GROUP NUMBER: 3 G3 CALCULATION GROUP Control Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5696607 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3245133 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3245133 2 0.3245133 1.0000000 GROUP NUMBER: 4 G4 CALCULATION GROUP Clinical Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5696607 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3245133 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3245133 2 0.3245133 1.0000000 GROUP NUMBER: 5 Epiphenomena model - fit model to Control A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3362271 2 0.3362271 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3245133 2 0.3245133 1.0000000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z))] 1 1 0.9726345 MATRIX D This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W))] 1 1 0.0200498 MATRIX E This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z))] 1 1 0.0064230 MATRIX F This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W))] 1 1 0.0008928 MATRIX G This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z))] 1 1 0.9494143 MATRIX H This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W))] 1 1 0.0279298 MATRIX I This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z))] 1 1 0.0198429 MATRIX J This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W))] 1 1 0.0028130 MATRIX K This is a computed FULL matrix of order 16 by 1 [=C.G_C.I.(Z-R)_C.H.(Z-R)_E.G.(Z-P)_D.G.(Z-P)_E.I+E.G.P+C.I.R_D.H+D.G.P+C.H.R_C.J.(Z-R).(Z-R)_E.H.(Z-P).(Z-R)_D.I.(Z-P).(Z-R)_E.J.(Z-R)+E.H.P.(Z-R)+C.J.R.(Z-R)_D.J.(Z-R)+D.I.P.(Z-R)+C.J.R.(Z-R)_F.G.(Z-P).(Z-P)_F.I.(Z-P)+F.G.P.(Z-P)+D.I.R.(Z-P)_F.H.(Z-P)+F.G.P.(Z-P)+E.H.R.(Z-P)_F.J+F.H.P+F.I.P+F.G.P.P+D.J.R+D.I.P.R+E.J.R+E.H.P.R+C.J.R.R] 1 1 0.9234331 2 0.0161026 3 0.0226652 4 0.0032498 5 0.0101446 6 0.0061729 7 0.0139512 8 0.0019046 9 0.0000798 10 0.0001769 11 0.0004632 12 0.0005803 13 0.0002407 14 0.0002556 15 0.0002401 16 0.0003394 MATRIX L This is a FULL matrix of order 4 by 1 1 1 1.0000000 2 1.0000000 3 1.0000000 4 1.0000000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.0000000 1.0000000 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 12.8900000 13.9100000 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670693 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1232404 -0.1232404 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.1656628 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0799201 -0.0799201 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 4.0299438 4.0299438 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 3.0319187 3.0319187 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1591549 MATRIX W This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 -0.2810097 -0.2810097 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 -0.0502427 -0.0502427 MATRIX Z This is a FULL matrix of order 1 by 1 1 1 1.0000000 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1591549 Function value of this group: 1226.7470466 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 6 Epiphenomena model - fit model to Clinical A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3362271 2 0.3362271 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3245133 2 0.3245133 1.0000000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z))] 1 1 0.9379064 MATRIX D This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W))] 1 1 0.0327459 MATRIX E This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z))] 1 1 0.0252218 MATRIX F This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W))] 1 1 0.0041259 MATRIX G This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z))] 1 1 0.9117238 MATRIX H This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W))] 1 1 0.0504494 MATRIX I This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z))] 1 1 0.0310132 MATRIX J This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W))] 1 1 0.0068136 MATRIX K This is a computed FULL matrix of order 12 by 1 [=(Z).(C.I.(Z-R))_(N).(E.G.(Z-P))_(Z+N).(E.I+E.G.P+C.I.R)_(Z).(C.J.(Z-R).(Z-R))_(N).(E.H.(Z-P).(Z-R))_(Z).(D.I.(Z-P).(Z-R))_(Z+N).(E.J.(Z-R)+E.H.P.(Z-R)+C.J.R.(Z-R))_(Z).(D.J.(Z-R)+D.I.P.(Z-R)+C.J.R.(Z-R))_(N).(F.G.(Z-P).(Z-P))_(Z+N).(F.I.(Z-P)+F.G.P.(Z-P)+D.I.R.(Z-P))_(N).(F.H.(Z-P)+F.G.P.(Z-P)+E.H.R.(Z-P))_(Z+N).(F.J+F.H.P+F.I.P+F.G.P.P+D.J.R+D.I.P.R+E.J.R+E.H.P.R+C.J.R.R)] 1 1 0.0242688 2 0.0024922 3 0.0196646 4 0.0044486 5 0.0001151 6 0.0004516 7 0.0018322 8 0.0014652 9 0.0002173 10 0.0013167 11 0.0002358 12 0.0017131 MATRIX L This is a FULL matrix of order 4 by 1 1 1 1.0000000 2 1.0000000 3 1.0000000 4 1.0000000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 1.0000000 0.0000000 MATRIX N This is a FULL matrix of order 1 by 1 1 1 0.2033665 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 15.0800000 18.1900000 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670693 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1232404 -0.1232404 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.1656628 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0799201 -0.0799201 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 4.0299438 4.0299438 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 3.0319187 3.0319187 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1591549 MATRIX W This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 -0.2810097 -0.2810097 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 -0.0502427 -0.0502427 MATRIX Z This is a FULL matrix of order 1 by 1 1 1 1.0000000 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1591549 Function value of this group: 1156.5559482 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE RED - Hessian/precision problem Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2383.303 Degrees of freedom >>>>>>>>>>>>>>>> 664 Akaike's Information Criterion >>>> 1055.303 Bayesian Information Criterion >>>> -971.233 Sample size Adjusted BIC >>>> 82.899 Deviance Information Criterion >>>> -361.058 ** Convergence has not been achieved in MAXIT steps: MAXIT = 500 ** ABNORMAL EXIT from NAG Library routine G01GCF: IFAIL = 3 ** NAG soft failure - control returned ** Convergence has not been achieved in MAXIT steps: MAXIT = 500 ** ABNORMAL EXIT from NAG Library routine G01GCF: IFAIL = 3 ** NAG soft failure - control returned I had some trouble calculating them, but perhaps 90% confidence intervals on Chi-squared could be: from 0.66400E+03 to 0.14620E+04 90% confidence intervals on AIC: from ******* to 134.039 90% confidence intervals on RMSEA: from 0.000 to 0.000 Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 11 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5798371 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3362111 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5798371 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3362111 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 GROUP NUMBER: 3 G3 CALCULATION GROUP Control Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5696663 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3245197 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 GROUP NUMBER: 4 G4 CALCULATION GROUP Clinical Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5696663 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3245197 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 GROUP NUMBER: 5 Epiphenomena model - fit model to Control A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z))] 1 1 0.9728097 MATRIX D This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W))] 1 1 0.0199244 MATRIX E This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z))] 1 1 0.0063822 MATRIX F This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W))] 1 1 0.0008836 MATRIX G This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z))] 1 1 0.9499207 MATRIX H This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W))] 1 1 0.0276537 MATRIX I This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z))] 1 1 0.0196576 MATRIX J This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W))] 1 1 0.0027680 MATRIX K This is a computed FULL matrix of order 16 by 1 [=C.G_C.I.(Z-R)_C.H.(Z-R)_E.G.(Z-P)_D.G.(Z-P)_E.I+E.G.P+C.I.R_D.H+D.G.P+C.H.R_C.J.(Z-R).(Z-R)_E.H.(Z-P).(Z-R)_D.I.(Z-P).(Z-R)_E.J.(Z-R)+E.H.P.(Z-R)+C.J.R.(Z-R)_D.J.(Z-R)+D.I.P.(Z-R)+C.J.R.(Z-R)_F.G.(Z-P).(Z-P)_F.I.(Z-P)+F.G.P.(Z-P)+D.I.R.(Z-P)_F.H.(Z-P)+F.G.P.(Z-P)+E.H.R.(Z-P)_F.J+F.H.P+F.I.P+F.G.P.P+D.J.R+D.I.P.R+E.J.R+E.H.P.R+C.J.R.R] 1 1 0.9240921 2 0.0159550 3 0.0224451 4 0.0032309 5 0.0100865 6 0.0061252 7 0.0138478 8 0.0018745 9 0.0000785 10 0.0001741 11 0.0004557 12 0.0005708 13 0.0002384 14 0.0002528 15 0.0002375 16 0.0003350 MATRIX L This is a FULL matrix of order 4 by 1 1 1 1.0000000 2 1.0000000 3 1.0000000 4 1.0000000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.0000000 1.0000000 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 12.8900000 13.9100000 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670750 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1230482 -0.1230482 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.1656655 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0795865 -0.0795865 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 4.0299333 4.0299333 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 3.0319172 3.0319172 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1591549 MATRIX W This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 -0.2809952 -0.2809952 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 -0.0502377 -0.0502377 MATRIX Z This is a FULL matrix of order 1 by 1 1 1 1.0000000 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1591549 Function value of this group: 1226.6425917 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 6 Epiphenomena model - fit model to Clinical A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z))] 1 1 0.9383347 MATRIX D This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W))] 1 1 0.0325111 MATRIX E This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z))] 1 1 0.0250738 MATRIX F This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W))] 1 1 0.0040804 MATRIX G This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z))] 1 1 0.9127078 MATRIX H This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W))] 1 1 0.0498781 MATRIX I This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z))] 1 1 0.0307226 MATRIX J This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W))] 1 1 0.0066915 MATRIX K This is a computed FULL matrix of order 12 by 1 [=(Z).(C.I.(Z-R))_(N).(E.G.(Z-P))_(Z+N).(E.I+E.G.P+C.I.R)_(Z).(C.J.(Z-R).(Z-R))_(N).(E.H.(Z-P).(Z-R))_(Z).(D.I.(Z-P).(Z-R))_(Z+N).(E.J.(Z-R)+E.H.P.(Z-R)+C.J.R.(Z-R))_(Z).(D.J.(Z-R)+D.I.P.(Z-R)+C.J.R.(Z-R))_(N).(F.G.(Z-P).(Z-P))_(Z+N).(F.I.(Z-P)+F.G.P.(Z-P)+D.I.R.(Z-P))_(N).(F.H.(Z-P)+F.G.P.(Z-P)+E.H.R.(Z-P))_(Z+N).(F.J+F.H.P+F.I.P+F.G.P.P+D.J.R+D.I.P.R+E.J.R+E.H.P.R+C.J.R.R)] 1 1 0.0240522 2 0.0024803 3 0.0195369 4 0.0043708 5 0.0001131 6 0.0004441 7 0.0017993 8 0.0014386 9 0.0002151 10 0.0013021 11 0.0002330 12 0.0016891 MATRIX L This is a FULL matrix of order 4 by 1 1 1 1.0000000 2 1.0000000 3 1.0000000 4 1.0000000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 1.0000000 0.0000000 MATRIX N This is a FULL matrix of order 1 by 1 1 1 0.2033701 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 15.0800000 18.1900000 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670750 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1230482 -0.1230482 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.1656655 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0795865 -0.0795865 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 4.0299333 4.0299333 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 3.0319172 3.0319172 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1591549 MATRIX W This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 -0.2809952 -0.2809952 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 -0.0502377 -0.0502377 MATRIX Z This is a FULL matrix of order 1 by 1 1 1 1.0000000 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1591549 Function value of this group: 1156.6495305 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2383.292 Degrees of freedom >>>>>>>>>>>>>>>> 664 Akaike's Information Criterion >>>> 1055.292 Bayesian Information Criterion >>>> -971.239 Sample size Adjusted BIC >>>> 82.893 Deviance Information Criterion >>>> -361.063 ** Convergence has not been achieved in MAXIT steps: MAXIT = 500 ** ABNORMAL EXIT from NAG Library routine G01GCF: IFAIL = 3 ** NAG soft failure - control returned ** Convergence has not been achieved in MAXIT steps: MAXIT = 500 ** ABNORMAL EXIT from NAG Library routine G01GCF: IFAIL = 3 ** NAG soft failure - control returned 90% confidence intervals on Chi-squared: from 0.66400E+03 to 0.14620E+04 90% confidence intervals on AIC: from ******* to 134.039 90% confidence intervals on RMSEA: from 0.000 to 0.000 Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 11 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5798371 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3362111 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5798371 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3362111 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 GROUP NUMBER: 3 G3 CALCULATION GROUP Control Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5696663 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3245197 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 GROUP NUMBER: 4 G4 CALCULATION GROUP Clinical Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5696663 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3245197 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 GROUP NUMBER: 5 Epiphenomena model - fit model to Control A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z))] 1 1 0.9728097 MATRIX D This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W))] 1 1 0.0199244 MATRIX E This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z))] 1 1 0.0063822 MATRIX F This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W))] 1 1 0.0008836 MATRIX G This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z))] 1 1 0.9499207 MATRIX H This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W))] 1 1 0.0276537 MATRIX I This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z))] 1 1 0.0196576 MATRIX J This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W))] 1 1 0.0027680 MATRIX K This is a computed FULL matrix of order 16 by 1 [=C.G_C.I.(Z-R)_C.H.(Z-R)_E.G.(Z-P)_D.G.(Z-P)_E.I+E.G.P+C.I.R_D.H+D.G.P+C.H.R_C.J.(Z-R).(Z-R)_E.H.(Z-P).(Z-R)_D.I.(Z-P).(Z-R)_E.J.(Z-R)+E.H.P.(Z-R)+C.J.R.(Z-R)_D.J.(Z-R)+D.I.P.(Z-R)+C.J.R.(Z-R)_F.G.(Z-P).(Z-P)_F.I.(Z-P)+F.G.P.(Z-P)+D.I.R.(Z-P)_F.H.(Z-P)+F.G.P.(Z-P)+E.H.R.(Z-P)_F.J+F.H.P+F.I.P+F.G.P.P+D.J.R+D.I.P.R+E.J.R+E.H.P.R+C.J.R.R] 1 1 0.9240921 2 0.0159550 3 0.0224451 4 0.0032309 5 0.0100865 6 0.0061252 7 0.0138478 8 0.0018745 9 0.0000785 10 0.0001741 11 0.0004557 12 0.0005708 13 0.0002384 14 0.0002528 15 0.0002375 16 0.0003350 MATRIX L This is a FULL matrix of order 4 by 1 1 1 1.0000000 2 1.0000000 3 1.0000000 4 1.0000000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.0000000 1.0000000 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 12.8900000 13.9100000 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670750 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1230482 -0.1230482 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.1656655 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0795865 -0.0795865 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 4.0299333 4.0299333 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 3.0319172 3.0319172 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1591549 MATRIX W This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 -0.2809952 -0.2809952 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 -0.0502377 -0.0502377 MATRIX Z This is a FULL matrix of order 1 by 1 1 1 1.0000000 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1591549 Function value of this group: 1226.6425917 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 6 Epiphenomena model - fit model to Clinical A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z))] 1 1 0.9383347 MATRIX D This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W))] 1 1 0.0325111 MATRIX E This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z))] 1 1 0.0250738 MATRIX F This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W))] 1 1 0.0040804 MATRIX G This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z))] 1 1 0.9127078 MATRIX H This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W))] 1 1 0.0498781 MATRIX I This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z))] 1 1 0.0307226 MATRIX J This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W))] 1 1 0.0066915 MATRIX K This is a computed FULL matrix of order 12 by 1 [=(Z).(C.I.(Z-R))_(N).(E.G.(Z-P))_(Z+N).(E.I+E.G.P+C.I.R)_(Z).(C.J.(Z-R).(Z-R))_(N).(E.H.(Z-P).(Z-R))_(Z).(D.I.(Z-P).(Z-R))_(Z+N).(E.J.(Z-R)+E.H.P.(Z-R)+C.J.R.(Z-R))_(Z).(D.J.(Z-R)+D.I.P.(Z-R)+C.J.R.(Z-R))_(N).(F.G.(Z-P).(Z-P))_(Z+N).(F.I.(Z-P)+F.G.P.(Z-P)+D.I.R.(Z-P))_(N).(F.H.(Z-P)+F.G.P.(Z-P)+E.H.R.(Z-P))_(Z+N).(F.J+F.H.P+F.I.P+F.G.P.P+D.J.R+D.I.P.R+E.J.R+E.H.P.R+C.J.R.R)] 1 1 0.0240522 2 0.0024803 3 0.0195369 4 0.0043708 5 0.0001131 6 0.0004441 7 0.0017993 8 0.0014386 9 0.0002151 10 0.0013021 11 0.0002330 12 0.0016891 MATRIX L This is a FULL matrix of order 4 by 1 1 1 1.0000000 2 1.0000000 3 1.0000000 4 1.0000000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 1.0000000 0.0000000 MATRIX N This is a FULL matrix of order 1 by 1 1 1 0.2033701 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 15.0800000 18.1900000 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670750 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1230482 -0.1230482 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.1656655 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0795865 -0.0795865 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 4.0299333 4.0299333 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 3.0319172 3.0319172 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1591549 MATRIX W This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 -0.2809952 -0.2809952 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 -0.0502377 -0.0502377 MATRIX Z This is a FULL matrix of order 1 by 1 1 1 1.0000000 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1591549 Function value of this group: 1156.6495305 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2383.292 Degrees of freedom >>>>>>>>>>>>>>>> 664 Akaike's Information Criterion >>>> 1055.292 Bayesian Information Criterion >>>> -971.239 Sample size Adjusted BIC >>>> 82.893 Deviance Information Criterion >>>> -361.063 ** Convergence has not been achieved in MAXIT steps: MAXIT = 500 ** ABNORMAL EXIT from NAG Library routine G01GCF: IFAIL = 3 ** NAG soft failure - control returned ** Convergence has not been achieved in MAXIT steps: MAXIT = 500 ** ABNORMAL EXIT from NAG Library routine G01GCF: IFAIL = 3 ** NAG soft failure - control returned 90% confidence intervals on Chi-squared: from 0.66400E+03 to 0.14620E+04 90% confidence intervals on AIC: from ******* to 134.039 90% confidence intervals on RMSEA: from 0.000 to 0.000 Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 11 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5798371 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3362111 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5798371 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3362111 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 GROUP NUMBER: 3 G3 CALCULATION GROUP Control Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5696663 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3245197 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 GROUP NUMBER: 4 G4 CALCULATION GROUP Clinical Correlation matrix B factor MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000000 MATRIX C This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5696663 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000000 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX X This is a computed FULL matrix of order 1 by 1 [=A*A'] 1 1 0.0000000 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=C*C'] 1 1 0.3245197 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 GROUP NUMBER: 5 Epiphenomena model - fit model to Control A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z))] 1 1 0.9728097 MATRIX D This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W))] 1 1 0.0199244 MATRIX E This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z))] 1 1 0.0063822 MATRIX F This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W))] 1 1 0.0008836 MATRIX G This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z))] 1 1 0.9499207 MATRIX H This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W))] 1 1 0.0276537 MATRIX I This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z))] 1 1 0.0196576 MATRIX J This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W))] 1 1 0.0027680 MATRIX K This is a computed FULL matrix of order 16 by 1 [=C.G_C.I.(Z-R)_C.H.(Z-R)_E.G.(Z-P)_D.G.(Z-P)_E.I+E.G.P+C.I.R_D.H+D.G.P+C.H.R_C.J.(Z-R).(Z-R)_E.H.(Z-P).(Z-R)_D.I.(Z-P).(Z-R)_E.J.(Z-R)+E.H.P.(Z-R)+C.J.R.(Z-R)_D.J.(Z-R)+D.I.P.(Z-R)+C.J.R.(Z-R)_F.G.(Z-P).(Z-P)_F.I.(Z-P)+F.G.P.(Z-P)+D.I.R.(Z-P)_F.H.(Z-P)+F.G.P.(Z-P)+E.H.R.(Z-P)_F.J+F.H.P+F.I.P+F.G.P.P+D.J.R+D.I.P.R+E.J.R+E.H.P.R+C.J.R.R] 1 1 0.9240921 2 0.0159550 3 0.0224451 4 0.0032309 5 0.0100865 6 0.0061252 7 0.0138478 8 0.0018745 9 0.0000785 10 0.0001741 11 0.0004557 12 0.0005708 13 0.0002384 14 0.0002528 15 0.0002375 16 0.0003350 MATRIX L This is a FULL matrix of order 4 by 1 1 1 1.0000000 2 1.0000000 3 1.0000000 4 1.0000000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.0000000 1.0000000 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 12.8900000 13.9100000 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670750 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1230482 -0.1230482 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.1656655 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0795865 -0.0795865 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 4.0299333 4.0299333 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 3.0319172 3.0319172 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1591549 MATRIX W This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 -0.2809952 -0.2809952 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 -0.0502377 -0.0502377 MATRIX Z This is a FULL matrix of order 1 by 1 1 1 1.0000000 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1591549 Function value of this group: 1226.6425917 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 6 Epiphenomena model - fit model to Clinical A/D data MATRIX A This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3362111 2 0.3362111 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.3245197 2 0.3245197 1.0000000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|Z))] 1 1 0.9383347 MATRIX D This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_Z|W))] 1 1 0.0325111 MATRIX E This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|Z))] 1 1 0.0250738 MATRIX F This is a computed FULL matrix of order 1 by 1 [=\MULN((A_(T+(O.Q)+(M.X))_(T+(O.Q)+(M.X))_W|W))] 1 1 0.0040804 MATRIX G This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|Z))] 1 1 0.9127078 MATRIX H This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_Z|W))] 1 1 0.0498781 MATRIX I This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|Z))] 1 1 0.0307226 MATRIX J This is a computed FULL matrix of order 1 by 1 [=\MULN((B_(U+(O.S)+(M.Y))_(U+(O.S)+(M.Y))_W|W))] 1 1 0.0066915 MATRIX K This is a computed FULL matrix of order 12 by 1 [=(Z).(C.I.(Z-R))_(N).(E.G.(Z-P))_(Z+N).(E.I+E.G.P+C.I.R)_(Z).(C.J.(Z-R).(Z-R))_(N).(E.H.(Z-P).(Z-R))_(Z).(D.I.(Z-P).(Z-R))_(Z+N).(E.J.(Z-R)+E.H.P.(Z-R)+C.J.R.(Z-R))_(Z).(D.J.(Z-R)+D.I.P.(Z-R)+C.J.R.(Z-R))_(N).(F.G.(Z-P).(Z-P))_(Z+N).(F.I.(Z-P)+F.G.P.(Z-P)+D.I.R.(Z-P))_(N).(F.H.(Z-P)+F.G.P.(Z-P)+E.H.R.(Z-P))_(Z+N).(F.J+F.H.P+F.I.P+F.G.P.P+D.J.R+D.I.P.R+E.J.R+E.H.P.R+C.J.R.R)] 1 1 0.0240522 2 0.0024803 3 0.0195369 4 0.0043708 5 0.0001131 6 0.0004441 7 0.0017993 8 0.0014386 9 0.0002151 10 0.0013021 11 0.0002330 12 0.0016891 MATRIX L This is a FULL matrix of order 4 by 1 1 1 1.0000000 2 1.0000000 3 1.0000000 4 1.0000000 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 1.0000000 0.0000000 MATRIX N This is a FULL matrix of order 1 by 1 1 1 0.2033701 MATRIX O This is a FULL matrix of order 1 by 2 1 2 1 15.0800000 18.1900000 MATRIX P This is a FULL matrix of order 1 by 1 1 1 0.4670750 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1230482 -0.1230482 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.1656655 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0795865 -0.0795865 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 4.0299333 4.0299333 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 3.0319172 3.0319172 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1591549 MATRIX W This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX X This is a FULL matrix of order 1 by 2 1 2 1 -0.2809952 -0.2809952 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 -0.0502377 -0.0502377 MATRIX Z This is a FULL matrix of order 1 by 1 1 1 1.0000000 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1591549 Function value of this group: 1156.6495305 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2383.292 Degrees of freedom >>>>>>>>>>>>>>>> 664 Akaike's Information Criterion >>>> 1055.292 Bayesian Information Criterion >>>> -971.239 Sample size Adjusted BIC >>>> 82.893 Deviance Information Criterion >>>> -361.063 ** Convergence has not been achieved in MAXIT steps: MAXIT = 500 ** ABNORMAL EXIT from NAG Library routine G01GCF: IFAIL = 3 ** NAG soft failure - control returned ** Convergence has not been achieved in MAXIT steps: MAXIT = 500 ** ABNORMAL EXIT from NAG Library routine G01GCF: IFAIL = 3 ** NAG soft failure - control returned 90% confidence intervals on Chi-squared: from 0.66400E+03 to 0.14620E+04 90% confidence intervals on AIC: from ******* to 134.039 90% confidence intervals on RMSEA: from 0.000 to 0.000 This problem used 7.2% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 2.20 Execution 0: 0: 1: 5.09 TOTAL 0: 0: 1: 7.29 Total number of warnings issued: 2 ______________________________________________________________________________