Perhaps I didn't explain my problem clearly. Sorry about the confusion. What I want is one model for 3 Ys, all 3 Ys are binary, and the error terms are correlated. And, the independent variables in the 3 equations are not the same. the model could be described as follow: Y1 = X1B1 + e1; Y2 = X2*B2 +e2; Y3=X3*B3+e3; and (e1,e2,e3|x) ~ N((0,0,0), (1,rho12, rho13, rho12, 1, rho23 rho13, rho23, 1) Basically, I want to use a similar function as proc qlim (for 2 binary outcomes), but I want to work for at least 3 binary outcomes. Thank you
... View more