跟着Cell学作图 | 2.柱状图+误差棒+散点+差异显著性检验

简介: 这篇2020年发表在cell上关于新冠的组学文章里面有大量的生信内容。今天带大家复现其中的一个Supplemental Figure:带散点的柱状图。

跟着 Cell 学作图 | 2.柱状图+误差棒+散点


数据格式


表一用于制作散点图


image.png


2021.4.14_2

表二用于制作柱状图。(对表一进行统计)


image.png


2021.4.14_3

注:以上数据纯属虚构,无实际意义!


思路


柱状图+显著性检验+误差棒+散点


绘制


---
title: "barplot"
author: "MZBJ"
date: "2020/4/14"
---
# 导入所需的包
library(ggplot2)
library(ggsignif)
library(ggpubr)
setwd("F:/HJH/mzbj/cell") 
# 导入数据
plot_data1 <- read.csv(file = "NPT_sod.csv")
plot_data2 <- read.csv(file = "NPT.csv")
#-----------------------
p<- ggplot()+ 
  geom_bar(data=plot_data1,mapping=aes(x=group,y=mean), #柱状图
           fill = "white",
           size = 1.5,
           color = c("#1E90FF","gold","red"),
           position="dodge", # 柱状图格式
           stat="identity", # 数据格式
           width = 0.6)+  # 柱状图尺寸
  geom_jitter(data=plot_data2, #散点
              mapping=aes(x=group,y=Retive_Abundance,fill = group,colour = group,shape = group),
              size = 2,
              height = 0.05,#散点纵向抖动
              width = 0.1)+ #散点横向抖动
  scale_color_manual(values = c("#1E90FF","gold","red"))+ #散点颜色
  geom_signif(data=plot_data2,mapping=aes(x=group,y=Retive_Abundance), # 不同组别的显著性
              comparisons = list(c("N", "P"), # 哪些组进行比较
                                 c("N", "T"),
                                 c("P", "T")),
              annotation=c("**"), # 显著性差异做标记
              map_signif_level=T, # T为显著性,F为p value
              tip_length=c(0,0,0,0,0,0), # 修改显著性线两端的长短
              y_position = c(46,54,49), # 设置显著性线的位置高度
              size=1, # 修改线的粗细
              textsize = 7, # 修改*标记的大小
              test = "t.test")+ # 检验的类型
  geom_errorbar(data=plot_data1,mapping=aes(x = group,ymin = mean-sd, ymax = mean+sd), # 误差线添加
                width = 0.3, #误差线的宽度
                color = c("#1E90FF","gold","red"), #颜色
                size=0.8)+ #粗细
  scale_y_continuous(limits =c(0, 60) ,expand = c(0,0))+ # y轴的范围
  theme_classic(  # 主题设置,这个是无线条主题
    base_line_size = 1 # 坐标轴的粗细
  )+
  labs(title="White blood cell(WBC)",x="",y="Retive_Abundance")+ # 添加标题,x轴,y轴内容
  theme(plot.title = element_text(size = 15,
                                  colour = "black",
                                  hjust = 0.5),
        axis.title.y = element_text(size = 15, 
                                    # family = "myFont", 
                                    color = "black",
                                    face = "bold", 
                                    vjust = 1.9, 
                                    hjust = 0.5, 
                                    angle = 90),
        legend.title = element_text(color="black", # 修改图例的标题
                                    size=15, 
                                    face="bold"),
        legend.text = element_text(color="black", # 设置图例标签文字
                                   size = 10, 
                                   face = "bold"),
        axis.text.x = element_text(size = 13,  # 修改X轴上字体大小,
                                   color = "black", # 颜色
                                   face = "bold", #  face取值:plain普通,bold加粗,italic斜体,bold.italic斜体加粗
                                   vjust = 0.5, # 位置
                                   hjust = 0.5, 
                                   angle = 0), #角度
        axis.text.y = element_text(size = 13,  
                                   color = "black",
                                   face = "bold", 
                                   vjust = 0.5, 
                                   hjust = 0.5, 
                                   angle = 0) 
  ) 
p

image.png

大功告成!

相关文章
R实战 | 对称云雨图 + 箱线图 + 配对散点 + 误差棒图 +均值连线
R实战 | 对称云雨图 + 箱线图 + 配对散点 + 误差棒图 +均值连线
2026 1
R实战 | 对称云雨图 + 箱线图 + 配对散点 + 误差棒图 +均值连线
|
数据可视化
R语言绘图教程丨Nature论文都在用的多组比较箱线图,自动计算显著性并标注,附带误差线
R语言绘图教程丨Nature论文都在用的多组比较箱线图,自动计算显著性并标注,附带误差线
sklearn中分类模型评估指标(一):准确率、Top准确率、平衡准确率
accuracy_score函数计算准确率分数,即预测正确的分数(默认)或计数(当normalize=False时)。 在多标签分类中,该函数返回子集准确率(subset accuracy)。 如果样本的整个预测标签集与真实标签集严格匹配,则子集准确率为 1.0; 否则为 0.0。
VS2019编译VTK-9.1.0+Qt5.15.2
本文介绍了在VS2019和Qt 5.15.2环境下编译VTK 9.1.0的过程,包括成功编译无警告和错误的截图、修改的pro文件内容,以及在QML中加入VTK的失败尝试和参考链接。
658 4
|
机器学习/深度学习 监控
在进行多任务学习时,如何确保模型不会过度拟合单一任务而忽视其他任务?
在进行多任务学习时,如何确保模型不会过度拟合单一任务而忽视其他任务?
251 2
|
Oracle 关系型数据库 Java
java操作多数据源将oracle数据同步达梦数据库
java操作多数据源将oracle数据同步达梦数据库
|
数据挖掘
跟着 NC 学作图 | 多组散点图+配对连线+差异分析
跟着 NC 学作图 | 多组散点图+配对连线+差异分析
481 0
【SPSS】相关分析和偏相关分析详细操作过程(附案例实战)
【SPSS】相关分析和偏相关分析详细操作过程(附案例实战)
1574 0
跟着 Cell 学作图 | 分组蜂群图+均值线+差异分析(组间+组内)
跟着 Cell 学作图 | 分组蜂群图+均值线+差异分析(组间+组内)
381 0
|
关系型数据库 数据库 MySQL