#尋找通路和基因的關(guān)系
library(clusterProfiler)
mmu_kegg <- clusterProfiler::download_KEGG("mmu")
names(mmu_kegg)
head(mmu_kegg$KEGGPATHID2NAME)
head(mmu_kegg$KEGGPATHID2EXTID)
PATH2ID <- mmu_kegg$KEGGPATHID2EXTID
PATH2NAME <- mmu_kegg$KEGGPATHID2NAME
PATH_ID_NAME <- merge(PATH2ID, PATH2NAME, by="from")
colnames(PATH_ID_NAME) <- c("KEGGID", "ENTREZID", "DESCRPTION")
head(PATH_ID_NAME)
#轉(zhuǎn)換EntrezID為Gene Symbol
PATH_ID_NAME$geneid <- PATH_ID_NAME$ENTREZID
head(PATH_ID_NAME)
View(map_dt)
# transform id
map_dt <- bitr(PATH_ID_NAME$ENTREZID, fromType = "ENTREZID",
toType = c( "SYMBOL"),
OrgDb = org.Mm.eg.db)
dt_merge <- merge(map_dt, PATH_ID_NAME, by.y = "geneid", by.x = "ENTREZID")
View(dt_merge)
【單細胞轉(zhuǎn)錄組】KEGG通路反推基因列表
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
- 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來疚俱,“玉大人劝术,你說我怎么就攤上這事〈艮龋” “怎么了夯尽?”我有些...
- 文/不壞的土叔 我叫張陵,是天一觀的道長登馒。 經(jīng)常有香客問我匙握,道長,這世上最難降的妖魔是什么陈轿? 我笑而不...
- 正文 為了忘掉前任圈纺,我火速辦了婚禮,結(jié)果婚禮上麦射,老公的妹妹穿的比我還像新娘蛾娶。我一直安慰自己,他們只是感情好潜秋,可當我...
- 文/花漫 我一把揭開白布蛔琅。 她就那樣靜靜地躺著,像睡著了一般峻呛。 火紅的嫁衣襯著肌膚如雪罗售。 梳的紋絲不亂的頭發(fā)上辜窑,一...
- 文/蒼蘭香墨 我猛地睜開眼放钦,長吁一口氣:“原來是場噩夢啊……” “哼色徘!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起操禀,我...
- 正文 年R本政府宣布,位于F島的核電站护戳,受9級特大地震影響翎冲,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜媳荒,卻給世界環(huán)境...
- 文/蒙蒙 一抗悍、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧钳枕,春花似錦缴渊、人聲如沸。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至,卻和暖如春指蚁,著一層夾襖步出監(jiān)牢的瞬間菩佑,已是汗流浹背。 一陣腳步聲響...
推薦閱讀更多精彩內(nèi)容
- 我們知道在單細胞轉(zhuǎn)錄組數(shù)據(jù)分析的下游,在我們走過了fig1(通常是個圖譜)之后牲剃,往往要沿著基因這條路深挖一下生物學(xué)...
- 差異分析完以后凿傅,就有了基因列表和差異倍數(shù)缠犀,有了這兩個東西,就可以用clusterprofiler做GSEA聪舒,然后有...
- 通過單細胞轉(zhuǎn)錄組測序剖析人類胎兒消化道發(fā)育過程中自噬相關(guān)基因的動態(tài)表達 利用單細胞轉(zhuǎn)錄組揭示人類胚胎腸道發(fā)育過程中...
- 轉(zhuǎn)載來自:劉小澤 鏈接:http://www.reibang.com/p/e659a2f164f7[https:...