upsetplot大家见得多,首先来个富集分析的实例:

library(DOSE)
data(geneList)
de <- names(geneList)[abs(geneList) > 2]
edo <- enrichDGN(de)

library(enrichplot)
upsetplot(edo)

在v >= 1.5.2版本的enrichplot中,我加入了对GSEA结果的支持,于是你现在也可以用upsetplot来画,但风格不一样哦。

再来个实例:

library(clusterProfiler)
kk2 <- gseKEGG(geneList     = geneList,
               organism     = 'hsa',
               nPerm        = 1000,
               minGSSize    = 120,
               pvalueCutoff = 0.05,
               verbose      = FALSE)
upsetplot(kk2) 

Gene set variation analysis (GSVA). Pathway analyses were predominantly performed on the 50 hallmark pathways described in the molecular signature database 16 , exported using the GSEABase package (version 1.36.0). We also assessed metabolic pathway activities using a described curated dataset 48 . To reduce pathway overlaps and pathway redundancies, each gene set associated with a pathway was trimmed to only contain unique genes, and all genes associated to two or more pathways were removed. Most gene sets retained >70% of their associated genes. Next, to assign pathway activity estimates to individual cells, we applied GSVA 49 using standard settings, as implemented in the GSVA package (version 1.22.4).

针对上面这篇文章所采用的策略,似乎我们现在所呈现的图,更好,也能给出更多的信息,你可能还可以发现一些有趣的cross-talk。PS:画图的集合,也是可选的。