其他
跟着Microbiome学绘图(1) 个性化组合进化树与热图
欢迎关注R语言数据分析指南
❝今天来介绍如何对进化树枝进行折叠，将属于同一分支的折叠，从下面的论文中挑出一份数据通过一个小例子来展示
❞
❝近期系统整理了一下公众号所写过的一下经典文档，如果需要获取全部代码的欢迎加入小编的VIP交流群，「99元即可获取上百篇文档」的代码，公众号右下角添加微信咨询即可
❞
加载R包
library(tidyverse)
library(ggtree)
library(ggdendro)
library(patchwork)
加载数据
df <- read_csv("Figure5-metabolicfunctions_logtpm.csv") %>% dplyr::rename(tax=`...1`) %>%
column_to_rownames(var="tax")
绘制进化树
此处先选择需要的点，后对其进行折叠，具体参数请参考 ggtree官方文档，https://yulab-smu.top/treedata-book/chapter6.html
tree <- hclust(dist(df)) %>% ggtree() %>%
collapse(174, 'min', fill="#E41A1C") %>%
collapse(173, 'max', fill="#1E90FF") %>%
collapse(146, 'mix', fill="#FF8C00") %>%
collapse(141, 'mixed', fill="#4C005C") %>%
collapse(155, 'min', fill="#003380") %>%
collapse(154, 'min', fill="#FF8C00") %>%
collapse(155, 'min', fill="#FF8C00") %>%
collapse(152, 'mixed', fill="#740AFF") %>%
collapse(151, 'min', fill="#FF8C00") %>%
collapse(150, 'min', fill="#8F7C00") %>%
collapse(143, 'min', fill="#E41A1C") %>%
collapse(144, 'min', fill="#94FFB5") %>%
collapse(159, 'min', fill="#0075DC") %>%
collapse(158, 'min', fill="#00FFFF") %>%
collapse(156, 'min', fill="#E41A1C") %>%
collapse(148, 'min', fill="#8F7C00") %>%
collapse(140, 'min', fill="#E41A1C")
提取聚类标签
hrdata <- hclust(dist(df)) %>% dendro_data(.,type="rectangle")
hcdata <- hclust(dist(t(df))) %>% dendro_data(.,type="rectangle")
定义热图数据
df2 <- read_csv("Figure5-metabolicfunctions_logtpm.csv") %>% dplyr::rename(tax=`...1`) %>%
pivot_longer(-tax)
自定义因子
df2$tax <- factor(df2$tax,levels = hrdata$labels %>% select(label) %>% pull())
df2$name <- factor(df2$name,levels = hcdata$labels %>% select(label) %>% pull())
绘制热图
heatmap <- df2 %>% ggplot(aes(name,tax,fill=value))+
geom_tile()+
labs(x=NULL,y=NULL)+
scale_y_discrete(expand=c(0,0),position="left")+
scale_x_discrete(expand=c(0,0))+
scale_fill_gradientn(colours = rev(RColorBrewer::brewer.pal(11,"RdBu")))+
theme(axis.text.x=element_text(color="black",angle=90,vjust=0.5,size=8,hjust=1),
axis.text.y=element_blank(),
axis.ticks = element_blank(),
legend.title = element_blank(),
legend.text = element_text(size=8,color="black"))+
guides(fill=guide_colorbar(direction = "vertical",reverse = F,barwidth = unit(.6, "cm"),
barheight = unit(18.5,"cm")))
拼图
tree + heatmap+ plot_layout(widths = c(.3,1))、
数据代码获取
❝本节介绍到此结束，也许某一天你将用到这个小技能。如需要获取相关数据可以公众号右下角添加小编微信 加入小编的
❞2022年VIP交流群，付费99元，群内会同步上传公众号文档代码；目前已上传2021-2022约150篇文档代码；也可对该篇文档进行打赏获取数据，添加小编微信时，请备注单位-方向-姓名以便高效处理
免费交流群
欢迎大家扫描下方二维码加入「QQ交流群」,与全国各地上千位小伙伴交流