R/gisco-get-healthcare.R

In giscoR: Download Map Data from GISCO API - Eurostat

#' Healthcare services in Europe
#'
#' @description
#' The dataset contains information on main healthcare services considered to
#' be 'hospitals' by Member States. The definition varies slightly from country
#' to country, but roughly includes the following:
#'
#' *"'Hospitals' comprises licensed establishments primarily engaged in
#' providing medical, diagnostic, and treatment services that include physician,
#' nursing, and other health services to in-patients and the specialised
#' accommodation services required by inpatients*.
#'
#' @family services
#' @inherit gisco_get_education return source
#' @inheritParams gisco_get_countries
#' @encoding UTF-8
#' @export
#'
#' @param year character string or number. Release year of the file. One of
#'   `2023`, `2020`.
#'
#'
#' @details
#' Files are distributed [EPSG:4326](https://epsg.io/4326).
#'
#' ```{r child = "man/chunks/healthcare_meta.Rmd"}
#' ```
#'
#' @examplesIf gisco_check_access()
#' health_benelux <- gisco_get_healthcare(
#'   country = c("BE", "NL", "LU"),
#'   year = 2023
#' )
#'
#' # Plot if downloaded
#' if (!is.null(health_benelux)) {
#'   benelux <- gisco_get_countries(country = c("BE", "NL", "LU"))
#'
#'   library(ggplot2)
#'   ggplot(benelux) +
#'     geom_sf(fill = "grey10", color = "grey20") +
#'     geom_sf(
#'       data = health_benelux, color = "red",
#'       size = 0.2, alpha = 0.25
#'     ) +
#'     theme_void() +
#'     theme(
#'       plot.background = element_rect(fill = "black"),
#'       text = element_text(color = "white"),
#'       panel.grid = element_blank(),
#'       plot.title = element_text(face = "bold", hjust = 0.5),
#'       plot.subtitle = element_text(face = "italic", hjust = 0.5)
#'     ) +
#'     labs(
#'       title = "Healthcare services", subtitle = "Benelux (2023)",
#'       caption = "Source: Eurostat, Healthcare 2023 dataset."
#'     ) +
#'     coord_sf(crs = 3035)
#' }
gisco_get_healthcare <- function(
  year = c(2023, 2020),
  cache = TRUE,
  update_cache = FALSE,
  cache_dir = NULL,
  verbose = FALSE,
  country = NULL
) {
  # Given vars
  year <- match_arg_pretty(year)

  api_entry <- paste0(
    "https://gisco-services.ec.europa.eu/pub/healthcare/",
    year,
    "/gpkg/EU.gpkg"
  )
  filename <- paste0("health_", year, "_", basename(api_entry))

  if (cache) {
    # Guess source to load
    namefileload <- download_url(
      api_entry,
      filename,
      cache_dir,
      "health",
      update_cache,
      verbose
    )
  } else {
    namefileload <- api_entry
  }

  if (is.null(namefileload)) {
    return(NULL)
  }

  data_sf <- read_geo_file_sf(namefileload)

  if (!is.null(country) && "cntr_id" %in% names(data_sf)) {
    country <- convert_country_code(country)
    data_sf <- data_sf[data_sf$cntr_id %in% country, ]
  }

  data_sf
}