Multivariate probit model