Я хотел бы видеть результаты моих данных на региональной карте Турции. Но я застреваю, когда пытаюсь увидеть распределение своих данных и результатов модели на карте.
Я знаю, что мне нужно написать еще несколько кодов о слиянии реальных данных и шейп-файла, я попробовал вот так нижеcode.
what I wrote until now as below;
temp <- tempfile(fileext = ".zip")
download.file("http://ec.europa.eu/eurostat/cache/GISCO/geodatafiles/NUTS_2013_01M_SH.zip", temp)
unzip(temp)
#load the data and filter it to Turkey and NUTS2 level
EU_NUTS = readOGR(dsn = "./NUTS_2013_01M_SH/data", layer = "NUTS_RG_01M_2013")
map_nuts2 <- subset(EU_NUTS, STAT_LEVL_ == 2) # set NUTS level
country <- substring(as.character(map_nuts2$NUTS_ID),1, 2)
map <- c("TR") # limit it to Turkey
map_nuts2a <- map_nuts2[country %in% map,]
fulldata = join (x = pisa, y = map_nuts2a, by = "region", type = "full") "ОНО ДАЕТ НАСТРОЕНИЕ ФИТИНГА"
#plot it
plot(map_nuts2a, col = colorRampPalette(c("white", "red","blue","green","black","orange",
"navy","yellow","pink","brown","gray","violet"))(nrow(map_nuts2a@data)))