Construção de indicadores de acessibilidade espacial à saúde com o pacote asha

Exemplo a partir do caso da cidade de São Paulo

Bruno Pinheiro1

Escola de Artes, Ciências e Humanidades da Universidade de São Paulo
bruno.pinheiro.oliveira@usp.br

Alexandre Ribeiro Leichsenring2

Escola de Artes, Ciências e Humanidades da Universidade de São Paulo
alexandre.leichsenring@usp.br

2019-04-05

Introdução

bla bla bla…

Cada projeto de dados enfrenta desafios particulates. Fatores como a existência ou não dos dados, se o acesso é público, se os dados são estruturados ou não.

O problema abordado pela pesquisa foi o da relação entre o sistema municipal de atenção básica de saúde, sua organização territorial e a distribuição espacial da população. A pergunta norteadora central foi “como otimizar a acessibilidade espacial da atenção básica de saúde?”. Desta derivaram perguntas secundárias: “como avaliar a acessibilidade espacial da atenção básica de saúde?”, “quais os dados necessários, quais os dados existentes e que modelos poderiam ser utilizados?”, “dada a definição da metodologia, como viabilizar replicabilidade?”, “como lidar com a realidade de que a escassez de recursos imporá barreiras à redistribuição espacial de unidades de saúde?”, “como tornar o processo o menos custoso possível na interação com as API’s?”

Além das consultas à bibliografia, a resposta a estas perguntas envolveu um conjunto de procedimentos computacionais e estatísticos de complexidades variadas, no caso desta pesquisa realizados no R: manipulação de bases de dados espaciais, rotina de captura de informações espaciais usando as API’s do Google, criação de mapas e gráficos, geração de estatísticas descritivas e cálculo de indicadores de acessibilidade espacial.

Ao fim a metodologia foi encapsulada no pacote do R chamado asha. A proposta deste artigo é, por meio do relato do processo de investigação, de construção da metodologia e dos indicadores, apresentar de forma prática um exemplo de uso de dados públicos para a construção de uma ferramenta propositiva para o processo de produção de políticas públicas, incluindo os desafios enfrentados e as soluções que permitiram superá-los.

“O Brasil já possui diversos casos práticos de publicação de dados na web: IBGE, Datasus, IPEA, INEP, Portal da Transparência, e muitos outros. Só que esses casos não estão aderentes às premissas dos dados abertos e da reusabilidade dos dados publicados, além de cada um fazer suas publicações em formatos e padrões diferentes.” (Miranda, p. 64)

Objetivos

Estrutura do pacote e funcionalidade

A instalação do pacote e carregado no ambiente do R conforme indicado a seguir:

Funcionalidades centrais do asha

Funcionalidades do pacote asha
Family Function
Obtenção de dados get_zip(), get_censo()
Estimação de origem-destino asha_nn(), asha_intersect()

Passo 1: Obtenção dos dados

Dados fundamentais

A construção das medidas de acessibilidade disponíveis no pacote asha depende de um conjunto específico de dados a respeito do território, da população e dos serviços de saúde pública.

