Abstract
Count data are increasingly ubiquitous in genetic association studies, where it is possible to observe excess zero counts as compared to what is expected based on standard assumptions. For instance, in rheumatology, data are usually collected in multiple joints within a person or multiple sub-regions of a joint, and it is not uncommon that the phenotypes contain enormous number of zeroes due to the presence of excessive zero counts in majority of patients. Most existing statistical methods assume that the count phenotypes follow one of these four distributions with appropriate dispersion-handling mechanisms: Poisson, Zero-inflated Poisson (ZIP), Negative Binomial, and Zero-inflated Negative Binomial (ZINB). However, little is known about their implications in genetic association studies. Also, there is a relative paucity of literature on their usefulness with respect to model misspecification and variable selection. In this article, we have investigated the performance of several state-of-the-art approaches for handling zero-inflated count data along with a novel penalized regression approach with an adaptive LASSO penalty, by simulating data under a variety of disease models and linkage disequilibrium patterns. By taking into account data-adaptive weights in the estimation procedure, the proposed method provides greater flexibility in multi-SNP modeling of zero-inflated count phenotypes. A fast coordinate descent algorithm nested within an EM (expectation-maximization) algorithm is implemented for estimating the model parameters and conducting variable selection simultaneously. Results show that the proposed method has optimal performance in the presence of multicollinearity, as measured by both prediction accuracy and empirical power, which is especially apparent as the sample size increases. Moreover, the Type I error rates become more or less uncontrollable for the competing methods when a model is misspecified, a phenomenon routinely encountered in practice.
Highlights
In genetic association studies, phenotypes are often measured in counts
In Section Discussion, we provide guidelines to researchers and practitioners alike for using zero-inflated count models arising in genetic association studies including some concluding remarks as well as directions for future research
We propose the EM Adaptive Least Absolute Shrinkage And Selection Operator (LASSO) (AL) estimator by modifying the EM algorithm developed by Wang et al (2015) for their LASSO estimator, which is implemented in the R (R Core Team, 2015) package mpath (Wang, 2015)
Summary
Phenotypes are often measured in counts. For such count data, the standard methods for modeling the genotype–phenotype relationship include the Poisson and Negative Binomial regression models. Some of the phenotypes of interest in arthritis research are calculated based on a summary of individual integervalued count measures (Teare et al, 2013), and it is highly plausible that they contain enormous number of zeroes due to the preponderance of zero counts in majority of patients. Some examples of such count phenotypes include the maximum number of pain sites (Holliday et al, 2010), total number of pain sites (Nicholl et al, 2011), and modified Larsen score (Teare et al, 2013), among others. These two standard approaches [viz. the Poisson regression (PR) and the Negative Binomial (NB) regression] fail to take into account the added variability associated with the extraneous zero observations
Sign-in/Register to view highlights & summary Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
