這篇文章給大家分享的是有關如何利用R語言的ggplot2包繪制GO富集柱形圖的內容。小編覺得挺實用的,因此分享給大家做個參考,一起跟隨小編過來看看吧。
為通化等地區用戶提供了全套網頁設計制作服務,及通化網站建設行業解決方案。主營業務為成都網站制作、網站建設、外貿網站建設、通化網站設計,以傳統方式定制建設網站,并提供域名空間備案等一條龍服務,秉承以專業、用心的態度為用戶提供真誠的服務。我們深信只要達到每一位用戶的要求,就會得到認可,從而選擇與我們長期合作。這樣,我們也可以走得更遠!
一 載入數據集和R包
利用各種生信工具得到富集分析結果,數據列可能不一致,但關鍵幾列都有。
library(ggplot2)data <- read.csv("GO_enrichment_significant.csv",header=TRUE)head(data)
二 對上述GO結果繪制基礎bar圖
參照之前ggplot2使用方法,更改geom即可繪制簡單的bar圖,按照GO_category分組顏色
ggplot(data=data, aes(x=GO_term,y=Num_of_symbols_in_list_in_GO, fill=GO_category)) + geom_bar(stat="identity", width=0.8)
可看出和文獻中的差距較大,體現在以下幾個方面:
A:標題,坐標軸“業余”;
B:GO_category順序未按照輸入文件,相同GO_category沒在一起;
C:橫坐標label太長,重疊在一起。
三 “細節”調整GO結果bar圖
3.1 坐標軸調整策略
#將GO_term設定為factor即可按照順序輸出GO_term_order=factor(as.integer(rownames(data)),labels=data$GO_term)ggplot(data=data, aes(x=GO_term_order,y=Num_of_symbols_in_list_in_GO, fill=GO_category)) + geom_bar(stat="identity", width=0.8) + coord_flip() + xlab("GO term") + ylab("Num of Genes") + theme_bw()
好像有一點能看了,嘗試其他策略。
3.2 調整橫坐標label策略
將label調整成一定角度傾斜
COLS <- c("#66C3A5", "#8DA1CB", "#FD8D62")ggplot(data=data, aes(x=GO_term_order,y=Num_of_symbols_in_list_in_GO, fill=GO_category)) + geom_bar(stat="identity", width=0.8) + scale_fill_manual(values = COLS) + theme_bw() + xlab("GO term") + ylab("Num of Genes") + labs(title = "The Most Enriched GO Terms")+ theme(axis.text.x=element_text(face = "bold", color="gray50",angle = 70,vjust = 1, hjust = 1 ))
嗯 ,標簽太長溢出,采取保留GO-term的前三個單詞(可以其他策略)后面...代替,可以excel或者R function 解決。
3.3 調整label長度后繪圖
GO_term_order=factor(as.integer(rownames(data)),labels=labels)
COLS <- c("#66C3A5", "#8DA1CB", "#FD8D62")
ggplot(data=data, aes(x=GO_term_order,y=Num_of_symbols_in_list_in_GO, fill=GO_category)) +
geom_bar(stat="identity", width=0.8) +
scale_fill_manual(values = COLS) + theme_bw() +
xlab("GO term") + ylab("Num of Genes") + labs(title = "The Most Enriched GO Terms")+ theme(axis.text.x=element_text(face = "bold", color="gray50",angle = 70,vjust = 1, hjust = 1 ))
感謝各位的閱讀!關于“如何利用R語言的ggplot2包繪制GO富集柱形圖”這篇文章就分享到這里了,希望以上內容可以對大家有一定的幫助,讓大家可以學到更多知識,如果覺得文章不錯,可以把它分享出去讓更多的人看到吧!
網頁標題:如何利用R語言的ggplot2包繪制GO富集柱形圖
網站鏈接:http://vcdvsql.cn/article30/jhgdso.html
成都網站建設公司_創新互聯,為您提供網站營銷、搜索引擎優化、標簽優化、、靜態網站、軟件開發
聲明:本網站發布的內容(圖片、視頻和文字)以用戶投稿、用戶轉載內容為主,如果涉及侵權請盡快告知,我們將會在第一時間刪除。文章觀點不代表本網站立場,如需處理請聯系客服。電話:028-86922220;郵箱:631063699@qq.com。內容未經允許不得轉載,或轉載時需注明來源: 創新互聯