多個基因集富集結果展示

通常我們會同時對多個基因集分別進行富集分析，結果放在一起展示。這時我們需要在富集結果后面加一列，標記該結果是哪個基因集的富集，在Excel中可以很方便地操作。如下面動圖所示，分組的名字自己根據實際取名即可。

有了這個多組基因富集后整合起來的數據，就可以用BIC繪圖了。數據粘貼就不展示了，直接看參數選擇。

與單組富集結果相比，最大的改動就在：

• 新增的Group列而非 log_odds_ratio列作為橫軸(X-axis)信息

提交后獲得結果。圖中每個點代表一個富集的條目，在Y軸有對應標記。每一列是一組基因的富集結果。三組共有的富集在最上面，2組共有的富集在中間，每組特有的富集在底部。每個點的大小代表用于分析的基因集中匹配到該通路的基因數目，顏色代表富集程度。

如果希望在顯示多組時，依然保留log_odds_ratio的信息，也可以。這里換一套數據更好展示(因為Group2、Group3是模擬數據，直接從Group1中抽取出來的，所以繪制出來會存在重疊)

GOID Ontology Term Level q m t k log_odds_ratio p Group GO:0006730 biological_process one-carbon metabolic process 4 340 57 45240 13378 1.012309306 0.001481151 Group1 GO:0007154 biological_process cell communication 2 2169 6843 45240 13378 0.100137585 0.007326261 Group1 GO:0007165 biological_process signal transduction 5 1955 6136 45240 13378 0.107606604 0.006325629 Group1 GO:0023052 biological_process signaling 1 2100 6613 45240 13378 0.102820905 0.006590727 Group1 GO:0044700 biological_process single organism signaling 2 2100 6613 45240 13378 0.102820905 0.006590727 Group1 GO:0050896 biological_process response to stimulus 1 3251 10438 45240 13378 0.074846633 0.012472089 Group1 GO:0006730 biological_process one-carbon metabolic process 4 340 57 45240 13378 1.212309306 0.001481151 Group2 GO:0007154 biological_process cell communication 2 2169 6843 45240 13378 0.200137585 0.007326261 Group2 GO:0007165 biological_process signal transduction 5 1955 6136 45240 13378 0.207606604 0.006325629 Group2 GO:0023052 biological_process signaling 1 2100 6613 45240 13378 0.302820905 0.006590727 Group2 GO:0007165 biological_process signal transduction 5 1955 6136 45240 13378 0.307606604 0.006325629 Group3 GO:0023052 biological_process signaling 1 2100 6613 45240 13378 0.202820905 0.006590727 Group3

• log_odds_ratio列依然作為橫軸(X-axis)信息

• 新增的Group列作為Shape variable，用不同的形狀表示不同的組

• Shape variable order是可選項，調節組的順序，默認不填寫或按需設置都可

• 修改下顏色，用colorPicker設置，前面工具有介紹

提交后獲得結果。圖中每個點代表一個富集的條目，在Y軸有對應標記。這些條目按其log_odds_ratio的值排序后展示，log_odds_ratio高的條目在Y軸上方展示；每個點的大小代表用于分析的基因集中匹配到該通路的基因數目，顏色代表富集程度。點的形狀則代表其所屬的組信息。

但是這個圖出現了一個問題，圖例顯示不全。最簡單的解決辦法就是把圖的寬度和高度調大。

結果就正常了，可以下載PDF版、PPT版（如果選了參數）和對應的R代碼

往期精品(點擊圖片直達文字對應教程)

機器學習

后臺回復“生信寶典福利第一波”或點擊閱讀原文獲取教程合集

創作挑戰賽新人創作獎勵來咯，堅持創作打卡瓜分現金大獎

總結

以上是生活随笔為你收集整理的多个基因集富集结果泡泡图绘制展示的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。