Bayesian estimation in animal breeding using the Dirichlet process prior for correlated random effects

In the case of the mixed linear model the random effects are usually assumed to be normally distributed in both the Bayesian and classical frameworks. In this paper, the Dirichlet process prior was used to provide nonparametric Bayesian estimates for correlated random effects. This goal was achieved by providing a Gibbs sampler algorithm that allows these correlated random effects to have a nonparametric prior distribution. A sampling based method is illustrated. This method which is employed by transforming the genetic covariance matrix to an identity matrix so that the random effects are uncorrelated, is an extension of the theory and the results of previous researchers. Also by using Gibbs sampling and data augmentation a simulation procedure was derived for estimating the precision parameter M associated with the Dirichlet process prior. All needed conditional posterior distributions are given. To illustrate the application, data from the Elsenburg Dormer sheep stud were analysed. A total of 3325 weaning weight records from the progeny of 101 sires were used.

[1]  Hani Doss.  Bayesian Nonparametric Estimation for Incomplete Data Via Successive Substitution Sampling , 1994 .

[2]  P. Müller,et al.  Bayesian curve fitting using multivariate normal mixtures , 1996 .

[3]  D. Gianola,et al.  Mixed effects linear models with t-distributions for quantitative genetic analysis: a Bayesian approach , 1999, Genetics Selection Evolution.

[4]  Joseph G. Ibrahim,et al.  Semiparametric Bayesian Methods for Random Effects Models , 1998 .

[5]  Adrian F. M. Smith,et al.  Bayesian Analysis of Constrained Parameter and Truncated Data Problems , 1991 .

[6]  M. Escobar.  Estimating Normal Means with a Dirichlet Process Prior , 1994 .

[7]  M. West,et al.  Hyperparameter estimation in Dirichlet process mixture models , 1992 .

[8]  Jun S. Liu.  Nonparametric hierarchical Bayes via sequential imputations , 1996 .

[9]  S. MacEachern,et al.  Estimating mixture of dirichlet process models , 1998 .

[10]  S. MacEachern,et al.  A semiparametric Bayesian model for randomised block designs , 1996 .

[11]  C. Antoniak.  Mixtures of Dirichlet Processes with Applications to Bayesian Nonparametric Problems , 1974 .

[12]  A. Rukhin.  Bayes and Empirical Bayes Methods for Data Analysis , 1997 .

[13]  S. MacEachern.  Estimating normal means with a conjugate style dirichlet process prior , 1994 .

[14]  T. Ferguson.  A Bayesian Analysis of Some Nonparametric Problems , 1973 .

[15]  Adrian F. M. Smith,et al.  Gibbs Sampling for Marginal Posterior Expectations , 1991 .

[16]  M. Escobar,et al.  Bayesian Density Estimation and Inference Using Mixtures , 1995 .

[17]  D. Gianola,et al.  Attenuating effects of preferential treatment with Student-t mixed linear models: a simulation study , 1998, Genetics Selection Evolution.

[18]  J G Ibrahim,et al.  A semiparametric Bayesian approach to the random effects model. , 1998, Biometrics.

[19]  Bradley P. Carlin,et al.  BAYES AND EMPIRICAL BAYES METHODS FOR DATA ANALYSIS , 1996, Stat. Comput..