From 50bb6d42c7af8451cb642b68448dc756a93156dc Mon Sep 17 00:00:00 2001 From: Hel1vs Date: Thu, 14 Jan 2021 10:06:45 +0100 Subject: [PATCH] figure improvements upon reviewer feedback --- src/Neutrality/neutrality_paper.R | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/src/Neutrality/neutrality_paper.R b/src/Neutrality/neutrality_paper.R index e522312..88597f7 100644 --- a/src/Neutrality/neutrality_paper.R +++ b/src/Neutrality/neutrality_paper.R @@ -116,7 +116,7 @@ change[Category=="2deg"]$Category<-"2 °C" S = ggplot() S = S + geom_pointrange(data=poyrange[Category%in%c("2 °C","1.5 °C")&!region%in%c("SAF [2]","MEX [2]")], aes(ymin=min,ymax=showyear,y=median, x=region, colour=variable),alpha=0.5,size=5,fatten=1,show.legend = F) #,position=position_dodge(width=0.66) -S = S + geom_point(data=poy[Category%in%c("2 °C","1.5 °C")&!region%in%c("SAF [2]","MEX [2]")],aes(x=region,y=showyear,shape=model,colour=variable),size=3) #,position=position_dodge(width=0.66) +S = S + geom_point(data=poy[Category%in%c("2 °C","1.5 °C")&!region%in%c("SAF [2]","MEX [2]")],aes(x=region,y=showyear,shape=model,colour=variable),size=4) #,position=position_dodge(width=0.66) S = S + geom_point(data=change[change=="yes"],aes(x=region,y=2030),shape=5,size=2) S = S + guides(colour=F) S = S + geom_text(data=poy[Category%in%c("2 °C","1.5 °C")],stat="identity",aes(x=region,y=showyear,label=label),size=10) @@ -367,7 +367,7 @@ S2a = ggplot() S2a = S2a + geom_pointrange(data=poyrange1[Category%in%c("2 °C","1.5 °C")&!region%in%c("World [6 models]")&variable=="diff"], aes(ymin=min,ymax=max,y=median, x=region),alpha=0.5,size=5,fatten=1,show.legend = F,colour="#0080ff") S2a = S2a + geom_point(data=poy1[Category%in%c("2 °C","1.5 °C")&!region%in%c("World [6 models]")&variable=="diff"], - aes(x=region,y=value,shape=model),size=3) #,colour=model + aes(x=region,y=value,shape=model),size=4) #,colour=model S2a=S2a+geom_text(data=poyrange1[Category%in%c("2 °C")&!region%in%c("World [6 models]")&variable=="diff"],aes(x="USA [6 models]",y=-70),label ="a)",size=10) S2a = S2a + guides(colour=F) S2a = S2a + geom_hline(yintercept=0) @@ -523,7 +523,7 @@ a1a = ggplot() a1a = a1a + geom_pointrange(data=poyrange1[Category%in%c("2 °C","1.5 °C")&!region%in%c("World [5 models]")&variable=="diff"], aes(ymin=min,ymax=max,y=median, x=region),alpha=0.5,size=5,fatten=1,show.legend = F,colour="#fe868e") a1a = a1a + geom_point(data=poy1[Category%in%c("2 °C","1.5 °C")&!region%in%c("World [5 models]")&variable=="diff"], - aes(x=region,y=value,shape=model),size=3) #,colour=model + aes(x=region,y=value,shape=model),size=4) #,colour=model a1a=a1a+geom_text(data=poyrange1[Category%in%c("2 °C")&!region%in%c("World [5 models]")&variable=="diff"],aes(x="TUR [2 models]",y=-70),label ="b)",size=10) a1a = a1a + guides(colour=F) a1a = a1a + geom_hline(yintercept=0) @@ -612,7 +612,7 @@ g1 = ggplot() g1 = g1 + geom_pointrange(data=poyrange1[Category%in%c("2 °C","1.5 °C")&!region%in%c("World [6 models]","MEX [2 models]","SAF [2 models]")&gwp=="diff"], aes(ymin=min,ymax=max,y=median, x=region),alpha=0.5,size=5,fatten=1,show.legend = F,colour="#66b2ff") g1 = g1 + geom_point(data=poy1[Category%in%c("2 °C","1.5 °C")&!region%in%c("World [6 models]","MEX [2 models]","SAF [2 models]")&gwp=="diff"], - aes(x=region,y=value,shape=model),size=3) #,colour=model + aes(x=region,y=value,shape=model),size=4) #,colour=model g1=g1+geom_text(data=poyrange1[Category%in%c("2 °C")&!region%in%c("World [6 models]","MEX [2 models]","SAF [2 models]")&gwp=="diff"],aes(x="USA [6 models]",y=-70),label ="c)",size=10) g1 = g1 + guides(colour=F) g1 = g1 + geom_hline(yintercept=0) @@ -677,7 +677,7 @@ S2a=S2a+theme(legend.position = "none")+theme(axis.text.y=element_text(size=20), g1=g1+theme(legend.position = "none")+theme(axis.text.y=element_text(size=20),axis.text.x=element_text(size=22),strip.text=element_text(size=24)) a1a=a1a+theme(legend.position = "none")+theme(axis.text.y=element_text(size=20),axis.text.x=element_text(size=22),strip.text=element_text(size=24)) e1 = e1+theme(legend.position = "none")+theme(axis.text.y=element_text(size=20),axis.text.x=element_text(size=22),strip.text=element_text(size=24)) -text<-textGrob("Difference in phase-out year due to \n 1) inventory vs. model LULUCF (top-left): <0 earlier if based on inventory, \n 2) BECCS allocation (top-right): <0 earlier if based on biomass production, \n 3) GWP (bottom-left): <0 earlier if based on AR5 instead of AR4, \n 4) Equity (bottom-right): <0 earlier if based on equity range instead of cost-optimal", +text<-textGrob("Difference in phase-out year due to \n a) inventory vs. model LULUCF (top-left): <0 earlier if based on inventory, \n b) BECCS allocation (top-right): <0 earlier if based on biomass production, \n c) GWP (bottom-left): <0 earlier if based on AR5 instead of AR4, \n d) Equity (bottom-right): <0 earlier if based on equity range instead of cost-optimal", gp=gpar(fontsize=24)) S2a=S2a+scale_y_continuous(limits=c(-80,90),breaks=c(-80,-60,-40,-20,0,20,40,60,80)) #+ylab("") g1=g1+scale_y_continuous(limits=c(-80,90),breaks=c(-80,-60,-40,-20,0,20,40,60,80)) #+ylab("") @@ -687,6 +687,7 @@ g1=g1+geom_text(aes(x="BRA [3 models]",y=-40,label="Earlier"),stat="identity",da g1=g1+geom_text(aes(x="BRA [3 models]",y=40,label="Later"),stat="identity",data=poyrange1,size=10) S2a=S2a+geom_text(aes(x="BRA [3 models]",y=-40,label="Earlier"),stat="identity",data=poyrange1,size=10) S2a=S2a+geom_text(aes(x="BRA [3 models]",y=40,label="Later"),stat="identity",data=poyrange1,size=10) +#S2a=S2a+geom_text(aes(x="USA [6 models]",y=40,label="No difference between default and sensitivity"),stat="identity",data=poyrange1,size=4) # data=poyrange1[Category%in%c("2 °C")], a1a=a1a+geom_text(aes(x="USA [5 models]",y=-60,label="Earlier"),stat="identity",data=poyrange1,size=10) a1a=a1a+geom_text(aes(x="USA [5 models]",y=60,label="Later"),stat="identity",data=poyrange1,size=10) e1=e1+geom_text(aes(x="BRA [median]",y=-20,label="Earlier"),stat="identity",data=poyrange1,size=10)