【R語言】解決GO富集分析繪圖，標簽重疊問題

前面我給大家詳細介紹過

創新互聯公司專注于企業網絡營銷推廣、網站重做改版、嘉興網站定制設計、自適應品牌網站建設、HTML5、商城網站建設、集團公司官網建設、成都外貿網站建設公司、高端網站制作、響應式網頁設計等建站業務，價格優惠性價比高，為嘉興等各大城市提供網站開發制作服務。

?GO簡介及GO富集結果解讀

?四種GO富集柱形圖、氣泡圖解讀

?GO富集分析四種風格展示結果—柱形圖，氣泡圖

?KEGG富集分析—柱形圖，氣泡圖，通路圖

? DAVID GO和KEGG富集分析及結果可視化

也用視頻給大家介紹過

? GO和KEGG富集分析視頻講解

最近有粉絲反映說，利用clusterProfiler這個包繪制GO富集分析氣泡圖和柱形圖的時候，發現GO條目的名字都重疊在一起了。

氣泡圖

柱形圖

這個圖別說美觀了，簡直不忍直視。經過我的認真研究，發現跟R版本有關。前面我給大家展示的基本都是R 3.6.3做出來的圖。很多粉絲可能用的都是最新版本的R 4.1.2。

我們知道R的版本在不停的更新，相應的R包也在不停的更新。我把繪制氣泡圖和柱形圖相關的函數拿出來認真的研究了一下，終于發現的癥結碰侍所在。

dotplot這個函數，多了個 label_format 參數

我們來看看這個參數究竟是干什么用的，看看參數說明

label_format :

a numeric value sets wrap length, alternatively a custom function to format axis labels. by default wraps names longer that 30 characters

原來這個參數默認值是30，當標簽的長度大于30個字符就會被折疊，用多行來展示。既然問題找到了，我們就來調節一下笑咐吵這個參數，把他設置成100，讓我們的標簽可以一行展示。

是不是還是原來的配方，還是熟悉的味道

同樣的柱形圖，我們也能讓他恢復原來的容貌。

關于如何使用R做GO和簡閉KEGG富集分析，可參考下文

GO和KEGG富集分析視頻講解

[R語言] GO富集分析可視化 GOplot::GOCircle

查看GOplot內示例數據的格式，對自己的數據做處理

觀察結論：

觀察自己的兩個數據表：

table.legend 設置為T時會顯示表格

本圓族圖中表格和隱擾圖例是出圖后灶腔旦剪切拼合而成，沒有用R中的拼圖包

轉錄組不求人系列(十三): GO、KEGG富集個性化作圖

當富集分析完成，拿到如下的分析結果后，就可以進行作圖了。

富集分析結果的可視化無非就是柱狀圖和氣泡圖，但是公司默認出圖實在是太丑，所以還是自己動手修改修改。

一、常規柱狀圖（ggplot2）

橫軸為gene counts，或者用-logP也行，填充相應的用P值或者gene counts。ggplot畫圖的好處就是可以進行很多調整。

二、常規氣泡圖（ggplot2）

氣泡圖與柱狀圖碰春如出一轍，只是在展示方式上出現了差別。一個用geom_bar()函數，氣泡圖類似于散點圖用geom_point()函數。

三、上下調同時展示（ggplot2）

很多時候研究者拿到差異基因后，上下調基因是分別富集的，在展示上需要同時體現二者，我們之前提到metascape可以做到： 轉錄組不求人系列(十二): Cell文章最喜歡用的差異基因GO、KEGG富集分析工具 ，除此之外，之前講過的氣泡圖也可以展現多組的結果： 復現《nature communications》圖表（四）：ggplot畫多組富集氣泡圖 。這里我們繼續提供一種bar圖的展顫吵租示方式。將茄兆down的數值調整為負，做一列分組，就可以展示了。

當然了，以上所說的可視化還是比較常規，在基礎上可以自己做調整。也有一些文章總是標新立意，有很多奇特的展現方式，我們會在之后的系列中講解。

R語言：clusterProfiler進行GO富集分析和Gene_ID轉換

ID轉換用到豎蔽的是 bitr() 函數，bitr()的使用方法：

org.Hs.eg.db包含有多種gene_name的類型

keytypes() ：keytypes(x)，查看注釋包中可以使用的類型耐拿

columns() ：類似于keytypes()，針對org.Hs.eg.db兩個函數返回值一致

select() ：select(x, keys, columns, keytype, ...) eg.

函數enrichGO()進行GO富集分析余畝州，enrichGO()的使用方法：

舉例：

當前標題：go分析作圖r語言 r語言igraph畫圖
當前URL：http://m.newbst.com/article28/ddpdocp.html

成都網站建設公司_創新互聯，為您提供品牌網站制作、定制開發、網站策劃、移動網站建設、、標簽優化

廣告

聲明：本網站發布的內容（圖片、視頻和文字）以用戶投稿、用戶轉載內容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網站立場，如需處理請聯系客服。電話：028-86922220；郵箱：631063699@qq.com。內容未經允許不得轉載，或轉載時需注明來源： 創新互聯