其他
跟着Nature学绘图——双向柱状图！
点击上方
“科研后花园”
关注我们
复现图片如下：
1、设置工作环境并加载所需R包：
#设置工作环境rm(list=ls())setwd("D:/nature文章原图复现系列/代码/双向柱状图")
#加载包library(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphics2、加载绘图数据（根据原文图片随机编写数据，无实际意义）：
#加载数据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))温馨提示
如果你喜欢本文，请分享到朋友圈，想要获得更多信息，请关注我。