calc stats by commune

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)

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'))