Мне нужно сопоставить базу данных, составленную из 81 страны, с картой мира на т.
В настоящее время я использую этот пакет
library("rnaturalearth")
library("rnaturalearthdata")
theme_set(theme_bw())
world <- ne_countries(scale = "medium", type="countries", returnclass = c("sf"))
class(world)
worldmap<-as(world,"Spatial")
class(worldmap)
plot(worldmap,cex=.01)
worldmap@data <- data.frame(database,worldmap@data
[match(database[, "code"],worldmap@data[, "gu_a3"]),])
, но когда я пытаюсь сопоставить данные с карта, поскольку в моей базе данных нет полного списка стран мира, она не совпадает корректно.
База данных имеет следующую форму
# A tibble: 6 x 12
code country sales2010 gdp gdppc population Export Import tradecost skills investmentcost
<chr> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
1 ARG Argent… 44287 4.24e11 10386. 40788453 8.11e10 6.88e10 2.83 9.48 0.509
2 ARM Armenia 4 9.26e 9 3218. 2877319 2.21e 9 4.54e 9 1.37 10.9 0.333
3 AUS Austra… 185329 1.15e12 52022. 22031750 2.60e11 2.54e11 2.23 11.8 0.917
4 AUT Austria 20222 3.92e11 46858. 8363404 1.98e11 1.86e11 1.02 9.89 0.944
5 AZE Azerba… 2019 5.29e10 5843. 9054332 2.82e10 1.02e10 1.38 10.5 0.417
6 BEL Belgium 130775 4.81e11 44142. 10895586 3.69e11 3.64e11 0.656 10.8 0.833
# … with 1 more variable: distance <dbl>
Я также пытался создать базу данных со всеми странами из списка на карте и оставить НС там, где у меня нет данных, но это создает много проблем с пространственной эконометрией c оценка.
Знаете ли вы какие-либо лучшее решение? Я могу использовать любой другой пакет карт, если могу преобразовать данные в формат "sf"
Спасибо!