这图怎么画| 富集分析之双向柱状图

# 创建示例数据
dt = data.frame(pathway = paste('pathway',1:20),
                Classification = rep(paste('Class',1:5),4),
                pvalue  = runif(20,0.0001,0.05),
                count = floor(runif(20,5,20)))
head(dt)  
dt = dt[order(dt$Classification),]
dt$pathway = factor(dt$pathway, levels =  dt$pathway)   
library(ggplot2)
# 左图
p1 <- ggplot(dt, aes( x = pathway,
                y = log10(pvalue), #为了方向向左，故不加负号
                fill = Classification))+
  scale_fill_manual(values = rev(c('#e4e242','#8ab92b','#80c5d8','#437eb8','#c4a2c5')))+ # rev是因为颜色写反了
  scale_y_continuous(expand = c(0,0),labels = c(3,2,1,0)) +#标签为负号，自行修改
  geom_bar(stat = 'identity')+
  coord_flip()+
  ylab('-log10(P-value)')+
  xlab('')
p1
# 右图
p2 <- ggplot(dt, aes( x = pathway,
                      y = count),fill = '#595758')+
  scale_y_continuous(expand = c(0,0)) +
  geom_bar(stat = 'identity')+
  coord_flip()+
  ylab('No. of genes')+
  xlab('') +
  theme(axis.text.y = element_blank(),
        axis.ticks.y = element_blank(),
        axis.title.y = element_blank())
p2
library(patchwork)
p1+p2+plot_layout(guides = 'collect')+ plot_layout(widths = c(2, 1.5))
# 保存
ggsave('enrich_barplot.pdf',width = 7,height = 5)