globeMap <- read.csv("D:/DOCUMENTOS FRANCISCO/Esp. Analitica Aplicada a los Negocios/2do Semestre/Comunicacion de los Analisis/Clase 1/Actividad RStudio/Data/globeMap.csv")
str(globeMap)
## 'data.frame':    104410 obs. of  13 variables:
##  $ X                        : int  1 2 3 4 5 6 7 8 9 10 ...
##  $ long                     : num  -69.9 -69.9 -69.9 -70 -70.1 ...
##  $ lat                      : num  12.5 12.4 12.4 12.5 12.5 ...
##  $ group                    : int  1 1 1 1 1 1 1 1 1 1 ...
##  $ order                    : int  1 2 3 4 5 6 7 8 9 10 ...
##  $ region                   : chr  "Aruba" "Aruba" "Aruba" "Aruba" ...
##  $ subregion                : chr  NA NA NA NA ...
##  $ Three_Letter_Country_Code: chr  NA NA NA NA ...
##  $ Continent_Name           : chr  NA NA NA NA ...
##  $ Continent_Code           : chr  NA NA NA NA ...
##  $ Country_Name             : chr  NA NA NA NA ...
##  $ Two_Letter_Country_Code  : chr  NA NA NA NA ...
##  $ Country_Number           : int  NA NA NA NA NA NA NA NA NA NA ...

MAPA DEL MUNDO

library(ggplot2)
myMap <- ggplot(globeMap, aes(x=long, y=lat, group = as.factor(group)))
myMap <- myMap + geom_polygon(aes(fill=Continent_Name, color=Continent_Name))
myMap <- myMap + scale_fill_viridis_d()
myMap <- myMap + scale_colour_manual(values = c("#999999", "#E69F00", "#56B4E9", "#009E73", "#F0E442", "#0072B2", "#D55E00", "#CC79A7"))
myMap <- myMap + labs(title = "Countries outlined by region", subtitle = "Using a mercator projection", caption = "Dataset retrieved from native ggplot map", x = "Longitude", y= "Latitude")
myMap <- myMap + theme_bw()
myMap

MAPA MUNDO CON IHDI

IHDI <- read.csv2("D:/DOCUMENTOS FRANCISCO/Esp. Analitica Aplicada a los Negocios/2do Semestre/Comunicacion de los Analisis/Clase 1/Actividad RStudio/Data/Coefficient_of_human_inequality.csv", na.strings="")
IHDI_2017 <- subset(IHDI, IHDI$year == "X2017")
colnames(IHDI_2017)[5] <- "Three_Letter_Country_Code"
library(dplyr)
mapData <- left_join(globeMap,IHDI_2017, by =("Three_Letter_Country_Code"))
myMap <- ggplot(mapData, aes(x=globeMap$long, y=globeMap$lat, group = as.factor(globeMap$group)))
myMap <- myMap + geom_polygon(aes(fill = value))
myMap <- myMap + scale_fill_viridis_c(option="magma")
myMap <- myMap + labs(title = "World Inequality-Adjusted Development Index (IHDI) for 2017", subtitle = "Average of education, health and income inequality indexes", caption = "Grey territories have no index. Source UN", x = "Longitude", y= "Latitude")
myMap <- myMap + theme_bw()
myMap