跟着PC学作图 | 小提琴图+Tufte箱形图+折线图的绘制教程

2023-11-10 18:10

本文主要是介绍跟着PC学作图 | 小提琴图+Tufte箱形图+折线图的绘制教程,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、 数据

本期教程是基于ggplot()绘制,数据类型也符合ggplot()函数类型即可。

本期,直接使用的文章中的数据。

二、绘图

导入相关包

library(ggplot2)
library(tidyverse)
library(reshape)
library(agricolae)library(ggthemes)

导入数据

diam <- read.csv("inputdata.csv", header = TRUE, sep = ',', stringsAsFactors=FALSE)
head(diam)

数据整理

此步更具自己数据需求操作

### Add groups to the data
diam$Groups <- paste(diam$Species,diam$CO2,diam$Day,sep="_")
diam$Groups <- factor(diam$Groups,levels = unique(gtools::mixedsort(diam$Groups)))

差异比较

maxDiam = diam %>% group_by(Groups) %>% summarize(MaxDiameter=max(diameter))### Do tukey tests
hsd_diameter=HSD.test(aov(diameter~Groups,data=diam), "Groups", group=T)
hsd_diameter$groups <- hsd_diameter$groups[levels(diam$Groups),]
tukeyDiameter <- data.frame(rownames(hsd_diameter$groups),do.call("rbind",strsplit(rownames(hsd_diameter$groups),"_")),hsd_diameter$groups,stringsAsFactors = T)
colnames(tukeyDiameter) <- c("Groups","Species","CO2","Day","median_diameter","groups")tukeyDiameter <- left_join(tukeyDiameter,meanDiam,by="Groups")
head(tukeyDiameter)
> head(tukeyDiameter)Groups Species      CO2 Day median_diameter groups
M82_ambient_7     M82_ambient_7     M82  ambient   7        321.9840     bc
M82_ambient_10   M82_ambient_10     M82  ambient  10        390.9108     ab
M82_ambient_13   M82_ambient_13     M82  ambient  13        430.4615      a
M82_elevated_7   M82_elevated_7     M82 elevated   7        316.3562     bc
M82_elevated_10 M82_elevated_10     M82 elevated  10        448.0713      a
M82_elevated_13 M82_elevated_13     M82 elevated  13        443.9278      a

ggplot绘图

  1. 箱线图
plotGG <- ggplot(diam,aes(x=Groups,y=diameter)) +scale_fill_manual(values=c("skyblue1", "black")) + scale_color_manual(values=c("skyblue1", "black")) +geom_violin(aes(color=paste(CO2)),position=position_dodge(1),trim = F,width=1)

Tufte箱形图

plotGG + geom_tufteboxplot(position=position_dodge(1),outlier.colour="transparent", width = 1,color="black")

line图

plotGG02 + geom_line(data=tukeyDiameter, aes(x=Groups, y=median_diameter,group = paste(CO2,Species)),alpha=0.5)

添加显著性差异

plotGG03 <- plotGG02 + geom_text(data=tukeyDiameter, aes(x=Groups, y=(maxDiam$MaxDiameter)*1.25,label=groups,group=NULL),vjust=0,color="black")

进行美化

