Mixed Logistic Regression in Genome-Wide Association Studies

Motivation Mixed linear models (MLM) have been widely used to account for population structure in case-control genome-wide association studies, the status being analyzed as a quantitative phenotype. Chen et al. proved that this method is inappropriate and proposed a score test for the mixed logistic regression (MLR). However this test does not allow an estimation of the variants’ effects. Results We propose two computationally efficient methods to estimate the variants’ effects. Their properties are evaluated on two simulations sets, and compared with other methods (MLM, logistic regression). MLR performs the best in all circumstances. The variants’ effects are well evaluated by our methods, with a moderate bias when the effect sizes are large. Additionally, we propose a stratified QQ-plot, enhancing the diagnosis of p -values inflation or deflation, when population strata are not clearly identified in the sample. Availability All methods are implemented in the R package milorGWAS available at https://github.com/genostats/milorGWAS . Contact [email protected] Supplementary information Supplementary data are available at Bioinformatics online.