edit figures

analysis/paper/paper.Rmd

@@ -791,6 +791,35 @@ Our results show that in 2015, higher-income people in higher-income countries h
 
 ```{r}
 
+pdat_int_energy = get_sector_summary_by_eu_ntile_direct(eu_q_count) %>%
+  ungroup() %>%
+  filter(year==2015) %>%
+  left_join(read_csv(here("analysis/data/derived/sectors_agg_method1_ixi.csv")), 
+            by="sector_agg_id") %>% 
+  mutate(intensity_energy = (total_energy_use_tj)/(total_fd_me)) %>%
+  select(five_sectors, eu_q_rank, 
+         intensity_energy) %>%
+  filter(eu_q_rank == 10) %>%
+  select(-eu_q_rank)
+
+pdat_final_demand = get_sector_summary_by_eu_ntile_direct(eu_q_count) %>%
+  ungroup() %>%
+  filter(year==2015) %>%
+  left_join(read_csv(here("analysis/data/derived/sectors_agg_method1_ixi.csv")), 
+            by="sector_agg_id") %>%
+  left_join(pdat_int_energy, by="five_sectors") %>%
+  mutate(total_energy_use_tj_new = (total_fd_me)*intensity_energy) %>%
+  mutate(total_energy_use_tj_diff = total_energy_use_tj-total_energy_use_tj_new) %>%
+  select(eu_q_rank,total_energy_use_tj_new) %>%
+  group_by(eu_q_rank) %>%
+  summarise(total_energy_use_tj_new = sum(total_energy_use_tj_new)) %>%
+  mutate(total_blah = total_energy_use_tj_new/33417583,
+         total_blah = total_blah*1000)
+
+ineq_curr_check = pdat_final_demand$total_blah[10]/pdat_final_demand$total_blah[1]
+
+
+
 df_energy_deciles = read_csv(here("analysis/data/raw/energy_figure_data.csv")) 
 
 ineq_curr = df_energy_deciles$pae_energy_use_gj[10]/df_energy_deciles$pae_energy_use_gj[1]