Bayesian linear mixed models with polygenic effects
Accepted version
Peer-reviewed
Repository URI
Repository DOI
Change log
Authors
Abstract
We considered Bayesian estimation of polygenic effects, in particular heritability in relation to a class of linear mixed models implemented in R. Our approach is applicable to both family-based and population-based studies in human genetics with which a genetic relationship matrix can be derived either from family structure or genome-wide data. Using a simulated and a real data, we demonstrate our implementation of the models in the generic statistical software systems JAGS and Stan as well as several R packages. In doing so, we have not only provided facilities in R linking standalone programs such as GCTA and other packages in R but also addressed some technical issues in the analysis. Our experience with a host of general and special software systems will facilitate investigation into more complex models for both human and nonhuman genetics.
Description
Keywords
Journal Title
Conference Name
Journal ISSN
1548-7660