diff --git a/analysis/preprocessing/full_code.Rmd b/analysis/preprocessing/full_code.Rmd
index ed175d2f12d96797a5a5020023f3368119ce4384..4ff4d0044b52ed345e9553a9566430e3113f6360 100644
--- a/analysis/preprocessing/full_code.Rmd
+++ b/analysis/preprocessing/full_code.Rmd
@@ -5617,15 +5617,15 @@ for (i in years_exb_ixi){
mutate(V1 = dplyr::recode(V1,"GR" = "EL","GB" = "UK"))
# TIVs
-
- # CO2 - combustion - air
-
+ ## CO2 - combustion - air
Exiobase_TIV_co2_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_co2_combustion_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_co2_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_combustion_air_", year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CO2_combustion_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_co2_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_combustion_air_", year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5648,14 +5648,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CO2_combustion_europe,TIV_CO2_combustion_not_europe)
- # CO2 - noncombustion - cement - air
-
+ ## CO2 - noncombustion - cement - air
Exiobase_TIV_co2_noncombustion_cement_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_co2_noncombustion_cement_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_co2_noncombustion_cement_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_noncombustion_cement_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CO2_noncombustion_cement_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_co2_noncombustion_cement_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_noncombustion_cement_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5678,14 +5679,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CO2_noncombustion_cement_europe,TIV_CO2_noncombustion_cement_not_europe)
- # CO2 - noncombustion - lime - air
-
+ ## CO2 - noncombustion - lime - air
Exiobase_TIV_co2_noncombustion_lime_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_co2_noncombustion_lime_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_co2_noncombustion_lime_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_noncombustion_lime_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CO2_noncombustion_lime_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_co2_noncombustion_lime_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_noncombustion_lime_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5708,14 +5710,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CO2_noncombustion_lime_europe,TIV_CO2_noncombustion_lime_not_europe)
- # CO2 - agriculture - peat decay - air
-
+ ## CO2 - agriculture - peat decay - air
Exiobase_TIV_co2_agriculture_peatdecay_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_co2_agriculture_peatdecay_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_co2_agriculture_peatdecay_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_agriculture_peatdecay_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CO2_agriculture_peatdecay_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_co2_agriculture_peatdecay_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_agriculture_peatdecay_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5738,14 +5741,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CO2_agriculture_peatdecay_europe,TIV_CO2_agriculture_peatdecay_not_europe)
- # CO2 - waste - biogenic - air
-
+ ## CO2 - waste - biogenic - air
Exiobase_TIV_co2_waste_biogenic_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_co2_biogenic_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_co2_waste_biogenic_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_biogenic_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CO2_waste_biogenic_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_co2_waste_biogenic_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_biogenic_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5768,14 +5772,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CO2_waste_biogenic_europe,TIV_CO2_waste_biogenic_not_europe)
- # CO2 - waste - fossil - air
-
+ ## CO2 - waste - fossil - air
Exiobase_TIV_co2_waste_fossil_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_co2_waste_fossil_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_co2_waste_fossil_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_waste_fossil_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CO2_waste_fossil_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_co2_waste_fossil_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_co2_waste_fossil_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5798,15 +5803,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CO2_waste_fossil_europe,TIV_CO2_waste_fossil_not_europe)
-
- # CH4 - combustion -air
-
+ ## CH4 - combustion -air
Exiobase_TIV_ch4_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_combustion_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_combustion_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_combustion_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_combustion_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5829,14 +5834,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_combustion_europe,TIV_CH4_combustion_not_europe)
- # CH4 - noncombustion - gas - air
-
+ ## CH4 - noncombustion - gas - air
Exiobase_TIV_ch4_noncombustion_gas_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_noncombustion_gas_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_noncombustion_gas_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_gas_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_noncombustion_gas_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_noncombustion_gas_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_gas_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5859,14 +5865,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_noncombustion_gas_europe,TIV_CH4_noncombustion_gas_not_europe)
- # CH4 - noncombustion - oil - air
-
+ ## CH4 - noncombustion - oil - air
Exiobase_TIV_ch4_noncombustion_oil_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_noncombustion_oil_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_noncombustion_oil_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_oil_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_noncombustion_oil_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_noncombustion_oil_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_oil_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5889,14 +5896,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_noncombustion_oil_europe,TIV_CH4_noncombustion_oil_not_europe)
- # CH4 - noncombustion - anthracite - air
-
+ ## CH4 - noncombustion - anthracite - air
Exiobase_TIV_ch4_noncombustion_anthracite_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_noncombustion_anthracite_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_noncombustion_anthracite_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_anthracite_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_noncombustion_anthracite_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_noncombustion_anthracite_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_anthracite_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5919,14 +5927,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_noncombustion_anthracite_europe,TIV_CH4_noncombustion_anthracite_not_europe)
- # CH4 - noncombustion - bituminouscoal - air
-
+ ## CH4 - noncombustion - bituminouscoal - air
Exiobase_TIV_ch4_noncombustion_bituminouscoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_noncombustion_bituminouscoal_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_noncombustion_bituminouscoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_bituminouscoal_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_noncombustion_bituminouscoal_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_noncombustion_bituminouscoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_bituminouscoal_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5949,14 +5958,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_noncombustion_bituminouscoal_europe,TIV_CH4_noncombustion_bituminouscoal_not_europe)
- # CH4 - noncombustion - cokingcoal - air
-
+ ## CH4 - noncombustion - cokingcoal - air
Exiobase_TIV_ch4_noncombustion_cokingcoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_noncombustion_cokingcoal_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_noncombustion_cokingcoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_cokingcoal_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_noncombustion_cokingcoal_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_noncombustion_cokingcoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_cokingcoal_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -5979,14 +5989,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_noncombustion_cokingcoal_europe,TIV_CH4_noncombustion_cokingcoal_not_europe)
- # CH4 - noncombustion - lignite - air
-
+ ## CH4 - noncombustion - lignite - air
Exiobase_TIV_ch4_noncombustion_lignite_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_noncombustion_lignite_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_noncombustion_lignite_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_lignite_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_noncombustion_lignite_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_noncombustion_lignite_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_lignite_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6009,14 +6020,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_noncombustion_lignite_europe,TIV_CH4_noncombustion_lignite_not_europe)
- # CH4 - noncombustion - subbituminouscoal - air
-
+ ## CH4 - noncombustion - subbituminouscoal - air
Exiobase_TIV_ch4_noncombustion_subbituminouscoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_noncombustion_subbituminouscoal_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_noncombustion_subbituminouscoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_subbituminouscoal_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_noncombustion_subbituminouscoal_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_noncombustion_subbituminouscoal_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_subbituminouscoal_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6039,14 +6051,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_noncombustion_subbituminouscoal_europe,TIV_CH4_noncombustion_subbituminouscoal_not_europe)
- # CH4 - noncombustion - oilrefinery - air
-
+ ## CH4 - noncombustion - oilrefinery - air
Exiobase_TIV_ch4_noncombustion_oilrefinery_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_noncombustion_oilrefinery_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_noncombustion_oilrefinery_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_oilrefinery_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_noncombustion_oilrefinery_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_noncombustion_oilrefinery_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_noncombustion_oilrefinery_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6069,14 +6082,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_noncombustion_oilrefinery_europe,TIV_CH4_noncombustion_oilrefinery_not_europe)
- # CH4 - agriculture - air
-
+ ## CH4 - agriculture - air
Exiobase_TIV_ch4_agriculture_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_agriculture_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_agriculture_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_agriculture_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_agriculture_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_agriculture_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_agriculture_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6099,14 +6113,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_agriculture_europe,TIV_CH4_agriculture_not_europe)
- # CH4 - waste - air
-
+ ## CH4 - waste - air
Exiobase_TIV_ch4_waste_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_ch4_CO2eq_waste_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_ch4_waste_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_waste_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_CH4_waste_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_ch4_waste_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_ch4_CO2eq_waste_air_", year_current, "_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6129,15 +6144,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_CH4_waste_europe,TIV_CH4_waste_not_europe)
-
- # N2O - combustion - air
-
+ ## N2O - combustion - air
Exiobase_TIV_n2o_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_n2o_CO2eq_combustion_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_n2o_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_n2o_CO2eq_combustion_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_N2O_combustion_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_n2o_combustion_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_n2o_CO2eq_combustion_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6160,14 +6175,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_N2O_combustion_europe,TIV_N2O_combustion_not_europe)
- # N2O - agriculture - air
-
+ ## N2O - agriculture - air
Exiobase_TIV_n2o_agriculture_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_n2o_CO2eq_agriculture_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_n2o_agriculture_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_n2o_CO2eq_agriculture_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_N2O_agriculture_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_n2o_agriculture_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_n2o_CO2eq_agriculture_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6190,14 +6206,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_N2O_agriculture_europe,TIV_N2O_agriculture_not_europe)
- # SF6 - air
-
+ ## SF6 - air
Exiobase_TIV_sf6_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_sf6_CO2eq_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_sf6_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_sf6_CO2eq_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_SF6_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_sf6_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_sf6_CO2eq_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6220,14 +6237,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_SF6_europe,TIV_SF6_not_europe)
- # HFC - air
-
+ ## HFC - air
Exiobase_TIV_hfc_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_hfc_CO2eq_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_hfc_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_hfc_CO2eq_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_HFC_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_hfc_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_hfc_CO2eq_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6250,14 +6268,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_HFC_europe,TIV_HFC_not_europe)
- # PFC - air
-
+ ## PFC - air
Exiobase_TIV_pfc_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_pfc_CO2eq_air_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_pfc_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_pfc_CO2eq_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_PFC_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_pfc_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_pfc_CO2eq_air_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%
@@ -6280,14 +6299,15 @@ for (i in years_exb_ixi){
ZA) %>%
select(TIV_PFC_europe,TIV_PFC_not_europe)
- # Energy use
-
+ ## energy use
Exiobase_TIV_energy_use_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_energy_carrier_use_",year_current,"_ixi.csv"))[,-1]
+ ### domestic part
Exiobase_TIV_country_breakdown_energy_use_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_energy_carrier_use_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
gather(country, TIV_energy_domestic)
+ ### non-domestic part
Exiobase_TIV_europe_breakdown_energy_use_bp = read.csv(paste0(data_dir_exiobase, "/IOT_",year_current,"_ixi/TIV_country_breakdown_energy_carrier_use_",year_current,"_ixi.csv"))[,-1] %>%
row_to_names(row_number = 1) %>%
mutate_at(vars(AT:ZA), funs(as.numeric(as.character(.)))) %>%