This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Too many fails converting gene SYMBOL to ENTREZID

Hello, my dataset is GSE37133. I use GPL information to annotate gene SYMBOL(No bioconductor package for this chip). Main code:

gpl <- getGEO("GPL2996", destdir=".")
# dataframe that can map probes id and gene SYMBOL
annot <- Table(gpl)[,c(1, 2, 5, 10)]
# expression data annotation
expr_2 <- merge(x = annot, y = diff_expr_tenmg, by = "ID")

I want to call enrichKEGG function in clusterProfiler which do not support gene SYMBOL. So I need to run bitr to convert gene SYMBOL to ENTREZID. Main code

# load genowide annotation database
library(org.Rn.eg.db)
library(clusterProfiler)

# convertion
gene_list <- bitr(expr_2$Gene_Symbol, fromType="SYMBOL", toType="ENTREZID", OrgDb=org.Rn.eg.db, drop=TRUE)

Here, about 30% gene SYMBOL fails at converting. Where am I wrong? Thanks.

microarray clusterprofiler geo

Your code is not reproducible because diff_expr_tenmg was never instantiated (?). Have you triied biomaRt to convert the gene symbols?

Sorry I didn't paste all my code here. I will try biomaRt later. Thanks Kevin.

0 answers

No answers yet.

Log in to answer this question.