** 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 P FULL 4 4 I IDEN 4 4 H FU 1 1 BEGIN ALGEBRA; X = A*A' ; Y = C*C' ; Z = E ; END ALGEBRA; COMPUTE \STND ((I-P)~ *( X+Y+Z| H@X+Y _ H@X+Y| X+Y+Z)*((I-P)~)') / SPECIFY P 0 0 0 0 5 0 0 0 0 0 0 0 0 0 5 0 MATRIX H .5 MATRIX A 0 0 0 MATRIX C 1 0 1 MATRIX E 0 SPECIFY C 1 0 2 SPECIFY E 0 BOUND 0 1 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 P FULL 4 4 = P1 I IDEN 4 4 = I1 H FU 1 1 ! .5 BEGIN ALGEBRA; X = A*A' ; Y = C*C' ; Z = E ; END ALGEBRA; COMPUTE \STND ((I-P)~ *( X+Y+Z| H@X+Y _ H@X+Y| X+Y+Z)*((I-P)~)') / 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 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 END ALGEBRA ; MATRIX J 1 1 1 1 MATRIX M 0 MATRIX V .159154943 MATRIX L 2.0 2.0 2.0 2.0 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 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 3 and 1 It is supplied as 1.000000000000000 which has absolute value > .99 Problem with correlation between variables 4 and 2 It is supplied as 1.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 RS OPTION TH=-5 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 0 2 MATRIX E This is a STANDARDISED matrix of order 2 by 2 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0 0 0 0 2 5 0 0 0 3 0 0 0 0 4 0 0 5 0 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 0 2 MATRIX E This is a STANDARDISED matrix of order 2 by 2 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0 0 0 0 2 5 0 0 0 3 0 0 0 0 4 0 0 5 0 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 = 10 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.3012 -0.1362 -0.3012 -0.1362 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.1309E-01 2 1.7890E-02 3 2.6822E-02 4 4.5682E-03 5 1.4580E-02 6 4.3937E-03 7 1.2789E-02 8 1.6052E-03 9 2.6463E-04 10 6.1045E-04 11 4.3471E-04 12 9.3205E-04 13 4.2144E-04 14 3.7570E-04 15 3.5753E-04 16 4.9068E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.2799 1.9312 1.8727 1.7093 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.7849 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.1814 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.3012 -0.1362 -0.3012 -0.1362 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.1814 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.8178E-02 2 2.9829E-03 3 1.7249E-02 4 3.3455E-03 5 2.3381E-04 6 8.6583E-04 7 2.2624E-03 8 1.7169E-03 9 2.3834E-04 10 1.2975E-03 11 2.6249E-04 12 2.2003E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.4632 2 0.0490 3 0.2836 4 0.0550 5 0.0038 6 0.0142 7 0.0372 8 0.0282 9 0.0039 10 0.0213 11 0.0043 12 0.0362 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7510 1.6110 1.7289 1.4859 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.0608 2 0.0608 3 0.0608 4 0.0608 5 0.0608 6 0.0608 7 0.0608 8 0.0608 9 0.0608 10 0.0608 11 0.0608 12 0.0608 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.0608 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: 1158.6849 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.470 Degrees of freedom >>>>>>>>>>>>>>>> 665 Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 10 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.3012 -0.1362 -0.3012 -0.1362 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.1309E-01 2 1.7890E-02 3 2.6822E-02 4 4.5682E-03 5 1.4580E-02 6 4.3937E-03 7 1.2789E-02 8 1.6052E-03 9 2.6463E-04 10 6.1045E-04 11 4.3471E-04 12 9.3205E-04 13 4.2144E-04 14 3.7570E-04 15 3.5753E-04 16 4.9068E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.2799 1.9312 1.8727 1.7093 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.7849 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.1814 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.3012 -0.1362 -0.3012 -0.1362 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.1814 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.8178E-02 2 2.9829E-03 3 1.7249E-02 4 3.3455E-03 5 2.3381E-04 6 8.6583E-04 7 2.2624E-03 8 1.7169E-03 9 2.3834E-04 10 1.2975E-03 11 2.6249E-04 12 2.2003E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.4632 2 0.0490 3 0.2836 4 0.0550 5 0.0038 6 0.0142 7 0.0372 8 0.0282 9 0.0039 10 0.0213 11 0.0043 12 0.0362 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7510 1.6110 1.7289 1.4859 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.0608 2 0.0608 3 0.0608 4 0.0608 5 0.0608 6 0.0608 7 0.0608 8 0.0608 9 0.0608 10 0.0608 11 0.0608 12 0.0608 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.0608 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: 1158.6849 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.470 Degrees of freedom >>>>>>>>>>>>>>>> 665 Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 10 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.3012 -0.1362 -0.3012 -0.1362 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.1309E-01 2 1.7890E-02 3 2.6822E-02 4 4.5682E-03 5 1.4580E-02 6 4.3937E-03 7 1.2789E-02 8 1.6052E-03 9 2.6463E-04 10 6.1045E-04 11 4.3471E-04 12 9.3205E-04 13 4.2144E-04 14 3.7570E-04 15 3.5753E-04 16 4.9068E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.2799 1.9312 1.8727 1.7093 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.7849 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.1814 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.3012 -0.1362 -0.3012 -0.1362 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.1814 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.8178E-02 2 2.9829E-03 3 1.7249E-02 4 3.3455E-03 5 2.3381E-04 6 8.6583E-04 7 2.2624E-03 8 1.7169E-03 9 2.3834E-04 10 1.2975E-03 11 2.6249E-04 12 2.2003E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.4632 2 0.0490 3 0.2836 4 0.0550 5 0.0038 6 0.0142 7 0.0372 8 0.0282 9 0.0039 10 0.0213 11 0.0043 12 0.0362 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7510 1.6110 1.7289 1.4859 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.0608 2 0.0608 3 0.0608 4 0.0608 5 0.0608 6 0.0608 7 0.0608 8 0.0608 9 0.0608 10 0.0608 11 0.0608 12 0.0608 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.0608 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: 1158.6849 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.470 Degrees of freedom >>>>>>>>>>>>>>>> 665 Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 10 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 CALCULATION GROUP Control 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.3012 -0.1362 -0.3012 -0.1362 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.1309E-01 2 1.7890E-02 3 2.6822E-02 4 4.5682E-03 5 1.4580E-02 6 4.3937E-03 7 1.2789E-02 8 1.6052E-03 9 2.6463E-04 10 6.1045E-04 11 4.3471E-04 12 9.3205E-04 13 4.2144E-04 14 3.7570E-04 15 3.5753E-04 16 4.9068E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.2799 1.9312 1.8727 1.7093 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.7849 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.1814 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.3012 -0.1362 -0.3012 -0.1362 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.1814 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.8178E-02 2 2.9829E-03 3 1.7249E-02 4 3.3455E-03 5 2.3381E-04 6 8.6583E-04 7 2.2624E-03 8 1.7169E-03 9 2.3834E-04 10 1.2975E-03 11 2.6249E-04 12 2.2003E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.4632 2 0.0490 3 0.2836 4 0.0550 5 0.0038 6 0.0142 7 0.0372 8 0.0282 9 0.0039 10 0.0213 11 0.0043 12 0.0362 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7510 1.6110 1.7289 1.4859 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.0608 2 0.0608 3 0.0608 4 0.0608 5 0.0608 6 0.0608 7 0.0608 8 0.0608 9 0.0608 10 0.0608 11 0.0608 12 0.0608 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.0608 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: 1158.6849 Where the fit function is -2 * Log-likelihood of raw ordinal Your model has 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.470 Degrees of freedom >>>>>>>>>>>>>>>> 665 Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 10 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.3012 -0.1362 -0.3012 -0.1362 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.1309E-01 2 1.7890E-02 3 2.6822E-02 4 4.5682E-03 5 1.4580E-02 6 4.3937E-03 7 1.2789E-02 8 1.6052E-03 9 2.6463E-04 10 6.1045E-04 11 4.3471E-04 12 9.3205E-04 13 4.2144E-04 14 3.7570E-04 15 3.5753E-04 16 4.9068E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.2799 1.9312 1.8727 1.7093 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.7849 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.1814 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.3012 -0.1362 -0.3012 -0.1362 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.1814 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.8178E-02 2 2.9829E-03 3 1.7249E-02 4 3.3455E-03 5 2.3381E-04 6 8.6583E-04 7 2.2624E-03 8 1.7169E-03 9 2.3834E-04 10 1.2975E-03 11 2.6249E-04 12 2.2003E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.4632 2 0.0490 3 0.2836 4 0.0550 5 0.0038 6 0.0142 7 0.0372 8 0.0282 9 0.0039 10 0.0213 11 0.0043 12 0.0362 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7510 1.6110 1.7289 1.4859 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.0608 2 0.0608 3 0.0608 4 0.0608 5 0.0608 6 0.0608 7 0.0608 8 0.0608 9 0.0608 10 0.0608 11 0.0608 12 0.0608 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.0608 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: 1158.6849 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 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.470 Degrees of freedom >>>>>>>>>>>>>>>> 665 Now I'm trying to improve on the current solution for you... Mx starting optimization; number of parameters = 10 *** WARNING! *** I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7050 2 0.0000 0.4797 MATRIX E This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.0000 1.0000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX I This is an IDENTITY matrix of order 4 by 4 MATRIX P This is a FULL matrix of order 4 by 4 1 2 3 4 1 0.0000 0.0000 0.0000 0.0000 2 1.0000 0.0000 0.0000 0.0000 3 0.0000 0.0000 0.0000 0.0000 4 0.0000 0.0000 1.0000 0.0000 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.4971 0.0000 2 0.0000 0.2301 MATRIX Z This is a computed FULL matrix of order 2 by 2 [=E] 1 2 1 1.0000 0.0000 2 0.0000 1.0000 EXPECTED MATRIX of this CALCULATION group 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.3012 -0.1362 -0.3012 -0.1362 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.1309E-01 2 1.7890E-02 3 2.6822E-02 4 4.5682E-03 5 1.4580E-02 6 4.3937E-03 7 1.2789E-02 8 1.6052E-03 9 2.6463E-04 10 6.1045E-04 11 4.3471E-04 12 9.3205E-04 13 4.2144E-04 14 3.7570E-04 15 3.5753E-04 16 4.9068E-04 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.2799 1.9312 1.8727 1.7093 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.7849 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.1814 MATRIX D This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.3012 -0.1362 -0.3012 -0.1362 MATRIX G This is a computed FULL matrix of order 1 by 1 [=B+C] 1 1 1.1814 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.8178E-02 2 2.9829E-03 3 1.7249E-02 4 3.3455E-03 5 2.3381E-04 6 8.6583E-04 7 2.2624E-03 8 1.7169E-03 9 2.3834E-04 10 1.2975E-03 11 2.6249E-04 12 2.2003E-03 MATRIX L This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.6200 3.0140 3.6200 3.0140 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.4632 2 0.0490 3 0.2836 4 0.0550 5 0.0038 6 0.0142 7 0.0372 8 0.0282 9 0.0039 10 0.0213 11 0.0043 12 0.0362 MATRIX O This is a FULL matrix of order 1 by 4 1 2 3 4 1 -0.1040 -0.0840 -0.1040 -0.0840 MATRIX R This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 1.0000 0.7409 0.3320 0.2460 2 0.7409 1.0000 0.2460 0.2666 3 0.3320 0.2460 1.0000 0.7409 4 0.2460 0.2666 0.7409 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.7510 1.6110 1.7289 1.4859 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.0608 2 0.0608 3 0.0608 4 0.0608 5 0.0608 6 0.0608 7 0.0608 8 0.0608 9 0.0608 10 0.0608 11 0.0608 12 0.0608 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.0608 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: 1158.6849 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 10 estimated parameters and 675 Observed statistics -2 times log-likelihood of data >>> 2385.470 Degrees of freedom >>>>>>>>>>>>>>>> 665 This problem used 0.8% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.13 Execution 0: 0: 8:38.10 TOTAL 0: 0: 8:38.23 Total number of warnings issued: 4 ______________________________________________________________________________ !save corrace.mxs ! no correlation model !drop 2 4 !end ______________________________________________________________________________