Multivariate probit model for binary/ordinal response