跟着Nature学作图：R语言ggplot2散点图和添加辅助线完整示例

library(tidyverse)
table(dat01$type)

dat01 %>% 
  filter(type=="cis_gwas") -> dat01.gwas

dim(dat01.gwas)

dat01 %>% 
  filter(type == "cis_local") -> dat01.local

dim(dat01.local)

left_join(dat01.gwas,dat01.local,by="ID") -> dat

library(ggplot2)
library(paletteer)
library(latex2exp)

help(package="latex2exp")

ggplot(data=dat,aes(x=h2.x,y=h2.y))+
  geom_point(aes(color=h2.x))+
  scale_color_paletteer_c(palette="ggthemes::Blue-Green Sequential",
                          direction = -1,
                          name="n neighbours")+
  theme_bw()+
  theme(panel.grid = element_blank(),
        legend.position = c(0.8,0.2))+
  geom_abline(slope = 1,intercept = 0,lty="dashed")+
  scale_x_continuous(limits = c(0,1),
                     breaks = seq(0,1,by=0.25))+
  scale_y_continuous(limits = c(0,1),
                     breaks = seq(0,1,by=0.25))+
  labs(x=TeX(r"(\textit{h}$^2$ {(\textit{cis})} { GWAS})"),
       y=TeX(r"(\textit{h}$^2$ {(\textit{cis})} { local})"))

library(patchwork)
p+
  scale_color_paletteer_c("ggthemes::Red-Green Diverging",
                          direction = -1) +
  p