跟着Nature Communications学作图：R语言ggtern包画三元相图（Ternary plots ）

跟着Nature Communications学作图：R语言ggtern包画三元相图（Ternary plots ）

论文

A highly conserved core bacterial microbiota with nitrogen-fixation capacity inhabits the xylem sap in maize plants

​​​​s41467-022-31113-w.pdf​​

数据代码链接

​​plots of OTUs in xylem sap across three fertilisation treatments. Size of each point represents relative abundance of OTU. Position is determined by the contribution of three fertilisation treatments to the total relative abundance, proximity to that vertex indicates enrichment of that OTU in that fertilisation treatment. Colours of circles correspond to different genera. Grey circles indicate OTUs with no significant differences in abundance.

作图数据部分截图

image.png

ggtern 这个R包是我第一次使用，首先需要安装，运行代码​​install.packages("ggtern")​​

读取数据集

plot_data <- read.table("data/20220612/Tern_data.txt", header=T, row.names= 1, sep="\t", comment.char = "")

作图代码

library(ggtern)p <- ggtern(data=plot_data, aes(x=CK, y=NPK, z=NPKM)) + geom_mask() + geom_point(aes(size=size, color=Genus),alpha=0.8) + scale_size(range = c(0, 10)) + scale_color_manual(values = c('#E31A1C','#228B22','#1F78B4', '#FDB462', '#8B658B', '#4876FF', '#00BFFF', '#EE82EE','#8B8682','#CDC9C9'), limits = c('Klebsiella','Pseudomonas','Enterobacteriaceae_unclassified','Rosenbergiella','Oxalobacteraceae_unclassified','Sphingobacterium','Lactococcus','Erwinia','Others','NotSig')) + guides(size="none") +theme_bw() + theme(axis.text=element_blank(), axis.ticks=element_blank())p

论文中的代码出图

image.png

调整图例的位置，可以用ggplot2中调整图例的方法

p+ theme(legend.position = "bottom")+ guides(color=guide_legend(nrow = 3, override.aes = list(size = 5), title.position = "top"))

image.png

试了一下和其他ggplot2的图进行组合，patchwork没有成功，找找有没有其他办法可以实现组合图

小明的数据分析笔记本

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。