Fitting genetic models using Markov Chain Monte Carlo algorithms with BUGS

S.M. van den Berg, A.L. Beem, D.I. Boomsma

Research output: Contribution to JournalArticleAcademicpeer-review

214 Downloads (Pure)


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.
Original languageEnglish
Pages (from-to)334-342
JournalTwin Research and Human Genetics
Issue number3
Publication statusPublished - 2006

Cohort Studies

  • Netherlands Twin Register (NTR)


Dive into the research topics of 'Fitting genetic models using Markov Chain Monte Carlo algorithms with BUGS'. Together they form a unique fingerprint.

Cite this