gmt = clusterProfiler::read.gmt(gmtfile = "C:\\文档\\liu\\生信\\GSEA资料\\MSigDB参考基因集\\c3.all.v6.2.symbols.gmt")
geneList = as.numeric(log2_foldchange.filter_sign)
names(geneList) = as.character(gene_name.filter_sign)
geneList.sort = sort(geneList, decreasing = T)
GSEA.res = GSEA(geneList = geneList.sort,TERM2GENE = gmt,pvalueCutoff = 0.05,maxGSSize = 1000, minGSSize = 10)
gseaplot(GSEA.res,geneSetID = 1,na.rm=TRUE)
一直到sort都没有问题:
之后的GSEA:
求大神~~
https://www.csdn.net/gather_24/Mtjagg1sNzIzMjgtYmxvZwO0O0OO0O0O.html
同求助……楼上的链接已经挂了
同求助大神
初始数据整理太重要啦!!!