Assorted repo updates

.gitignore

@@ -28,3 +28,4 @@ reminder_on_the_conditions_of_use.txt
 infuse_dist_lyr_2011.tgz
 cbs neighbourhoods.zip
 load_files
+minidocs/rsconnect/

analysis/rnet.R

@@ -37,6 +37,8 @@ line_widths <- line_widths + 0.3
 summary(line_widths)
 line_widths[line_widths > 8] <- 8
 
+library(leaflet)
+rg <
 leaflet() %>% addTiles() %>% addPolylines(data = rg, weight = line_widths, popup = rg@data$gov_target)
 
 plot(rg, lwd = rg@data$gov_target / 500)

buildmaster.R

@@ -1,6 +1,7 @@
 library(knitr)
 
-la_all <- c("Barnet", "Barnsley", "Bath and North East Somerset", "Bedford")
+la_all <- c("Barking-and-Dagenham", "Barnet",
+            "Barnsley", "Bath-and-North-East-Somerset")
 
 for(i in la_all){
   la <- i
@@ -17,9 +18,9 @@ for(i in la_all){
 }
 
 
-# las <- readOGR(dsn = "pct-bigdata/national/cuas-mf.geojson", layer = "OGRGeoJSON")
-# las_names <- las$CTYUA12NM
-# las_names <- las_names[order(las_names)]
-# las_names <- as.character(las_names)
-# head(las_names)
-# dput(las_names[2:5])
+las <- readOGR(dsn = "pct-bigdata/national/cuas-mf.geojson", layer = "OGRGeoJSON")
+las_names <- las$CTYUA12NM
+las_names <- las_names[order(las_names)]
+las_names <- as.character(las_names)
+head(las_names)
+dput(las_names[1:4])

cover-letter.Rmd

@@ -0,0 +1,31 @@
+---
+date: "August 23, 2015"
+output: word_document
+---
+
+<div style="float: right">
+
+Dr Robin Lovelace
+
+Leeds Institute for Data Analytics
+
+University of Leeds
+
+Leeds, LS2 9JT
+
+0113 343 0691
+</div>
+
+Dear Professor Ortúzar,
+
+Please consider the manuscript, entitled 'The Propensity to Cycle Tool: An open source online system for sustainable transport planning', for publication in Transportation Research Part A: Policy and Practice. The manuscript consists of approximately 6,000 words excluding references, eight figures and two tables.
+
+We think this paper is ideally suited to Transportation Part A because of high policy relevance and its framing from a policy perspective. Transportation Part A has a long history of methodological advance in transport policy evaluation which this paper seeks to extend in the field of strategic planning for new cycling infrastructure.
+
+The policy relevance of the research presented stems from its source of funding. The research was funded by the UK's Department for Transport (DfT), which intends to use the Propensity to Cycle Tool (PCT) to inform its investment strategies. This policy context, of increased interest and funding for cycling, is not unique to the UK. Indeed, we believe the PCT method will be of interest and potential use for transport planners and researchers in many countries. We have emphasised the potential re-use of the methods and code underlying the PCT internationally. The international reach and high impact of Transportation Part A makes it an ideal platform to communicate this research and make it available to researchers worldwide.
+
+We look forward to feedback from the reviewers and hope for a fast turn-around due to the policy-relevant nature of this research.
+
+Yours sincerely
+
+Robin

documents/national.Rmd

@@ -16,7 +16,7 @@ bibliography: Transport.bib
 
 ```{r, include=FALSE}
 pkgs <- c("grid", "png", "knitr", "tmap", "dplyr", "ggmap", "pander")
-lapply(pkgs, library, character.only = T)
+pload <- lapply(pkgs, require, character.only = T)
 opts_knit$set(root.dir = "../")
 # source("analysis/la-analysis.R", print.eval = T)
 ```

loading-data/load-flowsex.R

@@ -63,3 +63,24 @@ abline(a = 0, b = 1)
 
 plot(flow$Bicycle, flow$gendereq)
 abline(a = 0, b = 1)
+
+# joining flow with flowsex:
+flow <- readRDS("pct-bigdata/national/flow_eng_avlslope.Rds")
+flowsex <- read_csv("private-data/wu01bew_msoa_v1.csv", col_names=F, skip=12)
+names(flowsex) <- c("Area.of.residence", "Area.of.workplace", "All.fs", "Male", "Female")
+flowsex$id <- paste(flowsex$Area.of.residence, flowsex$Area.of.workplace)
+# Subset by zones in the study area
+o <- flow$Area.of.residence %in% cents$geo_code
+d <- flow$Area.of.workplace %in% cents$geo_code
+flow <- flow[o & d, ] # subset OD pairs with o and d in study area
+sel <- flow$All > mflow # subset OD pairs by n. people using it
+flow <- readRDS("pct-bigdata/national/flow_eng_avlslope.Rds")
+names(flowsex) <- c("Area.of.residence", "Area.of.workplace", "All.fs", "Male", "Female")
+flowsex$id <- paste(flowsex$Area.of.residence, flowsex$Area.of.workplace)
+# Subset by zones in the study area
+o <- flow$Area.of.residence %in% cents$geo_code
+d <- flow$Area.of.workplace %in% cents$geo_code
+flow$id <- paste(flow$Area.of.residence, flow$Area.of.workplace)
+flowsex <- dplyr::select(flowsex, Male, Female, id)
+flow <- left_join(flow, flowsex, by = "id")
+saveRDS(flow, "private-data/flowsex-merged.Rds")

set-up.R

@@ -2,7 +2,6 @@
 # NB: devtools allows installation of the latest packages
 if(!require(devtools)) install.packages("devtools")
 if(!require(stplanr)) install_github("robinlovelace/stplanr")
-if(!require(geojsonio)) install_github("ropensci/geojsonio")
 pkgs <- c(
   "ggmap",
   "tmap",
@@ -16,14 +15,15 @@ pkgs <- c(
   "geojsonio", # loads geojsons
   "tidyr", # tidies up your data!
   "readr", # reads your data fast
-  "knitr" # for knitting it all together
+  "knitr", # for knitting it all together
+  "geojsonio"
   )
 # Which packages do we require?
 # lapply(pkgs, library, character.only = T)
 reqs <- as.numeric(lapply(pkgs, require, character.only = TRUE))
 # Install packages we require
 if(sum(!reqs) > 0) install.packages(pkgs[!reqs])
-# Load publicly available test data
+ # Load publicly available test data
 
 # Option 1: clone the repository directly - if you have git installed
 # system2("git", args=c("clone", "[email protected]:Robinlovelace/pct-data.git", "--depth=1"))