跟着Nature学绘图——双向柱状图！

 

点击上方

“科研后花园”

关注我们

复现图片如下：

1、设置工作环境并加载所需R包：

#设置工作环境
rm(list=ls())
setwd("D:/nature文章原图复现系列/代码/双向柱状图")

#加载包
library(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphics

2、加载绘图数据（根据原文图片随机编写数据，无实际意义）：

#加载数据
df <- read.table("data.txt",header = T,sep='\t')
df$group <- factor(df$group,levels = c("HPRC.EAS","CPC"))
df$sample <- factor(df$sample,levels = c("OR4K2","OR4K1","GRPIN2","OR9G1","SULT1A3","OR4Q3",
                                         "PDPR","CYP2D6","PRAMEF18"))
#数据处理——将其中一组的数据转换为负值
df$value <- ifelse(df$group=="HPRC.EAS",-df$value,df$value)

3、绘图：

#颜色
col <- c("#f89f68","#4b84b3")
#绘图
ggplot(df,aes(sample,value,fill=group))+
  #柱状图
  geom_col(width = 0.8)+
  #坐标转换
  coord_flip()+
  #自定义颜色
  scale_fill_manual(values = col)+
  #主题设置
  theme_bw()+
  theme(panel.grid = element_blank(),
        axis.text.x = element_text(color = "black", size = 13),
        axis.text.y = element_text(color = "black",size = 14,face = "italic"),
        axis.title = element_text(color = "black",size = 16),
        legend.position = "top",
        legend.title = element_blank(),
        legend.text = element_text(color = "black",size = 15))+
  #标题设置
  labs(x=NULL,y="CNV Frequency")+
  #y轴范围设置
  scale_y_continuous(breaks = seq(-1, 1, 0.5), 
                     labels = as.character(abs(seq(-1, 1, 0.5))),
                     limits = c(-1, 1))

PS: 以上内容是小编个人学习代码笔记分享，仅供参考学习，欢迎大家一起交流学习。

代码及源数据获取：https://github.com/wzsBio/code

温馨提示

如果你喜欢本文，请分享到朋友圈，想要获得更多信息，请关注我。

更多推荐