clusterProfiler won’t read gene list

clusterProfiler won’t read gene list

0

So I have a list of DE genes that I would like to analyse for enriched GO and KEGG terms. I was going to use clusterProfiler for this, but I can’t seem to get past constructing the gene list. I have followed the vignette and I believe that I have constructed my gene list correctly as a sorted named vector of Entrez IDs, but it still doesn’t work. Below is the output of the sample gene list in the vignette and my list. There appears to be no difference between my named vector and the example, yet the example works and mine does not.

> str(geneList) # example gene list in vignette
 Named num [1:12495] 4.57 4.51 4.42 4.14 3.88 ...
 - attr(*, "names")= chr [1:12495] "4312" "8318" "10874" "55143" ...

> str(gene_list) # my gene list that I have constructed
Named num [1:18907] 4.64 4.25 3.89 3.64 3.5 ...
 - attr(*, "names")= chr [1:18907] "13048" "100043108" "13099" "23850" ...

I get the following error when I run

de <- names(geneList)[abs(gene_list) > 2]
edo <- enrichDGN(de)

--> No gene can be mapped....
--> Expected input gene ID: 229,151393,5076,2006,3622,8302
--> return NULL...


gsea


software error


clusterProfiler

• 454 views

Read more here: Source link