You need to sign in or sign up before continuing.
Newer
Older
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
setwd("D:/public/yang/CAPRIR_Project/caprir")
map_data <- caprir::prep_mapdata()
map_data_use <- map_data %>%
filter(!grepl("TR.*", CAPRI_NUTS_ID)) %>% # TR: Tuerkiye
filter(!grepl("UK.*", CAPRI_NUTS_ID)) %>%# UK:
filter(!grepl("PT20.*", CAPRI_NUTS_ID)) %>% # PT20/PT30
filter(!grepl("PT30.*", CAPRI_NUTS_ID)) %>%# PT20/PT30
filter(!grepl("ES70.*", CAPRI_NUTS_ID)) # ES70
setwd("D:/public/yang/CAPRIR_Project/capriv")
file1 <- "res_2_1230cap_after_2014_refdefaultA" # refdefautA
file2 <- "res_2_1230fta_import_bandefaultA" # bandefaultA
gamsPath <- "D://gams//win64//24.9"
igdx(gamsPath)
# CERE RAPE SOYA PULS
prod1 = c("CERE","RAPE","SOYA","PULS")
# MAIF OFAR GRAI GRAE
prod2 = c("MAIF", "OFAR","GRAI","GRAE")
# products %>% filter(key %in% prod_list)
outdata.dir <- paste0(getwd(), "/inst/outdata")
map_prod_data <- function(benchmark,scenario,prod, saveName){
benchmark <- capri_data(filename = paste0(gdx.dir,"/",benchmark),
selregion = "all",
seldim5 = "",
selcols = prod,
selrows = "LEVL" ,
simyear = "2030",
scenarioname = "benchmark")
scenario <- capri_data(filename = paste0(gdx.dir,"/",scenario),
selregion = "all",
seldim5 = "",
selcols = prod,
selrows = "LEVL" ,
simyear = "2030",
scenarioname = "scenario")
save(benchmark, scenario, file = paste0(outdata.dir,"/",saveName, ".rda"))
# map_capri(baseline = benchmark, scenario = scenario, comparison = FALSE)
# map_capri(baseline = benchmark, scenario = scenario, comparison = TRUE, percent_change = FALSE)
# map_capri(baseline = benchmark, scenario = scenario, comparison = TRUE, percent_change = TRUE)
}
map_prod_data(benchmark = file1, scenario = file2,prod = prod1,saveName = "prod1")
map_prod_data(benchmark = file1, scenario = file2,prod = prod2,saveName = "prod2")
t1 = Sys.time()
rm(benchmark, scenario)
load(file = paste0(outdata.dir,"/","prod1", ".rda"))
p1 <- map_capri(baseline = benchmark, scenario = scenario, comparison = TRUE, percent_change = TRUE, 11)
p1
ggsave(paste0( outdata.dir,"/png/cereals.png"), plot = p1,width = 16, height = 9, dpi = 100)
print(Sys.time()-t1)
rm(benchmark, scenario)
load(file = paste0(outdata.dir,"/","prod2", ".rda"))
p2 <- map_capri(baseline = benchmark, scenario = scenario, comparison = TRUE, percent_change = FALSE)
p2
ggsave(paste0( outdata.dir,"/png/MAIF.png"), plot = p2,width = 16, height = 9, dpi = 100)