clusterProfiler (KEGG Enrichment Analysis) error
0
0
Entering edit mode
7.1 years ago
ta_awwad ▴ 340

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.1k views
ADD COMMENT
1
Entering edit mode

Maybe it is this issue?

ADD REPLY
0
Entering edit mode

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"
ADD REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 2703 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