跟着Nature Genetics学作图：R语言ggplot2做簇状的堆积柱形图

library(readxl)
library(tidyverse)

dat<-read_excel("2024.data/20240109/41588_2023_1516_MOESM9_ESM.xlsx",
           sheet = "Sheet1",
           col_names = FALSE)
dat

colnames(dat)<-paste0("X",1:4)
dat %>% 
  mutate(X4=factor(X4,levels = c("ratio of down-relulated genes",
                                 "ratio of up-relulated genes",
                                 "ratio of None-differentially expressed genes"))) -> dat

library(ggplot2)
dat
ggplot(data=dat,aes(x=X1,y=X3))+
  geom_bar(aes(fill=X4),stat = "identity",
           position = "fill")+
  theme_bw(base_size = 15)+
  theme(panel.grid = element_blank(),
        panel.border = element_blank(),
        axis.line = element_line(),
        legend.position = "top",
        legend.title = element_blank())+
  guides(fill=guide_legend(ncol = 1))+
  scale_x_continuous(breaks = c(1,2,3,5,6,7,9,10,11,13,14,15),
                     labels = rep(c("D","I","N"),4))+
  scale_y_continuous(expand = expansion(mult = c(0,0)),
                     breaks = c(0,0.5,1),
                     labels = c(0,50,100))+
  labs(x=NULL,y="Percentage (%)")+
  scale_fill_manual(values = c("#f47720","#459943","#bdc3d2"))-> p1

p1
data.frame(x=c(2,6,10,14),
           y=1,
           label=c("Leaf","Fruit","Seed","Embryo")) %>% 
  ggplot(aes(x=x,y=y))+
  geom_text(aes(label=label),
            hjust=0.5)+
  theme_void()+
  scale_x_continuous(limits = c(1,15)) -> p2

library(patchwork)

p1+p2+
  plot_layout(ncol = 1,heights = c(10,1))