跟着 Cell 学作图 | 3.箱线图+散点+差异显著性检验
跟着 Cell 学作图 | 3.箱线图+散点+差异显著性检验
“实践是检验真理的唯一标准。”
“复现是学习R语言的最好办法。”
DOI: 10.1016/j.cell.2020.05.032
这篇2020年发表在cell
上关于新冠的组学文章里面有大量的生信内容。今天带大家复现其中的一个Supplemental Figure
:带散点的箱线图。
本文代码及示例数据领取:后台回复“20210416”
比较常用的展示差异的可视化方法。
数据格式
注:以上数据纯属虚构,无实际意义!
思路
箱线图+显著性检验+散点
均值比较的常见方法:
方法 | R实现函数 | 描述 |
---|---|---|
T-test | t.test() | 比较两组(参数检验) |
Wilcoxon test | wilcox.test() | 比较两组(非参数检验) |
ANOVA | aov()或anova() | 比较多组(参数检验) |
Kruskal-Wallis | kruskal.test() | 比较多组(非参数检验) |
绘制
#------
title: "boxplot"
author: "MZBJ"
date: "2020/4/16"
#-----# 导入所需的包
library(ggplot2)
library(ggsignif)
library(ggpubr)
library(RColorBrewer)
setwd("F:/HJH/mzbj/cell")
# 导入数据
plot_data <- read.csv(file = "ABCD.csv")
#-----------------------
p<- ggplot(data=plot_data)+ geom_boxplot(mapping=aes(x=group,y=Retive_Abundance,colour = group ), #箱线图alpha = 0.5,size=1.5,width = 0.6)+ geom_jitter(mapping=aes(x=group,y=Retive_Abundance,colour = group), #散点alpha = 0.3,size=3)+scale_color_manual(limits=c("A","B","C","D"), values=c("#85B22E","#5F80B4","#E29827","#922927"))+ #颜色geom_signif(mapping=aes(x=group,y=Retive_Abundance), # 不同组别的显著性comparisons = list(c("A", "B"), # 哪些组进行比较c("A", "C"),c("A", "D"),c("B", "C"),c("B", "D"),c("C", "D")),map_signif_level=T, # T显示显著性,F显示p valuetip_length=c(0,0,0,0,0,0,0,0,0,0,0,0), # 修改显著性线两端的长短y_position = c(40,41,42,39,38,40), # 设置显著性线的位置高度size=1, # 修改线的粗细textsize = 4, # 修改显著性标记的大小test = "t.test")+ # 检验的类型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
大功告成!
备注:有个瑕疵就是y轴过长。但是我设定y轴范围后,差异标记也被截掉了。目前能想到的解决办法就是AI简单向下拉一下Y轴。
往期内容:
跟着CELL学作图|1.火山图
跟着Cell学作图 | 2.柱状图+误差棒+散点+差异显著性检验
跟着 Cell 学作图 | 3.箱线图+散点+差异显著性检验相关推荐
- 跟着Cell学作图 | 12.韦恩图(Vennerable包)
"实践是检验真理的唯一标准." "复现是学习生信的最好办法." 2021.4.12_1 DOI: 10.1016/j.cell.2020.05.032 这篇20 ...
- 跟着 Cell 学作图 | 4.小提琴图
跟着 Cell 学作图 | 4.小提琴图 "实践是检验真理的唯一标准." "复现是学习R语言的最好办法." DOI: 10.1016/j.cell.2020.0 ...
- 跟着Cell学作图|9.PPI分析(GeNets数据库)
9.PPI分析(GeNets数据库) "实践是检验真理的唯一标准." "复现是学习R语言的最好办法." DOI: 10.1016/j.cell.2020.05. ...
- 跟着Cell学作图 | 5.UMAP降维分析
跟着 Cell 学作图 | 5.UMAP降维分析 "实践是检验真理的唯一标准." "复现是学习R语言的最好办法." 2021.4.12_1 DOI: 10.10 ...
- 跟着CELL学作图|1.火山图
跟着CELL学作图之火山图 "实践是检验真理的唯一标准." "复现是学习R语言的最好办法." DOI: 10.1016/j.cell.2020.05.032 这 ...
- 跟着 Cell 学作图 | 桑葚图(ggalluvial)
桑葚图 今天我们复现一幅2021年Cell上Graphical abstract的图. Title:Human oral mucosa cell atlas reveals a stromal-neu ...
- 跟着Cell学作图 | 2.柱状图+误差棒+散点+差异显著性检验
跟着 Cell 学作图 | 2.柱状图+误差棒+散点 "实践是检验真理的唯一标准." "复现是学习R语言的最好办法." 2021.4.12_1 DOI: 10. ...
- 跟着Nat Commun学作图 | 4.配对箱线图+差异分析
跟着Nat Commun学作图 | 4.配对箱线图+差异分析 今天要学习的图来自2021年10月29号发表在的Nature Communication上的一篇文章,题目是[新冠肺炎患者呼吸道菌群组成及 ...
- 跟着 NC 学作图 | 多组散点图+配对连线+差异分析
scatter_paired_line 今天我们复现一幅2022年7月发表在nature communications上的配对连线散点图. ❝ Title:Evaluation of transpla ...
最新文章
- 女装、退学、离职谷歌,他活成了别人不敢想的样子
- 这个勒索软件也太菜了!
- 今日头条的面试题(部分)
- 解决goldengate复制进程应用缓慢一例
- 三维重建1-位姿追踪:单应矩阵、本质矩阵和基本矩阵
- SAP系统和微信集成的系列教程之五:如何将SAP UI5应用嵌入到微信公众号菜单中
- 代码挑战“ Vrolijke Framboos”事后验尸
- 2021年中国电竞运动行业发展报告
- 液化气调压阀爆炸及使用问题分析
- linux中pad命令,1.1Linux初级命令(附例子)
- Spring-ConfigurationClassPostProcessor类
- linux运维脚本编写,Linux运维基础技能: 脚本编程与Linux命令
- 4.1 HTML5 音频
- Mysql 单表、连表 in 子查询条件拼接
- MAYA制作女猎人角色毛发-3D建模场景模型教程
- 企业应如何改善设备功率因数,有效规避罚款
- Linux服务器使用网络代理
- 关于键盘(总论8042)
- python之urlencode(),quote()及unquote()
- 如何去远程控制公司电脑 这三样必备一样也不能少