! ordinal2_bivariate_general.mx ! Fitting two-stage model to U.S. Veterans Twins data ! #define nvar 2 #define nvar2 4 #define maxthres 1 Analysis of regular smoking and smoking persistence data: MZM data NI=5 NG=3 la wtmz init persist xinit xersist Ordinal fi=smkmz.rec definition_variables wtmz / Begin matrices; T DI 1 1 ! Used to read weight variable M FU maxthres nvar fr ! Threshold values W LO nvar nvar fr ! additive genetic path X LO nvar nvar fr ! shared environmental path Y LO nvar nvar fr ! non-shared environment path end matrices; SP T -1 FI Y(2,1) MATRIX M -0.453 0.000 mat w 0.73 0.34 0.68 mat x 0.52 0.24 0.0 mat y 0.44 0.0 0.61 Begin algebra; A=W*W'; Q=\stnd(A); C=X*X'; R=\stnd(C); E=Y*Y'; S=\stnd(E); P=A+C+E; end algebra; TH M|M; FREQ T; CO ( P | A + C _ A' + C'| P ) ; bo 0.001 0.8 y(1,1) y(2,2) bo 0.0001 0.9 x(1,1) x(2,2) w(1,1) w(2,2) bo -0.9 0.9 x(2,1) y(2,1) bo -5.0 5.0 m(1,1) m(1,2) ! interval a(1,1) a(2,2) q(1,2) ! c(1,1) c(2,2) ! r(1,2) interval e(1,1) e(2,2) ! s(1,2) OPT func=1.E-12 OPT RS END Analysis of ordinal alcohol tolerance and dependence data: DZF data NI=5 la wtdz init persist xinit xersist OR fi=smkdz.rec definition_variables wtdz/ Begin matrices = group 1; N FU maxthres nvar fr g fu 1 1 U DI 1 1 ! Used to read weight variable end matrices; sp U -1 MATRIX N -0.55 -0.04 mat g 0.5 TH N|N; FREQ U; CO ( P | g@A + C _ g@A' + C'| P ) ' ; bo -5.0 5.0 n(1,1) n(1,2) OPT RS END Constraint function - constrain phenotypic variances to unity CO NI=1 Begin matrices = group 1; U unit 1 nvar end matrices; CO \d2v(P) = u; end