[ad_1]
[Esteartigofoipublicadopelaprimeiravezem[Thisarticlewasfirstpublishedon r.iresmi.net, e gentilmente contribuiu para R-blogueiros]. (Você pode relatar um problema sobre o conteúdo desta página aqui)
Deseja compartilhar seu conteúdo com R-blogueiros? clique aqui se você tiver um blog ou aqui se não tiver.
O pacote {geofacet} permite «organizar uma sequência de gráficos de dados para diferentes entidades geográficas em uma grade que se esforça para preservar parte da orientação geográfica original das entidades».
Como na postagem anterior, é interessante se você visualizar cada entidade como uma unidade e não se importar com seu tamanho ou peso real, e não desejar gastar muito tempo encontrando manualmente a melhor grade.
Usaremos novamente o mesmo conjunto de dados COVID-19. Adicionamos manualmente os departamentos no exterior assim que encontramos a grade correta (tentando sementes diferentes) e ajustamos a posição da Córsega.

# packages ---------------------------------------------------------------- library(tidyverse) library(httr) library(fs) library(sf) library(readxl) library(janitor) library(glue) library(geofacet) # also install ragg # sources ----------------------------------------------------------------- # https://www.data.gouv.fr/fr/datasets/donnees-hospitalieres-relatives-a-lepidemie-de-covid-19/ fichier_covid % stop_for_status() } if (!file_exists(fichier_pop)) { GET(url_donnees_pop, progress(), write_disk(fichier_pop)) %>% stop_for_status() } covid % clean_names() # adminexpress prétéléchargé dep % clean_names() %>% st_set_crs(2154) # construction de la grille ---------------------------------------- grid_fr % select(insee_dep, nom_dep) %>% grid_auto(names = "nom_dep", codes = "insee_dep", seed = 4) %>% add_row(row = 8, col = 1, name_nom_dep = "Guadeloupe", code_insee_dep = "971") %>% add_row(row = 9, col = 1, name_nom_dep = "Martinique", code_insee_dep = "972") %>% add_row(row = 10, col = 1, name_nom_dep = "Guyane", code_insee_dep = "973") %>% add_row(row = 7, col = 13, name_nom_dep = "Mayotte", code_insee_dep = "976") %>% add_row(row = 8, col = 13, name_nom_dep = "La Réunion", code_insee_dep = "974") grid_fr[grid_fr$code_insee_dep %in% c("2A", "2B"), "col"] % filter(sexe == 0) %>% rename(deces = dc, reanim = rea, hospit = hosp) %>% left_join(pop, by = c("dep" = "x1")) %>% mutate(incidence = deces / x2020_p * 100000) %>% rename(insee_dep = dep) %>% left_join(grid_fr %>% select(nom_dep = name_nom_dep, insee_dep = code_insee_dep)) %>% drop_na(insee_dep) %>% ggplot(aes(jour, incidence)) + geom_area() + facet_geo(~ nom_dep, grid = grid_fr) + labs(title = "Mortalité", subtitle = "COVID-19 - France", x = "date", y = "décès pourn100 000 hab.", caption = glue("http://r.iresmi.net/ndonnées SPF {Sys.Date()}")) + theme_minimal() + theme(strip.text = element_text(hjust = 0, size = 7)) ggsave(glue("resultats/covid_fr_mortalite_geofacette_{Sys.Date()}.png"), width = 25, height = 20, units = "cm", scaling = .8, res = 300, device = ragg::agg_png)
Relacionado
Se você chegou até aqui, por que não inscreva-se para atualizações do site? Escolha o seu sabor: e-mail, twitter, RSS ou facebook …
[ad_2]