Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2014 Jun;101(2):253-268.
doi: 10.1093/biomet/asu009.

Direct estimation of differential networks

Affiliations

Direct estimation of differential networks

Sihai Dave Zhao et al. Biometrika. 2014 Jun.

Abstract

It is often of interest to understand how the structure of a genetic network differs between two conditions. In this paper, each condition-specific network is modeled using the precision matrix of a multivariate normal random vector, and a method is proposed to directly estimate the difference of the precision matrices. In contrast to other approaches, such as separate or joint estimation of the individual matrices, direct estimation does not require those matrices to be sparse, and thus can allow the individual networks to contain hub nodes. Under the assumption that the true differential network is sparse, the direct estimator is shown to be consistent in support recovery and estimation. It is also shown to outperform existing methods in simulations, and its properties are illustrated on gene expression data from late-stage ovarian cancer patients.

Keywords: Differential network; Graphical model; High dimensionality; Precision matrix.

PubMed Disclaimer

Figures

Figure 1
Figure 1
Receiver operating characteristic curves for support recovery of Δ0=ΣY1ΣX1; solid line: thresholded direct estimator; dashed line: thresholded fused graphical lasso estimator with λ1 = 0; dotted line: fused graphical lasso estimator with λ1 = 0 · 1
Figure 2
Figure 2
Estimates of the differential networks between ovarian cancer subtypes. The direct and fused graphical lasso estimators were thresholded and the separate estimator was further sparsified; see text. Black edges show increase in conditional dependency from ovarian cancer subtype C1 to subtypes C2–C6, gray edges show decrease. (a)–(c): KEGG 04350, TGF-β pathway, (d)–(f): KEGG 04210, Apoptosis pathway. (a)–(e) tuned using L with (5), (f) tuned with λ1 = 0 and λ2 to give the same number of edges as (d); see text. Separate estimator not shown for apoptosis pathway.

References

    1. Bandyopadhyay S, Mehta M, Kuo D, Sung M-K, Chuang R, Jaehnig EJ, Bodenmiller B, Licon K, Copeland W, Shales M, Fiedler D, Dutkowski J, Guénolé A, van Attikum H, Shokat KM, Kolodner RD, Huh W-K, Aebersold R, Keogh M-C, Krogan NJ, Ideker T. Rewiring of genetic networks in response to DNA damage. Sci. Signal. 2010;330(6009):1385–1389. - PMC - PubMed
    1. Barabási A-L, Oltvai ZN. Network biology: understanding the cell’s functional organization. Nat. Rev. Genet. 2004;5(2):101–113. - PubMed
    1. Barabási A-L, Gulbahce N, Loscalzo J. Network medicine: a network-based approach to human disease. Nat. Rev. Genet. 2011;12(1):56–68. - PMC - PubMed
    1. Bellail AC, Qi L, Mulligan P, Chhabra V, Hao C. TRAIL agonists on clinical trials for cancer therapy: the promises and the challenges. Rev. Recent Clin. Trials. 2009;4(1):34–41. - PubMed
    1. Boyd S, Parikh N, Chu E, Peleato B, Eckstein J. Distributed optimization and statistical learning via the alternating direction method of multipliers. Foundat. Trends Mach. Learn. 2011;3(1):1–122.