Fitting Genetic Models Using Markov Chain Monte Carlo Algorithms With BUGS

Maximum likelihood estimation techniques are widely used in twin and family studies, but soon reach computational boundaries when applied to highly complex models (e.g., models including gene-by-environment interaction and gene-environment correlation, item response theory measurement models, repeated measures, longitudinal structures, extended pedigrees). Markov Chain Monte Carlo (MCMC) algorithms are very well suited to fit complex models with hierarchically structured data. This article introduces the key concepts of Bayesian inference and MCMC parameter estimation and provides a number of scripts describing relatively simple models to be estimated by the freely obtainable BUGS software. In addition, inference using BUGS is illustrated using a data set on follicle-stimulating hormone and luteinizing hormone levels with repeated measures. The examples provided can serve as stepping stones for more complicated models, tailored to the specific needs of the individual researcher.

[1]  Alaattin Erkanli,et al.  Markov Chain Monte Carlo Approaches to Analysis of Genetic and Environmental Components of Human Developmental Change and G × E Interaction , 2003 .

[2]  Dorret I. Boomsma,et al.  Variance Decomposition Using an IRT Measurement Model , 2007, Behavior genetics.

[3]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[4]  Lindon Eaves,et al.  Application of Bayesian Inference using Gibbs Sampling to Item-Response Theory Modeling of Multi-Symptom Genetic Data , 2005, Behavior genetics.

[5]  G. C. Tiao,et al.  Bayesian inference in statistical analysis , 1973 .

[6]  A. Brix.  Bayesian Data Analysis, 2nd edn , 2005 .

[7]  N. Horton.  Multilevel and Longitudinal Modeling Using Stata , 2006 .

[8]  Stan Lipovetsky,et al.  Generalized Latent Variable Modeling: Multilevel,Longitudinal, and Structural Equation Models , 2005, Technometrics.

[9]  Dorret I. Boomsma,et al.  Individual Differences in Puberty Onset in Girls: Bayesian Estimation of Heritabilities and Genetic Correlations , 2006, Behavior genetics.

[10]  L. Crittenden.  AN INTERPRETATION OF FAMILIAL AGGREGATION BASED ON MULTIPLE GENETIC AND ENVIRONMENTAL FACTORS , 1961, Annals of the New York Academy of Sciences.

[11]  Katrina J Scurrah,et al.  Covariance components models for longitudinal family data. , 2005, International journal of epidemiology.

[12]  A. Gelfand,et al.  Efficient parametrisations for normal linear mixed models , 1995 .

[13]  D I Boomsma,et al.  Using LISREL to analyze genetic and environmental covariance structure , 1986, Behavior genetics.

[14]  A. Erkanli,et al.  Genetic and environmental influences on the relative timing of pubertal change. , 2004, Twin research : the official journal of the International Society for Twin Studies.

[15]  Hoon Kim,et al.  Monte Carlo Statistical Methods , 2000, Technometrics.

[16]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[17]  Mark Von Tress,et al.  Generalized, Linear, and Mixed Models , 2003, Technometrics.

[18]  D W Fulker,et al.  Comparison of the biometrical genetical, MAVA, and classical approaches to the analysis of human behavior. , 1970, Psychological bulletin.

[19]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[20]  G. Willemsen,et al.  Netherlands Twin Register: A Focus on Longitudinal Research , 2002, Twin Research.

[21]  D. Falconer.  The inheritance of liability to certain diseases, estimated from the incidence among relatives , 1965 .