Abstract:Abstract Motivation Sparse regularized regression methods are now widely used in genome-wide association studies (GWAS) to address the multiple testing burden that limits discovery of potentially important predictors. Linear mixed models (LMMs) have become an attractive alternative to principal components (PCs) adjustment to account for population structure and relatedness in high-dimensional penalized models. However, their use in binary trait GWAS rely on the invalid assumption that the residual variance does not depend on the estimated regression coefficients. Moreover, LMMs use a single spectral decomposition of the covariance matrix of the responses, which is no longer possible in generalized linear mixed models (GLMMs). Results We introduce a new method called pglmm, a penalized GLMM that allows to simultaneously select genetic markers and estimate their effects, accounting for between-individual correlations and binary nature of the trait. We develop a computationally efficient algorithm based on penalized quasi-likelihood estimation that allows to scale regularized mixed models on high-dimensional binary trait GWAS. We show through simulations that when the dimensionality of the relatedness matrix is high, penalized LMM and logistic regression with PC adjustment fail to select important predictors, and have inferior prediction accuracy compared to pglmm. Further, we demonstrate through the analysis of two polygenic binary traits in a subset of 6731 related individuals from the UK Biobank data with 320K SNPs that our method can achieve higher predictive performance, while also selecting fewer predictors than a sparse regularized logistic lasso with PC adjustment. Availability and implementation Our Julia package PenalizedGLMM.jl is publicly available on github: https://github.com/julstpierre/PenalizedGLMM. Supplementary information Supplementary data are available at Bioinformatics online.

Variable Selection for High-dimensional Generalized Linear Models using an Iterated Conditional Modes/Medians Algorithm

Variable Selection in Logistic Regression Model with Genetic Algorithm.

Bayesian variable selection in linear regression models with instrumental variables

Variable selection in Bayesian generalized linear-mixed models: an illustration using candidate gene case-control association studies

Variable Selection Using Nonlocal Priors in High-Dimensional Generalized Linear Models With Application to fMRI Data Analysis

Variable Selection in Generalized Estimating Equations Via Empirical Likelihood and Gaussian Pseudo-Likelihood

Variable selection for partially linear models via Bayesian subset modeling with diffusing prior

Efficient penalized generalized linear mixed models for variable selection and genetic risk prediction in high-dimensional data

Advances in Bayesian model selection consistency for high-dimensional generalized linear models

High-dimensional generalized median adaptive lasso with application to omics data

Variable Selection and Minimax Prediction in High-dimensional Functional Linear Model

Bayesian Variable Selection for Single Index Logistic Model

A robust and efficient estimation and variable selection method for partially linear models with large-dimensional covariates

Regularization Methods for High-Dimensional Instrumental Variables Regression With an Application to Genetical Genomics

Statistical Inference for High-Dimensional Generalized Linear Models With Binary Outcomes

Projection Predictive Inference for Generalized Linear and Additive Multilevel Models

Efficient computation of high-dimensional penalized generalized linear mixed models by latent factor modeling of the random effects

Scalable Bayesian bi-level variable selection in generalized linear models

Broken Adaptive Ridge Method for Variable Selection in Generalized Partly Linear Models with Application to the Coronary Artery Disease Data

High-dimensional variable selection for Cox's proportional hazards model

Variable selection method for the identification of epistatic models.