(Acest articol a fost publicat pentru prima dată pe r.iresmi.netși cu amabilitate a contribuit la R-bloggeri). (Puteți raporta problema legată de conținutul acestei pagini aici)
Doriți să vă distribuiți conținutul pe R-bloggeri? dați clic aici dacă aveți un blog, sau aici dacă nu aveți.
Ziua 3 din 30DayMapChallenge: « Poligoane » (anterior).
O hartă clasică de coropletă a lumii. Vom folosi Baza de date comparativă privind pedeapsa cu moartea (Anckar și Denk 2024) disponibil pe https://dataverse.harvard.edu/dataset.xhtml?persistentId=doi:10.7910/DVN/LI3WYK și harta de bază de la {rnaturalearth}.
library(dplyr) library(tibble) library(readr) library(glue) library(readxl) library(janitor) library(ggplot2) library(scales) library(ggspatial) library(sf) library(rnaturalearth) library(countrycode)
Date
CDPD nu folosește coduri de țară ISO A3, așa că trebuie să le traducem din Corelații de război coduri cu {countrycode}.
codes <- tribble(
~deathpenalty, ~lib_deathpenalty, ~color,
0, "Abolished", "lightskyblue",
1, "Abolished for ordinary crimes only", "khaki1",
2, "Abolished for ordinary crimes only, but used during the last 10 years.", "gold1",
3, "Abolished in practice", "indianred1",
4, "Retained", "indianred4")
if (!file.exists("cdpd.rds")) {
download.file("https://dataverse.harvard.edu/api/access/datafile/10251764",
"cdpd.xlsx")
cdpd <- read_xlsx("cdpd.xlsx") |>
clean_names() |>
filter(year == 2022,
country != "MICRONESIA") |>
mutate(iso = countrycode(cowcode, "cown", "iso3c",
custom_match = c("342" = "SRB",
"348" = "MNE",
"818" = "VNM"))) |>
left_join(codes, join_by(deathpenalty)) |>
select(country, iso, lib_deathpenalty) |>
write_rds("cdpd.rds")
} else {
cdpd <- read_rds("cdpd.rds")
}
abolished <- nrow(filter(cdpd, lib_deathpenalty == "Abolished")) / nrow(cdpd)
Date de hartă
country_map <- ne_countries(scale = 110, returnclass = "sf") |> select(sovereignt, sov_a3, admin, adm0_a3)
Hartă
pal_dp <- codes |>
select(lib_deathpenalty, color) |>
deframe()
cdpd_map <- country_map |>
left_join(cdpd, join_by(adm0_a3 == iso))
cdpd_map |>
ggplot() +
geom_sf(aes(fill = lib_deathpenalty)) +
scale_fill_manual(values = pal_dp, na.value = "grey") +
coord_sf(crs = "+proj=eqearth") +
labs(title = glue("Death penalty is still retained in {percent(1 - abolished)} of the countries"),
subtitle = "2022",
fill = "Death penaltynstatus",
caption = glue("data: Comparative Death Penalty Database
map data: NaturalErth
r.iresmi.net - {Sys.Date()}")) +
theme_void() +
theme(plot.background = element_rect(color = NA,
fill = "white"),
legend.position = "bottom",
legend.text = element_text(size = 8),
plot.caption = element_text(size = 5,
color = "darkgrey"))