As bases de dados utilizadas na pesquisa realizada com dados da cidade de São Paulo estão incorporadas ao pacote e podem servir de exemplo para os dados que deverão ser obtidos para viabilizar a modelagem e construção dos indicadores. Para ver todos os dados disponíveis no pacote digite data(package = "asha") no console`.

Especificamente, a pesquisa baseou-se nas seguintes bases de dados:

Characterization of study data, 2010-2017
Information Source Period
Census tracts mesh Census 2010 (IBGE) 2010
Habitants and income Census 1010 aggregated by tracts (IBGE) 2010
Location of primary care units GEOSAMPA/PMSP; CNES; Google Geocode API 2015/2016
Health professionals CNES sep/2017
Primary care units coverage mesh Information request at SMS/PMSP 2017
Distances in length and travel time Google Distance Matrix API 2017

Após o download dos dados e da realização do pré-processamento foram definidas as seguintes bases de dados chave:

  • setores: malha de setores censitários
  • populacao: habitantes e renda
  • ubs_pontos: geolocalização das UBS
  • ubs_malhas: malha dos territórios de UBS
  • prof_saude: número de médicos e enfermeiros nas unidades de saúde de São Paulo

A base de dados final, resultado da modelagem da estrutura de origem e destino e da construção dos indicadores é a od_indicadores, cuja estrutura é apresentada abaixo:

Baixar os dados

Problema
  • Garantir replicabilidade da obtenção dos dados para outras unidades territoriais, dada a diversidade de fontes e tipos de dados
  • Dados de distâncias e tempos de viagem demandam interação com API’s do Google
Limitações
  • Dificuldade de obtenção programática de dados do DATASUS/CNES
  • Demanda por dados que ou não existem ou não estão ativamente disponíveis para acesso na internet
Soluções
  • Efetuar pedidos de acesso à informação
  • Obter dados do DATASUS/CNES manualmente
  • Estabelecer funções para obter dados do IBGE, para interagir com a Google Distance Matrix API, e para baixar qualquer arquivo comprimido em formato zip

Baixar dados com get_zip()

A função mais básica para obtenção de dados do pacote asha é a get_zip(). Ela é uma simples implementação conjunta das funções download.file() e unzip() do pacote utils, um dos core R packages. Portanto, além de baixar qualquer arquivo .zip da internet, ela também extrai os arquivos compactados. Além disso, a função dispõe de dois métodos para gerenciamento dos dados baixados.

O método básico é o temporário, caso em que os arquivos são locados no diretório temporário do sistema. O workflow deste método envolve o uso de duas outras funções, também disponibilizadas no pacote asha: list_tmp() facilita a listagem dos arquivos na pasta temporária, aceitando padrões de busca, inclusive expressões regulares e clean_tmp() exclui o diretório temporário, caso seja preciso limpá-lo.

O segundo método é o arquivamento local dos dados. Neste caso, os arquivos são guardados no diretório indicado na função: get_zip(savedir = "path")

No caso particular da investigação, esta função foi usada para obter os pontos geolocalizados das Unidades Básicas de Saúde, shapefile disponível no portal Geosampa.

O resultado mostra que os dados foram corretamente baixados. Como não indiquei nenhum diretório de salvamento os dados estão na pasta temporária. Aqui é que entra a função list_tmp(). Com ela eu posso listar os arquivos da pasta temporária.

list_tmp()
#>  [1] "/tmp/RtmpXLaA8V/file54e21ebab702.zip"                                            
#>  [2] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_AMBULATORIOS_ESPECIALIZADOS.dbf"       
#>  [3] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_AMBULATORIOS_ESPECIALIZADOS.shp"       
#>  [4] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_AMBULATORIOS_ESPECIALIZADOS.shx"       
#>  [5] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_HOSPITAL.dbf"                          
#>  [6] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_HOSPITAL.shp"                          
#>  [7] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_HOSPITAL.shx"                          
#>  [8] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_SAUDE_MENTAL.dbf"                      
#>  [9] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_SAUDE_MENTAL.shp"                      
#> [10] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_SAUDE_MENTAL.shx"                      
#> [11] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_UBS-POSTO_DE_SAUDE-CENTRO_DE_SAUDE.dbf"
#> [12] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_UBS-POSTO_DE_SAUDE-CENTRO_DE_SAUDE.shp"
#> [13] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_UBS-POSTO_DE_SAUDE-CENTRO_DE_SAUDE.shx"
#> [14] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_UNIDADES_DST-AIDS.dbf"                 
#> [15] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_UNIDADES_DST-AIDS.shp"                 
#> [16] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_UNIDADES_DST-AIDS.shx"                 
#> [17] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_URGENCIA-EMERGENCIA.dbf"               
#> [18] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_URGENCIA-EMERGENCIA.shp"               
#> [19] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_URGENCIA-EMERGENCIA.shx"               
#> [20] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_VIGILANCIA_EM_SAUDE.dbf"               
#> [21] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_VIGILANCIA_EM_SAUDE.shp"               
#> [22] "/tmp/RtmpXLaA8V/SAD69-96_SHP_TEMA_-_SAUDE_VIGILANCIA_EM_SAUDE.shx"               
#> [23] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_AMBULATORIOS_ESPECIALIZADOS.dbf"         
#> [24] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_AMBULATORIOS_ESPECIALIZADOS.shp"         
#> [25] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_AMBULATORIOS_ESPECIALIZADOS.shx"         
#> [26] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_HOSPITAL.dbf"                            
#> [27] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_HOSPITAL.shp"                            
#> [28] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_HOSPITAL.shx"                            
#> [29] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_SAUDE_MENTAL.dbf"                        
#> [30] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_SAUDE_MENTAL.shp"                        
#> [31] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_SAUDE_MENTAL.shx"                        
#> [32] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_UBS-POSTO_DE_SAUDE-CENTRO_DE_SAUDE.dbf"  
#> [33] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_UBS-POSTO_DE_SAUDE-CENTRO_DE_SAUDE.shp"  
#> [34] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_UBS-POSTO_DE_SAUDE-CENTRO_DE_SAUDE.shx"  
#> [35] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_UNIDADES_DST-AIDS.dbf"                   
#> [36] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_UNIDADES_DST-AIDS.shp"                   
#> [37] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_UNIDADES_DST-AIDS.shx"                   
#> [38] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_URGENCIA-EMERGENCIA.dbf"                 
#> [39] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_URGENCIA-EMERGENCIA.shp"                 
#> [40] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_URGENCIA-EMERGENCIA.shx"                 
#> [41] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_VIGILANCIA_EM_SAUDE.dbf"                 
#> [42] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_VIGILANCIA_EM_SAUDE.shp"                 
#> [43] "/tmp/RtmpXLaA8V/SIRGAS_SHP_TEMA_-_SAUDE_VIGILANCIA_EM_SAUDE.shx"

Além do arquivo baixado, que ganhou o nome temporário /tmp/RtmpXLaA8V/file54e21ebab702.zip, todos os arquivos que estavam compactados nele estão presentes, já descompactados.

É possível fazer buscas mais específicas. Por exemplo, para ver os arquivos .shp:

Ou para listar apenas o nome do arquivo e não o seu caminho completo:

Como o meu interesse é apenas no shapefile de UBS, e quero os dados projetados em Sirgas, posso usar expressões regulares e criar um padrão mais específico para extrair apenas o arquivo que desejo importar.

E então importar:

library(sf)
library(dplyr)
ubs_pontos <- read_sf(arquivo)
glimpse(ubs_pontos)
#> Observations: 455
#> Variables: 23
#> $ eq_cnes    <chr> "2788527", "3913007", "2091348", "2091321", "2752166"…
#> $ eq_enderec <chr> "R BANDEIRA DE ARACAMBI, 704", "AV DR ABRAAO RIBEIRO,…
#> $ eq_nome    <chr> "JD SAO FRANCISCO II - AMA/UBS INTEGRADA", "ALEXANDRE…
#> $ eq_bairro  <chr> "JD  SAO FRANCISCO", "BARRA FUNDA", "BUTANTA", "PACAE…
#> $ eq_cep     <chr> "08310010\n", "01133020\n", "05503000\n", "01246904\n…
#> $ eq_setcens <chr> "355030875000151", "355030806000026", "35503081200002…
#> $ eq_ar_pond <chr> "3550308005277", "3550308005029", "3550308005105", "3…
#> $ eq_cd_dist <chr> "75", "6", "12", "45", "1", "1", "1", "2", "3", "3", …
#> $ eq_distr   <chr> "SAO RAFAEL", "BARRA FUNDA", "BUTANTA", "JARDIM PAULI…
#> $ eq_cd_subp <chr> "30", "08", "10", "11", "25", "25", "25", "11", "01",…
#> $ eq_subpref <chr> "SAO MATEUS", "LAPA", "BUTANTA", "PINHEIROS", "MOOCA"…
#> $ eq_regiao5 <chr> "Leste", "Oeste", "Oeste", "Oeste", "Leste", "Leste",…
#> $ eq_regiao8 <chr> "Leste 2", "Oeste", "Oeste", "Oeste", "Leste 1", "Les…
#> $ eq_telefon <chr> "27516712", "34662502", "37267492", "30851134", "2605…
#> $ eq_cd_esfe <chr> "01", "04", "02", "02", "01", "01", "01", "01", "01",…
#> $ eq_esfera  <chr> "Municipal", "Privado", "Estadual", "Estadual", "Muni…
#> $ eq_cd_tipo <chr> "045", "018", "018", "018", "045", "045", "045", "045…
#> $ eq_tipo    <chr> "UNIDADE BASICA DE SAUDE", "CENTRO DE SAUDE", "CENTRO…
#> $ eq_classe  <chr> "UBS-POSTO_DE_SAUDE-CENTRO_DE_SAUDE", "UBS-POSTO_DE_S…
#> $ eq_leitos  <chr> "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"…
#> $ eq_cd_clas <chr> "5", "5", "5", "5", "5", "5", "5", "5", "5", "5", "5"…
#> $ eq_id      <chr> "206", "216", "217", "218", "445", "446", "447", "448…
#> $ geometry   <POINT> POINT (351300 7385828), POINT (330613.5 7397760), P…

Em resumo, o processo de importação deste arquivo, usando o método temporário, teria sido o seguinte:

Já usando o método de armazenamento local o processo, seria:

Baixar dados com get_censo()

Já a função get_censo() é específica. Ela aplica a função get_zip() para baixar expecificamente os seguintes dados do Censo 2010:

  • Resultados do universo agregados por setor censitário
  • Malhas digitais intramunicipais

Os dados são baixados de ftp://ftp.ibge.gob.br ou ftp://geoftp.ibge.gov.br. O primeiro disponibiliza os dados tabulares das pesquisas realizadas pelo IBGE. O segundo os dados espaciais produzidos pelo órgão, como as malhas digitais dos municípios, dos estados e dos setores censitários.

Dada a dimensão do município de São Paulo, seus arquivos são pesados. Por isto, exemplificarei o uso de get_censo() obtendo dados do estado de Roraima.

Baixar dados tabulares

Baixar dados espaciais

Para baixar os dados espaciais (a malha de setores censitários), basta adicionar o argumento mode na função getìbge().

Em qualquer um dos casos basta indicar o código de outra Unidade da Federação para baixar os dados de outro estado. O exemplo abaixo, por exemplo, baixaria os dados tabulares e as malhas de setores censitários de Santa Catarina

Por fim, a função get_censo pode ser usada para baixar outras malhas intramunicipais. Para isto é preciso indicar no argumento mesh, aceito quando datatype = "geo", qual delas será obtida. Por padrão a função baixa os setores censitários, mas também aceita os argumentos “distritos”, “municipios” e “sub_distritos”. Por exemplo, para baixar os distritos de Santa Catarina:

get_censo("sc", "geo", mesh = "distritos")

O script de importação dos dados utilizados na pesquisa em questão pode servir de exemplo para a obtenção dos dados.

Obtenção de dados de outras fontes

Uma informação essencial para a avaliação da acessibilidade do modelo territorial do sistema municipal de atenção básica em saúde é a delimitação territorial da cobertura das UBS. No caso da cidade de São Paulo, a Secretaria Municipal de Saúde dispunha de shapefiles utilizados na gestão, porém não disponibilizados à população. Por isto foi necessário impetrar um pedido de acesso à informação no e-Sic da Prefeitura de São Paulo. Além destes, os dados de saúde também foram obtidos manualmente junto ao DATASUS/CNES.

Obtenção de dados da Google Distance Matrix API

Além destes conjuntos de dados já relacionados acima, a metodologia depende da obtenção de dados de rotas (tempos de viagem), para a modelagem da estrutra OD do sistema de atenção básica. Esta é uma das vantagens centrais adotadas na pesquisa e implementada no pacote asha por meio da função asha_dists(), pois baseia a modelagem espacial no tempo efetivo da rota mais próxima entre os setores censitários e as UBS e não em distâncias euclidianas como muitos trabalhos na área adotam. Isto confere maior consitência aos indicadores espaciais.

Ainda que seja também uma função de obtenção de dados, ela será apresenta mais adiante. Como o objetivo desta vinheta é apresentar a aplicação da metodologia utilizando as funções do pacote asha de forma aplicada, esta função será apresentada na explicação sobre a modelagem da estrutura de sistema de atenção básica de saúde, mantendo-a no contexto da sua aplicação.

Passo 2: Preparação dos dados

Problema
  • Existência de variáveis desnecessárias
  • Classes de dados não padronizadas
  • Nomes das variáveis codificados, dificultando o uso
  • Malhas de UBS mais atual que geolocalização de pontos
  • Informações desatualizadas e inconsistens do CNES e da Prefeitura a respeito das UBS existentes no município
  • Dados mal documentados (metadados e variáveis)
Solução
  • Limpeza das variáveis não utilizadas
  • Recodificação das variáveis
  • Renomeação das variáveis
  • Geocodificação de endereços

Após o processo de importação foi realizada a preparação dos dados, antes de proceder com a modelagem da acessibilidade espacial. Como em todo projeto de dados, esta é das mais importantes etapas, pois torna todo o trabalho posterior mais fácil e organizado. Mas é a que mais consome tempo, pois implica conhecer os dados, suas variáveis, a interrelação entre os conjuntos de dados, e definir meios de solucionar problemas que se apresentem no processo cíclico de pre-processamento.

Particularmente, as maiores dificuldades encontradas foram relacionadas à definição das UBS existentes, em função de informações desatualizadas e erros de geolocalização. Isto demandou uma série de explorações nos dados, consulta ao órgão municipal de saúde e cruzamento de informações da Prefeitura e do CNES.

Após a identificação das UBS em atividade aquelas que não existiam na base de pontos, bastou geocodificar As coordenadas espaciais foram obtidas com a Google Geocoding API por meio da função geocode do pacote ggmap, que estabelece interface com a API.

A parte computacional deste procedimento é simples:

Depois de finalizado, quando estabelecidas as linhas OD da malha territorial vigente, percebi que havia algumas UBS (15 ao todo) com erros graves de geolocalização. Foi então realizada mais uma rodada de consulta aos endereços destas UBS nos sites da Prefeitura e do CNES para conferência de endereços. Após obtenção dos endereços corretos eles foram geocodificados e os pontos incosistentes corrigidos.

Todo do processo de tratament dos dados está nos scripts presentes podem ser visto aqui. Ao fim, foram geradas as seguintes bases de dados.

rm(list = ls())

Dados pré-processados

Modelagem OD

OD euclidiano

Isto sigfica identificar as “n” UBS mais proximas de cada setor censitario. No caso estamos trabalhando com as 5 UBS mais proximas.

Levantar Rotas

O levantamento de dados de rotas consiste na consulta à Google Distance Matrix API para obtenção de dados de rotas (distancia e tempo de viagem) entre as 94765 conexões de origem e destino criadas em od_euclidiano.

Gerar os indicadores

## od_indicadores da malha vigente -------------------------------------------------------------
od_indicadores <- list(vigente = asha::asha_intersect(ubs_malha, centroides, "cnes", "cd_geocodi"))

#' O resultado tem 18948 setores, de 18953. Os cinco nao incluidos
#' nao estao dentro de nenhuma area de ubs, por diferencas na geometria.
#' Abaixo incluo diretamente as linhas com os setores ausentes. Busquei
#' no mapa as UBS correspondentes.

cd_setores <- setores$cd_geocodi[!(setores$cd_geocodi %in% od_indicadores$vigente$cd_geocodi)]
cd_ubs <- c("3121135", "4049934", "2788039", "2788217", "2788500")

od_indicadores$vigente <- od_indicadores$vigente %>%
  rbind(tibble::tibble(cd_geocodi = cd_setores, cnes = cd_ubs)) %>%
  dplyr::mutate(malha = "vigente") %>%
  tibble::as_tibble()

## od_indicadores da malha experimental --------------------------------------------------------
od_indicadores$prox <- od_viagens %>%
  dplyr::group_by(cd_geocodi) %>% 
  # ha menor tempo com distancia igual
  dplyr::filter(tempo == min(tempo)) %>%
  dplyr::filter(distancias == min(distancias)) %>%
  dplyr::ungroup() %>%
  dplyr::select(cd_geocodi, cnes) %>%
  dplyr::mutate(malha = "experimento") %>%
  tibble::as_tibble()

# Extrair setores ausentes
cd_setores <- setores$cd_geocodi[!(setores$cd_geocodi %in% od_indicadores$prox$cd_geocodi)]
od_euclidiano <- asha::asha_nn(ubs_pontos, centroides, "cnes", "cd_geocodi", 5)

# terminar od_indicadores da malha experimental
od_indicadores$prox <- od_indicadores$prox %>%
  rbind(od_euclidiano %>%
          dplyr::filter(proximidade == 1, cd_geocodi %in% cd_setores) %>%
          dplyr::select(cd_geocodi, cnes) %>%
          dplyr::mutate(malha = "experimento"))


# UNIR BASES od_indicadores ---------------------------------------------------------------
od_indicadores <- do.call(rbind, od_indicadores) %>%
  dplyr::mutate(malha = factor(malha, levels = c("vigente", "experimento")))

od_indicadores <- od_indicadores %>%
  dplyr::left_join(populacao, by = "cd_geocodi") %>%
  dplyr::left_join(prof_saude, by = "cnes") %>%
  dplyr::mutate(oportunidades = enfermeiros + medicos,
         habitantes = tidyr::replace_na(habitantes, 0)) %>%
  dplyr::group_by(cnes) %>%
  dplyr::mutate(demanda = sum(habitantes, na.rm = TRUE)) %>%
  dplyr::ungroup() %>%
  dplyr::left_join(od_viagens[, c("cnes", "cd_geocodi", "distancias", "tempo")],
            by = c("cd_geocodi", "cnes")) %>%
  dplyr::rename(distancia = distancias)

# CALCULAR INDICADORES --------------------------------------------------------

## Acessibilidade viável -------------
od_indicadores <- asha::asha_av(od_indicadores, cnes, tempo, habitantes, malha, 15)

## Acessibilidade competitiva -------------
od_indicadores <- asha::asha_ac(od_indicadores, habitantes, cnes, malha, 1000)

  1. Estudante do curso de Gestão de Políticas Públicas da EACH/USP↩

  2. Professor do curso de Gestão de Políticas Públicas da EACH/USP↩