Como buscar a localização de estabelecimentos de saúde? Usando o pacote geobr em R.

geobr-read_health_facilities.R

# Dúvida da Carla Graciane: Como buscar a localização de estabelecimentos de saúde?
# A pergunta foi feita após a apresentação sobre o geobr no Meetup da R-Ladies SP
# Meetup: https://r-ladies-sao-paulo.github.io/2023-11-pacotes-br/slides.html
# Apresentação geobr: https://r-ladies-sao-paulo.github.io/2023-11-pacotes-br/slide-geobr.html

# Encontrei a função read_health_facilities()
# RECOMENDO ler essa página de documentação para entender melhor de onde vem os dados!
# https://ipeagit.github.io/geobr/reference/read_health_facilities.html


# Carregando pacotes -------
library(geobr)
library(ggplot2)
library(dplyr)


# Importando os dados ----------
locais_saude <- read_health_facilities()

# Quais são as colunas disponíveis? -------

glimpse(locais_saude)
# Rows: 360,177
# Columns: 8
# $ code_cnes    <int> 2684403, 2684411, 2684438, 2576694, 258337…
# $ code_muni    <int> 412510, 412510, 412510, 412520, 412520, 41…
# $ code_state   <dbl> 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41…
# $ abbrev_state <chr> "PR", "PR", "PR", "PR", "PR", "PR", "PR", …
# $ date_update  <chr> "2015-05-27", "2015-05-27", "2015-05-27", …
# $ year_update  <dbl> 2015, 2015, 2015, 2015, 2015, 2015, 2015, …
# $ data_source  <chr> "PMAQ", "PMAQ", "PMAQ", "PMAQ", "PMAQ", "P…
# $ geom         <POINT [°]> POINT (-50.297 -25.683), POINT (-50.…


# Seria interessante conseguir uma base de dados
# com a coluna do código CNES, para poder fazer join
# parece que tem aqui:
# https://cnes.datasus.gov.br/pages/downloads/arquivosBaseDados.jsp


# Atualização dos dados ------------

locais_saude |>
  count(year_update) |>
  arrange(desc(year_update))
#  year_update      n
# 1        2018   8925
# 2        2017  21667
# 3        2016  19477
# 4        2015 310108


# Atenção: parece que não tem sido atualizado desde 2018:
# É possível conversar com a equipe do pacote
# Perguntei aqui:
# https://github.com/ipeaGIT/geobr/issues/334


# Visualização ----------
# visualizando a localização dos pontos dos estabelecimentos no município de São Paulo
locais_saude |>
  # municipio de sao paulo
  filter(code_muni == 355030) |> # ATENÇÃO: Usar código IBGE com 6 dígitos
  ggplot() +
  geom_sf()