Prediction of patient reported outcome measures via multivariate ordered probit models