data-raw/create_presaved_steps_for_vignettes.R

In antaldaniel/retroharmonize: Ex Post Survey Data Harmonization

selected_eb_waves <- readRDS( file.path ('data-raw', 'trust_eb_9waves.rds'))
saveRDS(selected_eb_waves, file.path("inst", "eurob", "selected_eb_waves.rds") )


save ( ab_waves, documented_ab_waves, 
       ab_metadata, merged_ab, documenteded_merged_ab, R6, 
       harmonized_ab_waves, h_ab_structure, 
       file = file.path( "not_included", "afrob.rda"))

save (ab_metadata, to_harmonize, R6, 
      documented_ab_waves,documenteded_merged_ab, h_ab_structure,
      harmonized_ab_waves, file =file.path("inst", "afrob", "afrob_vignette.rda") , 
      compress = T)

R6 <- pull_survey ( merged_ab, id = "Afrobarometer_R6")
R6 <- select(R6, all_of(c("unique_id", "trust_president" )))


names ( test_trust)
test_trust <- select ( test_trust, all_of(
  c("rowid", "trust_european_commission")))
test_trust <- test_trust [1:100,]

names ( harmonized_ab_waves)

View ( harmonized_ab_waves )

as_factor ( harmonized_ab_waves$country)



load(system.file( file.path("afrob", "afrob_vignette.rda"),
             package = "retroharmonize"))

documented_eb_waves <- document_waves(eb_waves)
test_trust <- pull_survey(eb_waves, filename = "ZA4414_trust.rds")

save ( eb_trust_metadata, 
       test_trust, harmonized_eb_waves, documented_eb_waves, 
       file = file.path("inst", "eurob", "eurob_vignette.rda" ))
      
eb_plot <- numeric_harmonization %>%
  tidyr::pivot_longer( cols = contains("trust")) %>%
  filter ( !is.na(value))
names ( eb_plot )

ggplot (data = eb_plot, 
        aes ( 
          x = name,
        y = value , 
        color = name) ) +
  geom_point() +
  facet_wrap ( facets = "country")


names( numeric_harmonization)
library(ggplot2)

ab_plot <- harmonized_ab_waves %>%
  mutate_at ( vars(starts_with("trust")), 
              ~as_numeric(.)*within_country_weighting_factor) %>%
  select ( -all_of("within_country_weighting_factor") ) %>%
  group_by ( country, year ) %>%
  summarize_if ( is.numeric, mean, na.rm=TRUE ) %>%
  select ( all_of ( c("country", "year")), 
           contains("parliament"), 
           contains("president"), 
           contains("religious"),
           contains("traditional")) %>%
  pivot_longer ( cols = contains("trust")) %>%
  filter ( !is.nan(value)) %>%
  mutate ( name = case_when (
    grepl("president", name ) ~ "president", 
    grepl("parliament", name) ~ "parliament",
    grepl("religious", name) ~ "religious leaders",
    grepl("traditional", name ) ~ "traditional leaders", 
    TRUE ~ name 
  )) 

ab_plot_ordering <- ab_plot %>%
  select ( -all_of("name")) %>%
  group_by ( country ) %>%
  summarize ( avg_trust_level = mean(value)) %>%
  ungroup()

ab_plot2<- ab_plot %>%
  left_join ( ab_plot_ordering , by = 'country') %>%
  mutate ( country_code = countrycode::countrycode(
    country, "country.name", "iso2c"))  %>%
  ungroup()
names ( ab_plot2)

ab_plot_ordered  <- ab_plot2 %>%
  mutate ( country_code = forcats::fct_reorder(
    country_code, avg_trust_level), 
    country = forcats::fct_reorder(
      country, avg_trust_level)) %>%
  mutate ( year = as.factor(year))

levels ( ab_plot_ordered$country_code)


ab_plot1 <- ggplot ( data = ab_plot_ordered,
         aes ( x = year, y = value, 
               color = name)) +
  scale_y_continuous( breaks = c(0,1,2,3), 
                      n.breaks = 4,
                      labels = c("not", "little", "somewhat", "a lot")) +
  geom_point () + 
  scale_color_manual( values = c("#BAC615", "#4EC0E4", 
                                 "#E88500" , "#5C2320")) +
  facet_wrap ( facets = "country")+
  theme ( axis.text.x = element_text ( angle = 90), 
          legend.position = "bottom") +
  labs ( y = "trust level", x = "", color = "", 
         title = "Trust In African Leadership 2011-2018", 
         subtitle = "Based on Afrobarometer Rounds 5-7", 
         caption = "http://doi.org/10.5281/zenodo.3937746, R package website: retroharmonize.satellitereport.com \ua9 2020.")


ggsave(filename = file.path ("vignettes", "ab_plot1.png"), 
       plot = ab_plot1,
       units = "cm", width = 18, height = 13, dpi = 350)

here ("vignette")

file.exists(file.path ("vignette", "ab_plot1.png"))