.Rhistory

max_PTJA_values
# 4. Identify the Area with the Highest PTJA Value
# Identify the area with the highest PTJA value
highest_PTJA_area <- max_PTJA_values %>% filter(max_PTJA == max(max_PTJA))
highest_PTJA_area
LCR_max_PTJA_D <- LCR_dataset_interactive %>% as.data.frame() %>%
summarise(max_PTJA_D = max(PT_Job_Access_Index_Demand, na.rm = TRUE))
GMCA_max_PTJA_D <- GMCA_dataset_interactive %>% as.data.frame() %>%
summarise(max_PTJA_D = max(PT_Job_Access_Index_Demand, na.rm = TRUE))
WYCA_max_PTJA_D <- WYCA_dataset_interactive %>% as.data.frame() %>%
summarise(max_PTJA_D = max(PT_Job_Access_Index_Demand, na.rm = TRUE))
# Combine the results into a single data frame
max_PTJA_D_values <- bind_rows(
LCR_max_PTJA_D %>% mutate(area = "LCR"),
GMCA_max_PTJA_D %>% mutate(area = "GMCA"),
WYCA_max_PTJA_D %>% mutate(area = "WYCA"))
# Identify the area with the highest PTJA-D value
highest_PTJA_D_area <- max_PTJA_D_values %>% filter(max_PTJA_D == max(max_PTJA_D))
highest_PTJA_D_area
max_PTJA_D_values
# Interactive Map
(map <- LCR_dataset_interactive %>%
leaflet() %>%
addProviderTiles("OpenStreetMap") %>%
#  setView(lat = 53.407438, lng = -2.981652, zoom = 10) %>%
addPolygons(data = LCR_boundary,
color = "red",
weight = 1,
fill = FALSE) %>%
addPolygons(data = LCR_dataset_interactive,
color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
fillOpacity = .8,
weight = 1,
group = "PTJA-D",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index Demand: </strong>", PT_Job_Access_Index_Demand ),
htmltools::HTML)) %>%
addPolygons(data = GMCA_dataset_interactive,
color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
fillOpacity = .8,
weight = 1,
group = "PTJA-D",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index Demand: </strong>", PT_Job_Access_Index_Demand ),
htmltools::HTML)) %>%
addPolygons(data = WYCA_dataset_interactive,
color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
fillOpacity = .8,
weight = 1,
group = "PTJA-D",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index Demand: </strong>", PT_Job_Access_Index_Demand ),
htmltools::HTML)) %>%
addPolygons(data = LCR_dataset_interactive,
color = ~PTJA_pal(PT_Job_Access_Index),
fillOpacity = .8,
weight = 1,
group = "PTJA",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index: </strong>", PT_Job_Access_Index ),
htmltools::HTML)) %>%
addPolygons(data = GMCA_dataset_interactive,
color = ~PTJA_pal(PT_Job_Access_Index),
fillOpacity = .8,
weight = 1,
group = "PTJA",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index: </strong>", PT_Job_Access_Index ),
htmltools::HTML)) %>%
addPolygons(data = WYCA_dataset_interactive,
color = ~PTJA_pal(PT_Job_Access_Index),
fillOpacity = .8,
weight = 1,
group = "PTJA",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index: </strong>", PT_Job_Access_Index ),
htmltools::HTML)) %>%
addLegend("topright",
pal = PTJA_D_pal,
values = ~PT_Job_Access_Index_Demand,
title = "PTJA-D",
opacity = 1,
group = "PTJA-D") %>%
addLegend("topright",
pal = PTJA_pal,
values = ~PT_Job_Access_Index,
title = "PTJA",
opacity = 1,
group = "PTJA") %>%
addLayersControl("Layers",
position = "topleft",
baseGroups = c("PTJA", "PTJA-D"),
options = layersControlOptions(collapsed = FALSE)) %>%
addScaleBar((position = "bottomleft"))
)
PTJA_pal
# Combine interactive datasets
interactive_datasets <- rbind(LCR_dataset_interactive, GMCA_dataset_interactive, WYCA_dataset_interactive)
LCR_dataset_interactive
GMCA_dataset_interactive
WYCA_dataset_interactive
## Read GMCA datasets
# Load Manchester Geo data from file
GMCA_dataset_interactive <- read_sf("../Data/MANCH_dataset.shp") %>%
rename("LSOA_Code" = "LSOA21C",
"LSOA_Name" = "LSOA21N",
"TravelTime_Jobcentre" = "Tr__J_C",
"PT_Job_Access_Index" = "PT_Jb_A_I",
"PT_Job_Access_Index_bus" = "PT_J_A_I_B",
"PT_Job_Access_Index_tram" = "PT_J_A_I_T",
"PT_Job_Access_Index_walk" = "PT_J_A_I_W",
"Jobs" = "Employd",
"LSOA_Area" = "LSOA__2",
"Traveltime_empcent" = "Tr__E_C",
"PT_Job_Access_Index_Demand" = "PT_Jb_A_I_") %>%
st_transform(4326) %>%
mutate(PT_Job_Access_Index_bus = PT_Job_Access_Index_bus - PT_Job_Access_Index_walk,
PT_Job_Access_Index_tram = PT_Job_Access_Index_tram - PT_Job_Access_Index_walk) %>%
dplyr::select(LSOA_Code,
LSOA_Name,
Jobs,
LSOA_Area,
Traveltime_empcent,
PT_Job_Access_Index,
PT_Job_Access_Index_bus,
PT_Job_Access_Index_train,
PT_Job_Access_Index_walk,
PT_Job_Access_Index_Demand) %>%
as("Spatial")
## Read GMCA datasets
# Load Manchester Geo data from file
GMCA_dataset_interactive <- read_sf("../Data/MANCH_dataset.shp") %>%
rename("LSOA_Code" = "LSOA21C",
"LSOA_Name" = "LSOA21N",
"TravelTime_Jobcentre" = "Tr__J_C",
"PT_Job_Access_Index" = "PT_Jb_A_I",
"PT_Job_Access_Index_bus" = "PT_J_A_I_B",
"PT_Job_Access_Index_tram" = "PT_J_A_I_T",
"PT_Job_Access_Index_walk" = "PT_J_A_I_W",
"Jobs" = "Employd",
"LSOA_Area" = "LSOA__2",
"Traveltime_empcent" = "Tr__E_C",
"PT_Job_Access_Index_Demand" = "PT_Jb_A_I_") %>%
st_transform(4326) %>%
mutate(PT_Job_Access_Index_bus = PT_Job_Access_Index_bus - PT_Job_Access_Index_walk,
PT_Job_Access_Index_tram = PT_Job_Access_Index_tram - PT_Job_Access_Index_walk) %>%
dplyr::select(LSOA_Code,
LSOA_Name,
Jobs,
LSOA_Area,
Traveltime_empcent,
PT_Job_Access_Index,
PT_Job_Access_Index_bus,
PT_Job_Access_Index_walk,
PT_Job_Access_Index_Demand) %>%
as("Spatial")
# 2. Load Data ------------------------------------------------------------
# Load Liverpool Geo data from file
LCR_dataset_interactive <- read_sf("../Data/LCR_dataset.shp") %>%
rename("LSOA_Code" = "LSOA21C",
"LSOA_Name" = "LSOA21N",
"PT_Job_Access_Index" = "PT_Jb_A_I",
"PT_Job_Access_Index_bus" = "PT_J_A_I_B",
"PT_Job_Access_Index_train" = "PT_J_A_I_T",
"PT_Job_Access_Index_walk" = "PT_J_A_I_W",
"Jobs" = "Employd",
"LSOA_Area" = "LSOA__2",
"Traveltime_empcent" = "Tr__E_C",
"PT_Job_Access_Index_Demand" = "PT_Jb_A_I_") %>%
st_transform(4326) %>%
mutate(PT_Job_Access_Index_bus = PT_Job_Access_Index_bus - PT_Job_Access_Index_walk,
PT_Job_Access_Index_train = PT_Job_Access_Index_train - PT_Job_Access_Index_walk,
LSOA_Code = as.character(LSOA_Code),
LSOA_Area = round(LSOA_Area,2),
PT_Job_Access_Index = round(PT_Job_Access_Index,2),
) %>%
dplyr::select(LSOA_Code,
LSOA_Name,
Jobs,
LSOA_Area,
Traveltime_empcent,
PT_Job_Access_Index,
PT_Job_Access_Index_bus,
PT_Job_Access_Index_walk,
PT_Job_Access_Index_Demand) %>%
as("Spatial") # translates the shapefile into a spatial polygons dataframe
## Read WYCA datasets
WYCA_dataset_interactive <- read_sf("../Data/WYCA_dataset.shp") %>%
rename("LSOA_Code" = "LSOA21C",
"LSOA_Name" = "LSOA21N",
"PT_Job_Access_Index" = "PT_Jb_A_I",
"PT_Job_Access_Index_bus" = "PT_J_A_I_B",
"PT_Job_Access_Index_train" = "PT_J_A_I_T",
"PT_Job_Access_Index_walk" = "PT_J_A_I_W",
"Jobs" = "Employd",
"LSOA_Area" = "LSOA__2",
"Traveltime_empcent" = "Tr__E_C",
"PT_Job_Access_Index_Demand" = "PT_Jb_A_I_") %>%
st_transform(4326) %>%
mutate(PT_Job_Access_Index_bus = PT_Job_Access_Index_bus - PT_Job_Access_Index_walk,
PT_Job_Access_Index_train = PT_Job_Access_Index_train - PT_Job_Access_Index_walk,
LSOA_Code = as.character(LSOA_Code),
LSOA_Area = round(LSOA_Area,2),
PT_Job_Access_Index = round(PT_Job_Access_Index,2),
) %>%
dplyr::select(LSOA_Code,
LSOA_Name,
Jobs,
LSOA_Area,
Traveltime_empcent,
PT_Job_Access_Index,
PT_Job_Access_Index_bus,
PT_Job_Access_Index_walk,
PT_Job_Access_Index_Demand) %>%
as("Spatial") # translates the shapefile into a spatial polygons dataframe
# Combine interactive datasets
interactive_datasets <- rbind(LCR_dataset_interactive, GMCA_dataset_interactive, WYCA_dataset_interactive)
interactive_datasets
# Colour palettes based on the PTJA of with the highest max value
PTJA_pal <- colorNumeric(palette = "Spectral",
domain = interactive_datasets@data$PT_Job_Access_Index,
n=5,
reverse = TRUE)
PTJA_D_pal <- colorNumeric(palette = "Spectral",
domain = interactive_datasets@data$PT_Job_Access_Index_Demand,
n=5,
reverse = TRUE)
# Interactive Map
(map <- interactive_datasets %>%
leaflet() %>%
addProviderTiles("OpenStreetMap") %>%
#  setView(lat = 53.407438, lng = -2.981652, zoom = 10) %>%
addPolygons(data = interactive_datasets,
color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
fillOpacity = .8,
weight = 1,
group = "PTJA-D",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index Demand: </strong>", PT_Job_Access_Index_Demand ),
htmltools::HTML)) %>%
#    addPolygons(data = GMCA_dataset_interactive,
#                color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
#                fillOpacity = .8,
#               weight = 1,
#                group = "PTJA-D",
#                label = ~lapply(paste0(
#                  "<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
#                  "<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
#                ),  htmltools::HTML),
#                highlight = highlightOptions(weight =3, color ="red",
#                                             bringToFront =TRUE),
#                popup = ~lapply(paste0(
#                              "<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
#                              "<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                              #                "<strong>Population: </strong>", Population, "<br>",
#                              "<strong>Jobs: </strong>", Jobs, "<br>",
#                              "<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
3                              "<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
# Interactive Map
(map <- interactive_datasets %>%
leaflet() %>%
addProviderTiles("OpenStreetMap") %>%
#  setView(lat = 53.407438, lng = -2.981652, zoom = 10) %>%
addPolygons(data = interactive_datasets,
color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
fillOpacity = .8,
weight = 1,
group = "PTJA-D",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index Demand: </strong>", PT_Job_Access_Index_Demand ),
htmltools::HTML)) %>%
#    addPolygons(data = GMCA_dataset_interactive,
#                color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
#                fillOpacity = .8,
#               weight = 1,
#                group = "PTJA-D",
#                label = ~lapply(paste0(
#                  "<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
#                  "<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
#                ),  htmltools::HTML),
#                highlight = highlightOptions(weight =3, color ="red",
#                                             bringToFront =TRUE),
#                popup = ~lapply(paste0(
#                              "<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
#                              "<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                              #                "<strong>Population: </strong>", Population, "<br>",
#                              "<strong>Jobs: </strong>", Jobs, "<br>",
#                              "<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
#                              "<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
#                              "<strong>PT Job Access Index Demand: </strong>", PT_Job_Access_Index_Demand ),
#                  htmltools::HTML)) %>%
addPolygons(data = interactive_datasets,
color = ~PTJA_pal(PT_Job_Access_Index),
fillOpacity = .8,
weight = 1,
group = "PTJA",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index: </strong>", PT_Job_Access_Index ),
htmltools::HTML)) %>%
addLegend("topright",
pal = PTJA_D_pal,
values = ~PT_Job_Access_Index_Demand,
title = "PTJA-D",
opacity = 1,
group = "PTJA-D") %>%
addLegend("topright",
pal = PTJA_pal,
values = ~PT_Job_Access_Index,
title = "PTJA",
opacity = 1,
group = "PTJA") %>%
addLayersControl("Layers",
position = "topleft",
baseGroups = c("PTJA", "PTJA-D"),
options = layersControlOptions(collapsed = FALSE)) %>%
addScaleBar((position = "bottomleft"))
)
# Save Map
saveWidget(map, file="docs/index.html", selfcontained = TRUE)
# Interactive Map
(map <- interactive_datasets %>%
leaflet() %>%
addProviderTiles("OpenStreetMap") %>%
#  setView(lat = 53.407438, lng = -2.981652, zoom = 10) %>%
addPolygons(data = interactive_datasets,
color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
fillOpacity = .8,
weight = 1,
group = "PTJA-D",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index Demand: </strong>", PT_Job_Access_Index_Demand ),
htmltools::HTML)) %>%
#    addPolygons(data = GMCA_dataset_interactive,
#                color = ~PTJA_D_pal(PT_Job_Access_Index_Demand),
#                fillOpacity = .8,
#               weight = 1,
#                group = "PTJA-D",
#                label = ~lapply(paste0(
#                  "<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
#                  "<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
#                ),  htmltools::HTML),
#                highlight = highlightOptions(weight =3, color ="red",
#                                             bringToFront =TRUE),
#                popup = ~lapply(paste0(
#                              "<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
#                              "<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                              #                "<strong>Population: </strong>", Population, "<br>",
#                              "<strong>Jobs: </strong>", Jobs, "<br>",
#                              "<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
#                              "<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
#                              "<strong>PT Job Access Index Demand: </strong>", PT_Job_Access_Index_Demand ),
#                  htmltools::HTML)) %>%
addPolygons(data = interactive_datasets,
color = ~PTJA_pal(PT_Job_Access_Index),
fillOpacity = .8,
weight = 1,
group = "PTJA",
label = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>"
),  htmltools::HTML),
highlight = highlightOptions(weight =3, color ="red",
bringToFront =TRUE),
popup = ~lapply(paste0(
"<strong>LSOA Code: </strong>", as.character(LSOA_Code), "<br>",
"<strong>LSOA Name: </strong>", LSOA_Name, "<br>",
#                "<strong>Population: </strong>", Population, "<br>",
"<strong>Jobs: </strong>", Jobs, "<br>",
"<strong>LSOA Area: </strong>", LSOA_Area, "<br>",
"<strong>Time to nearest town centre: </strong>", Traveltime_empcent, "<br>",
"<strong>PT Job Access Index: </strong>", PT_Job_Access_Index ),
htmltools::HTML)) %>%
addLegend("topright",
pal = PTJA_D_pal,
values = ~PT_Job_Access_Index_Demand,
title = "PTJA-D",
opacity = 1,
group = "PTJA-D") %>%
addLegend("topright",
pal = PTJA_pal,
values = ~PT_Job_Access_Index,
title = "PTJA",
opacity = 1,
group = "PTJA") %>%
addLayersControl(
position = "topleft",
baseGroups = c("PTJA", "PTJA-D"),
options = layersControlOptions(collapsed = FALSE)) %>%
addScaleBar((position = "bottomleft"))
)
# Save Map
saveWidget(map, file="docs/index.html", selfcontained = TRUE)