大数跨境
0
0

R语言ggtree按照指定的节点旋转树

R语言ggtree按照指定的节点旋转树 小明的数据分析笔记本
2021-01-12
0
导读:R语言里的ggtree这个包可视化进化树有一个默认的顺序,如果想要改变枝的相对位置应该如何实现呢?

R语言里的ggtree这个包可视化进化树有一个默认的顺序,如果想要改变枝的相对位置应该如何实现呢?通过查找ggtree作者写的帮助文档找到了对应的办法,可以使用rotate()函数

ggtree的帮助文档链接

http://yulab-smu.top/treedata-book/index.html

首先我们使用R语言内置的数据集mtcars做一个层次聚类
mtcars
df<-mtcars[,c(1,3:6)]
df
df.hclust<-hclust(dist(df))
接下来使用ggtree对层次聚类的结果进行展示
library(ggtree)
ggtree(df.hclust)+
  geom_tiplab(offset = 2)+
  xlim(NA,280)+
  geom_highlight(node = 34,fill="blue")
image.png

我们看到图上标记蓝色的一个分支默认是在最底下的,如果想要把这个分支放到顶上应该如何修改呢?可以直接用ggtree中的rotate()函数。rotate()接受两个参数,一个是需要旋转的节点。另外一个就是树

先通过geom_text()函数给每个节点添加上文字标签
ggtree(df.hclust)+
  geom_tiplab(offset = 2)+
  xlim(NA,280)+
  geom_text(aes(label=node))
image.png

通过上图我们看到需要旋转的是33节点

p1<-ggtree(df.hclust)+
  geom_tiplab(offset = 2)+
  xlim(NA,280)
  #geom_highlight(node = 34,fill="blue")
p2<-ggtree::rotate(p1,33)
p2+
  geom_hilight(node=34,fill="blue")
image.png

这样就把这一大块放到上面去了。

细心的读者可能发现了,这里在使用rotate()这个函数的时候写法是ggtree::rotate(p1,33),这样是为了使用指定包里的某个函数,因为R语言里的函数很多,有可能会重名,有时候你用到的函数可能并不是想实现功能的那个函数,所以比较保险的做法还是加上包的命名然后用两个冒号链接函数

还有一个知识点是如果想要这个树的开口朝下(现在是开口朝又),可以加
ggtree(df.hclust)+
  geom_tiplab(angle=90,offset = -70)+
  layout_dendrogram()
image.png

欢迎大家关注我的公众号

小明的数据分析笔记本

ggtree的其他推文

R语言读入比对好的fasta文件然后做NJ树并做boostrap检验


【声明】内容源于网络
0
0
小明的数据分析笔记本
分享R语言和python在生物信息领域做数据分析和数据可视化的简单小例子;偶尔会分享一些组学数据处理相关的内容
内容 971
粉丝 0
小明的数据分析笔记本 分享R语言和python在生物信息领域做数据分析和数据可视化的简单小例子;偶尔会分享一些组学数据处理相关的内容
总阅读437
粉丝0
内容971