A Bivariate Ordered Probit Estimator with Mixed Effects