R code to optimize probit, ordered probit, bivariate probit, and bivariate ordered probit models