MyNixOS website logo
Description

Collaborative Graphical Lasso - Multi-Omics Network Reconstruction.

Reconstruct networks from multi-omics data sets with the collaborative graphical lasso (coglasso) algorithm described in Albanese, A., Kohlen, W., and Behrouzi, P. (2024) <arXiv:2403.18602>. Build multiple networks using the coglasso() function, select the best one with stars_coglasso().

coglasso - Collaborative Graphical Lasso

R-CMD-check

codecov

Coglasso implements collaborative graphical lasso, an algorithm for network reconstruction from multi-omics data sets (Albanese, Kohlen and Behrouzi, 2024). Our algorithm joins the principles of the graphical lasso by Friedman, Hastie and Tibshirani (2008) and collaborative regression by Gross and Tibshirani (2015).

Installing from CRAN (not released on CRAN yet)

You will be able to install the CRAN release of coglasso with:

install.packages("coglasso")

Installing the development version

To install the development version of coglasso from GitHub you need to make sure to install devtools with:

if (!require("devtools")) {
  install.packages("devtools")
}

You can then install the development version with:

devtools::install_github("DrQuestion/coglasso")

Usage

Here follows an example on how to reconstruct a multi-omics network with collaborative graphical lasso. For a more exhaustive example we refer the user to the vignette vignette("coglasso"). The package provides example multi-omics data sets of different dimensions, here we will use multi_omics_sd_small. Please notice that the current version of the coglasso package expects multi-omics data sets with two “omic” layers, where the single layers are grouped by column. For example, in multi_omics_sd_small the first 14 columns represent transcript abundances, and the other 5 columns represent metabolite abundances. To default usage of coglasso() only needs the input dataset and the dimension of the first “omic” layer.

library(coglasso)

cg <- coglasso(multi_omics_sd_small, pX = 14)

coglasso() explores several combinations of the hyperparameters characterizing collaborative graphical lasso. To select the combination yielding the most stable, yet sparse network, the package provides the function stars_coglasso(). This function implements a coglasso-adapted version of the StARS selection algorithm (Liu, Roeder and Wasserman, 2010).

sel_cg <- stars_coglasso(cg)

References

Albanese, A., Kohlen, W., & Behrouzi, P. (2024). Collaborative graphical lasso (arXiv:2403.18602). arXivhttps://doi.org/10.48550/arXiv.2403.18602

Friedman, J., Hastie, T., & Tibshirani, R. (2008). Sparse inverse covariance estimation with the graphical lasso. Biostatistics, 9(3), 432–441. https://doi.org/10.1093/biostatistics/kxm045

Gross, S. M., & Tibshirani, R. (2015). Collaborative regression. Biostatistics, 16(2), 326–338. https://doi.org/10.1093/biostatistics/kxu047

Liu, H., Roeder, K., & Wasserman, L. (2010). Stability Approach to Regularization Selection (StARS) for High Dimensional Graphical Models (arXiv:1006.3316). arXivhttps://doi.org/10.48550/arXiv.1006.3316

Metadata

Version

1.0.2

License

Unknown

Platforms (77)

    Darwin
    FreeBSD
    Genode
    GHCJS
    Linux
    MMIXware
    NetBSD
    none
    OpenBSD
    Redox
    Solaris
    WASI
    Windows
Show all
  • aarch64-darwin
  • aarch64-freebsd
  • aarch64-genode
  • aarch64-linux
  • aarch64-netbsd
  • aarch64-none
  • aarch64-windows
  • aarch64_be-none
  • arm-none
  • armv5tel-linux
  • armv6l-linux
  • armv6l-netbsd
  • armv6l-none
  • armv7a-darwin
  • armv7a-linux
  • armv7a-netbsd
  • armv7l-linux
  • armv7l-netbsd
  • avr-none
  • i686-cygwin
  • i686-darwin
  • i686-freebsd
  • i686-genode
  • i686-linux
  • i686-netbsd
  • i686-none
  • i686-openbsd
  • i686-windows
  • javascript-ghcjs
  • loongarch64-linux
  • m68k-linux
  • m68k-netbsd
  • m68k-none
  • microblaze-linux
  • microblaze-none
  • microblazeel-linux
  • microblazeel-none
  • mips-linux
  • mips-none
  • mips64-linux
  • mips64-none
  • mips64el-linux
  • mipsel-linux
  • mipsel-netbsd
  • mmix-mmixware
  • msp430-none
  • or1k-none
  • powerpc-netbsd
  • powerpc-none
  • powerpc64-linux
  • powerpc64le-linux
  • powerpcle-none
  • riscv32-linux
  • riscv32-netbsd
  • riscv32-none
  • riscv64-linux
  • riscv64-netbsd
  • riscv64-none
  • rx-none
  • s390-linux
  • s390-none
  • s390x-linux
  • s390x-none
  • vc4-none
  • wasm32-wasi
  • wasm64-wasi
  • x86_64-cygwin
  • x86_64-darwin
  • x86_64-freebsd
  • x86_64-genode
  • x86_64-linux
  • x86_64-netbsd
  • x86_64-none
  • x86_64-openbsd
  • x86_64-redox
  • x86_64-solaris
  • x86_64-windows