R语言ggplot2画分组堆积柱形图展示密码子偏向性的RSCU值

#rm(list=ls())
library(readr)
library(dplyr)
library(ggplot2)
library(ggstar)
library(aplot)
help(package="ggstar")
df<-read.csv("RSCU/aa.csv",header=F)
head(df)
# df %>% 
#   arrange(V2,V4,V1) -> df
# dim(df)
read.csv("RSCU/CBS1_rscu.txt",
               sep="\t",
               na.strings = "*",
         header = F) %>% 
  select(c("V2","V3","V5","V6","V7")) -> cbs1

cbs1 %>% 
  group_by(V7) %>% 
  summarise(y=sum(V5)) -> cbs1.1

read.csv("RSCU/DF2_rscu.txt",
         sep="\t",
         na.strings = "*",
         header = F) %>% 
  select(c("V2","V3","V5","V6","V7")) -> df2

df2 %>% 
  group_by(V7) %>% 
  summarise(y=sum(V5)) -> df2.1


read.csv("RSCU/GS1_rscu.txt",
         sep="\t",
         na.strings = "*",
         header = F) %>% 
  select(c("V2","V3","V5","V6","V7")) -> gs1

gs1 %>% 
  group_by(V7) %>% 
  summarise(y=sum(V5)) -> gs1.1


read.csv("RSCU/XZXS1_rscu.txt",
         sep="\t",
         na.strings = "*",
         header = F) %>% 
  select(c("V2","V3","V5","V6","V7")) -> xzxs1

xzxs1 %>% 
  group_by(V7) %>% 
  summarise(y=sum(V5)) -> xzxs1.1

p1<-ggplot()+
  geom_bar(data=cbs1,
           aes(x=V7,
               y=V5,
               fill=as.character(V6)),
           position = "stack",
           stat = "identity",
           width = 0.2,
           color="black")+
  geom_star(data=cbs1.1,aes(x=V7,y=y+0.2),
            starshape=11,size=3,
            fill="#ff61cc",color="#ff61cc")+
  geom_bar(data=df2,
           aes(x=V7+0.23,
               y=V5,
               fill=as.character(V6)),
           position = "stack",
           stat = "identity",
           width = 0.2,
           color="black")+
  geom_star(data=df2.1,aes(x=V7+0.23,y=y+0.2),
            starshape=12,size=3,
            fill="#ff61cc",color="#ff61cc")+
  geom_bar(data=gs1,
           aes(x=V7+0.47,
               y=V5,
               fill=as.character(V6)),
           position = "stack",
           stat = "identity",
           width = 0.2,
           color="black")+
  geom_star(data=gs1.1,aes(x=V7+0.47,y=y+0.2),
            starshape=13,size=3,
            fill="#ff61cc",color="#ff61cc")+
  geom_bar(data=xzxs1,
           aes(x=V7+0.71,
               y=V5,
               fill=as.character(V6)),
           position = "stack",
           stat = "identity",
           width = 0.2,
           color="black")+
  geom_star(data=xzxs1.1,aes(x=V7+0.71,y=y+0.2),
            starshape=14,size=3,
            fill="#ff61cc",color="#ff61cc")+
  theme_bw()+
  scale_y_continuous(expand=c(0,0),
                     limits = c(0,8.5))+
  scale_x_continuous(expand = c(0,0),
                     limits = c(0.5,21),
                     breaks = seq(1.35,20.35,1),
                     labels = df$V1)+
  labs(y="RSCU",x="")+
  theme(legend.position = "none",
        panel.grid = element_blank())+
  geom_star(aes(x=1.2,y=8),
            starshape=11,
            size=3,
            fill="#ff61cc",
            color="#ff61cc")+
  geom_text(aes(x=1.3,y=8),label="CBS1",hjust=0)+
  geom_star(aes(x=1.2,y=7.7),
            starshape=12,
            size=3,
            fill="#ff61cc",
            color="#ff61cc")+
  geom_text(aes(x=1.3,y=7.7),label="DF2",hjust=0)+
  geom_star(aes(x=1.2,y=7.4),
            starshape=13,
            size=3,
            fill="#ff61cc",
            color="#ff61cc")+
  geom_text(aes(x=1.3,y=7.4),label="GS1",hjust=0)+
  geom_star(aes(x=1.2,y=7.1),
            starshape=14,
            size=3,
            fill="#ff61cc",
            color="#ff61cc")+
  geom_text(aes(x=1.3,y=7.1),label="XZXS1",hjust=0)
ggsave(filename = "RSCU/p1-2.pdf",
       p1,
       width = 15,
       height = 8)

df2
p2<-ggplot(df2,aes(x=V7+0.35,y=V6))+
  geom_label(aes(label=V3,fill=as.character(V6)),
             size=2)+
  labs(x="",y="")+ylim(2.4,6.1)+
  theme_minimal()+
  theme(legend.position = "none",
        axis.text = element_blank(),
        axis.ticks = element_blank(),
        panel.grid = element_blank())
p3<-p1%>%
  insert_bottom(p2,height = 0.3)

ggsave(filename = "RSCU/p3-1.pdf",
       p3,
       width = 20,
       height = 8)