科研有捷徑，輸入代碼，一鍵獲取科研成果！就是這么省事，來具體看下有多方便！

搜索http://985.so/a9kb查看全部代碼（目前共計50+持續(xù)新增中），也可以點擊右側(cè)【目錄】，可以看到更多有趣的代碼；真香提示：文末可以知道如何獲取代碼~

場景需求：

GO功能富集圖我們常見的是點圖和柱狀圖，像這樣的：

或者像我們之前推文里介紹的網(wǎng)絡(luò)圖一樣，詳情查看推文“代碼分享│R可視化：基因與功能之間的關(guān)系--GO功能富集網(wǎng)絡(luò)圖繪制”（鏈接：https://mp.weixin.qq.com/s/TdTKHsKhZMl1Ukjlz4MjkQ）。

上面的點圖、柱狀圖和網(wǎng)絡(luò)圖體現(xiàn)的是總體差異基因的GO富集結(jié)果，并不能直觀的體現(xiàn)組間差異基因富集的結(jié)果，但是我們實際的需求可能是下圖這樣的，對差異上下調(diào)基因分別富集到的功能進行展示并繪制在一張圖上：

下面我們來具體看下怎么操作。代碼相關(guān)文件見如下3個文件夾，可在文末知道如何領(lǐng)取。

1.????? 加載相應(yīng)R依賴包：

library("clusterProfiler")

library("org.Mm.eg.db")

library("enrichplot")

library("ggplot2")

2.????? 數(shù)據(jù)導(dǎo)入

輸入數(shù)據(jù)格式如下所示，只需要2列即可，一列是分組信息，一列是基因ID。

3.????? 基因ID轉(zhuǎn)換

我們將基因symbol轉(zhuǎn)換成ENTREZID，OrgDb參數(shù)是選擇物種的，這里我們使用的是小鼠基因，所以是org.Mm.eg.db。

4.????? 常規(guī)GO富集圖

我們先得到常規(guī)的GO富集點圖和柱狀圖，ont這個參數(shù)還可以選擇"BP", "MF", "CC"這三個功能。

5.????? 組間差異比較GO富集圖

接下來就是最重要的組間比較結(jié)果圖了，輸入文件必須要有上面提到的2列，其中一列是分組信息，然后利用compareCluster這個函數(shù)進行分析，再利用simplify去除冗余的結(jié)果。

最后是出圖并輸出整體結(jié)果?？梢钥吹浇Y(jié)果展示了兩組間各自富集的GO功能，這樣就能很好的展示各自組特異的功能富集結(jié)果。

如需代碼及示例數(shù)據(jù)等文件，請掃碼聊天框回復(fù) “B30”領(lǐng)??！?