! ordinal_bivariate_simulated.mx #define nvar 2 #define nvar2 4 #define maxthres 1 Analysis of ordinal depression (0/1) and smoking initiation/persistence (0/1/2) data NI=5 NG=3 LA twina1 twina2 twinb1 twinb2 mzwt Ordinal fi=mzfull.rec definition_variables mzwt / 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 Z FU 1 1 ! weight matrix end matrices; SP Z -1 MAT L 1.0 MATRIX M 0.65 0.65 SP M 1 2 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; FREQ Z; 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 -5.0 5.0 m(1,1) m(1,2) 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=5 LA twina1 twina2 twinb1 twinb2 dzwt Ordinal fi=dzfull.rec definition_variables dzwt/ Begin matrices = group 1; N FU maxthres nvar fr U fu 1 1 ! weight matrix g fu 1 1 end matrices; SP U -1 MATRIX N 0.65 0.65 SP N 101 102 mat g 0.5 FREQ U; TH L*N | L*N ; CO P | g@A + C _ g@A' + C' | P ; bo -5.0 5.0 n(1,1) n(1,2) OPT RS END Data constraint CO NI=1 Begin matrices = group 1; U unit 1 nvar end matrices; CO \d2v(P) = u; end