** Mx startup successful ** **MX-Linux version 1.52b** 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 ! 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 C.I.(Z-R)_ ! 2 C.H.(Z-R)_ ! 3 E.G.(Z-P)_ ! 4 D.G.(Z-P)_ ! 5 E.I + E.G.P + C.I.R_ ! 6 D.H + D.G.P + C.H.R_ ! 7 C.J.(Z-R).(Z-R)_ ! 8 E.H.(Z-P).(Z-R)_ ! 9 D.I.(Z-P).(Z-R)_ ! 10 E.J.(Z-R) + E.H.P.(Z-R) + C.J.R.(Z-R)_ ! 11 D.J.(Z-R) + D.I.P.(Z-R) + C.J.R.(Z-R)_ ! 12 F.G.(Z-P).(Z-P)_ ! 13 F.I.(Z-P) + F.G.P.(Z-P) + D.I.R.(Z-P)_ ! 14 F.H.(Z-P) + F.G.P.(Z-P) + E.H.R.(Z-P)_ ! 15 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 END ALGEBRA; MATRIX L 1 1 1 1 MATRIX V .159154943 MATRIX P .5 MATRIX R 0 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 - - N FULL 1 1 FREE 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 (N).(E.G.(Z-P))_ ! 4 (Z+N).(E.I + E.G.P + C.I.R)_ ! 6 (Z).(C.J.(Z-R).(Z-R))_ ! 8 (N).(E.H.(Z-P).(Z-R))_ ! 9 (Z).(D.I.(Z-P).(Z-R))_ ! 10 (Z+N).(E.J.(Z-R) + E.H.P.(Z-R) + C.J.R.(Z-R))_ ! 11 (Z).(D.J.(Z-R) + D.I.P.(Z-R) + C.J.R.(Z-R))_ ! 12 (N).(F.G.(Z-P).(Z-P))_ ! 13 (Z+N).(F.I.(Z-P) + F.G.P.(Z-P) + D.I.R.(Z-P))_ ! 14 (N).(F.H.(Z-P) + F.G.P.(Z-P) + E.H.R.(Z-P))_ ! 15 (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 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 !Start .9 all 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 It has no free parameters specified 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 It has no free parameters specified 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 = 10 *** 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.5194860 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.2698657 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2698657 2 0.2698657 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.5194860 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.2698657 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2698657 2 0.2698657 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.5386772 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.2901732 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2901732 2 0.2901732 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.5386772 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.2901732 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2901732 2 0.2901732 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.2698657 2 0.2698657 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.2901732 2 0.2901732 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.9605342 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.0271558 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.0109985 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.0013114 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.9494823 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.0257854 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.0222317 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.0025006 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.9120102 2 0.0213543 3 0.0247677 4 0.0046487 5 0.0114778 6 0.0060387 7 0.0150064 8 0.0024019 9 0.0001262 10 0.0002687 11 0.0001849 12 0.0004029 13 0.0002467 14 0.0003205 15 0.0003226 16 0.0004216 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.5548462 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1120552 -0.1120552 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0702095 -0.0702095 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 3.6917062 3.6917062 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 2.8695648 2.8695648 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.2292056 -0.2292056 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 0.0136526 0.0136526 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: 1225.9554190 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.2698657 2 0.2698657 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.2901732 2 0.2901732 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.9178933 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.0439681 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.0329842 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.0051544 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.9158369 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.0501217 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.0285218 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.0055196 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.0261799 2 0.0027415 3 0.0213105 4 0.0050664 5 0.0001500 6 0.0005582 7 0.0013235 8 0.0009385 9 0.0001907 10 0.0014824 11 0.0002611 12 0.0020545 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.2038731 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.5548462 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1120552 -0.1120552 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0702095 -0.0702095 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 3.6917062 3.6917062 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 2.8695648 2.8695648 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.2292056 -0.2292056 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 0.0136526 0.0136526 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: 1159.2249244 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 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.180 Degrees of freedom >>>>>>>>>>>>>>>> 665 ** 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 ** 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.66500E+03 to 0.14630E+04 90% confidence intervals on AIC: from ******* to 133.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 = 10 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.5194588 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.2698374 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2698374 2 0.2698374 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.5194588 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.2698374 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2698374 2 0.2698374 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.5386544 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.2901486 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2901486 2 0.2901486 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.5386544 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.2901486 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2901486 2 0.2901486 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.2698374 2 0.2698374 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.2901486 2 0.2901486 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.9604290 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.0272265 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.0110274 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.0013171 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.9492901 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.0258843 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.0223096 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.0025160 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.9117257 2 0.0214268 3 0.0248600 4 0.0046600 5 0.0115055 6 0.0060542 7 0.0150451 8 0.0024164 9 0.0001271 10 0.0002704 11 0.0001861 12 0.0004055 13 0.0002478 14 0.0003219 15 0.0003240 16 0.0004234 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.5548409 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1121386 -0.1121386 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0703330 -0.0703330 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 3.6917009 3.6917009 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 2.8695484 2.8695484 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.2292107 -0.2292107 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 0.0136260 0.0136260 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.0232836 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.2698374 2 0.2698374 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.2901486 2 0.2901486 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.9176573 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.0440989 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.0330651 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.0051787 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.9154767 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.0503377 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.0286274 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.0055582 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.0262702 2 0.0027476 3 0.0213594 4 0.0051005 5 0.0001511 6 0.0005620 7 0.0013330 8 0.0009456 9 0.0001916 10 0.0014892 11 0.0002624 12 0.0020649 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.2039041 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.5548409 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1121386 -0.1121386 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0703330 -0.0703330 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 3.6917009 3.6917009 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 2.8695484 2.8695484 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.2292107 -0.2292107 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 0.0136260 0.0136260 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: 1159.1553892 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.179 Degrees of freedom >>>>>>>>>>>>>>>> 665 ** 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 ** 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.66500E+03 to 0.14630E+04 90% confidence intervals on AIC: from ******* to 133.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 = 10 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.5194588 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.2698374 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2698374 2 0.2698374 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.5194588 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.2698374 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2698374 2 0.2698374 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.5386544 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.2901486 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2901486 2 0.2901486 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.5386544 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.2901486 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2901486 2 0.2901486 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.2698374 2 0.2698374 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.2901486 2 0.2901486 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.9604290 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.0272265 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.0110274 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.0013171 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.9492901 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.0258843 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.0223096 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.0025160 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.9117257 2 0.0214268 3 0.0248600 4 0.0046600 5 0.0115055 6 0.0060542 7 0.0150451 8 0.0024164 9 0.0001271 10 0.0002704 11 0.0001861 12 0.0004055 13 0.0002478 14 0.0003219 15 0.0003240 16 0.0004234 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.5548409 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1121386 -0.1121386 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0703330 -0.0703330 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 3.6917009 3.6917009 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 2.8695484 2.8695484 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.2292107 -0.2292107 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 0.0136260 0.0136260 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.0232836 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.2698374 2 0.2698374 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.2901486 2 0.2901486 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.9176573 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.0440989 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.0330651 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.0051787 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.9154767 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.0503377 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.0286274 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.0055582 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.0262702 2 0.0027476 3 0.0213594 4 0.0051005 5 0.0001511 6 0.0005620 7 0.0013330 8 0.0009456 9 0.0001916 10 0.0014892 11 0.0002624 12 0.0020649 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.2039041 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.5548409 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1121386 -0.1121386 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0703330 -0.0703330 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 3.6917009 3.6917009 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 2.8695484 2.8695484 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.2292107 -0.2292107 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 0.0136260 0.0136260 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: 1159.1553892 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.179 Degrees of freedom >>>>>>>>>>>>>>>> 665 ** 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 ** 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.66500E+03 to 0.14630E+04 90% confidence intervals on AIC: from ******* to 133.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 = 10 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.5194588 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.2698374 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2698374 2 0.2698374 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.5194588 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.2698374 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2698374 2 0.2698374 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.5386544 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.2901486 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2901486 2 0.2901486 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.5386544 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.2901486 EXPECTED MATRIX of this CALCULATION group 1 2 1 1.0000000 0.2901486 2 0.2901486 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.2698374 2 0.2698374 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.2901486 2 0.2901486 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.9604290 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.0272265 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.0110274 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.0013171 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.9492901 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.0258843 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.0223096 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.0025160 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.9117257 2 0.0214268 3 0.0248600 4 0.0046600 5 0.0115055 6 0.0060542 7 0.0150451 8 0.0024164 9 0.0001271 10 0.0002704 11 0.0001861 12 0.0004055 13 0.0002478 14 0.0003219 15 0.0003240 16 0.0004234 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.5548409 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1121386 -0.1121386 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0703330 -0.0703330 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 3.6917009 3.6917009 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 2.8695484 2.8695484 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.2292107 -0.2292107 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 0.0136260 0.0136260 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.0232836 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.2698374 2 0.2698374 1.0000000 MATRIX B This is a constrained FULL matrix of order 2 by 2 1 2 1 1.0000000 0.2901486 2 0.2901486 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.9176573 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.0440989 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.0330651 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.0051787 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.9154767 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.0503377 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.0286274 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.0055582 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.0262702 2 0.0027476 3 0.0213594 4 0.0051005 5 0.0001511 6 0.0005620 7 0.0013330 8 0.0009456 9 0.0001916 10 0.0014892 11 0.0002624 12 0.0020649 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.2039041 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.5548409 MATRIX Q This is a FULL matrix of order 1 by 2 1 2 1 -0.1121386 -0.1121386 MATRIX R This is a FULL matrix of order 1 by 1 1 1 0.0000000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 -0.0703330 -0.0703330 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 3.6917009 3.6917009 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 2.8695484 2.8695484 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.2292107 -0.2292107 MATRIX Y This is a FULL matrix of order 1 by 2 1 2 1 0.0136260 0.0136260 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: 1159.1553892 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.179 Degrees of freedom >>>>>>>>>>>>>>>> 665 ** 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 ** 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.66500E+03 to 0.14630E+04 90% confidence intervals on AIC: from ******* to 133.039 90% confidence intervals on RMSEA: from 0.000 to 0.000 This problem used 0.7% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.12 Execution 0: 0: 1:27.67 TOTAL 0: 0: 1:27.79 Total number of warnings issued: 2 ______________________________________________________________________________ !Save biepiace.mxs !Drop 6 !End !Get biepiace.mxs !Drop 5 !End ______________________________________________________________________________