Check GDS R stack¶

library(arm)

library(feather)

library(geojsonio)

library(ggmap)

library(gstat)

library(hexbin)

library(igraph)

library(kableExtra)

library(knitr)

library(lidR)

library(lme4)

library(mapdata)

library(maptools)

library(mapview)

library(ncdf4)

library(nlme)

library(plyr)

library(proj4)

library(RColorBrewer)

library(RandomFields)

library(RNetCDF)

library(randomForest)

library(raster)

library(RCurl)

library(reshape2)

library(rgdal)

library(rgeos)

library(rmarkdown)

library(RODBC)

library(RSQLite)

library(sf)

library(shiny)

library(sp)

library(spacetime)

library(spatstat)

library(spdep)

library(splancs)

# Tidyverse
library(ggplot2)
library(tibble)
library(tidyr)
library(readr)
library(purrr)
library(stringr)
library(forcats)

library(tmap)

library(TraMineR)

library(tufte)

db = st_read('https://ndownloader.figshare.com/files/20232174')

plot(db$geom)

data(World, metro)
metro$growth <- (metro$pop2020 - metro$pop2010) / (metro$pop2010 * 10) * 100

mapWorld <- tm_shape(World) +
    tm_polygons("income_grp", palette="-Blues", contrast=.7, id="name", title="Income group") +
    tm_shape(metro) +
    tm_bubbles("pop2010", col = "growth", 
               border.col = "black", border.alpha = .5, 
               style="fixed", breaks=c(-Inf, seq(0, 6, by=2), Inf),
               palette="-RdYlBu", contrast=1, 
               title.size="Metro population", 
               title.col="Growth rate (%)", id="name") + 
    tm_style_gray() + tm_format_World()

mapWorld