The Annals of Applied Statistics
- Ann. Appl. Stat.
- Volume 6, Number 2 (2012), 432-456.
Simultaneous SNP identification in association studies with missing data
Association testing aims to discover the underlying relationship between genotypes (usually Single Nucleotide Polymorphisms, or SNPs) and phenotypes (attributes, or traits). The typically large data sets used in association testing often contain missing values. Standard statistical methods either impute the missing values using relatively simple assumptions, or delete them, or both, which can generate biased results. Here we describe the Bayesian hierarchical model BAMD (Bayesian Association with Missing Data). BAMD is a Gibbs sampler, in which missing values are multiply imputed based upon all of the available information in the data set. We estimate the parameters and prove that updating one SNP at each iteration preserves the ergodic property of the Markov chain, and at the same time improves computational speed. We also implement a model selection option in BAMD, which enables potential detection of SNP interactions. Simulations show that unbiased estimates of SNP effects are recovered with missing genotype data. Also, we validate associations between SNPs and a carbon isotope discrimination phenotype that were previously reported using a family based method, and discover an additional SNP associated with the trait. BAMD is available as an R-package from http://cran.r-project.org/package=BAMD.
Ann. Appl. Stat., Volume 6, Number 2 (2012), 432-456.
First available in Project Euclid: 11 June 2012
Permanent link to this document
Digital Object Identifier
Mathematical Reviews number (MathSciNet)
Zentralblatt MATH identifier
Li, Zhen; Gopal, Vikneswaran; Li, Xiaobo; Davis, John M.; Casella, George. Simultaneous SNP identification in association studies with missing data. Ann. Appl. Stat. 6 (2012), no. 2, 432--456. doi:10.1214/11-AOAS516. https://projecteuclid.org/euclid.aoas/1339419602
- Supplementary material: Theory and additional simulations. The Supplemental Information contains details on the variable selector, and the proof of convergence of the two Markov chains (the Gibbs sampler and the model search). In addition, there are further comparisons between BAMD and BIMBAM.