plotGG03 +theme_bw()+ggtitle("Root diameter change by species in response to eCO2 over time") + xlab('')+ylab(expression(paste("Root diameter (",mu,"m)"))+theme(#legend.position="bottom",legend.text = element_text(size=12),axis.text.y = element_text(size=12, colour="black", family = "sans", angle = 0), axis.text.x = element_text(size=12, colour="black", family = "sans", angle = 45, hjust = 1),axis.title= element_text(size=12),strip.text.x = element_text(size=12, angle = 0),strip.text.y = element_text(size=12, angle = 0),plot.title = element_text(size=12, angle = 0))

绘图全部代码

ggplot(diam,aes(x=Groups,y=diameter)) +scale_fill_manual(values=c("skyblue1", "black")) + scale_color_manual(values=c("skyblue1", "black")) +geom_violin(aes(color=paste(CO2) ),position=position_dodge(1),trim = F,width=1) + #geom_boxplot(aes(color=paste(CO2)),position=position_dodge(1),width=0.09,outlier.shape = NA) + geom_tufteboxplot(position=position_dodge(1),outlier.colour="transparent", width = 1,color="black") + #geom_dotplot(aes(fill=CO2),binaxis='y', stackdir='center',position=position_dodge(1),dotsize = 0.5)  + #Looks cleaner than with#geom_jitter(aes(color=paste(CO2)),width = 0.05,size=0.75) + geom_line(data=tukeyDiameter, aes(x=Groups, y=median_diameter,group = paste(CO2,Species)),alpha=0.5) +facet_grid(.~Species,scales = "free_x",as.table = T) +coord_cartesian(ylim = c(min(diam$diameter)*.6, max(diam$diameter)*1.5)) +theme(axis.text.x = element_text(angle = 45, hjust = 1)) +## 添加显著性差异geom_text(data=tukeyDiameter, aes(x=Groups, y=(maxDiam$MaxDiameter)*1.25,label=groups,group=NULL),vjust=0,color="black")+theme_bw()+ggtitle("Root diameter change by species in response to eCO2 over time") + xlab('')+ylab(expression(paste("Root diameter (",mu,"m)")))+theme(#legend.position="bottom",legend.text = element_text(size=12),axis.text.y = element_text(size=12, colour="black", family = "sans", angle = 0), axis.text.x = element_text(size=12, colour="black", family = "sans", angle = 45, hjust = 1),axis.title= element_text(size=12),strip.text.x = element_text(size=12, angle = 0),strip.text.y = element_text(size=12, angle = 0),plot.title = element_text(size=12, angle = 0))

往期文章:

1. 最全WGCNA教程(替换数据即可出全部结果与图形)

WGCNA分析 | 全流程分析代码 | 代码一

WGCNA分析 | 全流程分析代码 | 代码二

WGCNA分析 | 全流程代码分享 | 代码三


2. 精美图形绘制教程

精美图形绘制教程


小杜的生信筆記,主要发表或收录生物信息学的教程,以及基于R的分析和可视化(包括数据分析,图形绘制等);分享感兴趣的文献和学习资料!!

这篇关于跟着PC学作图 | 小提琴图+Tufte箱形图+折线图的绘制教程的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/384391

相关文章

使用Docker构建Python Flask程序的详细教程

《使用Docker构建PythonFlask程序的详细教程》在当今的软件开发领域,容器化技术正变得越来越流行,而Docker无疑是其中的佼佼者,本文我们就来聊聊如何使用Docker构建一个简单的Py... 目录引言一、准备工作二、创建 Flask 应用程序三、创建 dockerfile四、构建 Docker

使用Python绘制3D堆叠条形图全解析

《使用Python绘制3D堆叠条形图全解析》在数据可视化的工具箱里,3D图表总能带来眼前一亮的效果,本文就来和大家聊聊如何使用Python实现绘制3D堆叠条形图,感兴趣的小伙伴可以了解下... 目录为什么选择 3D 堆叠条形图代码实现:从数据到 3D 世界的搭建核心代码逐行解析细节优化应用场景:3D 堆叠图

深度解析Spring AOP @Aspect 原理、实战与最佳实践教程

《深度解析SpringAOP@Aspect原理、实战与最佳实践教程》文章系统讲解了SpringAOP核心概念、实现方式及原理,涵盖横切关注点分离、代理机制(JDK/CGLIB)、切入点类型、性能... 目录1. @ASPect 核心概念1.1 AOP 编程范式1.2 @Aspect 关键特性2. 完整代码实

Java Web实现类似Excel表格锁定功能实战教程

《JavaWeb实现类似Excel表格锁定功能实战教程》本文将详细介绍通过创建特定div元素并利用CSS布局和JavaScript事件监听来实现类似Excel的锁定行和列效果的方法,感兴趣的朋友跟随... 目录1. 模拟Excel表格锁定功能2. 创建3个div元素实现表格锁定2.1 div元素布局设计2.

SpringBoot连接Redis集群教程

《SpringBoot连接Redis集群教程》:本文主要介绍SpringBoot连接Redis集群教程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1. 依赖2. 修改配置文件3. 创建RedisClusterConfig4. 测试总结1. 依赖 <de

Nexus安装和启动的实现教程

《Nexus安装和启动的实现教程》:本文主要介绍Nexus安装和启动的实现教程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、Nexus下载二、Nexus安装和启动三、关闭Nexus总结一、Nexus下载官方下载链接:DownloadWindows系统根

CnPlugin是PL/SQL Developer工具插件使用教程

《CnPlugin是PL/SQLDeveloper工具插件使用教程》:本文主要介绍CnPlugin是PL/SQLDeveloper工具插件使用教程,具有很好的参考价值,希望对大家有所帮助,如有错... 目录PL/SQL Developer工具插件使用安装拷贝文件配置总结PL/SQL Developer工具插

Java中的登录技术保姆级详细教程

《Java中的登录技术保姆级详细教程》:本文主要介绍Java中登录技术保姆级详细教程的相关资料,在Java中我们可以使用各种技术和框架来实现这些功能,文中通过代码介绍的非常详细,需要的朋友可以参考... 目录1.登录思路2.登录标记1.会话技术2.会话跟踪1.Cookie技术2.Session技术3.令牌技

Python使用Code2flow将代码转化为流程图的操作教程

《Python使用Code2flow将代码转化为流程图的操作教程》Code2flow是一款开源工具,能够将代码自动转换为流程图,该工具对于代码审查、调试和理解大型代码库非常有用,在这篇博客中,我们将深... 目录引言1nVflRA、为什么选择 Code2flow?2、安装 Code2flow3、基本功能演示

Java Spring 中的监听器Listener详解与实战教程

《JavaSpring中的监听器Listener详解与实战教程》Spring提供了多种监听器机制,可以用于监听应用生命周期、会话生命周期和请求处理过程中的事件,:本文主要介绍JavaSprin... 目录一、监听器的作用1.1 应用生命周期管理1.2 会话管理1.3 请求处理监控二、创建监听器2.1 Ser