[nmicrobiol图表复现] ggtree绘制离散型环状热图

欢迎关注R语言数据分析指南

❝

本节来复现「nmicrobiol」上的一张图，此图有些小细节因此适用于「ggtree」绘制不需要其它扩展包；下面小编就通过一个详细的案例介绍如何绘制此图；「数据+代码已经上传VIP群，加群的观众老爷请自行下载」

❞

❝

有需要学习数据可视化的朋友欢迎加入小编2022年度VIP群，目前群内已经上传公众号文档「数据+代码约270篇」，VIP交流群(1)已经500人满员,随着内容不断增多，为了更好的创做内容现在进群需「付费149元」 双十一临近小编也做一波优惠活动 「从即日起至11月底无需转发文档到朋友圈积赞，可直接已129元的价格入群」

❞

淘宝店铺个性化绘图服务

❝

最近有部分观众老爷找小编做个性化绘图，为了给观众老爷提供更加优质的服务；小编开通了 「淘宝店铺(R语言数据分析指南)」 ，有特殊绘图需求的欢迎关注咨询 (「由于小编时间精力有限有需求的请及时下单」)

❞

参考文档

❝

https://yulab-smu.top/treedata-book/chapter7.html

❞

加载R包

library(tidyverse)
library(ggtree)
library(ggtreeExtra)
library(ggnewscale)

读入树文件

tree <- read.tree("RAxML_bipartitions.allUK_1000.nwk")

绘制进化树

cir <- ggtree(tree, layout="circular")+
  geom_tiplab(size=2,align=T,linesize=0,color="black",offset =0.4)

导入热图数据

热图中包含进化树名称，因此只需要将其转换为行名即可

df <- read_csv("metadata.csv") %>%
  column_to_rownames(var="label")

绘制heatmap

此处使用「color=NULL」来取消热图边框颜色,offset设置热图与进化树直接的距离

p1 <- gheatmap(cir,df %>% select(1),offset=0.6,width=.1,
               colnames_offset_y=0,colnames = F,color=NULL)+
  scale_fill_manual(values=c("#440154","#3b528b"))+
  new_scale_fill()

绘制外圈

gheatmap(p1,df %>% select(2),offset=0.696,width=.1,
           colnames_offset_y=0,colnames = F,color=NULL)+
  scale_fill_manual(values=c("#21918c","#5ec962","#fde725"))+
  theme(legend.title = element_blank(),
        legend.text=element_text(color="black"),
        legend.background = element_blank(),
        legend.key = element_blank(),
        legend.spacing.x = unit(0.1,'cm'),
        legend.key.width=unit(0.4,'cm'),
        legend.key.height=unit(0.4,'cm'))

❝

本节介绍到此结束，喜欢的观众老爷欢迎分享转发；「当然更推荐大家加入我的VIP交流群」扫描下方二维码加小编微信「添加小编微信请备注来意，以便高效处理」，无需求请勿扰

❞

小编微信

关注下方公众号下回更新不迷路