生信分析中经常会得到一些基因,然后做GO富集分析,达到对基因进行注释和分类的目的。
本文利用R语言的ggplot2包,从头带您绘制可发表级别的GO富集分析结果图。
一 载入数据集和R包
利用各种生信工具得到富集分析结果,数据列可能不一致,但关键几列都有。
library(ggplot2)data "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)