Updated attraction and mode choice

Author: CorinStaves

.gitignore

@@ -3,3 +3,4 @@
 .RData
 .Ruserdata
 /data
+/result/local

script/Manchester/destinationChoice/observedAttraction.R

@@ -40,33 +40,57 @@ test <- trips %>% count(EndPurpose, purpose) %>% arrange(desc(n))
 # Remove home/work/business destinations
 trips <- filter(trips, purpose != "work" & purpose != "home" & purpose != "education" & purpose != "business")
 
+# Add MSOA/LSOA
+OA_lookup <- readr::read_csv(paste0("data/Manchester/OA_lookup.csv"), col_select = c("OA11CD","LSOA11CD","MSOA11CD"))
+
+trips <- left_join(trips,OA_lookup,by = c("EndOutputArea" = "OA11CD")) %>% rename(EndMSOA = MSOA11CD, EndLSOA = LSOA11CD)
+
 # Get Attraction of each OA
-rawAttractions <- trips %>%
-  count(EndOutputArea,EndPurpose) %>%
-  pivot_wider(names_from = EndPurpose,
-              values_from = n) %>%
-  filter(EndOutputArea != "         ")
 
-zoneAttractions <- trips %>%
-  count(EndOutputArea,purpose) %>%
-  pivot_wider(names_from = purpose,
-              values_from = n) %>%
-  filter(EndOutputArea != "         ")
+getAttractions <- function(regionType) {
+  rawAttractions <- trips %>%
+    count(!!enquo(regionType),EndPurpose) %>%
+    pivot_wider(names_from = EndPurpose,
+                values_from = n) %>%
+    filter(!!enquo(regionType) != "         ")
+  
+  zoneAttractions <- trips %>%
+    count(!!enquo(regionType),purpose) %>%
+    pivot_wider(names_from = purpose,
+                values_from = n) %>%
+    filter(!!enquo(regionType) != "         ")
+  
+  return(full_join(rawAttractions,zoneAttractions))
+}
+
+attractions_OA <- getAttractions(EndOutputArea)
+attractions_LSOA <- getAttractions(EndLSOA)
+attractions_MSOA <- getAttractions(EndMSOA)
+
 
-attractions <- full_join(rawAttractions,zoneAttractions)
 
 # Attach attractions to OAs
 
-all_OAs <- sf::st_read("~/Documents/TfGM/destinationChoice/gm_oa.gpkg")
-OA_lookup <- readr::read_csv(paste0("data/Manchester/OA_lookup.csv"), col_select = c("OA11CD","LAD17CD","LAD17NM"))
+OAs <- sf::st_read("~/Documents/manchester/zones/gm_oa.shp")
+LSOAs <- sf::st_read("~/Documents/manchester/zones/LSOA_studyArea.shp")
+MSOAs <- sf::st_read("~/Documents/manchester/zones/MSOA_studyArea.shp")
 
 
-gm_districts <- c("Bolton","Bury","Manchester","Oldham","Rochdale","Salford","Stockport","Tameside","Trafford","Wigan")
-gm_oa_list <- OA_lookup %>% filter(LAD17NM %in% gm_districts) %>% pull(OA11CD)
+attr_OAs <- OAs %>% select(-fid) %>%
+  left_join(attractions_OA, by = c("OA11CD" = "EndOutputArea")) %>%
+  mutate(across(where(is.integer),function(x) replace_na(x,0)))
+
+attr_LSOAs <- LSOAs %>% 
+  left_join(attractions_LSOA, by = c("LSOA11CD" = "EndLSOA")) %>%
+  mutate(across(where(is.integer),function(x) replace_na(x,0)))
 
-gm_OAs <- all_OAs %>% 
-  filter(geo_code %in% gm_oa_list) %>%
-  left_join(attractions, by = c("geo_code" = "EndOutputArea")) %>%
+attr_MSOAs <- MSOAs %>% 
+  left_join(attractions_MSOA, by = c("MSOA11CD" = "EndMSOA")) %>%
   mutate(across(where(is.integer),function(x) replace_na(x,0)))
 
-sf::st_write(gm_OAs,"result/attractions.gpkg",delete_layer = TRUE)
+
+# Write
+sf::st_write(attr_OAs,"result/attractions_OA.gpkg",delete_layer = TRUE)
+sf::st_write(attr_LSOAs,"result/attractions_LSOA.gpkg",delete_layer = TRUE)
+sf::st_write(attr_MSOAs,"result/attractions_MSOA.gpkg",delete_layer = TRUE)
+

