# 效应图
rm(list=ls())
library(tidyverse)
# setwd("D:\\李远存文件（千万别打开！！！）\\课题-GBD相关\\APC模型\\全球2型糖尿病与糖网趋势分析\\APC结果绘图用")  # 已注释：请按本地环境设置工作目录
DR <- read.csv("DR.csv",header = T)
DM <- read.csv("DM.csv",header = T)
DR$group <- factor(DR$group,
                   levels = c("age","period","cohort"),
                   labels = c("age","period","cohort"))
DM$group <- factor(DM$group,
                   levels = c("age","period","cohort"),
                   labels = c("age","period","cohort"))
str(DR)
x <- c(1:21)
DM$x <- x
DR$x <- x
fig1 <- ggplot(data=DR)+
  geom_ribbon(aes(x=x,ymin=lower,ymax=upper,fill=group),alpha=0.6)+
  geom_point(aes(x=x,y=coef,color=group))+
  geom_line(aes(x=x,y=coef,color=group),size=1)+
  scale_x_continuous(breaks=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21),
                    labels=c(DR$variable))+
  scale_y_continuous(limits = c(-2,2))+
  scale_fill_discrete(name="Effect")+
  scale_color_discrete(name="Effect")+
  facet_grid(~group,scales = "free_x", space = "free_x")+
  xlab("")+
  ylab("Coefficients for APC model of T2DR")+
  theme(axis.text.x = element_text(angle = 45,vjust = 1.0,hjust = 1.0))
fig1

fig2 <- ggplot(data=DM)+
  geom_ribbon(aes(x=x,ymin=lower,ymax=upper,fill=group),alpha=0.6)+
  geom_point(aes(x=x,y=coef,color=group))+
  geom_line(aes(x=x,y=coef,color=group),size=1)+
  scale_x_continuous(breaks=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21),
                     labels=c(DR$variable))+
  scale_y_continuous(limits = c(-2,2))+
  scale_fill_discrete(name="Effect")+
  scale_color_discrete(name="Effect")+
  facet_grid(~group,scales = "free_x", space = "free_x")+
  xlab("")+
  ylab("Coefficients for APC model of T2DM")+
  theme(axis.text.x = element_text(angle = 45,vjust = 1.0,hjust = 1.0))
fig2

fig2/fig1+plot_annotation(tag_levels = "A")