Computationally scalable regression modeling for ultrahigh-dimensional omics data with ParProx

Abstract

Statistical analysis of ultrahigh-dimensional omics scale data has long depended on univariate hypothesis testing. With growing data features and samples, the obvious next step is to establish multivariable association analysis as a routine method for understanding genotype-phenotype associations. Here we present ParProx, a state-of-the-art implementation to optimize overlapping group lasso regression models for time-to-event and classification analysis, guided by biological priors through coordinated variable selection. ParProx not only enables model fitting for ultrahigh-dimensional data within the architecture for parallel or distributed computing, but also allows users to obtain interpretable regression models consistent with known biological relationships among the independent variables, a feature long neglected in statistical modeling of omics data. We demonstrate ParProx using three different omics data sets of moderate to large numbers of variables, where we use genomic regions and pathways to arrive at sparse regression models comprised of biologically related independent variables. ParProx is naturally applicable to a wide range of studies using ultrahigh-dimensional omics data, ranging from genome-wide association analysis to single cell sequencing studies where multivariable modeling is computationally intractable.