基因Id转化

library(biomaRt)
ensembl = useMart(biomart = "ENSEMBL_MART_ENSEMBL",
                  dataset="mmusculus_gene_ensembl", 
                  host = "www.ensembl.org")
values = c("NM_144958", "NM_144959", "NM_144960")
getBM(c("refseq_mrna", "external_gene_name"), "refseq_mrna", gset_expr_annot$symbol,mart=ensembl)

keytypes(org.Hs.eg.db)
values = c("NM_144958", "NM_144959", "NM_144960")
bitr(c("PTGS1"), fromType = "REFSEQ",
     toType = c("ENTREZID","SYMBOL","ALIAS","REFSEQ"),
     OrgDb = org.Hs.eg.db)

参考

https://stackoverflow.com/questions/34612838/convert-genbank-nucleotide-accession-to-gene-symbol