diff --git a/analysis/figures/figure2.pdf b/analysis/figures/figure2.pdf
index c7fa31a7a91bf3e0dad00f74fe535248fb772b55..d9e2b9320a6386a809203586d5906acace4b191d 100644
Binary files a/analysis/figures/figure2.pdf and b/analysis/figures/figure2.pdf differ
diff --git a/analysis/figures/figure3.pdf b/analysis/figures/figure3.pdf
index 5cc97ae647bc283f765fef7bd8bea3a1f73a909d..b2110d5e59442ba6dd8699f9e4ec9d8e5082ac51 100644
Binary files a/analysis/figures/figure3.pdf and b/analysis/figures/figure3.pdf differ
diff --git a/analysis/figures/figure4.pdf b/analysis/figures/figure4.pdf
index d9d4b1c9d030ffc3770b1452e385a23490688cd7..2d40e37d1782830f1a596c91ca677d46ccfe617b 100644
Binary files a/analysis/figures/figure4.pdf and b/analysis/figures/figure4.pdf differ
diff --git a/analysis/figures/figure5.pdf b/analysis/figures/figure5.pdf
index 25a8298706fa2410b37fe480cd1544a2ee0b826a..0f90c3ef5a29d8a33beb2f05886fca17715758c6 100644
Binary files a/analysis/figures/figure5.pdf and b/analysis/figures/figure5.pdf differ
diff --git a/analysis/paper/paper.Rmd b/analysis/paper/paper.Rmd
index 8a78475ba784f26fda74ecacdabf5336d98d840d..a07009db904276c5460bb0f027fe6d83ed8a16c1 100644
--- a/analysis/paper/paper.Rmd
+++ b/analysis/paper/paper.Rmd
@@ -425,7 +425,7 @@ knitr::include_graphics(here::here("analysis", "figures", "figure1.pdf"))
 
 Figures 1d-e show that energy and carbon intensities of consumption are particularly high in the lower four deciles, while the higher deciles do not show large differences in weighted average energy and carbon intensity. The different intensities of household consumption across European expenditure deciles can be attributed to a combination of two plausible causes: first, if the composition of consumption baskets systematically differs according to the level of household expenditure. Second, if energy and carbon intensity within individual consumption sectors systematically differs according to the level of household expenditure.
 
-```{r , fig.width=8, fig.height=2}
+```{r , fig.width=8, fig.height=2.5}
 
 pdat_basket = get_sector_summary_by_eu_ntile_direct(eu_q_count) %>%
   ungroup() %>%
@@ -462,20 +462,35 @@ pdat = pdat_int_co2eq %>%
 
 pdat$facet = factor(pdat$indicator, levels = c("Expenditure share (%)", "Energy intensity (MJ/€)", "Carbon intensity (kgCO2eq/€)"))
 
-a = ggplot(pdat, aes(x=value, y=five_sectors, color=factor(eu_q_rank))) + #I added factor on eu_q_rank, to go back, remove factor call and remove 'discrete'
-  #geom_point(shape=18, alpha=0.75) +
+#a = ggplot(pdat, aes(x=value, y=five_sectors, color=eu_q_rank)) + 
+#  #geom_point(shape=18, alpha=0.75) +
+#  geom_text(label="I", alpha=0.75, fontface="bold") +
+#  scale_color_viridis(option = "A", end = 0.8,
+#                      direction = -1, 
+#                      breaks = c(1, 5, 10), 
+#                      name="European\ndecile") +
+#  facet_wrap(~facet, scales="free_x") +
+#  theme_minimal() +
+#  labs(x="",y="") +
+#  theme(text=element_text(size=13, family="Liberation Sans Narrow"),
+#        panel.spacing = unit(2.5, "lines")) 
+
+a = ggplot(pdat, aes(x=value, y=five_sectors, color=eu_q_rank)) +
   geom_text(label="I", alpha=0.75, fontface="bold") +
-  scale_color_viridis(option = "A", end = 0.8,
-                      direction = -1, 
-                      discrete = TRUE, # I added here, remove to go back to original
-                      breaks = c(1, 5, 10), 
-                      name="European\ndecile") +
-  facet_wrap(~facet, scales="free_x") +
+  scale_colour_viridis_b(option = "A", end = 0.8,
+                         direction = -1,
+                         breaks = c(2,3,4,5,6,7,8,9),
+                         #guide_legend(),
+                         guide = guide_colorsteps(show.limits=TRUE),
+                         name="European\ndecile") +
+  facet_wrap(~facet, scales="free_x",
+             strip.position="bottom") +
   theme_minimal() +
   labs(x="",y="") +
   theme(text=element_text(size=13, family="Liberation Sans Narrow"),
-        panel.spacing = unit(2.5, "lines")) 
-  
+        panel.spacing = unit(2.5, "lines"),
+        strip.placement="outside") 
+
 ggsave(here("analysis", "figures", "figure2.pdf"), device=cairo_pdf)
 ```
 
diff --git a/analysis/paper/paper.docx b/analysis/paper/paper.docx
index 1e98f6eae7d22a4bc64ac86ffa2457efca7827f2..65770c141803e9cffe5c16fb1421072565ba8c79 100644
Binary files a/analysis/paper/paper.docx and b/analysis/paper/paper.docx differ