基因ID转换---R语言

这里使用的是Y叔的R包clustProfilter，里面有个函数bitr()

library(org.Hs.eg.db)
library("clusterProfiler")
> keytypes(org.Hs.eg.db)
 [1] "ACCNUM"       "ALIAS"        "ENSEMBL"      "ENSEMBLPROT"  "ENSEMBLTRANS" "ENTREZID"    
 [7] "ENZYME"       "EVIDENCE"     "EVIDENCEALL"  "GENENAME"     "GO"           "GOALL"       
[13] "IPI"          "MAP"          "OMIM"         "ONTOLOGY"     "ONTOLOGYALL"  "PATH"        
[19] "PFAM"         "PMID"         "PROSITE"      "REFSEQ"       "SYMBOL"       "UCSCKG"      
[25] "UNIGENE"      "UNIPROT"
#这里需要确认你的输入的gene ID是什么类型的
> gene.df <- bitr("TP53", fromType = "ALIAS" , 
+                 toType = c("ENSEMBL", "SYMBOL"), 
+                 OrgDb = org.Hs.eg.db)
'select()' returned 1:1 mapping between keys and columns
> gene.df
  ALIAS         ENSEMBL SYMBOL
1  TP53 ENSG00000141510   TP53

但是使用这种方法总是会有一些基因比对不上，就会有类似的warning

warning

如果有更好的方法，欢迎大家一起探讨！