install.packages(c("cowplot", "googleway", "ggplot2", "ggrepel", 
                   "ggspatial", "libwgeom", "sf", "rnaturalearth", "rnaturalearthdata", "rgeos", "haven", "readstata13", "dplyr"))
                 
                 library("haven")
                 library("readstata13")
                 library("dplyr")
                 library("tidyr")
                 library("ggplot2")
                 theme_set(theme_bw())
                 library("sf")
                 library("rnaturalearth")
                 library("rnaturalearthdata")
                 library("rgeos")
                 data = read.dta13(file.choose(),nonint.factors=T)
                 t = read.dta13(file.choose(),nonint.factors=T)
                 data = merge(data, t, by = 'location_name')
                 data$Location = data$location_name
                 data$QCI = data$pca
                 data$iso = data$iso3_countries
                 data2 = filter(data, year == 2017, sex_name %in% c('Male', 'Female'), age_name == "Age-standardized", type %in% c("Country"))
                 data3 = select(data2, location_name, iso, sex_name, QCI)
                 data4 = spread(data3, sex_name, QCI)
                 data4$GDR = data4$Female / data4$Male
                 world <- ne_countries(scale = "medium", returnclass = "sf")
                 world$iso = world$iso_a3
                 world = merge(world, data2, by = 'iso')
                 ggplot(data = world) +
                   geom_sf(aes(fill = GDR)) +
                   scale_fill_gradient2(high = 'blue', mid = 'grey24', low = 'red')
                 data3 = select(data2, Location, QCI)