This is a test version of Biostars. For the public version, visit https://www.biostars.org.
clusterProfiler (KEGG Enrichment Analysis) error

Hello everybody, I am currently running RNA_seq pipeline and it was absolutely fine until I reached the KEGG step. first, I have 400 DE genes . I created matrix of ENTREZID and log2foldchange. here is the head of my matrix:

> head(cp_matrix)
     35526      41092      39976      40831      36778      33636 
-0.4144379 -0.6321093 -0.3219794 -2.0373447  2.3645869  0.6998755

then I have run:

kegg_enrich <- enrichKEGG(gene = names(cp_matrix),
                          organism = "fly",
                          pvalueCutoff = 0.05, 
                          qvalueCutoff = 0.10)

and I have the same error message:

No gene can be mapped....
--> return NULL...

any idea how to solve this??

thanks much

rna-seq chip-seq kegg pathway
• 4,887 views • link
written by ta_awwad ta_awwad

Maybe it is this issue?

• 0 views • link
written by h.mon h.mon

My IDs are already character

class (names(cp_matrix))
[1] "character"
names(cp_matrix)
  [1] "36544"    "39244"    "36420"    "31904"    "34189"    "44149"    "35526"    "41092"   
  [9] "31130"    "32471"    "37815"    "43183"    "45398"    "40850"    "53425"    "40831"
• 0 views • link
written by ta_awwad ta_awwad

Did you convert genes IDs according to ENTREZID using clusterProfiler::bitr?

• 0 views • link
written by H.Hasani H.Hasani

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • clusterProfiler Results not showing
    written by geoffkerr22 • geoffkerr22
    Hello, I am having some issues with clusterProfiler and the enrichKEGG function. Here is my code: library(clusterProfiler) kegg_enrich <- enrichKEGG(gene = results_sig_entrez$entrez, organism = 'hsa', …