R Code to Call WinBugs and Estimate Bayesian Ordered Probit and Bivariate Probit Models