! ordinal_bivariate.mx The following MX script lines were read for group 1 #DEFINE NVAR 2 #DEFINE NVAR2 4 #DEFINE MAXTHRES 2 ANALYSIS OF ORDINAL DEPRESSION (0/1) AND SMOKING INITIATION/PERSISTENCE (0/1/2) DATA NI=NVAR2 NG=3 ORDINAL FI=DEPSMKMF.REC Ordinal data read initiated NOTE: Rectangular file contained 1317 records with data that contained a total of 4535 observations BEGIN MATRICES; M FU MAXTHRES NVAR FR L LO MAXTHRES MAXTHRES W LO NVAR NVAR FR X LO NVAR NVAR FI Y LO NVAR NVAR FR END MATRICES; MAT L 1.0 1.0 1.0 MATRIX M 0.5294 0.7191 0.0 0.5 SP M 1 2 0 3 ST 0.7 Y(1,1) Y(2,2) W(1,1) W(2,2) ST 0.0001 X(1,1) X(2,2) X(2,1) ST 0.2 W(2,1) Y(2,1) BEGIN ALGEBRA; A=W*W'; O=\STND(A); C=X*X'; R=\STND(C); E=Y*Y'; Q=\STND(E); P=A+C+E; END ALGEBRA; TH L*M|L*M; CO P | A + C _ A' + C' | P ; BO 0.001 1.0 Y(1,1) Y(2,2) BO 0.0001 0.999 X(1,1) X(2,2) W(1,1) W(2,2) BO -0.999 0.999 X(2,1) Y(2,1) W(2,1) BO 0.001 3.0 M(2,2) BO -5.0 5.0 M(1,1) ! interval a(1,1) a(2,2) c(1,1) c(2,2) e(1,1) e(2,2) o(1,2) r(1,2) q(1,2) OPT FUNC=1.E-12 OPT RS END The following MX script lines were read for group 2 ANALYSIS OF ORDINAL DEPRESSION AND SMOKING DATA: DZF DATA NI=NVAR2 ORDINAL FI=DEPSMKDF.REC Ordinal data read initiated NOTE: Rectangular file contained 877 records with data that contained a total of 2796 observations BEGIN MATRICES = GROUP 1; N FU MAXTHRES NVAR FR G FU 1 1 END MATRICES; MATRIX N 0.5781 0.6884 0 0.72 SP N 101 102 0 103 MAT G 0.5 TH L*N | L*N ; CO P | G@A + C _ G@A' + C' | P ; BO 0.001 3.0 N(2,2) BO -5.0 5.0 N(1,1) OPT RS END The following MX script lines were read for group 3 DATA CONSTRAINT CO NI=1 BEGIN MATRICES = GROUP 1; U UNIT 1 NVAR END MATRICES; CO \D2V(P) = U; END Summary of VL file data for group 1 Code 1.0000 2.0000 3.0000 4.0000 Number 1013.0000 1286.0000 982.0000 1254.0000 Mean 0.1925 0.6454 0.2169 0.6555 Variance 0.1554 0.7234 0.1699 0.7426 Minimum 0.0000 0.0000 0.0000 0.0000 Maximum 1.0000 2.0000 1.0000 2.0000 Summary of VL file data for group 2 Code 1.0000 2.0000 3.0000 4.0000 Number 598.0000 826.0000 586.0000 786.0000 Mean 0.1940 0.6525 0.2526 0.7468 Variance 0.1564 0.7086 0.1888 0.7820 Minimum 0.0000 0.0000 0.0000 0.0000 Maximum 1.0000 2.0000 1.0000 2.0000 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Analysis of ordinal depression (0/1) and smoking initiation/persistence (0/1/2) MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 2 1 2 1 1 2 2 0 3 MATRIX O This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Q This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX R This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX W This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 5 2 6 7 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 8 2 9 10 GROUP NUMBER: 2 Analysis of ordinal depression and smoking data: DZF MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX G This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX L This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 2 1 2 1 1 2 2 0 3 MATRIX N This is a FULL matrix of order 2 by 2 1 2 1 101 102 2 0 103 MATRIX O This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Q This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX R This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX W This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 5 2 6 7 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 8 2 9 10 GROUP NUMBER: 3 Data constraint MATRIX A This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX L This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified MATRIX M This is a FULL matrix of order 2 by 2 1 2 1 1 2 2 0 3 MATRIX O This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX P This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX Q This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX R This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX U This is a UNIT matrix of order 1 by 2 MATRIX W This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 5 2 6 7 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 It has no free parameters specified MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 8 2 9 10 Mx starting optimization; number of parameters = 12 MX PARAMETER ESTIMATES GROUP NUMBER: 1 Analysis of ordinal depression (0/1) and smoking initiation/persistence (0/1/2) MATRIX A This is a computed FULL matrix of order 2 by 2 [=W*W'] 1 2 1 0.3739 0.0779 2 0.0779 0.7728 MATRIX C This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 1.0000E-08 1.0000E-08 2 1.0000E-08 2.0000E-08 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.6261 0.0277 2 0.0277 0.2272 MATRIX L This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0000 2 1.0000 1.0000 MATRIX M This is a FULL matrix of order 2 by 2 1 2 1 0.8246 0.2694 2 0.0000 0.3990 MATRIX O This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.1449 2 0.1449 1.0000 MATRIX P This is a computed FULL matrix of order 2 by 2 [=A+C+E] 1 2 1 1.0000 0.1056 2 0.1056 1.0000 MATRIX Q This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 0.0734 2 0.0734 1.0000 MATRIX R This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 0.7071 2 0.7071 1.0000 MATRIX W This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6115 2 0.1274 0.8698 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0000E-04 2 1.0000E-04 1.0000E-04 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.7913 2 0.0350 0.4754 Matrix of EXPECTED thresholds 1 2 3 4 Threshold 1 0.8246 0.2694 0.8246 0.2694 Threshold 2 0.8246 0.6684 0.8246 0.6684 EXPECTED COVARIANCE MATRIX 1 2 3 4 1 1.0000 2 0.1056 1.0000 3 0.3739 0.0779 1.0000 4 0.0779 0.7728 0.1056 1.0000 Function value of this group: 6219.2085 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 2 Analysis of ordinal depression and smoking data: DZF MATRIX A This is a computed FULL matrix of order 2 by 2 [=W*W'] 1 2 1 0.3739 0.0779 2 0.0779 0.7728 MATRIX C This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 1.0000E-08 1.0000E-08 2 1.0000E-08 2.0000E-08 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.6261 0.0277 2 0.0277 0.2272 MATRIX G This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX L This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0000 2 1.0000 1.0000 MATRIX M This is a FULL matrix of order 2 by 2 1 2 1 0.8246 0.2694 2 0.0000 0.3990 MATRIX N This is a FULL matrix of order 2 by 2 1 2 1 0.7618 0.1733 2 0.0000 0.4495 MATRIX O This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.1449 2 0.1449 1.0000 MATRIX P This is a computed FULL matrix of order 2 by 2 [=A+C+E] 1 2 1 1.0000 0.1056 2 0.1056 1.0000 MATRIX Q This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 0.0734 2 0.0734 1.0000 MATRIX R This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 0.7071 2 0.7071 1.0000 MATRIX W This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6115 2 0.1274 0.8698 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0000E-04 2 1.0000E-04 1.0000E-04 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.7913 2 0.0350 0.4754 Matrix of EXPECTED thresholds 1 2 3 4 Threshold 1 0.7618 0.1733 0.7618 0.1733 Threshold 2 0.7618 0.6228 0.7618 0.6228 EXPECTED COVARIANCE MATRIX 1 2 3 4 1 1.0000 2 0.1056 1.0000 3 0.1869 0.0390 1.0000 4 0.0390 0.3864 0.1056 1.0000 Function value of this group: 4295.6685 Where the fit function is -2 * Log-likelihood of raw ordinal GROUP NUMBER: 3 Data constraint MATRIX A This is a computed FULL matrix of order 2 by 2 [=W*W'] 1 2 1 0.3739 0.0779 2 0.0779 0.7728 MATRIX C This is a computed FULL matrix of order 2 by 2 [=X*X'] 1 2 1 1.0000E-08 1.0000E-08 2 1.0000E-08 2.0000E-08 MATRIX E This is a computed FULL matrix of order 2 by 2 [=Y*Y'] 1 2 1 0.6261 0.0277 2 0.0277 0.2272 MATRIX L This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0000 2 1.0000 1.0000 MATRIX M This is a FULL matrix of order 2 by 2 1 2 1 0.8246 0.2694 2 0.0000 0.3990 MATRIX O This is a computed FULL matrix of order 2 by 2 [=\STND(A)] 1 2 1 1.0000 0.1449 2 0.1449 1.0000 MATRIX P This is a computed FULL matrix of order 2 by 2 [=A+C+E] 1 2 1 1.0000 0.1056 2 0.1056 1.0000 MATRIX Q This is a computed FULL matrix of order 2 by 2 [=\STND(E)] 1 2 1 1.0000 0.0734 2 0.0734 1.0000 MATRIX R This is a computed FULL matrix of order 2 by 2 [=\STND(C)] 1 2 1 1.0000 0.7071 2 0.7071 1.0000 MATRIX U This is a UNIT matrix of order 1 by 2 MATRIX W This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.6115 2 0.1274 0.8698 MATRIX X This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0000E-04 2 1.0000E-04 1.0000E-04 MATRIX Y This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.7913 2 0.0350 0.4754 Your model has 12 estimated parameters and 7333 Observed statistics Observed statistics include 2 constraints. -2 times log-likelihood of data >>> 10514.877 Degrees of freedom >>>>>>>>>>>>>>>> 7321 This problem used 11.7% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 1.42 Execution 0: 0: 0: 8.58 TOTAL 0: 0: 0:10.00 Total number of warnings issued: 0 ______________________________________________________________________________