在 R 语言中跟随 Nature Communications 学习绘图,涉及树状图、分支颜色调整、分组注释、热图注释与柱状堆积图注释的具体操作如下:
环境准备与数据读取:
清除全局环境,并设定工作路径为桌面。加载相关 R 包,包括 ggtree、ape、vegan、ggtreeExtra、ggnewscale 和 ggstar。读取数据文件,计算距离矩阵,并使用类平均法进行层次聚类。将聚类结果转换为系统发育格式,并将分组信息与树结构组合。
初步展示进化树:
使用 ggtree 对转换后的系统发育树进行初步展示。
整理注释数据:
整理分组注释信息,包括 Diet、Diet 、Habitat、Captivewild 和 Sample type 等。整理物种丰度与分类信息,以及处理重复数据。
添加注释信息:
自定义分支颜色,根据分组信息对分支进行着色。添加第一层至第六层的注释信息,覆盖上述分组注释内容。使用 ggtree 的相关函数,如 geom_tiplab、geom_text 或 annotate 等,添加注释。
绘制进化树并添加注释:
利用 ggtree 绘制最终的进化树。在进化树上添加详细的注释信息,包括分组注释、热图注释、柱状堆积图注释等。通过调整注释的位置、大小和颜色等参数,使图形更加美观和易于理解。
可视化与展示:
完成上述步骤后,使用 ggsave 函数将绘制好的图形保存为文件,如 PNG、PDF 或 SVG 格式。在 RStudio 或其他图形查看器中查看和展示绘制好的进化树及其注释信息。
注意:以上步骤为一般性指导,具体实现时可能需要根据数据特点和需求进行调整和优化。同时,热图注释和柱状堆积图注释的具体实现可能需要借助其他 R 包或函数,如 pheatmap、ggplot2 等。