HighDimMixedModels.jl: Robust High Dimensional Mixed Models across Omics Data

High dimensional mixed-effect models are an increasingly important form of regression in modern biology, in which the number of variables often matches or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate gradient descent (CGD) algorithm that lacks guarantees of convergence to a global optimum. Here, we study empirically the behavior of the algorithm across a number of common study types in modern omics datatypes. In particular, we study the empirical performance of high dimensional mixed-effect models fit to data simulated to mimic the features of transcriptome, genome-wide association, and microbiome data. In addition, we study the performance of the model on real data from each of these study types. To facilitate these simulations, we implement the algorithm in an open source Julia package HighDimMixedModels.jl. We compare the performance of two commonly used penalties, namely LASSO and SCAD, within the HighDimMixedModels.jl framework. Our results demonstrate that the SCAD penalty consistently outperforms LASSO in terms of both variable selection and estimation accuracy across omics data. Through our comprehensive analysis, we illuminate the intricate relationship between algorithmic behavior, penalty selection, and dataset properties such as the correlation structure among features, providing valuable insights for researchers employing high dimensional mixed-effect models in biological investigations. Author SummaryHigh dimensional mixed-effect models are increasingly indispensable in modern biology, particularly in omics studies, where the number of variables often equals or surpasses the number of samples, and data are collected in clusters or groups. In our research, we concentrate on the penalized likelihood approach to fitting these models, employing a coordinate gradient descent (CGD) algorithm. While CGD is a widely used optimization technique, its convergence to a global optimum lack guarantees, prompting our empirical investigation into its behavior across various study types common in modern omics datasets. Our study provides insights into the performance of high dimensional mixed-effect models fitted to data simulating transcriptome, genome-wide association, and microbiome datasets. Additionally, we evaluate the models performance on real datasets from each of these study types. To facilitate reproducibility and further research, we have implemented the algorithm in an open-source Julia package, HighDimMixedModels.jl. Notably, HighDimMixedModels.jl stands out as the first package capable of seamlessly handling various omics datasets without errors, offering a user-friendly solution for researchers across disciplines. While numerous software packages are available for implementing high dimensional mixed-effects models on omics data, there is currently no comprehensive review source summarizing all methods. We provide a table summarizing existing methods, available in the Supplementary Material.