Entering edit mode
2.9 years ago
bookorg
▴
20
Hello, I have analysed three dataset and found 3 common DEG among them.But when i go for the GSEA analysis and try to draw the plot ,there I faced a trouble .When i run it,in return i get the bellow messages in R.I am very much confused about that.
preparing geneSet collections... GSEA analysis... no term enriched under specific pvalueCutoff... Warning message: In preparePathwaysAndStats(pathways, stats, minSize, maxSize, gseaParam, : All values in the stats vector are greater than zero and scoreType is "std", maybe you should switch to scoreType = "pos".
Here is my code:
library(readr)
sign_DEG_ <- read_csv("For GSEA,ORA,Pathway analysis/sign_DEG.csv")
View(sign_DEG)
library(dplyr)
library(stringr)
library(ggplot2)
library(ggpubr)
library(clusterProfiler)
library(ReactomePA)
library(org.Hs.eg.db)
library(EnsDb.Hsapiens.v86)
library(enrichplot)
geneid <- sign_DEG$SYMBOL
genes <- ensembldb::select(org.Hs.eg.db,
keys=geneid,
columns=c("ENTREZID"),
keytype="SYMBOL")
genes <- na.omit(genes)
sign_DEG <- merge(genes, sign_DEG, by = 'SYMBOL')
write.csv(sign_DEG_,"./Sign DEG after omit.csv", row.names = TRUE)
sign_DEGLogFC <- sign_DEG$logFC
names(sign_DEGLogFC) <-sign_DEG$ENTREZID
sign_DEG <- sort(sign_DEG, decreasing = TRUE)
bpSars <- gseGO(sign_DEGLogFC,
ont = 'BP',
org.Hs.eg.db)
bpSars <- setReadable(bpSars, org.Hs.eg.db)
bpSars@result$Description <- str_to_title(str_wrap(bpSars@result$Description, 50))
write.csv(bpSars, './Result/gsea-bp-sars.csv')
dotplot(bpSars, showCategory=10, split=".sign") + facet_grid(.~.sign)
ggsave('./Result/gsea-bp-sarscov.tiff', height = 8, width = 8)
Can you include your code?
it means there are no significantly differentially expressed genes (nothing statistically significant). Most immediately, I'd do two things:
1) set my p-value cut-off very liberally to see if this fixes the problem
2) set the scoreType argument to pos as they are suggesting.
If you still get errors after that, then maybe post those
according to your two points i got the 1st one but I am confused about the second point.Can you plz explain the second point in deeply plz.