Metropolis sampling in pedigree analysis

This paper reviews and develops the applications of Markov chain Monte Carlo methods in pedigree analysis, with particular stress on the Metropolis algorithm. In likelihood based genetic analyses, standard deterministic algorithms often fail because of the computational complexity of the observed pedigree data under a proposed genetic model. The new Monte Carlo methods permit approximate maximum likelihood estimation in the presence of such complexity. Monte Carlo implementation of the EM algorithm is the key to successful maximum likelihood analysis. Gibbs sampling and the Metropolis algorithm are alternative ways of defining Markov chains for performing the E step of the EM algorithm. Two applications illustrate the power and simplicity of the Metropolis algorithm. One of these applications involves a discrete model for variance component analysis of quantitative traits; the other application involves a Monte Carlo version of location scores for multipoint linkage analysis.

[1]  E A Thompson,et al.  A Monte Carlo method for combined segregation and linkage analysis. , 1992, American journal of human genetics.

[2]  D. Thomas,et al.  A Gibbs sampling approach to linkage analysis. , 1992, Human heredity.

[3]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[4]  J. Ott,et al.  Maximum likelihood estimation by counting methods under polygenic and mixed models in human pedigrees. , 1979, American journal of human genetics.

[5]  N. Sheehan,et al.  On the irreducibility of a Markov chain defined on a space of genotype configurations by a sampling scheme. , 1993, Biometrics.

[6]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[7]  K. Lange,et al.  Programs for pedigree analysis: Mendel, Fisher, and dGene , 1988, Genetic epidemiology.

[8]  J. Ott.  Computer-simulation methods in human linkage analysis. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Charles J. Geyer,et al.  Reweighting Monte Carlo Mixtures , 1991 .

[10]  D E Weeks,et al.  Efficient computation of lod scores: genotype elimination, genotype redefinition, and hybrid maximum likelihood algorithms , 1989, Annals of human genetics.

[11]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[12]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[13]  M. Boehnke,et al.  Estimating the power of a proposed linkage study for a complex genetic trait. , 1989, American journal of human genetics.

[14]  E. Lander,et al.  Construction of multilocus genetic linkage maps in humans. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[15]  G. Lathrop,et al.  Construction of human linkage maps: Likelihood calculations for multilocus linkage analysis , 1986, Genetic epidemiology.

[16]  T. Goradia,et al.  Localization of an ataxia-telangiectasia gene to an ∼500-kb interval on chromosome 11q23.1 : linkage analysis of 176 families by an international consortium , 1991 .

[17]  H. Cann,et al.  Centre d'etude du polymorphisme humain (CEPH): collaborative genetic mapping of the human genome. , 1990, Genomics.

[18]  Probability functions on complex pedigrees , 1978 .

[19]  Emile H. L. Aarts,et al.  Simulated Annealing: Theory and Applications , 1987, Mathematics and Its Applications.

[20]  N. A. Barnicot.  Red hair in African Negroes; a preliminary study. , 1953, Annals of eugenics.

[21]  Calyampudi R. Rao,et al.  Linear statistical inference and its applications , 1965 .

[22]  E A Thompson,et al.  Evaluation of likelihood ratios for complex genetic models. , 1991, IMA journal of mathematics applied in medicine and biology.

[23]  T. Goradia,et al.  Localization of an ataxia-telangiectasia locus to a 3-cM interval on chromosome 11q23: linkage analysis of 111 families by an international consortium. , 1991, American journal of human genetics.

[24]  E. Thompson,et al.  Monte Carlo estimation of variance component models for large complex pedigrees. , 1991, IMA journal of mathematics applied in medicine and biology.

[25]  R. Elston,et al.  A general model for the genetic analysis of pedigree data. , 1971, Human heredity.

[26]  N E Morton,et al.  Analysis of family resemblance. 3. Complex segregation of quantitative traits. , 1974, American journal of human genetics.

[27]  J. Ott,et al.  Strategies for multilocus linkage analysis in humans. , 1984, Proceedings of the National Academy of Sciences of the United States of America.

[28]  K. Lange,et al.  Bayesian methods and optimal experimental design for gene mapping by radiation hybrids , 1992, Annals of human genetics.

[29]  S. B. Holt,et al.  Genetics of dermal ridges; bilateral asymmetry in finger ridge-counts. , 1954, Annals of eugenics.

[30]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[31]  K. Lange,et al.  Extensions to pedigree analysis. IV. Covariance components models for multivariate traits. , 1983, American journal of medical genetics.

[32]  K Lange,et al.  Simulation of pedigree genotypes by random walks. , 1989, American journal of human genetics.

[33]  K Lange,et al.  A random walk method for computing genetic location scores. , 1991, American journal of human genetics.

[34]  Andrew L. Rukhin,et al.  Tools for statistical inference , 1991 .

[35]  P M Nadkarni,et al.  Fast computation of genetic likelihoods on human pedigree data. , 1992, Human heredity.