Bayesian prediction of breeding values for multivariate binary and continuous traits in simulated horse populations using threshold-linear models with Gibbs sampling.