Lightweight Iterative Gene set Enrichment in R
R C++
Switch branches/tags
Nothing to show
Clone or download
Latest commit 0fdf27e Jun 27, 2018
Permalink
Failed to load latest commit information.
R modify ignore files Jun 9, 2018
data Make examples. Clean up code Sep 27, 2015
man modify ignore files Jun 9, 2018
src use different rand generator Jun 10, 2018
tools make images smaller; move to tools folder Jun 10, 2018
vignettes update docs Jun 27, 2018
.Rbuildignore make images smaller; move to tools folder Jun 10, 2018
.gitignore make images smaller; move to tools folder Jun 10, 2018
.travis.yml Add travis testing and codecov Sep 27, 2015
DESCRIPTION fix description based on feedback from cran Jun 12, 2018
LICENSE prepare for cran Jun 8, 2018
NAMESPACE modify ignore files Jun 9, 2018
NEWS.md
README.md update docs Jun 27, 2018
cran_build_instructions.md modify ignore files Jun 9, 2018

README.md

Lightweight Iterative Gene set Enrichment in R (LIGER)

Build Status DOI

Gene Set Enrichment Analysis (GSEA) is a computational method that determines whether an a priori defined set of genes shows statistically significant, concordant differences between two biological states. The original algorithm is detailed in Subramanian, Tamayo, et al. with Java implementations available through the Broad Institute.

The liger package provides a lightweight R implementation of this enrichment test on a list of values. Given a list of values, such as p-values or log-fold changes derived from differential expression analysis or other analyses comparing biological states, this package enables you to test a priori defined set of genes for enrichment to enable interpretability of highly significant or high fold-change genes.

Tutorials

Sample plots

Testing individual gene sets

> gsea(values=vals, geneset=gs, mc.cores=1, plot=TRUE)

Testing multiple gene sets

> ptm <- proc.time()
> bulk.gsea(vals, org.Hs.GO2Symbol.list[1:10])
                  p.val     q.val     sscore       edge
  GO:0000002 0.00009999 0.0000000  2.6139724  70.912194
  GO:0000003 0.25647435 0.4074000  0.5400972  13.170093
  GO:0000012 0.24347565 0.4074000  0.5550191   8.392397
  GO:0000014 0.13428657 0.3779000  0.6906745  -4.458762
  GO:0000018 0.14448555 0.7215200 -0.6498716  11.111976
  GO:0000022 0.33446655 0.8023667 -0.5014560 -11.015244
> proc.time() - ptm
   user  system elapsed
  6.200   0.065   6.414
> ptm <- proc.time()
> iterative.bulk.gsea(values=vals, set.list=org.Hs.GO2Symbol.list[1:10])
                  p.val   q.val     sscore       edge
  GO:0000002 0.00009999 0.00000  2.6054741  70.912194
  GO:0000003 0.25867413 0.41730  0.5352656  13.170093
  GO:0000012 0.24727527 0.41730  0.5545500   8.392397
  GO:0000014 0.14128587 0.40300  0.6903051  -4.458762
  GO:0000018 0.13988601 0.72684 -0.6485933  11.111976
  GO:0000022 0.31666833 0.81440 -0.4975180 -11.015244
> proc.time() - ptm
    user  system elapsed
  11.250   0.107   5.320

Install

require(devtools)
devtools::install_github("JEFworks/liger")

Contributing

We welcome any bug reports, enhancement requests, and other contributions. To submit a bug report or enhancement request, please use the liger GitHub issues tracker. For more substantial contributions, please fork this repo, push your changes to your fork, and submit a pull request with a good commit message.