Statistical Applications in Genetics and Molecular Biology
Editor-in-Chief: Stumpf, Michael P.H.
IMPACT FACTOR increased in 2015: 1.265
5-year IMPACT FACTOR: 1.423
Rank 42 out of 123 in category Statistics & Probability in the 2015 Thomson Reuters Journal Citation Report/Science Edition
SCImago Journal Rank (SJR) 2014: 0.740
Source Normalized Impact per Paper (SNIP) 2014: 0.470
Impact per Publication (IPP) 2014: 0.926
Mathematical Citation Quotient (MCQ) 2014: 0.17
Locating Multiple Interacting Quantitative Trait Loci with the Zero-Inflated Generalized Poisson Regression
1Technische Universität München
2Wroclaw University of Technology and Purdue University
3Technische Universität München
Citation Information: Statistical Applications in Genetics and Molecular Biology. Volume 9, Issue 1, ISSN (Online) 1544-6115, DOI: 10.2202/1544-6115.1545, June 2010
- Published Online:
We consider the problem of locating multiple interacting quantitative trait loci (QTL) influencing traits measured in counts. In many applications the distribution of the count variable has a spike at zero. Zero-inflated generalized Poisson regression (ZIGPR) allows for an additional probability mass at zero and hence an improvement in the detection of significant loci. Classical model selection criteria often overestimate the QTL number. Therefore, modified versions of the Bayesian Information Criterion (mBIC and EBIC) were successfully used for QTL mapping. We apply these criteria based on ZIGPR as well as simpler models. An extensive simulation study shows their good power detecting QTL while controlling the false discovery rate. We illustrate how the inability of the Poisson distribution to account for over-dispersion leads to an overestimation of the QTL number and hence strongly discourages its application for identifying factors influencing count data. The proposed method is used to analyze the mice gallstone data of Lyons et al. (2003). Our results suggest the existence of a novel QTL on chromosome 4 interacting with another QTL previously identified on chromosome 5. We provide the corresponding code in R.