In rfaridi/hiesBD: What the Package Does (One Line, Title Case)

Smoking

library(survey)
library(srvyr)
load(here("./data/HIES_2016/smoker_2016.rda"))

smokers_svy <- smoker_2016 %>%  
                  as_survey(id=psu,
                            strata=stratum,
                            weights=hhwgt)
class(smokers_svy)
summary(smokers_svy)

smokers_svy %>% 
   summarise(national=survey_mean(smoker)) %>%  # note the use of survey mean
   mutate(national=round(national,2)) %>% 
   pull(national)

Rural urban category wise household size

smokers_svy %>% 
    group_by(urbrural) %>% 
    summarize(smoker_perc=survey_mean(smoker)) %>% 
    ungroup()

smokers_svy %>% 
    group_by(division_name) %>% 
    summarize(smoker_perc=survey_mean(smoker)) %>% 
    ungroup()  %>% 
    rename(Division=division_name)  -> smokers_division

smokers_svy %>% 
    group_by(zila_name) %>% 
    summarize(smoker_perc=survey_mean(smoker)) %>% 
    ungroup()  %>% 
    rename(District=zila_name)  -> smokers_district
save(smokers_district, file=here("./data/HIES_2016/smokers_district.rda"))

bd_adm1 <- read_sf("~/Dropbox/miscStuff/spatialStuff/BGD_adm/BGD2015/ADM1/bgd_bnd_adm1/bgd_bnd_adm1.shp") # this needs to be changed too
bd_adm2 <- read_sf("~/Dropbox/miscStuff/spatialStuff/BGD_adm/BGD2015/ADM2/bgd_bnd_adm2.shp")

bd_adm1 %>% 
    mutate(division=Division,
       Division=toupper(Division)) -> bd_adm1
bd1_smokers <- bd_adm1 %>% 
                 left_join(smokers_division, by=c("Division"))

ggplot(bd2_smokers) +
    geom_sf(aes(fill=smoker_perc)) +
    geom_sf_text(aes(label=district),size=3, color="white")+
    scale_fill_gradient("smoker_perc", high = "darkred", low = "darkolivegreen1")+
    theme_void()

bd_adm2 %>% 
    mutate(district=District,
       District=toupper(District)) -> bd_adm2
bd2_smokers <- bd_adm2 %>% 
               left_join(smokers_district, by=c("District"))

ggplot(bd2_smokers) +
    geom_sf(aes(fill=smoker_perc)) +
    geom_sf_text(aes(label=district),size=3, color="white") +
    scale_fill_gradient("smoker_perc", high = "darkred", low = "darkolivegreen1") +
    theme_void()

ggplot(bd1_smokers) +
    geom_sf(aes(fill=smoker_perc)) +
    geom_sf_text(aes(label=division),size=3, color="white")+
    scale_fill_gradient("smoker_perc", high = "darkred", low = "darkolivegreen1") +
    theme_void()

rfaridi/hiesBD documentation built on Dec. 22, 2021, 3:05 p.m.