

❝给予长期支持我们的读者们一个特别待遇:购买小编VIP会员文档的读者,「将自动获得2024年及以后更新的绘图文档代码,无需额外付费」。目前会员文档(2023+2024)「已经更新上传了150+案例文档」,每个案例都附有相应的数据和代码,并配有对应的注释文档,方便大家学习和参考。 ❞
library(tidyverse)
library(magrittr)
library(ggtext)
library(grid)
library(cowplot)
col <- c("#ABDDDE","#FAD510","#C6CDF7","#F4B5BD","#0A9F9D","#FAEED1","#005295","#E6A0C4","#C52E19",
"#9C9D9D")
df <- read_tsv("data.txt") %>%
pivot_longer(-c("pathway","group")) %>%
set_colnames(c("pathway","group","Species<br>(contribution percentage>10%)","value")) %>%
mutate(value=case_when(group =="Enriched in stable low DDS" ~ -value,
TRUE ~ value))
df$pathway <- factor(df$pathway,levels = df$pathway %>% unique() %>% rev())
plot <- df %>% dplyr::rename("Percentage of species contributing<br>to each pathway,%"="value") %>%
ggplot(aes(`Percentage of species contributing<br>to each pathway,%`,
pathway,
fill=`Species<br>(contribution percentage>10%)`))+
geom_col(color="black",size=0.5)+
scale_fill_manual(values = c("#ABDDDE","#FAD510","#C6CDF7","#F4B5BD","#0A9F9D","#FAEED1","#005295","#E6A0C4","#C52E19",
"#9C9D9D"))+
coord_cartesian(clip = "off")+
theme_bw()+
theme(legend.title = element_markdown(color="black",face="bold"),
legend.text = element_text(color="black",face= "bold.italic"),
axis.ticks = element_blank(),
axis.text=element_text(color="black",face="bold"),
axis.title.x = element_markdown(color="black",face="bold"),
axis.title.y = element_blank(),
plot.margin=unit(c(2,0.3,0.3,0.3),unit="cm"))