Hi all, I'm using enricher to do GO and KEGG enrichment, I have seen that https://guangchuangyu.github.io/2015/05/use-clusterprofiler-as-an-universal-enrichment-analysis-tool/
My annotation file is hear anno <- read.delim("/path/GO_anno.txt") gene <- read.delim("/path/id.txt",header = T) Query KEGG description TRINITY_DN3944_c0_g1_i8.p1 ko:K00001,ko:K18857 alcohol dehydrogenase TRINITY_DN40699_c0_g1_i1.p1 ko:K00002,ko:K00011 Aldo-keto reductase family 4 member TRINITY_DN28700_c0_g1_i1.p1 ko:K00002,ko:K22374 Aldo-keto reductase family 4 member
I use this code to create TERM2GENE and TERM2NAME: TERM2GENE <- anno[,c('GO', 'Query')] TERM2NAME <- anno[,c('GO', 'description')]
When I run x = enricher(gene, TERM2GENE=TERM2GENE, TERM2NAME=TERM2NAME), I get this ouput: --> No gene can be mapped.... --> Expected input gene ID: TRINITY_DN1095_c0_g5_i2,...... --> return NULL...
What should I do? Thanks in advance.