今天我们做一下这篇《Cell》文章的富集结果图。

这是一个双侧柱状图，特点是，可以理解为上下调分开展示。当然，还可以将P值等信息添加在图上。其实我们之前好几篇文章做过双侧柱状图，但不足之处是没有像这篇文章一样，坐标轴位于中间，且通路分别展示在两侧，以及添加更多信息！

接下来盘它！

读入数据，可以将上下调提前标记分开。

setwd("E:/生物信息学/复现Cell双侧柱状图GO")

A <- read.csv('GO.csv', header = T)

library(ggplot2)

library(tidyverse)

A$group <- ''

A$group[which(A$ratio >0)]='up'

A$group[which(A$ratio <0)]='down'

ggplot作图：

ggplot(A,aes(reorder(Description, ratio),ratio,fill=group))+

geom_col()+

theme_bw()+

theme(panel.grid.major=element_blank(),

panel.grid.minor=element_blank(),

panel.border = element_blank(),

legend.title = element_blank(),

axis.text = element_text(color="black",size=10),

axis.line.x = element_line(color='black'),

axis.ticks.y = element_blank(),

axis.text.y = element_blank(),

legend.position = 'none')+

coord_flip()+

geom_segment(aes(y=0, yend=0,x=0,xend=18.5))+

geom_text(data = A[which(A$ratio>0),],aes(x=Description, y=-0.01, label=Description),

hjust=1, size=4)+

geom_text(data = A[which(A$ratio<0),],aes(x=Description, y=0.01, label=Description),

hjust=0, size=4)+

geom_text(data = A[which(A$ratio>0),],aes(label=Padj),

hjust=-0.1, size=4, color='red')+

geom_text(data = A[which(A$ratio<0),],aes(label=Padj),

hjust=1.1, size=4, color="red")+

scale_fill_manual(values = c("#1084A4",

"#8D4873"))+

scale_x_discrete(expand = expansion(mult = c(0,0)))+

ylim(-0.5, 0.5)+

labs(x='', y='Ratio')

转自：KS科研分享与服务

如有侵权，请联系本站删除！