-
Notifications
You must be signed in to change notification settings - Fork 1
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Showing
15 changed files
with
1,274 additions
and
48 deletions.
There are no files selected for viewing
File renamed without changes.
File renamed without changes.
File renamed without changes.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
196 changes: 196 additions & 0 deletions
196
...er-radicalisation-enviro/presentation_cpsa/figure3_descriptive_breakdown_responsibility.R
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,196 @@ | ||
| # Packages ---------------------------------------------------------------- | ||
| library(dplyr) | ||
| library(ggplot2) | ||
| source("_papier-radicalisation-enviro/article/environment.R") | ||
|
|
||
| # Data -------------------------------------------------------------------- | ||
|
|
||
| data <- readRDS("_SharedFolder_quorum-enviro/data/cleanData/data.rds") %>% | ||
| filter(ses_age >= 18) %>% | ||
| mutate(age_cat = get_age_category(ses_age), | ||
| ses_region = case_when( | ||
| ses_region_qc == 1 ~ "Quebec", | ||
| ses_region_ont == 1 ~ "Ontario", | ||
| ses_region_west == 1 ~ "Western provinces", | ||
| ses_region_mari == 1 ~ "Maritimes" | ||
| ), | ||
| educ_level = ifelse(ses_educBHS == 1, 0, NA), | ||
| educ_level = ifelse(ses_educCollege == 1, 0.5, | ||
| educ_level), | ||
| educ_level = ifelse(ses_educUniv == 1, 1, | ||
| educ_level), | ||
| responsibility_climateChange_Govt = (responsability_climateChangeFedGovt + responsability_climateChangeProvGovt) / 2) %>% | ||
| select(age_cat, educ_level, ses_region, | ||
| responsibility_climateChange_Govt, | ||
| responsibility_climateChange_Enterprise = responsability_climateChangeEnterprise, | ||
| responsibility_climateChange_Citizens = responsability_climateChangeCitizens) %>% | ||
| tidyr::pivot_longer(., cols = c(starts_with("responsibility_climateChange")), | ||
| names_to = "actor", | ||
| names_prefix = "responsibility_climateChange_", | ||
| values_to = "responsibility") %>% | ||
| mutate(actor = case_when( | ||
| actor == "Govt" ~ "Government\n", | ||
| actor == "Enterprise" ~ "Enterprises\n", | ||
| actor == "Citizens" ~ "Citizens\n" | ||
| ), | ||
| actor = factor(actor, levels = c("Citizens\n", "Government\n", "Enterprises\n")), | ||
| ) | ||
|
|
||
| # Age -------------------------------- | ||
|
|
||
| data_group_age <- data %>% | ||
| filter(age_cat < 70) %>% | ||
| group_by(age_cat, actor) %>% | ||
| summarise(mean_responsibility = mean(responsibility, na.rm = T), | ||
| sd = sd(responsibility, na.rm = T), | ||
| n = n()) %>% | ||
| mutate(me = margin_error(mean_responsibility, sd, n), | ||
| conf_low = mean_responsibility - me, | ||
| conf_high = mean_responsibility + me) | ||
|
|
||
| plot_age <- ggplot(data_group_age, aes(x = mean_responsibility, y = age_cat * -1)) + | ||
| facet_wrap(~actor, ncol = 1, | ||
| strip.position = "left") + | ||
| geom_rect(data = subset(data_group_age, actor == 'Enterprises\n'), | ||
| fill = "grey90", xmin = -2,xmax = 2, | ||
| ymin = -100,ymax = 100) + | ||
| geom_rect(data = subset(data_group_age, actor == 'Government\n'), | ||
| fill = "grey95", xmin = -1.5,xmax = 1.5, | ||
| ymin = -100,ymax = 100) + | ||
| geom_rect(data = subset(data_group_age, actor == 'Citizens\n'), | ||
| fill = "white", xmin = -1.5,xmax = 1.5, | ||
| ymin = -100,ymax = 100) + | ||
| geom_point() + | ||
| geom_linerange(aes(xmin = conf_low, xmax = conf_high)) + | ||
| scale_y_continuous(name = "\nAge (grouped by 5 years)\n", | ||
| labels = c(abs(seq(from = -18, to = -70, by = -5))), | ||
| breaks = seq(from = -18, to = -70, by = -5), | ||
| expand = c(0.1, 0.1)) + | ||
| scale_x_continuous(breaks = c(0.05, 0.35, 0.64, 0.945), | ||
| limits = c(0, 1), | ||
| labels = c("None at all", "A little bit", "Partially", "Completely")) + | ||
| clessnize::theme_clean_light(base_size = 15) + | ||
| theme(strip.background.x = element_rect(fill = "#F2F2F2", color = NA), | ||
| strip.background.y = element_blank(), | ||
| strip.text = element_text(size = 9), | ||
| axis.title.x = element_blank(), | ||
| axis.title.y = element_text(hjust = 0.5), | ||
| axis.text.y = element_text(size = 7), | ||
| strip.placement = "outside", | ||
| panel.grid.major.y = element_blank(), | ||
| axis.line.x = element_blank(), | ||
| axis.ticks.x = element_blank(), | ||
| legend.title = element_text(size = 7.5)) | ||
|
|
||
| # Educ -------------------------------------------------------------------- | ||
|
|
||
| data_group_educ <- data %>% | ||
| group_by(educ_level, actor) %>% | ||
| summarise(mean_responsibility = mean(responsibility, na.rm = T), | ||
| sd = sd(responsibility, na.rm = T), | ||
| n = n()) %>% | ||
| mutate(me = margin_error(mean_responsibility, sd, n), | ||
| conf_low = mean_responsibility - me, | ||
| conf_high = mean_responsibility + me) | ||
|
|
||
| plot_educ <- ggplot(data_group_educ, aes(x = mean_responsibility, y = educ_level)) + | ||
| facet_wrap(~actor, ncol = 1, | ||
| strip.position = "right") + | ||
| geom_rect(data = subset(data_group_educ, actor == 'Enterprises\n'), | ||
| fill = "grey90", xmin = -2,xmax = 2, | ||
| ymin = -100,ymax = 100) + | ||
| geom_rect(data = subset(data_group_educ, actor == 'Government\n'), | ||
| fill = "grey95", xmin = -1.5,xmax = 1.5, | ||
| ymin = -100,ymax = 100) + | ||
| geom_rect(data = subset(data_group_educ, actor == 'Citizens\n'), | ||
| fill = "white", xmin = -1.5,xmax = 1.5, | ||
| ymin = -100,ymax = 100) + | ||
| geom_point() + | ||
| geom_linerange(aes(xmin = conf_low, xmax = conf_high)) + | ||
| scale_y_continuous(name = "\nEducation level\n", | ||
| labels = c("None", "College", "University"), | ||
| breaks = c(0, 0.5, 1), | ||
| expand = c(0.1, 0.1)) + | ||
| scale_x_continuous(breaks = c(0.05, 0.35, 0.64, 0.945), | ||
| limits = c(0, 1), | ||
| labels = c("None at all", "A little bit", "Partially", "Completely")) + | ||
| clessnize::theme_clean_light(base_size = 15) + | ||
| theme(strip.background.x = element_rect(fill = "#F2F2F2", color = NA), | ||
| strip.background.y = element_blank(), | ||
| strip.text = element_text(size = 9), | ||
| axis.title.x = element_blank(), | ||
| axis.title.y = element_text(hjust = 0.5), | ||
| axis.text.y = element_text(size = 10), | ||
| strip.placement = "outside", | ||
| panel.grid.major.y = element_blank(), | ||
| axis.line.x = element_blank(), | ||
| axis.ticks.x = element_blank(), | ||
| legend.title = element_text(size = 7.5)) | ||
|
|
||
|
|
||
| # Region ------------------------------------------------------------------ | ||
|
|
||
| data_group_region <- data %>% | ||
| group_by(ses_region, actor) %>% | ||
| summarise(mean_responsibility = mean(responsibility, na.rm = T), | ||
| sd = sd(responsibility, na.rm = T), | ||
| n = n()) %>% | ||
| mutate(me = margin_error(mean_responsibility, sd, n), | ||
| conf_low = mean_responsibility - me, | ||
| conf_high = mean_responsibility + me) | ||
|
|
||
| plot_region <- ggplot(data_group_region, aes(x = mean_responsibility, y = ses_region)) + | ||
| facet_wrap(~actor, ncol = 1, | ||
| strip.position = "left") + | ||
| geom_rect(data = subset(data_group_region, actor == 'Enterprises\n'), | ||
| fill = "grey90", xmin = -2,xmax = 2, | ||
| ymin = -100,ymax = 100) + | ||
| geom_rect(data = subset(data_group_region, actor == 'Government\n'), | ||
| fill = "grey95", xmin = -1.5,xmax = 1.5, | ||
| ymin = -100,ymax = 100) + | ||
| geom_rect(data = subset(data_group_region, actor == 'Citizens\n'), | ||
| fill = "white", xmin = -1.5,xmax = 1.5, | ||
| ymin = -100,ymax = 100) + | ||
| geom_point() + | ||
| geom_linerange(aes(xmin = conf_low, xmax = conf_high)) + | ||
| scale_y_discrete(name = "\nRegion\n", | ||
| expand = c(0.1, 0.1)) + | ||
| scale_x_continuous(breaks = c(0.05, 0.35, 0.64, 0.945), | ||
| limits = c(0, 1), | ||
| labels = c("None at all", "A little bit", "Partially", "Completely")) + | ||
| clessnize::theme_clean_light(base_size = 15) + | ||
| theme(strip.background.x = element_rect(fill = "#F2F2F2", color = NA), | ||
| strip.background.y = element_blank(), | ||
| strip.text = element_text(size = 9), | ||
| axis.title.x = element_blank(), | ||
| axis.title.y = element_text(hjust = 0.5), | ||
| axis.text.y = element_text(size = 10), | ||
| strip.placement = "outside", | ||
| panel.grid.major.y = element_blank(), | ||
| axis.line.x = element_blank(), | ||
| axis.ticks.x = element_blank(), | ||
| legend.title = element_text(size = 7.5)) | ||
|
|
||
|
|
||
| # Merge graphs ------------------------------------------------------------ | ||
|
|
||
| top_row <- cowplot::plot_grid(plot_age, plot_educ, ncol = 2) | ||
| bottom_row <- cowplot::plot_grid(NULL, plot_region, NULL, ncol = 3, rel_widths = c(1, 2.5, 1), | ||
| axis = "tblr") | ||
|
|
||
| # Combinez les rangées en une seule figure | ||
| combined_plot <- cowplot::plot_grid(top_row, NULL, bottom_row, | ||
| ncol = 1, rel_heights = c(1, 0.07, 1)) | ||
|
|
||
| plot <- cowplot::ggdraw() + | ||
| cowplot::draw_plot(combined_plot, 0, 0.1, 1, 0.9) + | ||
| cowplot::draw_label("Mean responsibility", | ||
| x = 0.5, y = 0.075, hjust = 0.5, | ||
| size = 15) + | ||
| cowplot::draw_label("The error bars represent the 95% confidence interval.\nDue to large error margins, respondents over 70 years old were removed from the age visualisation.", | ||
| x = 0.05, y = 0.025, hjust = 0, size = 9, fontface = "italic") + | ||
| theme(plot.background = element_rect(fill = "white", color = NA)) | ||
|
|
||
| ggsave(plot, filename = "_SharedFolder_quorum-enviro/_papier-radicalisation-enviro/graphs/dans_article/figure3_descriptive_breakdown_responsibility.png", | ||
| width = 10, height = 9) | ||
|
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Oops, something went wrong.