** Mx startup successful ** **MX-Linux version 1.52b** The following MX script lines were read for group 1 G1 CALCULATION GROUP CONTROL GROUP CORRELATION MATRIX A FACTOR DA CALC NG=4 MATRICES A LO 2 2 FIXED C LO 2 2 FREE E ST 2 2 FREE H FU 1 1 BEGIN ALGEBRA; X = A*A' ; Y = C*C' ; Z = E ; END ALGEBRA; COMPUTE \STND( X+Y+Z| H@X+Y _ H@X+Y| X+Y+Z) ; MATRIX H .5 MATRIX A 0 0 0 MATRIX C .8 .8 .8 MATRIX E .5 OPTION RS SE STANDARD ERRORS withdrawn in this version Please use interval command instead END The following MX script lines were read for group 2 G2 CALCULATION GROUP CLINICAL GROUP CORRELATION MATRIX A FACTOR DA CALC MATRICES A LO 2 2 = A1 C LO 2 2 = C1 E ST 2 2 = E1 H FU 1 1 ! .5 BEGIN ALGEBRA; X = A*A' ; Y = C*C' ; Z = E ; END ALGEBRA; COMPUTE \STND( X+Y+Z| H@X+Y _ H@X+Y| X+Y+Z) ; MATRIX H .5 OPTION RS SE STANDARD ERRORS withdrawn in this version Please use interval command instead END The following MX script lines were read for group 3 CORRELATED LIABILITY MODEL: CONTROL GROUP 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 I IDEN 1 1 R FULL 4 4 =%E1 ! CORRELATION MATRIX A1B1A2B2 A FULL 1 4 ! MATRIX FOR AGE DEFINITION VARIABLE S FULL 1 4 ! MATRIX FOR SEX DEFINITION VARIABLE L FULL 1 4 ! ESTIMATED THRESHOLD FOR A O FULL 1 4 ! AGE DIFFERENCE IN THRESHOLD D FULL 1 4 ! SEX DIFFERENCE IN THRESHOLD M FULL 1 1 ! MEAN OF DUMMY VARIABLE V FULL 1 1 ! VARIANCE OF DUMMY VARIABLE J FULL 4 1 ! FOR EXTRACTING RIGHT PART OF K Z ZERO 1 1 BEGIN ALGEBRA ; T = L+(A.O)+(S.D); K = \MULN(R_T_T_(I|I|I|I)) _ ! 1 \MULN(R_T_T_(I|Z|I|I)) _ ! 2 \MULN(R_T_T_(I|I|I|Z)) _ ! 3 \MULN(R_T_T_(Z|I|I|I)) _ ! 4 \MULN(R_T_T_(I|I|Z|I)) _ ! 5 \MULN(R_T_T_(Z|Z|I|I)) _ ! 6 \MULN(R_T_T_(I|I|Z|Z)) _ ! 7 \MULN(R_T_T_(I|Z|I|Z)) _ ! 8 \MULN(R_T_T_(Z|I|I|Z)) _ ! 9 \MULN(R_T_T_(I|Z|Z|I)) _ ! 10 \MULN(R_T_T_(Z|Z|I|Z)) _ ! 11 \MULN(R_T_T_(I|Z|Z|Z)) _ ! 12 \MULN(R_T_T_(Z|I|Z|I)) _ ! 13 \MULN(R_T_T_(Z|Z|Z|I)) _ ! 14 \MULN(R_T_T_(Z|I|Z|Z)) _ ! 15 \MULN(R_T_T_(Z|Z|Z|Z)) ; ! 16 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 END ALGEBRA ; MATRIX J 1 1 1 1 MATRIX M 0 MATRIX V .159154943 MATRIX L 2.0 1.5 2.0 1.5 MATRIX O -.05 -.05 -.05 -.05 MATRIX D -.05 -.05 -.05 -.05 SPECIFY L 20 21 20 21 ! THRESHOLD FOR A SIB 1; THRESHOLD FOR B SIB 1; THRESHOLD FOR A SIB 2; THRESHOLD FOR B SIB 2; SPECIFY O 30 31 30 31 SPECIFY D 40 41 40 41 SP A AGE1 AGE1 AGE2 AGE2; ! AGE FOR SIB 1; AGE FOR SIB 1; AGE FOR SIB 2; AGE FOR SIB 2; SP S SEX1 SEX1 SEX2 SEX2; ! AGE FOR SIB 1; AGE FOR SIB 1; AGE FOR SIB 2; AGE FOR SIB 2; SP J -5 0 -5 0; MEANS M; COVARIANCE V; WEIGHT (\PART(K,J)) / OPTION RS END The following MX script lines were read for group 4 CORRELATED LIABILITY MODEL: INEFFICIENT APPROACH; CLINICAL GROUP 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 I IDEN 1 1 R FULL 4 4 =%E2 ! CORRELATION MATRIX A1B1A2B2 A FULL 1 4 ! MATRIX FOR AGE DEFINITION VARIABLE S FULL 1 4 ! MATRIX FOR SEX DEFINITION VARIABLE L FULL 1 4 =L3! ESTIMATED THRESHOLD FOR A O FULL 1 4 =O3! AGE DIFFERENCE IN THRESHOLD D FULL 1 4 =D3! SEX DIFFERENCE IN THRESHOLD M FULL 1 1 ! MEAN OF DUMMY VARIABLE V FULL 1 1 ! VARIANCE OF DUMMY VARIABLE J FULL 4 1 ! FOR EXTRACTING RIGHT PART OF K Z ZERO 1 1 B FULL 1 1 C FULL 1 1 FREE BEGIN ALGEBRA ; G = B + C; T = L+(A.O)+(S.D); K = B.(\MULN(R_T_T_(I|Z|I|I))) _ ! 2 C.(\MULN(R_T_T_(Z|I|I|I))) _ ! 4 G.(\MULN(R_T_T_(Z|Z|I|I))) _ ! 6 B.(\MULN(R_T_T_(I|Z|I|Z))) _ ! 8 C.(\MULN(R_T_T_(Z|I|I|Z))) _ ! 9 B.(\MULN(R_T_T_(I|Z|Z|I))) _ ! 10 G.(\MULN(R_T_T_(Z|Z|I|Z))) _ ! 11 B.(\MULN(R_T_T_(I|Z|Z|Z))) _ ! 12 C.(\MULN(R_T_T_(Z|I|Z|I))) _ ! 13 G.(\MULN(R_T_T_(Z|Z|Z|I))) _ ! 14 C.(\MULN(R_T_T_(Z|I|Z|Z))) _ ! 15 G.(\MULN(R_T_T_(Z|Z|Z|Z))) ; ! 16 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Problem with correlation between variables 2 and 1 It is supplied as 1.140000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 2 It is supplied as 2.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.280000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 3 It is supplied as 9.000000000000000 which has absolute value > .99 Y = \SUM(K); U = Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y; N = (K%U); END ALGEBRA ; MATRIX J 1 1 1 1 MATRIX M 0 MATRIX V .159154943 MATRIX B 1 MATRIX C .5 SP A AGE1 AGE1 AGE2 AGE2; ! AGE FOR SIB 1; AGE FOR SIB 1; AGE FOR SIB 2; AGE FOR SIB 2; SP S SEX1 SEX1 SEX2 SEX2; ! AGE FOR SIB 1; AGE FOR SIB 1; AGE FOR SIB 2; AGE FOR SIB 2; SP J -5 0 -5 0; MEANS M; COVARIANCE V; WEIGHT (\PART(N,J)) / OPTION FUNC=.000001 OPTION RS OPTION TH=-20 END Summary of VL file data for group 3 TYPE SEX2 SEX1 AGE2 AGE1 Code -5.0000E+00 -4.0000E+00 -3.0000E+00 -2.0000E+00 -1.0000E+00 Number 4.2200E+02 4.2200E+02 4.2200E+02 4.2200E+02 4.2200E+02 Mean 1.6303E+00 5.4976E-01 9.1943E-01 1.6690E+01 1.6478E+01 Variance 3.8017E+00 2.4752E-01 7.4077E-02 1.1376E+01 2.2599E+00 Minimum 1.0000E+00 0.0000E+00 0.0000E+00 1.0960E+01 1.1960E+01 Maximum 1.6000E+01 1.0000E+00 1.0000E+00 2.4950E+01 2.1050E+01 DUMMY Code 1.0000E+00 Number 4.2200E+02 Mean 1.0000E+00 Variance 2.6940E-16 Minimum 1.0000E+00 Maximum 1.0000E+00 Summary of VL file data for group 4 TYPE SEX2 SEX1 AGE2 AGE1 DUMMY Code -5.0000 -4.0000 -3.0000 -2.0000 -1.0000 1.0000 Number 253.0000 253.0000 253.0000 253.0000 253.0000 253.0000 Mean 3.1779 0.5375 0.9289 17.0966 16.4344 1.0000 Variance 7.7826 0.2486 0.0661 10.7495 1.4494 0.0000 Minimum 1.0000 0.0000 0.0000 11.1500 13.0900 1.0000 Maximum 12.0000 1.0000 1.0000 24.5500 19.7100 1.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 4 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Y This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Z This is a computed FULL matrix of order 2 by 2 It has no free parameters specified GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 4 0 MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Y This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Z This is a computed FULL matrix of order 2 by 2 It has no free parameters specified GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 -1 -2 -2 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 40 41 40 41 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 -5 2 0 3 -5 4 0 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 1 by 4 1 2 3 4 1 20 21 20 21 MATRIX M This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 30 31 30 31 MATRIX R This is a constrained FULL matrix of order 4 by 4 It has no free parameters specified MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 -3 -3 -4 -4 MATRIX T This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX V This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a NULL matrix of order 1 by 1 GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 -1 -1 -2 -2 MATRIX B This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX C This is a FULL matrix of order 1 by 1 1 1 42 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 40 41 40 41 MATRIX G This is a computed 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 J This is a FULL matrix of order 4 by 1 1 1 -5 2 0 3 -5 4 0 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 1 by 4 1 2 3 4 1 20 21 20 21 MATRIX M This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX N This is a computed FULL matrix of order 12 by 1 It has no free parameters specified MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 30 31 30 31 MATRIX R This is a constrained FULL matrix of order 4 by 4 It has no free parameters specified MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 -3 -3 -4 -4 MATRIX T This is a computed FULL matrix of order 1 by 4 It has no free parameters specified MATRIX U This is a computed FULL matrix of order 12 by 1 It has no free parameters specified MATRIX V This is a 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 MATRIX Z This is a NULL matrix of order 1 by 1 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 group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6561 2 0.5142 0.4217 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7821 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.4305 0.3373 2 0.3373 0.4422 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7821 2 0.7821 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7794 0.3009 0.2349 2 0.7794 1.0000 0.2349 0.3066 3 0.3009 0.2349 1.0000 0.7794 4 0.2349 0.3066 0.7794 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6561 2 0.5142 0.4217 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7821 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.4305 0.3373 2 0.3373 0.4422 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7821 2 0.7821 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7794 0.3009 0.2349 2 0.7794 1.0000 0.2349 0.3066 3 0.3009 0.2349 1.0000 0.7794 4 0.2349 0.3066 0.7794 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2477 -0.1046 -0.2477 -0.1046 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1359E-01 2 1.8907E-02 3 2.5737E-02 4 4.3314E-03 5 1.2138E-02 6 5.5369E-03 7 1.3486E-02 8 2.1060E-03 9 1.7983E-04 10 4.0979E-04 11 5.5294E-04 12 1.0809E-03 13 3.0940E-04 14 3.1233E-04 15 2.7533E-04 16 5.6686E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5854 2.8710 3.5854 2.8710 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1032 -0.0762 -0.1032 -0.0762 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7794 0.3009 0.2349 2 0.7794 1.0000 0.2349 0.3066 3 0.3009 0.2349 1.0000 0.7794 4 0.2349 0.3066 0.7794 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2551 1.8886 1.9021 1.7062 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1227.1680 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.1547 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2477 -0.1046 -0.2477 -0.1046 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.1547 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.6255E-02 2 2.1866E-03 3 1.7767E-02 4 3.7943E-03 5 1.2379E-04 6 6.0664E-04 7 2.3096E-03 8 1.9759E-03 9 1.6105E-04 10 1.0410E-03 11 1.7790E-04 12 2.3274E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5854 2.8710 3.5854 2.8710 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4471 2 0.0372 3 0.3025 4 0.0646 5 0.0021 6 0.0103 7 0.0393 8 0.0336 9 0.0027 10 0.0177 11 0.0030 12 0.0396 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1032 -0.0762 -0.1032 -0.0762 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7794 0.3009 0.2349 2 0.7794 1.0000 0.2349 0.3066 3 0.3009 0.2349 1.0000 0.7794 4 0.2349 0.3066 0.7794 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7813 1.6170 1.7081 1.4846 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0587 2 0.0587 3 0.0587 4 0.0587 5 0.0587 6 0.0587 7 0.0587 8 0.0587 9 0.0587 10 0.0587 11 0.0587 12 0.0587 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0587 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1156.4578 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.626 Degrees of freedom >>>>>>>>>>>>>>>> 664 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 group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 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 >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 Now I'm trying to improve on the current solution for you... 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 1 We probably have a minimum here, but you might consider trying different starting values. You can randomize these with TH=n on the OU line, where n is the number of times you wish to do this. I STRONGLY recommend BOundaries to be set if you use TH MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 2 G2 CALCULATION GROUP Clinical group Correlation matrix A factor MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.0000 2 0.0000 0.0000 MATRIX C This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6234 2 0.4963 0.4142 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.7808 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a computed FULL matrix of order 2 by 2 [=A*A'] 1 2 1 0.0000 0.0000 2 0.0000 0.0000 MATRIX Y This is a computed FULL matrix of order 2 by 2 [=C*C'] 1 2 1 0.3886 0.3094 2 0.3094 0.4179 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.7808 2 0.7808 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 GROUP NUMBER: 3 Correlated liability model: control group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 12.8900 12.8900 13.9100 13.9100 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 16 by 1 [=\MULN(R_T_T_(I|I|I|I))_\MULN(R_T_T_(I|Z|I|I))_\MULN(R_T_T_(I|I|I|Z))_\MULN(R_T_T_(Z|I|I|I))_\MULN(R_T_T_(I|I|Z|I))_\MULN(R_T_T_(Z|Z|I|I))_\MULN(R_T_T_(I|I|Z|Z))_\MULN(R_T_T_(I|Z|I|Z))_\MULN(R_T_T_(Z|I|I|Z))_\MULN(R_T_T_(I|Z|Z|I))_\MULN(R_T_T_(Z|Z|I|Z))_\MULN(R_T_T_(I|Z|Z|Z))_\MULN(R_T_T_(Z|I|Z|I))_\MULN(R_T_T_(Z|Z|Z|I))_\MULN(R_T_T_(Z|I|Z|Z))_\MULN(R_T_T_(Z|Z|Z|Z))] 1 1 9.1187E-01 2 1.8637E-02 3 2.6942E-02 4 4.5409E-03 5 1.2353E-02 6 5.5899E-03 7 1.3953E-02 8 2.0856E-03 9 1.9320E-04 10 3.9328E-04 11 5.5740E-04 12 1.0439E-03 13 2.9400E-04 14 2.8862E-04 15 2.7020E-04 16 5.3220E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 0.0000 0.0000 1.0000 1.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 2.2487 1.8937 1.8937 1.6893 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1226.9447 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 4 Correlated liability model: inefficient approach; clinical group data MATRIX A This is a FULL matrix of order 1 by 4 1 2 3 4 1 15.0800 15.0800 18.1900 18.1900 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.2007 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.2493 -0.1261 -0.2493 -0.1261 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.2007 MATRIX I This is an IDENTITY matrix of order 1 by 1 MATRIX J This is a FULL matrix of order 4 by 1 1 1 1.0000 2 1.0000 3 1.0000 4 1.0000 MATRIX K This is a computed FULL matrix of order 12 by 1 [=B.(\MULN(R_T_T_(I|Z|I|I)))_C.(\MULN(R_T_T_(Z|I|I|I)))_G.(\MULN(R_T_T_(Z|Z|I|I)))_B.(\MULN(R_T_T_(I|Z|I|Z)))_C.(\MULN(R_T_T_(Z|I|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|I)))_G.(\MULN(R_T_T_(Z|Z|I|Z)))_B.(\MULN(R_T_T_(I|Z|Z|Z)))_C.(\MULN(R_T_T_(Z|I|Z|I)))_G.(\MULN(R_T_T_(Z|Z|Z|I)))_C.(\MULN(R_T_T_(Z|I|Z|Z)))_G.(\MULN(R_T_T_(Z|Z|Z|Z)))] 1 1 2.7558E-02 2 2.8995E-03 3 1.9241E-02 4 3.8053E-03 5 1.5794E-04 6 6.5274E-04 7 2.3755E-03 8 2.0102E-03 9 2.0175E-04 10 1.0800E-03 11 2.1787E-04 12 2.3221E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.5843 2.8842 3.5843 2.8842 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a computed FULL matrix of order 12 by 1 [=(K%U)] 1 1 0.4408 2 0.0464 3 0.3077 4 0.0609 5 0.0025 6 0.0104 7 0.0380 8 0.0322 9 0.0032 10 0.0173 11 0.0035 12 0.0371 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1036 -0.0768 -0.1036 -0.0768 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7770 0.2799 0.2205 2 0.7770 1.0000 0.2205 0.2947 3 0.2799 0.2205 1.0000 0.7770 4 0.2205 0.2947 0.7770 1.0000 MATRIX S This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 0.0000 0.0000 MATRIX T This is a computed FULL matrix of order 1 by 4 [=L+(A.O)+(S.D)] 1 2 3 4 1 1.7725 1.5994 1.6996 1.4865 MATRIX U This is a computed FULL matrix of order 12 by 1 [=Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y_Y] 1 1 0.0625 2 0.0625 3 0.0625 4 0.0625 5 0.0625 6 0.0625 7 0.0625 8 0.0625 9 0.0625 10 0.0625 11 0.0625 12 0.0625 MATRIX V This is a FULL matrix of order 1 by 1 1 1 0.1592 MATRIX Y This is a computed FULL matrix of order 1 by 1 [=\SUM(K)] 1 1 0.0625 MATRIX Z This is a NULL matrix of order 1 by 1 Matrix of EXPECTED thresholds DUMMY Threshold 1 0.0000 (OBSERVED MATRIX is nonexistent for raw data) EXPECTED COVARIANCE MATRIX DUMMY DUMMY 0.1592 Function value of this group: 1155.8376 Where the fit function is -2 * Log-likelihood of raw ordinal *** WARNING! *** Minimization may not be successful. See above CODE GREEN - it probably was OK Your model has 11 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2382.782 Degrees of freedom >>>>>>>>>>>>>>>> 664 This problem used 0.8% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0:13:56.97 TOTAL 0: 0:13:57.11 Total number of warnings issued: 40 ______________________________________________________________________________ !save corrace.mxs ! no correlation model !drop 2 4 !end ______________________________________________________________________________