Skip to content

Instantly share code, notes, and snippets.

@beatrizmilz
Last active December 1, 2023 17:38
Show Gist options
  • Save beatrizmilz/191eddf764670a6884284a95cce51d37 to your computer and use it in GitHub Desktop.
Save beatrizmilz/191eddf764670a6884284a95cce51d37 to your computer and use it in GitHub Desktop.
Como buscar a localização de estabelecimentos de saúde? Usando o pacote geobr em 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()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment