Approximate Bayesian inference for cumulative probit regression models