Open Access

Multivariate Bayesian analysis of Gaussian, right censored Gaussian, ordered categorical and binary traits using Gibbs sampling

  • Inge Riis Korsgaard1Email author,
  • Mogens Sandø Lund1,
  • Daniel Sorensen1,
  • Daniel Gianola2,
  • Per Madsen1 and
  • Just Jensen1
Genetics Selection Evolution200335:159

Received: 5 October 2001

Accepted: 3 September 2002

Published: 15 March 2003


A fully Bayesian analysis using Gibbs sampling and data augmentation in a multivariate model of Gaussian, right censored, and grouped Gaussian traits is described. The grouped Gaussian traits are either ordered categorical traits (with more than two categories) or binary traits, where the grouping is determined via thresholds on the underlying Gaussian scale, the liability scale. Allowances are made for unequal models, unknown covariance matrices and missing data. Having outlined the theory, strategies for implementation are reviewed. These include joint sampling of location parameters; efficient sampling from the fully conditional posterior distribution of augmented data, a multivariate truncated normal distribution; and sampling from the conditional inverse Wishart distribution, the fully conditional posterior distribution of the residual covariance matrix. Finally, a simulated dataset was analysed to illustrate the methodology. This paper concentrates on a model where residuals associated with liabilities of the binary traits are assumed to be independent. A Bayesian analysis using Gibbs sampling is outlined for the model where this assumption is relaxed.


categoricalGaussianmultivariate Bayesian analysisright censored Gaussian

(To access the full article, please see PDF)

Authors’ Affiliations

Department of Animal Breeding and Genetics, Danish Institute of Agricultural Sciences
Department of Meat and Animal Sciences, University of Wisconsin-Madison


© INRA, EDP Sciences 2003