EDP Sciences Journals List
Free access article

Issue Genet. Sel. Evol.
Volume 39, Number 1, January-February 2007
Page(s) 27 - 38
DOI 10.1051/gse:2006032

Genet. Sel. Evol. 39 (2007) 27-38
DOI: 10.1051/gse:2006032

Improved techniques for sampling complex pedigrees with the Gibbs sampler

K. Joseph Abrahama, Liviu R. Totirb and Rohan L. Fernandob, c

a  1301 Agronomy Hall, Iowa State University, Ames, IA 50011, USA
b  Department of Animal Science, Iowa State University, Ames, IA 50011, USA
c  Lawrence H. Baker Center for Bioinformatics and Biological Statistics, Iowa State University, Ames, IA 50011, USA

(Received 5 January 2006; accepted 3 September 2006; published online 11 January 2007)

Abstract - Markov chain Monte Carlo (MCMC) methods have been widely used to overcome computational problems in linkage and segregation analyses. Many variants of this approach exist and are practiced; among the most popular is the Gibbs sampler. The Gibbs sampler is simple to implement but has (in its simplest form) mixing and reducibility problems; furthermore in order to initiate a Gibbs sampling chain we need a starting genotypic or allelic configuration which is consistent with the marker data in the pedigree and which has suitable weight in the joint distribution. We outline a procedure for finding such a configuration in pedigrees which have too many loci to allow for exact peeling. We also explain how this technique could be used to implement a blocking Gibbs sampler.


Key words: Gibbs sampler / Markov chain Monte Carlo / pedigree peeling / Elston Stewart algorithm

Correspondence and reprints: abraham@iastate.edu

© INRA, EDP Sciences 2006


What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.
  • If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
  • You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
  • You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.