! 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
______________________________________________________________________________