2019-08-26 biomaRt 進行基因轉(zhuǎn)換

library("biomaRt")

human = useMart("ensembl", dataset = "hsapiens_gene_ensembl")

mouse = useMart("ensembl", dataset = "mmusculus_gene_ensembl")

gene<- read.table("geneinrepeat1000kb.txt")

? gene<-gene[,1]

? aa<-getBM(attributes = c("refseq_mrna","external_gene_name"),filters = "refseq_mrna",values = gene,mart = mouse)

? write.table(aa,file="geneinreat1000kb-genename.txt")

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容