R_scripts/Sensitivity_1B.R at master · ashreeta/R_scripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
library(plotly)
library(Hmisc)
library(reshape2)
library(ggplot2)
library(plyr)
library(wesanderson)
library(RColorBrewer)

Production_electricity_cost <- read_excel("~/WORK_PROJECTS/Iwb/IWB/Results_Final/R_data/Sensitivity/Sensitivity_1B/Production_electricity_cost.xlsx")
Production_electricity_carbon <- read_excel("~/WORK_PROJECTS/Iwb/IWB/Results_Final/R_data/Sensitivity/Sensitivity_1B/Production_electricity_carbon.xlsx")


data<- Production_electricity_cost
data<- Production_electricity_carbon

df <-melt(as.data.frame(data),c("Years"))

df$cat <- ''
df[df$variable == 'Demand + HPs + E-Mobility',]$cat <- "Demand"
df[df$variable != 'Demand + HPs + E-Mobility',]$cat <- "Production"


cols <- c("Grid Import" = "#4393C3", "CHP" = "#D6604D", "Hydroelectricity CH" = "#E0E0E0",
          "Hydroelectricity BS" = "#878787","PV"="#4D4D4D", "micro-CHP"="coral4", "Demand + HPs + E-Mobility"="peachpuff")

p4<-ggplot(data=df, aes(x = cat, y = value, fill=factor(variable, levels=c("Grid Import","Hydroelectricity CH",
                                                                  "Hydroelectricity BS","CHP","micro-CHP", "PV", "Demand + HPs + E-Mobility")),width = .8)) +
  geom_bar(stat = 'identity', position = 'stack') + facet_grid(~ Years)+scale_fill_manual(values=cols)+
  labs(x="Years",y="GWh/year") +
  theme_minimal() + theme(legend.title=element_blank())+ ggtitle("Electricity production")+theme(text = element_text(size = 12))
p4


#production only graph

h<- subset(df,df$variable %in% c("2014", "Demand + HPs + E-Mobility"))
g<- subset(df,df$variable %in% c("PV", "Hydroelectricity BS","Hydroelectricity CH","micro-CHP","Grid Import"))

p4 <- ggplot(data=df, aes(x=Years, y=value)) +
  geom_bar(data=df,stat="identity", aes(fill=factor(variable, levels=c("Grid Import","Hydroelectricity CH",
                                                                      "Hydroelectricity BS","CHP","micro-CHP", "PV", "Demand + HPs + E-Mobility"))),width = 8,
           position='stack') +
  scale_fill_manual(values=cols)+
  labs(x="Years",y="GWh/year") +  scale_x_continuous(breaks=c(2020, 2035, 2050)) +
  theme_minimal() + theme(legend.title=element_blank())+ ggtitle("Electricity production")+theme(text = element_text(size = 12))+facet_grid( ~ Years)
p4

# heat
Production_heat_cost <- read_excel("~/WORK_PROJECTS/Iwb/IWB/Results_Final/R_data/Sensitivity/Sensitivity_1B/Production_heat_cost.xlsx")


data<- Production_heat_cost
data<- Production_heat_carbon


df <-melt(as.data.frame(data),c("Years"))

df$cat <- ''
df[df$variable == 'Heat demand',]$cat <- "Demand"
df[df$variable != 'Heat demand',]$cat <- "Production"


cols <- c("Oil boilers" = "#A6CEE3", "Gas boilers" = "#1F78B4", "Wood boilers" = "#B2DF8A",
          "Waste boilers" = "#33A02C","HP"="#FB9A99","ST"="#E31A1C","DH gas"="#FDBF6F",
          "Gas HP"="#FF7F00","micro-CHP"= "#CAB2D6", "DH Wood" = "#fee1ba", "DH Waste"="#dbdb73", "DH HP"="#fc9d24","Heat demand"= "sienna1")

p4<-ggplot(data=df, aes(x = cat, y = value, fill=factor(variable, levels=c("ST","micro-CHP","Oil boilers","Gas boilers","Wood boilers","Waste boilers",
                                                                           "HP", "Gas HP","DH HP","DH gas",
                                                                           "DH Wood", "DH Waste", "Heat demand")),width = .8)) +
  geom_bar(stat = 'identity', position = 'stack') + facet_grid(~ Years)+scale_fill_manual(values=cols)+
  labs(x="Years",y="GWh/year") +
  theme_minimal() + theme(legend.title=element_blank())+ ggtitle("Heat production")+theme(text = element_text(size = 12))
p4


#prodcution graph only


df2 <-subset(df,df$variable %in% c("Months","Oil boilers","Gas boilers","Wood boilers","Waste boilers","HP", "ST",
                                   "Gas HP", "micro-CHP","DH gas", "DH wood", "DH Waste", "DH HP"))

p4 <- ggplot(data=df, aes(x=Years, y=value)) +
  geom_bar(data=df2,stat="identity", aes(fill=factor(variable, levels=c("ST","micro-CHP","Oil boilers","Gas boilers","Wood boilers","Waste boilers",
                                                                      "HP", "Gas HP","DH HP","DH gas",
                                                                      "DH wood", "DH Waste" ))),width = 8) +
  scale_fill_manual(values=cols)+
  labs(x="Years",y="GWh/year") +  scale_x_continuous(breaks=c(2020, 2035, 2050)) +
  theme_minimal() + theme(legend.title=element_blank())+ ggtitle("Heat production")+theme(text = element_text(size = 12))
p4