查看原文
其他

如约 为大家带来了graphlan的物种分类树

文涛聊科研 微生信生物 2022-05-08


本次更新,我将整合两批代码并且即将与大家见面。这次我将物种分类树的笔记写给大家,方便理解。其实轮子我在进化树的时候就已经造的差不多了,本次重点解决其他问题,精修了轮子,期待我最后整合的《graphlan进化树和分类树一体化解决方案》

上次物种进化树可以在这里查看:

graphlan进一步学习笔记

物种分类树

写在前面

正如之前上一篇graphlan进化树所谈到的,graphlan可视化在本次物种树中又有着新的表现。其实我就是准备tree文件上下了些功夫,因为需要整合到phyloseq种,方便任何人在任何时候都可以随时重复。当然我做到了,phyloseq就像是一个核心,只要建立起来,我们就可以使用这个核心达到任何目的。其次设置环宇进化树之间的相对关系花费了一些时间,或许没有找到合适的方案吧。今天我需要记录的是一下几个点:

  1. graphlan物种树不同分类层次之间的分隔符号为英文状态下的“.”;

  2. 注释信息不可用特殊符号,例如:"()",如果出现特殊符号,相比连出图都不可能出现,并且没有任何报错信息,我想这个坑将会很难跨越。

  3. 不同于进化树的是物种分类树我设置了每一个分类分级的名称标签,标注在对应的环上,与外环的标签一同构成标签链。设置参数为:internal_label;

  4. 注意注释阴影部分默认颜色为灰色,注释信息默认排布为垂直环;

  5. 相对于进化树,物种分类树有些参数我们可以舍弃,例如:ignore_branch_len。

  6. 环宇进化树的相对关系是难点,这里参数如下设置达到内注释标签效果:annotation_background_separation -0.001 annotation_background_offset -0.001

外注释效果:annotation_background_separation -0.001 annotation_background_offset -0.4 

  1. 最后,其实我可以直接使用OTU来构建物种分类树,现在所看到的物种分类树只有六个层次,到达属等级,使用OTU等级都贱分类树可以达到八个分类层次,这是便可以直接使用OTU水平的分析展示在外环。方法都是一样的,大家有兴趣可以试试。


将分支宽度加粗 

这是我们原始尚未添加任何注释的的分类树


欢迎加入微生信生物讨论群,扫描下方二维码添加小编微信,小编带你入伙啦,大牛如云,让交流变得简单。

历史目录

R语言分析技术

  1. 《扩增子16s核心OTU挑选-基于otu_table的UpSet和韦恩图》

  2. 《分类堆叠柱状图顺序排列及其添加合适条块标签》

  3. 《R语言绘制带有显著性字母标记的柱状图》

  4. 《使用R实现批量方差分析(aov)和多重比较(LSD)》

  5. Rstudio切换挂载R版本及本地安装一些包

  6. 玩转R包

  7. 用ggplot画你们心中的女神

  8. ggplot版钢铁侠

  9. 想用ggplot做一张致谢ppt,但是碰到了520,画风就变了

扩增子专题

  1. 《16s分析之Qiime数据整理+基础多样性分析》

  2. 《16s分析之差异展示(热图)》

  3. 迅速提高序列拼接效率,得到后续分析友好型输入,依托qiime

  4. https://mp.weixin.qq.com/s/6zuB9JKYvDtlomtAlxSmGw》

  5. 16s分析之网络分析一(MENA)

  6. 16s分析之进化树+差异分析(一)

  7. 16s分析之进化树+差异分析(二)

  8. Qiime2学习笔记之Qiime2网站示例学习笔记

  9. PCA原理解读

  10. PCA实战

  11. 16s分析之LEfSe分析

基于phyloseq的微生物群落分析

  1. 开年工作第一天phyloseq介绍

  2. phyloseq入门

  3. R语言微生物群落数据分析:从原始数据到结果(No1)

  4. R语言微生物群落数据分析:从原始数据到结果(No2))

  5. phyloseq进化树可视化

  6. 基于phyloseq的排序分析

代谢组专题

  1. 非靶向代谢组学数据分析连载(第零篇引子)

  2. 非靶向代谢组学分析连载(第一篇:缺失数据处理、归一化、标准化)

当科研遇见python

1.当科研遇见python

科学知识图谱

1.citespace 的基本术语与下载安装

杂谈

  1. 我的生物信息之路


您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存