From 0b4910f87f7937cd4e3fcb42b2f7bf5ddb40c3e5 Mon Sep 17 00:00:00 2001 From: flaneuse Date: Tue, 13 Sep 2016 14:50:11 -0400 Subject: [PATCH] calc stats by commune --- R/HT_04_improvedToilets.R | 7 +++++++ R/HT_05_improvedWater.R | 9 +++++++++ 2 files changed, 16 insertions(+) diff --git a/R/HT_04_improvedToilets.R b/R/HT_04_improvedToilets.R index 1ebde23..218f7f1 100644 --- a/R/HT_04_improvedToilets.R +++ b/R/HT_04_improvedToilets.R @@ -143,6 +143,8 @@ toilet_admin2$admin1 = factor(toilet_admin2$admin1, levels = toilet_admin1$admin # Weights have no effect on EAs, but useful way to summarise by EA toilet_ea = calcPtEst('improved_toilet', by_var = 'cluster_id', design = DHSdesign, df = hh) +# -- By Admin3 (communes) -- +toilet_admin3 = calcPtEst('improved_toilet', by_var = 'admin3', design = DHSdesign, df = hh) # Admin1 map -------------------------------------------------------------- # NOTE!!!! Keep track of this, cuz it'll save you HOURS of debugging. @@ -213,6 +215,7 @@ plotMap(haiti_polygons, exportPlot = TRUE, fileName = '~/Creative Cloud Files/MAV/Projects/Haiti_WASH-PAD_2016-09/exported_R/HTI_imprtoilet_adm2.pdf') + # -- dot plot -- pairGrid(toilet_admin2, y_var = 'admin2', @@ -264,3 +267,7 @@ p = ur_pairGrid(urb_rural_toilet, size_SE = 2, file_name = '~/Creative Cloud Files/MAV/Projects/Haiti_WASH-PAD_2016-09/exported_R/HTI_UR_toilets.pdf') +# Admin 3 dot plot -------------------------------------------------------- + +pairGrid_admin3(toilet_admin3, savePlots = FALSE, y_var = 'admin3', + colorDot = colour_toilet, fill_limits = colour_limits) + geom_vline(xintercept = 0.27558) diff --git a/R/HT_05_improvedWater.R b/R/HT_05_improvedWater.R index 411cb1f..01918c7 100644 --- a/R/HT_05_improvedWater.R +++ b/R/HT_05_improvedWater.R @@ -136,11 +136,20 @@ water_admin2 = left_join(water_admin2, water_admin1_sum, by = 'admin1') # Resort order of admin1 based on the admin1 averages. Best (most improved toilet) == top. water_admin2$admin1 = factor(water_admin2$admin1, levels = water_admin1$admin1) +# -- By Admin3 -- +water_admin3 = calcPtEst('impr_water_under30min', by_var = 'admin3', design = DHSdesign, df = hh) + # Admin2 dot plot --------------------------------------------------------- pairGrid(water_admin2, savePlots = F, y_var = 'admin2') + +# Admin 3 dot plot -------------------------------------------------------- + +pairGrid_admin3(water_admin3, savePlots = F, y_var = 'admin3', + facet_var = NA, incl_comparison = FALSE, fill_limits = colour_limits) + # Admin1 map -------------------------------------------------------------- haiti_polygons = right_join(water_admin1_PaP, dhs_geo$df, by = c('region_name' = 'DHSREGFR'))