Skip to content

Instantly share code, notes, and snippets.

Embed
What would you like to do?
## devtools::install_github("stephenturner/msigdf")
library(msigdf)
library(dplyr)
library(clusterProfiler)
c2 <- msigdf.human %>%
filter(collection == "c2") %>% select(geneset, entrez) %>% as.data.frame
data(geneList)
de <- names(geneList)[1:100]
x <- enricher(de, TERM2GENE = c2)
y <- GSEA(geneList, TERM2GENE = c2)
head(x)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
You can’t perform that action at this time.