I wanted to perform the kegg enrichment method. Here's what I did: ran my fasta to kaas to get the KEGG Orthology IDs (K#####) then used these ID to perform enrichKEGG.
- An example of the command was performed to test the function.
keggID <- c("K15179", "K03377") x <- enrichKEGG(keggID, organism='ko', keyType='kegg')
--> No gene can be mapped.... --> Expected input gene ID: K10780,K16305,K08093,K14275,K18819,K22554 --> return NULL...
When I use other KO ID like K05164, the program ran without any errors. This could mean that the KEGG IDs (in my example above) could not be mapped because they are not enriched? Also for the IDs that were accepted (e.g. K05164), it didn't have qvalue. I wonder if what I am doing wrong.
- My other question is that, are the Kegg orthologous IDs invalid for bitr_kegg() function? I could not convert my K##### IDs. This shows up everytime I try different IDs (except for the example, ofc):
Warning message: In bitr_kegg("K15179", fromType = "kegg", toType = "uniprot", organism = "hsa") : 100% of input gene IDs are fail to map...
This is also true for ncbi protein and geneID. Am I using the wrong organism?
- Also, what is the proper way to perform KEGG enrichment for non-model organism? I am really confused as to whether I would just use the closest model organism with available db to perform KEGG or use all available databases?
Thank you for the help!