install.packages( "ggplot2") library(ggplot2) setwd('C:/Users/r07am15.UOA/Documents/My papers/SSP comparison/') #d data<-read.csv('R_carbonloss.csv', na.strings="") str(data) ggplot(data=data, aes(x=Model, y=Loss, fill=Carbon)) + geom_bar(stat="identity") data <- with(data, data[order(Region, SSP, Loss,]) ggplot(data=data, aes(x=Model, y=Loss, fill=Carbon)) + geom_bar(stat="identity") + labs(y="Carbon loss (tonnes)") + facet_grid(~SSP) +ggtitle("d) Model comparison") + theme_bw() + theme(legend.direction = "horizontal", panel.grid.minor = element_blank(), axis.title.x = element_blank(), panel.grid.major=element_blank(), panel.grid.major.y = element_line( size=0.5, color="grey" )) ggplot(data=data, aes(x=Model, y=Loss, fill=Carbon)) + geom_bar(stat="identity") + labs(y="Carbon loss (tonnes)") + facet_grid(~SSP) +ggtitle("(d) Global comparison") + theme_bw() + theme(panel.grid.minor = element_blank(), axis.title.x = element_blank(), panel.grid.major=element_blank(), panel.grid.major.y = element_line( size=0.5, color="grey" )) #a) par(mfrow=c(2,2)) data1<-read.csv('R_carbonloss_SSP1.csv', na.strings="") str(data1) par(mfrow=c(2,2)) data1 <- with(data1, data1[order(Region, SSP, Loss,]) ggplot(data=data1, aes(x=Model, y=Loss, fill=Carbon)) + geom_bar(stat="identity") + labs(y="Carbon loss (tonnes)") + ylim(0, 2.0e+10) + facet_grid(~Region) +ggtitle("(a) SSP1")+ theme_bw() + theme(panel.grid.minor = element_blank(), axis.title.x = element_blank(), panel.grid.major=element_blank(), panel.grid.major.y = element_line( size=0.5, color="grey" )) #B) data2<-read.csv('R_carbonloss_SSP2.csv', na.strings="") str(data2) data2 <- with(data2, data2[order(Region, SSP, Loss,]) ggplot(data=data2, aes(x=Model, y=Loss, fill=Carbon)) + geom_bar(stat="identity") + labs(y="Carbon loss (tonnes)") + ylim(0, 2.0e+10) + facet_grid(~Region) +ggtitle("(b) SSP2")+ theme_bw() + theme(panel.grid.minor = element_blank(), axis.title.x = element_blank(), panel.grid.major=element_blank(), panel.grid.major.y = element_line( size=0.5, color="grey" )) #C) data3<-read.csv('R_carbonloss_SSP3.csv', na.strings="") str(data2) data3 <- with(data3, data3[order(Region, SSP, Loss,]) ggplot(data=data3, aes(x=Model, y=Loss, fill=Carbon)) + geom_bar(stat="identity") + labs(y="Carbon loss (tonnes)") + facet_grid(~Region) +ggtitle("(c) SSP3")+ theme_bw() + theme(panel.grid.minor = element_blank(), axis.title.x = element_blank(), panel.grid.major=element_blank(), panel.grid.major.y = element_line( size=0.5, color="grey" ))