跟着Nature学作图：R语言ggplot2箱线图/散点图完整示例

library(readxl)
dat01<-read_excel("data/20220711/41586_2022_4808_MOESM8_ESM.xlsx",
                  sheet = "Fig4e",
                  skip = 1)
dat01

library(ggplot2)

dat01$Type<-factor(dat01$Type,
                   levels = c("SNP","InDel","SV","SV array"))

ggplot(data=dat01,aes(x=Type,y=Accuracy))+
  
  geom_line(aes(group=Metabolic),
            color="grey")+
  geom_boxplot(aes(color=Type),
               fill="transparent")+
  geom_point(aes(color=Type))+
  theme_bw()+
  theme(panel.grid = element_blank(),
        legend.position = "none")+
  scale_color_manual(values = c("#3288bd","#66c2a5",
                                "#f46d43","#f1226e"))+
  labs(x=NULL,y="Prediction accuracy") -> p1

p1

dat02<-read_excel("data/20220711/41586_2022_4808_MOESM8_ESM.xlsx",
                  sheet = "Fig4d",
                  skip = 1)
dat02

library(paletteer)
library(latex2exp)

ggplot(data=dat02,aes(x=accuracy_snp,y=accuracy_sv))+
  geom_point(aes(color=h2_snp))+
  scale_x_continuous(limits = c(0,0.5),
                     expand = expansion(mult = c(0,0)))+
  scale_y_continuous(limits = c(0,0.5),
                     expand = expansion(mult = c(0,0)))+
  geom_abline(slope = 1,intercept = 0,lty="dashed")+
  theme_bw()+
  theme(panel.grid = element_blank(),
        panel.border = element_blank(),
        axis.line = element_line(),
        legend.position = c(0.8,0.3))+
  scale_color_paletteer_c(palette="grDevices::heat.colors",
                          direction = -1,
                          breaks=c(0.01,0.5,0.99),
                          labels=c(0,"0.50","1.00"),
                          name=TeX(r"(\textit{h}${^2}$)"))+
  labs(x="Accuracy of SNP",y="Accuracy of SV") -> p2

p2

library(patchwork)

p2+p1