Only positive enrichment score in GSEA
0
0
Entering edit mode
8 months ago
enee ▴ 20

Hi everyone, I am a novice in the enrichment analyses and I have a problem with the gseGO function of the clusterprofiler package in R. Basically I'm analyzing mouse expression data and when I go to run the analysis with my code there seem to be only enriched pathways (with positive NES). I don't understand why if I use the human annotation org.Hs.eg.db my results and graphs come out as I expect them to. The first plot comes from the following code with mouse annotation (org.Mm.eg.db) and the second plot comes from the code with the human annotation.

GO_gsea <- function(df) { library(clusterProfiler) library(org.Mm.eg.db) geneList = df[,2] names(geneList) = as.character(df[,7]) geneList = sort(geneList, decreasing = TRUE) data_GO_GSEA <- gseGO(geneList = geneList, OrgDb = org.Mm.eg.db, ont = "BP", minGSSize = 3, maxGSSize = 800, pvalueCutoff = 0.05, verbose = TRUE, keyType = "SYMBOL", pAdjustMethod = "BH", eps = 0) }
CeO2vsDiff1_gsea <- GO_gsea(CeO2vsDiff1) 
dotplot(CeO2vsDiff1_gsea, showCategory=10, split=".sign") + facet_grid(.~.sign)

enter image description here enter image description here

enrichmentanalysis clusterprofiler gsea bioconductor rnaseq • 354 views
ADD COMMENT

Login before adding your answer.

Traffic: 1865 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6