Verify cycle count data at regional level by comparing with NTS data

[joeytalbot]

No description provided.

[Robinlovelace]

Links with #76

[Robinlovelace]

Next step I guess is to convert the numbers here of km cycled per person per year into bkm cycled/yr using mid year population estimates, e.g. from here: https://www.ons.gov.uk/peoplepopulationandcommunity/populationandmigration/populationestimates/datasets/middlesuperoutputareamidyearpopulationestimates

Suggested place to put that code: https://github.com/saferactive/saferactive/blob/master/code/tests/regional-uptake-nts.R that could just continue that reproducible script.

[Robinlovelace]

Just a thought: it may be worth uploading regions_simple and reading that in instead of the detailed raw geojson object here:

saferactive/code/tests/regional-uptake-nts.R

Lines 18 to 21 in d709f4a

	regions = read_sf("https://github.com/saferactive/saferactive/releases/download/0.1.4/Regions_.December_2020._EN_BGC.geojson")
	regions = sf::st_make_valid(regions)
	qtm(regions)
	regions_simple = rmapshaper::ms_simplify(regions, 0.01)

Was thinking of doing that but will hold my horses for now. Let me know Joey how it goes and happy to provide further input. That region object could be used to get regional mid year population estimates from the MSOA data, suggest making centroids of those MSOAs and then aggregating with aggregate() or st_join() and then group_by() + summarise() pipeline.

[Robinlovelace]

Link to GAM results at regional level: https://github.com/saferactive/saferactive/releases/download/0.1.4/gam-full-results-grid-national.Rds