Вы можете использовать пакеты sp или sf.
Скачать "gadm36_AFG_1_sp.rds" с gadm.org (уровень 1 = провинция).
AFG_adm1 <- readRDS("gadm36_AFG_1_sp.rds")
library(sp)
op <- par(mar=c(0,0,0,0))
plot(AFG_adm1)
Вы можете объединить эту карту со своими данными, используя названия провинций (NAME_1).
names(AFG_adm1@data)
#[1] "GID_0" "NAME_0" "GID_1" "NAME_1"
#[5] "VARNAME_1" "NL_NAME_1" "TYPE_1" "ENGTYPE_1"
#[9] "CC_1" "HASC_1"
Но будьте осторожны с написанием.
AFG_adm1$NAME_1
[1] "Badakhshan" "Badghis" "Baghlan" "Balkh"
[5] "Bamyan" "Daykundi" "Farah" "Faryab"
[9] "Ghazni" "Ghor" "Hilmand" "Hirat"
[13] "Jawzjan" "Kabul" "Kandahar" "Kapisa"
[17] "Khost" "Kunar" "Kunduz" "Laghman"
[21] "Logar" "Nangarhar" "Nimroz" "Nuristan"
[25] "Paktika" "Paktya" "Panjshir" "Parwan"
[29] "Samangan" "Sari Pul" "Takhar" "Uruzgan"
[33] "Wardak" "Zabul"