-
-
Save ArnulfoPerez/e315135bd34df2521c2472da4ae24aca to your computer and use it in GitHub Desktop.
Mapa de casos COVID-19
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| ## Auto-Install the following packages | |
| .packs <- c("tidyverse", "lubridate", "ggrepel", "viridis", "scales") | |
| .success <- suppressWarnings(sapply(.packs, require, character.only = TRUE)) | |
| if (length(names(.success)[!.success])) { | |
| install.packages(names(.success)[!.success]) | |
| sapply(names(.success)[!.success], require, character.only = TRUE) | |
| } | |
| if (!require(mxmaps)) | |
| devtools::install_github("diegovalle/mxmaps") | |
| library(mxmaps) | |
| options(stringsAsFactors = FALSE) | |
| temp <- tempfile() | |
| download.file("http://187.191.75.115/gobmx/salud/datos_abiertos/datos_abiertos_covid19.zip",temp) | |
| files <- unzip(temp, list = TRUE) | |
| df <- read_csv(unz(temp, arrange(files, Length)[1,]$Name)) | |
| unlink(temp) | |
| pos <- filter(df, RESULTADO == 1) | |
| pos$region <- str_mxmunicipio(pos$ENTIDAD_RES, pos$MUNICIPIO_RES) | |
| muns <- pos %>% | |
| group_by(region) %>% | |
| tally() | |
| muns <- left_join(muns, df_mxmunicipio, by = "region") | |
| muns$value2 <- muns$n / muns$pop * 10^5 | |
| muns$value <- if_else(muns$value2 > 400, 400, muns$value2) | |
| muns$name <- paste(muns$state_name, muns$municipio_name) | |
| cities <- filter(muns, name %in% c("Coahuila Monclova", | |
| "Baja California Sur Los Cabos", | |
| "Ciudad de México Cuajimalpa de Morelos", | |
| "Quintana Roo Benito Juárez", | |
| "Sinaloa Culiacán", | |
| "Tabasco Centro", | |
| "Baja California Mexicali", | |
| "Yucatán Mérida", | |
| "Puebla Puebla", | |
| "Sonora Sáric", | |
| "Guerrero Acapulco de Juárez", | |
| "Chihuahua Juárez", | |
| "Coahuila Piedras Negras", | |
| "Michoacán Lázaro Cárdenas", | |
| "Oaxaca Santa María Huatulco", | |
| "Colima Manzanillo", | |
| "Jalisco Puerto Vallarta")) | |
| cities$municipio_name <- str_replace(cities$municipio_name, | |
| "Cuajimalpa de Morelos", | |
| "Cuajimalpa") | |
| cities$municipio_name <- str_replace(cities$municipio_name, | |
| "Santa María Huatulco", | |
| "Huatulco") | |
| cities$municipio_name <- str_replace(cities$municipio_name, | |
| "Benito Juárez", | |
| "Cancún") | |
| cities$municipio_name <- str_replace(cities$municipio_name, | |
| "Juárez", | |
| "Ciudad Juárez") | |
| cities$municipio_name <- str_replace(cities$municipio_name, | |
| "Acapulco de Juárez", | |
| "Acapulco") | |
| cities$municipio_name <- str_replace(cities$municipio_name, | |
| "Centro", | |
| "Villahermosa") | |
| cities$group <- NA | |
| mxmunicipio_choropleth(muns, | |
| num_colors = 1, | |
| title = "Mapa de casos COVID-19 confirmados, por municipio de residencia (14 de abril del 2020)", | |
| legend = "tasa por\n100 mil\nhabitantes") + | |
| scale_fill_viridis("tasa por\n100 mil\nhabitantes", | |
| trans = scales::pseudo_log_trans(sigma = 0.001)) + | |
| geom_label_repel(data = cities, aes(long, lat, label = municipio_name), | |
| size = 3, | |
| force = .1, alpha = .8, | |
| box.padding = 3.3, label.padding = 0.18) + | |
| theme(legend.key.size = unit(2, "cm")) + | |
| theme(plot.title = element_text(size=32)) | |
| #ggsave("graphs/map_covid.png", dpi = 100, width = 19, height = 14) | |
| mxmunicipio_choropleth(muns, | |
| num_colors = 1, | |
| title = "Mapa de casos COVID-19 confirmados en el la región\ncentral de México, por municipio de residencia\n(14 de abril del 2020)", | |
| legend = "tasa por\n100 mil\nhabitantes", | |
| zoom = subset(df_mxmunicipio, state_name %in% c("Ciudad de México", | |
| "Puebla", | |
| "Morelos", | |
| "México", | |
| "Hidalgo", | |
| "Tlaxcala"))$region) + | |
| scale_fill_viridis("tasa por\n100 mil\nhabitantes", | |
| trans = scales::pseudo_log_trans(sigma = 0.001)) + | |
| theme(legend.key.size = unit(2, "cm")) + | |
| theme(plot.title = element_text(size=30)) | |
| #ggsave("graphs/map_centro_covid.png", dpi = 100, width = 16, height = 11) | |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment