Untitled

```{r}
Stations = read.csv(file = "Stations/Combined.csv")
LSOA21 = st_read("Shapefiles/21LSOA/LSOA_(Dec_2021)_Boundaries_Full_Clipped_EW_(BFC).shp")
LSOA0121 = st_read("Shapefiles/0121/0121.shp")
ETH21 = read.csv(file = "Ethnicity/Manchester/2021ETH.csv")
OSGB36 = "+init=epsg:27700 +towgs84=375,-111,431,0,0,0,0"
```

```{r}
MANLSNM21 =  read.csv(file = "Shapefiles/MANLSNM/MANLSNM21..shp")
```
```{r}
LSOA21MAN = filter(LSOA21, LSOA21NM == c("Bury", "Bolton", "Tameside", "Trafford", "Manchester", "Oldham", "Rochdale", "Salford"))
```
```{r}
LSOA21NM = 
LSOA21MAN = gather(LSOA21NM) %>% 
      separate(LSOA21NM, into = c("Area", "Short.Code"))
```
%>%
      group_by(key) %>%
      mutate(i1 = row_number()) %>% 
      ungroup() %>% 
      select(-key) %>% 
      spread(header, value) %>%
      select(-i1)
```

LSOA21MAN = filter(LSOA21MAN, Area == c("Bury", "Bolton", "Tameside", "Trafford", "Manchester", "Oldham", "Rochdale", "Salford"))


```

```{r}
LSOA21MAN= separate(LSOA21, LSOA21NM, into = c("Area", "Short.Code"))
```
LSOA21MAN = filter(LSOA21MAN, Area == c("Bury", "Bolton", "Tameside", "Trafford", "Manchester", "Oldham", "Rochdale", "Salford"))

21Code

21Code = read.csv(file = "LSOACode/21Code.csv")
separate(LSOA21NM, into = c("Area", "Short.Code"))
ETH21
LSOA21

LSOA21CD
Code



flights2 |>
  left_join(airlines)
#> Joining with `by = join_by(carrier)`

DETHNames = DETHNames |>
  left_join(LSOA21NM,join_by("Code"))


facet_wrap(~variety)