NC图表复现-时序分析预测图

library(tidyverse)
library(ggtext)

data <- read_tsv("data2.txt")

df <- read_tsv("data1.txt") %>% mutate(full_upper=GPPFULL+full_se,
                                       full_lower=GPPFULL-full_se,
                                       bgc_upper=GPPBGC+bgc_se,
                                       bgc_lower=GPPBGC-bgc_se)

ggplot()+
  annotate(geom="rect",xmin=2011,xmax=Inf,ymin=-Inf,ymax=0,alpha=.2)+
  geom_line(data=df %>% select(1,2,4) %>% pivot_longer(-id),
            aes(x = id, y = value,color=name),size=1)+
  geom_ribbon(data=df %>% select(1,2,6,7),
              aes(x=id,y=GPPFULL,ymin = full_upper, ymax = full_lower),
              fill = "#0B7751", alpha = 0.5)+
  geom_ribbon(data=df %>% select(1,4,8,9),
              aes(x=id,y=GPPBGC,ymin = bgc_upper, ymax = bgc_lower),
              fill = "#7294D4", alpha = 0.5)+
  geom_hline(yintercept = 100,linetype="dashed",color="#FF7F00",size=1)+
  geom_errorbar(data=data,aes(MedianYear,mean*500,
                              ymin=(mean-se)*500,
                              ymax=(mean+se)*500),
                width=1.8,size=1,color="#FF7F00")+
  geom_point(data=data,aes(MedianYear,mean*500),color="#FF7F00",
             shape=23,fill="#FDD262",size=4)+
  scale_color_manual(values = c(GPPBGC="#7294D4",GPPFULL="#0B7751"))+
  scale_x_continuous(expa=c(0,0))+
  theme_test() +
  theme(axis.title.y.right = element_markdown(size=10),
        axis.title.y.left = element_markdown(size=10,color="#FF7F00"),
        axis.text.x=element_text(color="black"),
        axis.text.y.left = element_text(color="#FF7F00"),
        plot.margin = unit(c(1,1,1,1),unit="cm"),
        legend.title = element_blank(),
        legend.background = element_blank(),
        legend.key = element_blank(),
        legend.position =c(0.5,0.9))