script/Manchester/modeChoice/cycleSpatialVariation.R

@@ -1,31 +1,50 @@
 setwd("~/Documents/JIBE/travelDiaryProcessing")
 source("script/Manchester/read.R")
 
-# Spatial variation based on origin
-cycleTripsByStartOA <- TRADS$trips %>% 
-  filter(t.startOA != "         ") %>%
-  group_by(t.startOA) %>% 
-  summarise(cycleStart = sum(t.m_main == "Bicycle"),
-            cycleStartShare = cycleStart / n())
-
-cycleTripsByEndOA <- TRADS$trips %>% 
-  filter(t.endOA != "         ") %>%
-  group_by(t.endOA) %>% 
-  summarise(cycleEnd = sum(t.m_main == "Bicycle"),
-            cycleEndShare = cycleEnd / n())
+OAs <- sf::st_read("~/Documents/manchester/zones/gm_oa.shp")
+LSOAs <- sf::st_read("~/Documents/manchester/zones/LSOA_studyArea.shp")
+MSOAs <- sf::st_read("~/Documents/manchester/zones/MSOA_studyArea.shp")
 
-cycleTripsByHomeOA <- TRADS$trips %>%
-  left_join(TRADS$households) %>%
-  group_by(hh.OA) %>%
-  summarise(cycleHome = sum(t.m_main == "Bicycle"),
-            cycleHomeShare = cycleHome / n())
+trips <- TRADS$trips %>% left_join(select(TRADS$households,hh.id,hh.OA))
 
-OAs <- sf::read_sf("~/Documents/manchester/JIBE/accessibility/OAs/gm_oa.shp")
+# Spatial variation based on origin
+getTripsAndShares <- function(var,suffix,areaType) {
+  lookup <- OA_lookup[,c("OA11CD",areaType)]
+  colnames(lookup)[1] <- var
+  
+  result <- trips %>% 
+    left_join(lookup) %>% 
+    group_by_(areaType) %>%
+    summarise(cycleCount = sum(t.m_main == "Bicycle"),
+              walkCount = sum(t.m_main == "Walk"),
+              trips = n(),
+              cycleShare = cycleCount / n(),
+              walkShare = walkCount / n())
+  
+  colnames(result)[2:6] <- paste(colnames(result)[2:6],suffix,sep = "_")
+  
+  return(result)
+}
 
 OAs <- OAs %>%
-  left_join(cycleTripsByStartOA,by = c("geo_code" = "t.startOA")) %>%
-  left_join(cycleTripsByEndOA,by = c("geo_code" = "t.endOA")) %>%
-  left_join(cycleTripsByHomeOA,by = c("geo_code" = "hh.OA")) %>%
-  mutate(across(starts_with("cycle"),function(x) replace_na(x,0)))
+  left_join(getTripsAndShares("t.startOA","start","OA11CD"),by = c("geo_code" = "OA11CD")) %>%
+  left_join(getTripsAndShares("t.endOA","end","OA11CD"),by = c("geo_code" = "OA11CD")) %>%
+  left_join(getTripsAndShares("hh.OA","hh","OA11CD"),by = c("geo_code" = "OA11CD")) %>%
+  mutate(across(starts_with("cycle") | starts_with("walk") | starts_with("trip"),function(x) replace_na(x,0)))
+  
+MSOAs <- MSOAs %>%
+  left_join(getTripsAndShares("t.startOA","start","MSOA11CD")) %>%
+  left_join(getTripsAndShares("t.endOA","end","MSOA11CD")) %>%
+  left_join(getTripsAndShares("hh.OA","hh","MSOA11CD")) %>%
+  mutate(across(starts_with("cycle") | starts_with("walk") | starts_with("trip"),function(x) replace_na(x,0)))
+
+LSOAs <- LSOAs %>%
+  left_join(getTripsAndShares("t.startOA","start","LSOA11CD")) %>%
+  left_join(getTripsAndShares("t.endOA","end","LSOA11CD")) %>%
+  left_join(getTripsAndShares("hh.OA","hh","LSOA11CD")) %>%
+  mutate(across(starts_with("cycle") | starts_with("walk") | starts_with("trip"),function(x) replace_na(x,0)))
+  
 
-sf::write_sf(OAs,"result/bikeShares.shp")
+sf::st_write(OAs,"result/activeSharesOA.gpkg",append = FALSE)
+sf::st_write(MSOAs,"result/activeSharesMSOA.gpkg",append = FALSE)
+sf::st_write(LSOAs,"result/activeSharesLSOA.gpkg",append = FALSE)