### Lab 06 

####################################################################
# Project 1: Aiming to create a choropleth map of German 
# election results of the party list vote for SPD by state 
####################################################################
library(dplyr)

## Step 0: Collect Election Results 
library(rvest)
url <- 'https://en.wikipedia.org/wiki/Results_of_the_2021_German_federal_election'

web_data <- read_html(url)

elec_results <- web_data %>% 
  html_table() %>%
  .[[10]]

## Step 0.5: Data Cleaning 
elec_results %>%
  View()

### Remove Non Data Rows 
names(elec_results) <- elec_results[1,]
elec_results <- elec_results %>% 
  slice(3:18)

### Convert data type to numeric 
elec_results$SPD <- as.numeric(elec_results$SPD)


## Step 1: Load Shape file (using sf package)
library(sf)
de_shp <- sf::read_sf('Germany_shapefile/de_admin1.shp')

de_shp %>% 
  glimpse()

de_shp$geometry[1]

## Step 1.5: Check the map
library(ggplot2)

ggplot(de_shp) + 
  geom_sf()

## Step 2: Join Data Sets
table(de_shp$NAME_1)
table(elec_results$State)

elec_results <- elec_results %>% 
  mutate(id = case_when(
    State == 'Bavaria' ~ 'Bayern',
    State == 'Hesse' ~ 'Hessen',
    State == 'Lower Saxony' ~ 'Niedersachsen', 
    State == 'North Rhine-Westphalia' ~ 'Nordrhein-Westfalen',
    State == 'Rhineland-Palatinate' ~ 'Rheinland-Pfalz',
    State == 'Saxony' ~ 'Sachsen', 
    State == 'Saxony-Anhalt' ~ 'Sachsen-Anhalt',
    State == 'Thuringia' ~ 'Thüringen',
    TRUE ~ State
  ))

de_shp <- de_shp %>% 
  left_join(elec_results, by = c('NAME_1' = 'id'))

de_shp %>% 
  glimpse()

## Step 3: Plot 
ggplot(de_shp) + 
  geom_sf(aes(fill = SPD))


## Step 4: Adjust 
ggplot(de_shp) + 
  geom_sf(aes(fill = SPD)) + 
  coord_sf() + 
  scale_fill_gradient(low = 'white', 
                      high = '#EB001F') + 
  theme_void() + 
  theme(legend.position = 'bottom')

####################################################################
# Project 2: Using sf, let's create a choropleth map of 
# disposable income of private households in Europe (NUTS-2)
####################################################################
library(sf)

#Load data
eu_data <- read.csv('https://raw.githubusercontent.com/apodkul/ppol6803_03/main/Data/eu_disp_income2016.csv')
  
#Load mapping file   
eu_shape <- read_sf('EU/nuts2.shp')

eu_shape %>% 
  left_join(eu_data, by = c('NUTS_ID' = 'geo')) %>%
  ggplot() + 
  geom_sf()

eu_shape %>% 
  left_join(eu_data, by = c('NUTS_ID' = 'geo')) %>%
  ggplot(aes(fill = values, color = values)) + 
  geom_sf() 

eu_shape2 <- eu_shape %>% 
  left_join(eu_data, by = c('NUTS_ID' = 'geo'))

eu_crop2 <- st_crop(eu_shape2, 
                   c(xmin = -20, 
                     xmax = 40, 
                     ymin = 30, 
                     ymax = 75)
                   )
eu_crop2 %>% 
  ggplot(aes(fill = values, color = values)) + 
  geom_sf() +
  coord_sf() + 
  theme_void()