website_examples.R

In flaneuse/llamar: Plotting Functions and Styles for Livelihoods Analysis Work

library(dplyr)
library(haven)
library(RColorBrewer)


# data sets ---------------------------------------------------------------

dhs_orig = read_dta('~/Documents/USAID/Rwanda/processeddata/DHS_2010_2015_analysis.dta')
dhs = removeAttributes(dhs_orig) %>% 
  filter(year == 2014)

dhs = factorize(dhs, dhs_orig, 'lvdzone', 'lz')
df2 = data.frame(avg = sample(-100:100, 10)/100, region = letters[1:10], ci = sample(1:100, 20)/1000) %>% mutate(lb = avg - ci, ub = avg + ci)


# plot_dot ----------------------------------------------------------------
plot_avg_dot(dhs, value_var = 'stunted2', by_var = 'lz', 
              percent_vals = TRUE, weight_var = 'cweight',
             dot_size = 9, include_n = FALSE, x_limits = c(0.15, 0.65), x_breaks = seq(0.2, 0.6, by = 0.2),
             dot_fill_cont = rev(brewer.pal(11, 'Spectral')[1:6]), sat_threshold = 0.65) +
  theme_stroke()

save_plot('~/GitHub/llamar/img/plot_dot1.png', width = 8, height = 8)

plot_dot(df2, by_var = 'region', value_var = 'avg', ref_line = 0, 
         ref_text = 'no change', label_ref = FALSE, lollipop = TRUE, value_label_offset = .075,
         dot_fill_cont = brewer.pal(10, 'RdYlBu'), percent_vals = TRUE) +
  theme_stroke()
save_plot('~/GitHub/llamar/img/plot_dot2.png', width = 5, height = 5)


# plot_avg_dot ------------------------------------------------------------

p = plot_avg_dot(dhs, value_var = 'stunted2', by_var = 'lz', 
             percent_vals = TRUE, weight_var = 'cweight',
             dot_size = 9, include_n = TRUE, x_breaks = c(0.2, 0.4, 0.6),
             dot_fill_cont = rev(brewer.pal(11, 'Spectral')[1:6]), sat_threshold = 0.65) 

save_plot('~/GitHub/llamar/img/plot_avg_dot.png',p, width = 8, height = 8)


# DHS musings -------------------------------------------------------------

plot_avg_dot(dhs, value_var = 'improvedWater', by_var = 'lz', 
             percent_vals = TRUE, weight_var = 'cweight',
             dot_size = 9, include_n = TRUE,
             dot_fill_cont = (brewer.pal(9, "Blues")), sat_threshold = 0.65) +
  theme_stroke()

w = read_dta('~/Documents/USAID/Rwanda/rawdata/RW_2014-15_DHS/rwir70dt/RWIR70FL.DTA')

w = w %>% mutate(contrac = ifelse(v313 == 3, 1, 
                                  ifelse(is.na(v313), NA, 0)))

w = factorize(w, w, 'sdistrict', 'dist')

plot_avg_dot(w, value_var = 'contrac', by_var = 'dist', 
             percent_vals = TRUE, 
             dot_size = 9, include_n = TRUE,
             dot_fill_cont = (brewer.pal(9, "Blues")), sat_threshold = 0.65) +
  theme_stroke()