GSEA富集分析图

绘制多条通路的富集图，如下

1.输入文件：差异基因富集表

（1）diff.symbol.csv

（2）pathway_TERM2SYMBOL.txt

（3）pathway_TERM2NAME.txt

2. 实现代码：

library(clusterProfiler)
library(enrichplot)
data <- read.csv(file="diff.SYMBOL.csv", header=T, sep=",", check.name=F)
names(data)[1] <- "gene"
names(data)[2] <- "logFC"
genelist <- data$logFC
names(genelist) <- data$gene
genelist <- sort(genelist, decreasing=T)
genelist <- genelist[genelist != 0]

##gsea富集
set.seed(1234)
res <- GSEA(
    genelist,
    nPerm = 1000,
    pvalueCutoff = 1,
    pAdjustMethod = "BH",
    TERM2GENE = term2gene,   ##两列，第一列keggID,第二列gene name
    TERM2NAME = term2name,   ##两列， 第一列KEGGID，第二列通路描述信息
    seed = T
)

##作图
pngout <- "gsea.pathway.png"
p <- gseaplot2(res, 1:10)
png(pngout, width=1780, height=1280, res=200)
print(p)
dev.off()

3.图片结果

见顶部！！！！