Location of Repository

Evaluating the performance of a block updating McMC sampler in a simple genetic application

By Nuala A. Sheehan, B. Guldbrandtsen and D.A. Sorensen

Abstract

This is the authors' final draft of the paper published as Journal of Agricultural, Biological, and Environmental Statistics, 2007, 12(2), pp.272-99. The definitive published version is available from http://lysander.asa.catchword.org/vl=437914/cl=24/nw=1/rpsv/cw/asa/10857117/v12n2/s8/p272Markov chain Monte Carlo (McMC) methods have provided an enormous breakthrough in the analysis of large complex problems such as those which frequently arise in genetic applications. The richness of the inference and the flexibility of an McMC Bayesian\ud approach in terms of design, data structure that can be analysed, and models that can be posed, is indisputable. However, despite the\ud strengths of the Bayesian approach, it is important to acknowledge that there are other, often easier, ways of tackling a problem. This is so, especially when simpler, qualitative answers are sought, such\ud as presence or absence of a quantitative trait locus. We critically evaluate the behaviour of a Bayesian McMC block sampler for the detection of a quantitative trait locus by linkage with marker data, and\ud compare it with a traditional least squares method. Some practical issues are illustrated by discussing the pros and cons of a Bayesian block updating sampling scheme versus the least squares method in\ud the context of a simple genetic mapping problem. Depending on the focus of analysis, we show that the McMC sampler does not always\ud outperform the simpler approach from a frequentist perspective, and, more to the point, may not always perform appropriately in any particular replication

Publisher: American Statistical Association and the International Biometric Society
Year: 2007
OAI identifier: oai:lra.le.ac.uk:2381/2663

Suggested articles

Preview

Citations

  1. (1996). A Monte Carlo method for Bayesian analysis of linkage between single markers and quantitative trait loci I. doi
  2. (1996). An approximation to the likelihood for a pedigree with loops, doi
  3. (1999). Analysis of Human Genetic Linkage, third edn, The Johns Hopkins doi
  4. (1998). Bayesian mapping of multiple quantitative trait loci from incomplete inbred line cross data,
  5. (2000). Bayesian mapping of quantitative trait loci under the identity-by-descent-based variance component model,
  6. (2002). Bayesian MCMC mapping of quantitative trait loci in a half-sib design: a graphical model perspective, doi
  7. (1996). Descent graphs in pedigree analysis: applications to haplotyping, location scores, and marker-sharing statistics,
  8. (1994). Empirical threshold values for quantitative trait mapping,
  9. (1953). Equations of state calculations by fast computing machines, doi
  10. (1992). Estimating the locations and the sizes of the effects of quantitative trait loci using flanking markers, doi
  11. (2000). Estimation of conditional multilocus gene identity among relatives, doi
  12. (1989). Hugin—a shell for building Bayesian belief universes for expert systems, doi
  13. (1994). Monte Carlo likelihood in genetic mapping, doi
  14. (2001). Monte Carlo methods on genetic structures, in doi
  15. (2000). Multilocus linkage analysis by blocked Gibbs sampling,
  16. (2002). On a misconception about irreducibility of the single-site gibbs sampler in a pedigree application,
  17. (1993). On the irreducibility of a Markov chain defined on a space of genotype configurations by a sampling scheme, doi
  18. (1999). Replication of linkage studies of complex traits: An examination of variation in location estimates, doi
  19. (1992). Sampling genotypes on complex pedigrees with phenotypic constraints: the origin of the B allele among the Polar Eskimos, doi
  20. (2000). Statistical Inference from Genetic Data on Pedigrees,
  21. (1998). Statistics in Human Genetics, Applications of Statistics Series, doi

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.