跟着NC学作图：R语言ggplot2分组散点图和拟合曲线

用户7010445

发布于 2024-04-30 12:30:31

1420

发布于 2024-04-30 12:30:31

文章被收录于专栏：小明的数据分析笔记本小明的数据分析笔记本

论文

Genetic gains underpinning a little-known strawberry Green Revolution

https://www.nature.com/articles/s41467-024-46421-6

（这个论文还挺有意思的，抽时间仔细看看）

完整的作图代码

library(tidyverse)
library(readxl)
library(patchwork)


fig1.dat<-read_excel("2024.data/20240429/Figure 1_UN-FAO_Observed_Predicted_Values.xlsx")

fig1.dat %>% colnames()

p1<-ggplot(data=fig1.dat,aes(x=Year,y=`Yield (kg/ha)`))+
  geom_point(aes(color=Soure),size=5)+
  geom_smooth(aes(color=Soure),method = "lm")+
  scale_color_manual(values=c("United_States"="#1c86ee","Euorpe"="#7f7f7f"))+
  theme_bw(base_size = 15)+
  theme(panel.grid = element_blank())+
  scale_y_continuous(labels = function(x){x/1000})+
  labs(y="Yield\n()kg/ha x 1,000)")+
  scale_x_continuous(breaks = seq(1960,2020,10))


p2<-ggplot(data=fig1.dat,aes(x=Year,y=`Area Harvested (ha)`))+
  geom_point(aes(color=Soure),size=5)+
  geom_smooth(aes(color=Soure),method = "lm")+
  scale_color_manual(values=c("United_States"="#1c86ee","Euorpe"="#7f7f7f"))+
  theme_bw(base_size = 15)+
  theme(panel.grid = element_blank())+
  scale_y_continuous(labels = function(x){x/1000})+
  labs(y="Acreage\n()ha x 1,000)")+
  scale_x_continuous(breaks = seq(1960,2020,10))


p3<-ggplot(data=fig1.dat,aes(x=Year,y=`Production (kg)`))+
  geom_point(aes(color=Soure),size=5)+
  geom_smooth(aes(color=Soure),method = "lm")+
  scale_color_manual(values=c("United_States"="#1c86ee","Euorpe"="#7f7f7f"))+
  theme_bw(base_size = 15)+
  theme(panel.grid = element_blank())+
  scale_y_continuous(labels = function(x){x/1000000000})+
  labs(y="Production\n()Kg x 1,000,000,000)")+
  scale_x_continuous(breaks = seq(1960,2020,10))

p1+p2+p3+
  plot_layout(nrow=1,guides = "collect")&theme(legend.position = "bottom")