跟着Nature Plants学作图：R语言ggplot2画分组折线图和置信区间
论文
The flying spider-monkey tree fern genome provides insights into fern evolution and arborescence
https://www.nature.com/articles/s41477-022-01146-6#Sec44
数据下载链接
https://doi.org/10.6084/m9.figshare.19125641
今天的推文重复一下论文中的Figure1d中左下角的小图
image.png
论文中提供的原始数据集如下
image.png
需要将其整理成3个单独的数据集
image.png
首先是做数据整理的代码
library(readxl) library(ggplot2) library(tidyverse) dat01<-read_excel("data/20220518/dat01.xlsx") dat01 dat01 %>% mutate(x=8:2, mean_value=`Duplication Per`, lower95=`Duplication Per`, upper95=`Duplication Per`, group="Experimental") %>% select(x,mean_value,lower95,upper95,group)-> dat01 dat02<-read_excel("data/20220518/dat02.xlsx") head(dat02) dat02 %>% group_by(MRCA) %>% summarise(mean_value = mean(`Duplication Per`), lower95 = Rmisc::CI(`Duplication Per`)[3], upper95 = Rmisc::CI(`Duplication Per`)[1]) %>% mutate(x=8:2, group="Null simulation") %>% select(x,mean_value,lower95,upper95,group) -> dat02 dat03<-read_excel("data/20220518/dat03.xlsx") head(dat03) dat03 %>% group_by(MRCA) %>% summarise(mean_value = mean(`Duplication Per`), lower95 = Rmisc::CI(`Duplication Per`)[3], upper95 = Rmisc::CI(`Duplication Per`)[1]) %>% mutate(x=8:2, group="Positive simulation") %>% select(x,mean_value,lower95,upper95,group) -> dat03 dat<-bind_rows(dat01,dat02,dat03) dat %>% mutate(group=fct_rev(group)) -> dat作图代码
ggplot(data = dat,aes(x=x,y=mean_value))+ geom_line(aes(color=group),size=1.5)+ geom_ribbon(aes(ymin=lower95-1,ymax=upper95+1, fill=group), alpha=0.5, show.legend = FALSE)+ theme_minimal()+ theme(panel.grid = element_blank(), axis.line = element_line(), axis.ticks = element_line(), legend.title = element_blank())+ scale_x_continuous(breaks = 2:8, labels = paste0("N",8:2))+ scale_y_continuous(breaks = c(0,25))+ labs(x="Subtree node",y="Percent subtree")最终结果
image.png
这里最终的结果和论文中的图还是有些差异的，因为我没有看到论文中用的是置信区间作为数据范围还是其他，我这里选择的是置信区间，然后对数值进行了加减1
示例数据可以到论文中去下载，代码可以在推文中复制
AI大模型学习福利
作为一名热心肠的互联网老兵，我决定把宝贵的AI知识分享给大家。 至于能学习到多少就看你的学习毅力和能力了 。我已将重要的AI大模型资料包括AI大模型入门学习思维导图、精品AI大模型学习书籍手册、视频教程、实战学习等录播视频免费分享出来。
一、全套AGI大模型学习路线
AI大模型时代的学习之旅：从基础到前沿，掌握人工智能的核心技能！
因篇幅有限，仅展示部分资料，需要点击文章最下方名片即可前往获取
二、640套AI大模型报告合集
这套包含640份报告的合集，涵盖了AI大模型的理论研究、技术实现、行业应用等多个方面。无论您是科研人员、工程师，还是对AI大模型感兴趣的爱好者，这套报告合集都将为您提供宝贵的信息和启示。
因篇幅有限，仅展示部分资料，需要点击文章最下方名片即可前往获
三、AI大模型经典PDF籍
随着人工智能技术的飞速发展，AI大模型已经成为了当今科技领域的一大热点。这些大型预训练模型，如GPT-3、BERT、XLNet等，以其强大的语言理解和生成能力，正在改变我们对人工智能的认识。 那以下这些PDF籍就是非常不错的学习资源。
因篇幅有限，仅展示部分资料，需要点击文章最下方名片即可前往获
四、AI大模型商业化落地方案
因篇幅有限，仅展示部分资料，需要点击文章最下方名片即可前往获
作为普通人，入局大模型时代需要持续学习和实践，不断提高自己的技能和认知水平，同时也需要有责任感和伦理意识，为人工智能的健康发展贡献力量