In j-hagedorn/TBDfun: TBD Functions

library(tidyverse); library(dplyr); library(sf); library(leaflet)

master <- st_read(system.file("extdata","mi_master_polygons.shp", package = "TBDfun"))
head(master)

tract <- master %>%
  group_by(NAME) %>%
  summarize(population = sum(estimate, na.rm = TRUE))

tract_map <- leaflet(tract) %>% addPolygons() %>% addTiles()
print(tract_map)

county <- master %>%
  group_by(county) %>%
  summarize(population = sum(estimate, na.rm = TRUE))

county_map <- leaflet(county) %>% addPolygons() %>% addTiles()
print(county_map)

cmhsp <- master %>%
  group_by(CMHSP) %>%
  summarize(population = sum(estimate, na.rm = TRUE))

cmhsp_map <- leaflet(cmhsp) %>% addPolygons() %>% addTiles()
print(cmhsp_map)

pihp <- master %>%
  group_by(PIHP) %>%
  summarize(
    population = sum(estimate, na.rm = TRUE))

pihp_map <- leaflet(pihp) %>% addPolygons() %>% addTiles()
print(pihp_map)

agency <- master %>%
  group_by(agency) %>%
  summarize(
    population = sum(estimate, na.rm = TRUE))

agency_map <- leaflet(agency) %>% addPolygons() %>% addTiles()
print(agency_map)