# tolerance trends, estimate plus raw data, selected countries
countries <- c("United States", "Ireland", "Sweden", "Russia", "Japan",
"Uruguay", "Spain", "South Africa", "India", "Australia")
c8_rd <- gm %>%
filter(country %in% countries) %>%
mutate(estimate = y_r/n,
country = factor(country, levels = countries)) %>%
left_join(g_res, by = c("rcode", "variable_cp")) %>%
rename(gamma_mean = mean)
c8_res <- a_res %>%
filter(country %in% countries) %>%
mutate(country = factor(country, levels = countries))
p_rd <- ggplot(data = c8_rd, aes(x = year, y = estimate)) +
geom_line(aes(color = variable_cp, alpha = 1/gamma_mean)) +
geom_point(aes(color = variable_cp, alpha = 1/gamma_mean)) + theme_bw() +
theme(legend.position="none") +
coord_cartesian(xlim = c(1975, 2016), ylim = c(0, 1)) +
labs(x = NULL, y = "Tolerance") +
geom_ribbon(data = c8_res, aes(ymin = lb, ymax = ub, linetype=NA), alpha = .25) +
geom_line(data = c8_res) +
facet_wrap(~country, nrow = 2) +
theme(axis.text.x = element_text(angle = 45, hjust = 1),
strip.background = element_rect(fill = "white", colour = "white"))
ggsave("paper/figures/sc_line.pdf", plot = p_rd, width = 6, height = 4)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.