Mapa de Pernambuco Segundo Grau de Instrução

Pacotes Utilizados:

library("Rcpp")
library(leaflet)
library(tidyverse)
library(stringr)
library(dplyr)
library(readr)
library(readxl)
library(rgdal)

Mapa de PE com percentual de pessoas sem instrução em cada município

setwd("C:\\Users\\User\\Desktop\\Mapas")

pessoas <- read_csv("C:\\Users\\User\\Desktop\\Mapas\\Amostra_Pessoas_26.txt", 
                    col_names = FALSE)

municipios <- read_excel("C:\\Users\\User\\Desktop\\Mapas\\municipios.xls", skip=2)



instrucao <- pessoas %>% transmute(COD_Municipio = str_sub(string = X1, start = 1, end = 7),
                                   COD_Instrucao = str_sub(string = X1, start = 158, end = 158),
                                   Instrucao = if_else(COD_Instrucao == "1", "Sem instrução e fundamental incompleto",
                                                       if_else(COD_Instrucao == "2","Fundamental completo e médio incompleto",
                                                               if_else(COD_Instrucao == "3", "Médio completo e superior incompleto",
                                                                       if_else(COD_Instrucao == "4", "Superior completo", "Não determinado")))))


dados <- instrucao %>% left_join(municipios, by = c("COD_Municipio" = "Município"))



shp <- readOGR("C:\\Users\\User\\Desktop\\Mapas", "26MUE250GC_SIR", stringsAsFactors=FALSE, encoding="UTF-8")

dados_aux = dados %>% group_by(COD_Municipio) %>% summarise(count = sum(COD_Instrucao==1)/(sum(COD_Instrucao==1)+sum(COD_Instrucao!=1)))
dados_aux

dados_completos <- merge(shp,dados_aux, by.x = c("CD_GEOCMU"), by.y = c("COD_Municipio"))



qpal <- colorQuantile("Reds", dados_completos$count, n = 4)
mapa <- leaflet() %>% addTiles() %>% addPolygons(data=dados_completos,stroke = FALSE, fillColor = ~qpal(dados_completos$count), smoothFactor = 0.2,fillOpacity = 1) %>%  addLegend("bottomright", pal = qpal, values = dados_completos$count,
                                                                                                                                                                                  title = "% sem instrução",
                                                                                                                                                                                  opacity = 1)

mapa

Mapa de Pernambuco Segundo Grau de Instrução