Entering edit mode
5.4 years ago
MatthewP
★
1.4k
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.
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.