【塊】Gene id格式轉(zhuǎn)換的多種方法

本文將逐步增加各種方法星爪,有時(shí)間再研究各個(gè)方法的區(qū)別。
每個(gè)人在轉(zhuǎn)換id格式時(shí)曾掂,可能會(huì)遇到很多與以下流程不同的情況惫谤。比如說gene id有小版本、去重后數(shù)據(jù)大量丟失等等珠洗,建議具體問題再去查詢解決方案溜歪,也可以把問題私信給我,我學(xué)習(xí)后整理在此文最后许蓖。

假設(shè)文件名文“genes_exp”蝴猪,其中有一列為需要轉(zhuǎn)換的“geneid”

1、使用gtf文件注釋

可在以下兩個(gè)鏈接下載gtf文件(分別是GENCODE和Ensembl):
https://www.gencodegenes.org/human/#
http://ftp.ensembl.org/pub/release-104/gtf/homo_sapiens/

library(rtracklayer)
gff<-readGFF(“”)#自己修改路徑

#提取所需要的部分注釋膊爪,假設(shè)是從“gene_id”轉(zhuǎn)換到”gene_name”
anno<-gff[gff$type == “gene”, c(“gene_id”,”gene_name”)]

#根據(jù)“gene_id”匹配轉(zhuǎn)換genes_exp與anno自阱,自己修改by.y參數(shù)
result<-merge(anno,genes_exp, by.x= “gene_id”, by.y= “geneid”)

#去重與檢查空值與NA
result<-result[!duplicated(result$gene_name),]
table(result$gene_name == “”)
table(is.na(result$gene_name))

2、使用bitr(clusterProfiler包)

有些人的數(shù)據(jù)有小數(shù)點(diǎn)米酬,可以用str_sub處理一下

library(clusterProfiler)
library(org.Hs.eg.db)

# 查看有哪些gene id格式
keytypes(org.Hs.eg.db)

#適用bitr()函數(shù)進(jìn)行轉(zhuǎn)換动壤,可同時(shí)轉(zhuǎn)換成多種格式
result$gene_name<-bitr(genes_exp$geneid, fromType= “ENSEMBL”, toType= c(“SYMBOL”)

#接下來去重,檢查空值與NA
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末淮逻,一起剝皮案震驚了整個(gè)濱河市琼懊,隨后出現(xiàn)的幾起案子阁簸,更是在濱河造成了極大的恐慌,老刑警劉巖哼丈,帶你破解...
    沈念sama閱讀 206,839評(píng)論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件启妹,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡醉旦,警方通過查閱死者的電腦和手機(jī)饶米,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來车胡,“玉大人檬输,你說我怎么就攤上這事⌒偌” “怎么了丧慈?”我有些...
    開封第一講書人閱讀 153,116評(píng)論 0 344
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)主卫。 經(jīng)常有香客問我逃默,道長(zhǎng),這世上最難降的妖魔是什么簇搅? 我笑而不...
    開封第一講書人閱讀 55,371評(píng)論 1 279
  • 正文 為了忘掉前任完域,我火速辦了婚禮,結(jié)果婚禮上瘩将,老公的妹妹穿的比我還像新娘吟税。我一直安慰自己,他們只是感情好姿现,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,384評(píng)論 5 374
  • 文/花漫 我一把揭開白布乌妙。 她就那樣靜靜地躺著,像睡著了一般建钥。 火紅的嫁衣襯著肌膚如雪藤韵。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,111評(píng)論 1 285
  • 那天熊经,我揣著相機(jī)與錄音泽艘,去河邊找鬼。 笑死镐依,一個(gè)胖子當(dāng)著我的面吹牛匹涮,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播槐壳,決...
    沈念sama閱讀 38,416評(píng)論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼然低,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起雳攘,我...
    開封第一講書人閱讀 37,053評(píng)論 0 259
  • 序言:老撾萬榮一對(duì)情侶失蹤带兜,失蹤者是張志新(化名)和其女友劉穎,沒想到半個(gè)月后吨灭,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體刚照,經(jīng)...
    沈念sama閱讀 43,558評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,007評(píng)論 2 325
  • 正文 我和宋清朗相戀三年喧兄,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了无畔。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 38,117評(píng)論 1 334
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡吠冤,死狀恐怖浑彰,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情拯辙,我是刑警寧澤郭变,帶...
    沈念sama閱讀 33,756評(píng)論 4 324
  • 正文 年R本政府宣布,位于F島的核電站薄风,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏拍嵌。R本人自食惡果不足惜遭赂,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,324評(píng)論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望横辆。 院中可真熱鬧撇他,春花似錦、人聲如沸狈蚤。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽脆侮。三九已至锌畸,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間靖避,已是汗流浹背潭枣。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評(píng)論 1 262
  • 我被黑心中介騙來泰國(guó)打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留幻捏,地道東北人盆犁。 一個(gè)月前我還...
    沈念sama閱讀 45,578評(píng)論 2 355
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像篡九,于是被迫代替她去往敵國(guó)和親谐岁。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,877評(píng)論 2 345

推薦閱讀更多精彩內(nèi)容