Despite the negative impacts of COVID-19, I am sure there have been positive developments in healthcare. For example, the pandemic has increased focus on public health and hygiene, which could prevent the spread of other diseases.
And we might not realise, despite the overloaded health system and staffing criss, it has also led to more investment and innovation in healthcare infrastructure and resources. It results in more covid test sites, health facilities, and medical professionals, which would provide opportunities to improvement and better supports in a long run.
Let’s see before and after COVID, the increased number of geospatial health units in NZ. The number here compares the number of units in 2019 and 2022
Please wait a few secs for the map to load ….. 🙂
Code
library(leaflet)
library(leaflegend) # Required for addLegendBin
# Calculate centroids for visualization
# ta - territorial authority
# class(ta)[1] "sf" "tbl_df" "tbl" "data.frame"
location <- st_centroid(ta)
# Filtering out certain values from 'ta' dataset
ta <- ta %>% filter(ta2018_v1_00 != 999)
# Define color bins for legend
color_bin <- colorBin("PRGn", domain = location$gain, bins = 4, pretty = FALSE, na.color = "#808080")
# Create leaflet map
leaflet() %>%
setView(lng = 174.768, lat = -36.852, zoom = 6) %>%
addProviderTiles("Stadia.StamenTonerBackground") %>% # this might not work if you aren't registered with Stadia
addCircleMarkers(data = location,
radius = ~sqrt(abs(gain)) * 3,
label = ~paste0(ta$ta2018_v1_00_name, ":", ta$gain),
weight = 10,
stroke = FALSE,
opacity = 0,
popup = ~paste0(ta$ta2018_v1_00_name, ": ", ta$gain),
fillOpacity = 0.8,
color = ~color_bin(gain)) %>%
addLegendBin(pal = color_bin,
values = location$gain,
position = "bottomright",
orientation = 'vertical',
title = "Increases in units in 2 years",
labelStyle = 'font-size: 10px; font-weight: bold;',
shape = 'circle'
)It seems there are more health units post COVID time. The rise is majorly significant in the urban areas.
