Analysis of Multiple Binary Responses Using a Threshold Model

Several discrete responses, such as health status, reproduction performance and meat quality, are routinely collected for several livestock species. These traits are often of binary or discrete nature. Genetic evaluation for these traits is frequently conducted using a single-trait threshold model, or they are considered continuous responses either in univariate or in multivariate context. Implementation of threshold models in the presence of several binary responses or a mixture of binary and continuous responses is far from simple. The complexity of such implementation is primarily due to the incomplete randomness of the residual (co)variance matrix. In the current study, a multiple binary trait simulation was carried out in order to implement and validate a new procedure for dealing with the consequences of the restrictions imposed to the residual variance using threshold models. Using three and eight binary responses, the proposed method was able to estimate all unknown parameters without any noticeable bias. In fact, for simulated residual correlations ranging from −0.8 to 0.8, the resulting HPD 95% intervals included the true values in all cases. The proposed procedure involved limited additional computational cost and is straightforward to implement independent of the number of binary responses involved in the analysis. Monitoring of the convergence of the procedure must be conducted at the identifiable scale, and special care must be placed on the selection of the prior of the non-identifiable model. The latter could have serious consequences on the final results due to potential truncation of the parameter space.