! ordinal_bivariate.mx #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=depsmkmm.rec Begin matrices; M FU maxthres nvar fr L LO maxthres maxthres W LO nvar nvar fr X LO nvar nvar fr 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.2 x(1,1) x(2,2) st 0.2 w(2,1) x(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 Analysis of ordinal depression and smoking data: DZF data NI=nvar2 Ordinal fi=depsmkdm.rec 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 Data constraint CO NI=1 Begin matrices = group 1; U unit 1 nvar end matrices; CO \d2v(P) = u; end