Scotland.rmd

---
title: "Scotland UK IMD Datazone Maps"
author: "Will Ball"
date: "16/05/2020"
output: 
  html_document:
    code_folding: hide
    df_print: paged
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

```{r, warning=FALSE, message=FALSE}
library(tidyverse)
library(curl)
library(sf)
library(rmapshaper)
library(viridis)
```


## IMD data

Read in straight from source and rename `area_code` to `dz` for joining

```{r, message=FALSE, warning=FALSE}
imd <- read_tsv("https://data.bris.ac.uk/datasets/1ef3q32gybk001v77c1ifmty7x/uk_imd_scores_data.txt") %>% 
  rename(dz = area_code)
imd
```

## Shapefile

Zip file this time so needs to be read in to a `tempfile()` and then unzipped into another.

Also rename `dz` 

```{r}
temp_1 <- tempfile()
temp_2 <- tempfile()
source <- "http://sedsh127.sedsh.gov.uk/Atom_data/ScotGov/ZippedShapefiles/SG_DataZoneBdry_2001.zip"
temp_1 <- curl_download(url = source, destfile = temp_1, quiet = FALSE)
unzip(temp_1, exdir = temp_2)
dz_shp <- read_sf(file.path(temp_2,"SG_DataZone_Bdry_2001.shp"))
dz_shp <- ms_simplify(dz_shp) %>% 
  select(-DZ_NAME, -DZ_GAELIC, -STDAREA_HA, -Shape_Leng, -Shape_Area) %>% 
  rename(dz = DZ_CODE)
dz_shp
```


## DZ lookup

Using the PHS opendata. 

```{r, warning=FALSE, message=FALSE}
lookup <- read_csv("https://www.opendata.nhs.scot/dataset/9f942fdb-e59e-44f5-b534-d6e17229cc7b/resource/e92d19d4-ced7-40c8-b628-e28e4528fc41/download/dz2001_codes_and_labels_21042020.csv") %>% 
  select(DataZone, CA, CAName, HB, HBName) %>% 
  rename(dz = DataZone)
lookup
```


## Join

Join and then coerce to `sf` object

```{r}
join <- left_join(lookup, imd) %>% 
  left_join(dz_shp) %>% 
  st_as_sf(.)
join
```


## Plot

Changes as required:
filter - by CAName/ctry/HBName
title - to specify area selected
ggsave - file name

```{r, fig.width=9, fig.height=6}
theme_set(theme_minimal())
join %>% 
  filter(ctry == "Scotland") %>% 
  ggplot() +                                      
  geom_sf(aes(fill = uk_imd_england_quintile),
          color = NA)+
  scale_fill_viridis(option = "plasma", discrete = T) +
  scale_color_viridis() +
  labs(x = NULL, y = NULL,                                                          
       title = "UK IMD Quintiles for Scotland by Data Zone (2001)",              
       subtitle = "Sources: Abel et al (2016) BMJ Open and Office for National Statistics under 
       Open Government License v3.0",                                                   
       caption = "Plot by @WillBall12 | Contains OS data © Crown copyright and database right (2020)") +  
  theme(axis.line=element_blank(), 
        axis.ticks=element_blank(), 
        axis.text=element_blank(),
        axis.title=element_blank(),
        panel.grid = element_blank(),
        plot.caption.position = 'plot',
        legend.title = element_blank()) 

ggsave("Scotland.png")
```