通路激活或者抑制
2020-03-11 本文已影响0人
一路向前_莫问前程_前程似锦
image.png
NES>0, 激活
NES <0, 抑制
facet_grid(~.sign)命令
data(geneList, package="DOSE")
gene <- names(geneList)[abs(geneList) > 2]
library(clusterProfiler)
library(clusterProfiler)
## 没有富集到任何数据
kk2 <- gseKEGG(geneList = geneList,
organism = 'hsa',
nPerm = 1000,
minGSSize = 120,
pvalueCutoff = 0.05,
verbose = FALSE)
head(kk2)
library(ggplot2)
dotplot(kk2,showCategory=12,split=".sign")+facet_grid(~.sign)
gseaKEGG_results <- kk2@result
