跟着NatureCommunications学作图:R语言ggtree根据分组给进化树上色
语言 分组 根据 进化 跟着 作图 树上 ggtree
2023-06-13 09:16:32 时间
论文
MiDAS 4: A global catalogue of full-length 16S rRNA gene sequences and taxonomy for studies of bacterial communities in wastewater treatment plants
https://www.nature.com/articles/s41467-022-29438-7
数据链接
https://figshare.com/articles/dataset/Dueholm2021a_data_zip/16566408/1
代码链接
https://github.com/msdueholm/MiDAS4
今天的推文我们重复一下论文中的Figure2
image.png
我试了下论文中提供的数据,这个数据量有点大,运行要好长时间,这里的示例数据我们不用论文中的数据了,用自己随便构造的
最基本的进化树可视化
library(ggtree)
library(treeio)
tree <- read.tree("data/20220829/fig2.nwk")
ggtree(tree)+
geom_tiplab()
image.png
准备表示分组的数据,这里需要用到一个列表格式的数据
groupInfo<-list(group1=c("A","B","C","D","E"),
group2=c("F","G","H","I","J","K","L","M"))
然后使用groupOTU()
函数将分组信息和进化树组合到一起
tree01<-groupOTU(tree,groupInfo)
对结果进行展示
ggtree(tree01,aes(color=group))+
geom_tiplab(show.legend=F,offset = 1)+
geom_tippoint(aes(shape=group),
show.legend = F,
size=5)
image.png
这里 color=group
的group是默认生成的,不用改
更改配色
ggtree(tree01,aes(color=group),size=2)+
geom_tiplab(show.legend=F,offset = 1)+
geom_tippoint(aes(shape=group),
show.legend = F,
size=5)+
scale_color_manual(values = c("#fb8d63","#69c0a2"))
image.png
封面图
ggtree(tree01,
aes(color=group),
size=2,
layout = "circular",
branch.length = "none")+
geom_tiplab(show.legend=F,offset = 1)+
geom_tippoint(aes(shape=group),
show.legend = F,
size=5)+
scale_color_manual(values = c("#fb8d63","#69c0a2")) -> p1
p1
ggtree(tree01,aes(color=group),size=2)+
geom_tiplab(show.legend=F,offset = 1)+
geom_tippoint(aes(shape=group),
show.legend = F,
size=5)+
scale_color_manual(values = c("#fb8d63","#69c0a2")) -> p2
p2
library(patchwork)
p1+theme(legend.position = "none")+
p2+theme(legend.position = "none")
相关文章
- C语言冒泡法_冒泡编程c语言
- c语言之选择排序法
- 【R语言】因子在临床分组中的应用
- 【R语言】rep函数生成差异表达分组变量
- GoLand 2022 for Mac(GO语言集成开发工具环境) v2022.2.3中文激活版
- R语言缺失数据变量选择LASSO回归:Bootstrap重(再)抽样插补和推算
- 跟着Nature Plants学作图:R语言ggplot2画分组折线图并对坐标轴添加一些额外注释
- 跟着NaturePlants学作图:R语言ggplot2分组折线图完整示例
- 跟着Nature Communications学作图:R语言ggplot2绘制带有条纹的分组柱形图
- Meta打造首个「蛋白质宇宙」全景图!用150亿参数语言模型,预测了6亿+蛋白质结构
- ICLR 2023 | PromptPG:当强化学习遇见大规模语言模型
- WordPress多语言WPMLv4.6.3插件 自动翻译多国语言
- Go语言复数
- Go语言和其它编程语言的对比
- 国人开发的编程语言-“明”语言
- Go语言中转换JSON数据简单例子