Gibbs sampling, adaptive rejection sampling and robustness to prior specification for a mixed linear model


Theobald C. M., Firat M. Z., Thompson R.

GENETICS SELECTION EVOLUTION, cilt.29, sa.1, ss.57-72, 1997 (SCI-Expanded) identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 29 Sayı: 1
  • Basım Tarihi: 1997
  • Doi Numarası: 10.1051/gse:19970104
  • Dergi Adı: GENETICS SELECTION EVOLUTION
  • Derginin Tarandığı İndeksler: Science Citation Index Expanded (SCI-EXPANDED), Scopus
  • Sayfa Sayıları: ss.57-72
  • Anahtar Kelimeler: adaptive rejection sampling, Bayesian inference, Gibbs sampling, mixed linear model, sire model, BAYESIAN-INFERENCE, PARAMETERS
  • Akdeniz Üniversitesi Adresli: Hayır

Özet

Markov chain Monte-Carlo methods are increasingly being applied to make inferences about the marginal posterior distributions of parameters in quantitative genetic models. This paper considers the application of one such method, Gibbs sampling, to Bayesian inferences about parameters in a normal mixed linear model when a restriction is imposed on the relative values of the variance components. Two prior distributions are proposed which incorporate this restriction. For one of them, the adaptive rejection sampling technique is used to sample from the conditional posterior distribution of a variance ratio. Simulated data from a balanced sire model are used to compare different implementations of the Gibbs sampler and also inferences based on the two prior distributions. Marginal posterior distributions of the parameters are illustrated. Numerical results suggest that it is not necessary to discard any iterates, and that similar inferences are made using the two prior specifications.

Markov chain Monte-Carlo methods are increasingly being applied to make inferences about the marginal posterior distributions of parameters in quantitative genetic models. This paper considers the application of one such method, Gibbs sampling, to Bayesian inferences about parameters in a normal mixed linear model when a restriction is imposed on the relative values of the variance components. Two prior distributions are proposed which incorporate this restriction. For one of them, the adaptive rejection sampling technique is used to sample from the conditional posterior distribution of a variance ratio. Simulated data from a balanced sire model are used to compare different implementations of the Gibbs sampler and also inferences based on the two prior distributions. Marginal posterior distributions of the parameters are illustrated. Numerical results suggest that it is not necessary to discard any iterates, and that similar inferences are made using the two prior specifications.