GO ENRICHMENT ANALYSIS- DESEQ

I keep getting this Error message:

preparing geneSet collections... --> Expected input gene ID: C9orf78,CNOT9,SPIDR,FAM9B,SDC1,CCDC182 Error in check_gene_id(geneList, geneSets) : --> No gene can be mapped....

Here is my code, i attempting to perform gene set enrichement analysis:

GroupA <- results(dds, name = "Cancer_Type_healthy_vs_MESO") original_gene_list <- GroupA$log2FoldChange names(original_gene_list) <- GroupA$X gene_list<-na.omit(original_gene_list) gene_list = sort(gene_list, decreasing = TRUE)

gse <- gseGO(geneList=gene_list, ont ="ALL", keyType = "SYMBOL", nPerm = 10000, minGSSize = 3, maxGSSize = 800, pvalueCutoff = 0.05, verbose = TRUE, OrgDb = org.Hs.eg.db, pAdjustMethod = "none")

Please can someone tell me what is wrong with my code?

• 958 views • link

updated 2.9 years ago by Basti •

• written 2.9 years ago by odi •

You need to show an example of your list

• 0 views • link

written 2.9 years ago by Basti •

0 answers

No answers yet.

More posts like this

Errors in Functional Enrichment Analysis with Clusterprofiler
written 2.8 years ago by Tanjim •

library(clusterProfiler) library(org.Hs.eg.db) library(tidyverse) library(DOSE) library(ReactomePA) library(enrichplot) library(fgsea) library(data.table) library(ggplot2) keytypes(org.Hs.eg.db) res = read.csv("coex.Csv") head(res) original_gene_list = res$correlation names(original_gene_list) <- res$gene gene_list<-na.omit(original_gene_list) gene_list = sort(gene_list, decreasing …
How to Load a Molecular Signature Database into clusterProfiler gseGO?
written 3.1 years ago by Manuel Sokolov Ravasqueira •

I am using gseGO this way: ```r gse <- gseGO(geneList=gene_list, ont ="ALL", minGSSize = 3, maxGSSize = 800, pvalueCutoff = 0.05, verbose = TRUE, OrgDb …
could not find function "gseGO" in ClusterProfiler
written 3.1 years ago by Nitin •

I am trying to run the code using the example data set, ```r gse <- gseGO(geneList=gene_list, ont ="ALL", keyType = "ENSEMBL", nPerm = 10000, minGSSize …
go analysis result error
written 3.3 years ago by 1_GOld •

Hello everyone in biostars. I'm trying to do a go, kegg analysis. But I'm having an error. I'm using the following code. reference from [github](https://github.com/Bio-protocol/GO_KEGG_enrichment_analysis_for_rice_genes/blob/master/workflow/5_run_enrichKEGG.R#L21) …
How to sort graph based on fold change and instead of activated and suppressed, how to add up and d…
written 3.4 years ago by sharmatina189059 •

BiocManager::install("clusterProfiler", version = "3.16") BiocManager::install("pathview") BiocManager::install("enrichplot") library(clusterProfiler) library(enrichplot) # we use ggplot2 to add x axis labels (ex: ridgeplot) library(ggplot2) # reading in data from …
do.call with gseGO not working
written 3.7 years ago by Pahi •

Hi! I'm working on a highly configurable R script. I have made a function in witch there are several more functions, like gseGO from clusterProfiler. …
extracting gene names from clusterProfiler following gseGO
written 3.8 years ago by Elizabeth Alice •

Hello All, I have run the gseGO function from clusterProfiler. gse <- gseGO(geneList=gene_list, ont ="MF", keyType = "ENSEMBL", minGSSize = 10, maxGSSize = 800, pvalueCutoff …
No gene can be mapped
written 4.6 years ago by yang.zhizhang •

Hi When I do my data's Gene Set Enrichment Analysis with ClusterProfiler using codes of Mohammed Khalfan from website, when I run the following code …
genes underlying gseGO categories
written 5.4 years ago by LacquerHed •

Wondering how I can get a list of the genes I submitted that fall within each category/term represented on the eventual dotplot. Thanks. library(clusterProfiler) library(enrichplot) …
Error in GSEA function of clusterProfiler on R
written 6.1 years ago by j.thomas93 •

I am trying to run GSEA on RStudio using this: library(clusterProfiler) library(enrichplot) # SET THE DESIRED ORGANISM HERE organism = "org.Hs.eg.db" library(organism, character.only = TRUE) …