Ridge estimation of network models from time-course omics data

Abstract

Time-course omics experiments enable the reconstruction of the dynamics of the cellular regulatory network. Here, we describe the means for this reconstruction and the downstream exploitation of the inferred network. It is assumed that one of the various vector-autoregressive models (VAR) models presented here serves as a reasonably accurate description of the time-course omics data. The models are estimated through ridge penalized likelihood maximization, accompanied by functionality for the determination of optimal penalty paramaters. Prior knowledge on the network topology is accommodated by the estimation procedures. Various routes that translate the fitted models into more tangible implications for the medical researcher are described. The network is inferred from the-nonsparse-ridge estimates through empirical Bayes probabilistic thresholding. The influence of a (trait of a) molecular entity at the current time on those at future time points is assessed by mutual information, impulse response analysis, and path decomposition of the covariance. The presented methodology is applied to the omics data from the p53 signaling pathway during HPV-induced cellular transformation. All methodology is implemented in the ragt2ridges package, freely available from the Comprehensive R Archive Network.

Keywords: cervical cancer; constrained estimation; maximum likelihood; mutual information; path analysis; penalized estimation; time series analysis; vector autoregressive process.