Regularized estimation of large-scale gene association networks using graphical Gaussian models

Abstract

Graphical Gaussian models are popular tools for the estimation of (undirected) gene association networks from microarray data. A key issue when the number of variables greatly exceeds the number of samples is the estimation of the matrix of partial correlations. Since the (Moore-Penrose) inverse of the sample covariance matrix leads to poor estimates in this scenario, standard methods are inappropriate and adequate regularization techniques are needed. In this article, we investigate a general framework for combining regularized regression methods with the estimation of Graphical Gaussian models. This framework includes various existing methods as well as two new approaches based on ridge regression and adaptive lasso, respectively. These methods are extensively compared both qualitatively and quantitatively within a simulation study and through an application to six diverse real data sets. In addition, all proposed algorithms are implemented in the R package "parcor", available from the R repository CRAN.

added additional experiments

References in corpus (5)

• High-dimensional graphs and variable selection with the Lasso
• Sparse permutation invariant covariance estimation
• Adaptive Lasso for High Dimensional Regression and Gaussian Graphical Modeling
• Sparse inverse covariance estimation with the lasso
• Structure Variability in Bayesian Networks

Cited by in corpus (19)

• A Tutorial on Regularized Partial Correlation Networks
• The Gaussian Graphical Model in Cross-sectional and Time-series Data
• Estimating Psychological Networks and their Accuracy: A Tutorial Paper
• Brief Report on Estimating Regularized Gaussian Networks from Continuous and Ordinal Data
• Efficient Learning of Optimal Markov Network Topology with k-Tree Modeling
• Bayesian Network Structure Learning with Permutation Tests
• Conic Optimization Theory: Convexification Techniques and Numerical Algorithms
• Graphical Lasso and Thresholding: Equivalence and Closed-form Solutions
• On generating random Gaussian graphical models
• Graph quilting: graphical model selection from partially observed covariances
• A partial orthogonalization method for simulating covariance and concentration graph matrices
• Comparisons of penalized least squares methods by simulations
• Hub discovery in partial correlation graphical models
• On the Computation and Applications of Large Dense Partial Correlation Networks
• Asymptotic control of FWER under Gaussian assumption: application to correlation tests
• Graphical models for circular variables
• Incorporating prior information and borrowing information in high-dimensional sparse regression using the horseshoe and variational Bayes
• On the use of cross-validation for the calibration of the adaptive lasso
• A Log-Linear Graphical Model for Inferring Genetic Networks from High-Throughput Sequencing Data