R/get_geography_data.R

In konradzdeb/SmarterScotland: Source Public Data from statistics.gov.scot

#' @title Get Geography Data
#'
#' @description The function fetches statistical data for a selected geography
#'   framework from: \href{http://statistics.gov.scot}{statistics.gov.scot}
#'   for a selected period for a chosen indicator.
#'
#' @param data_set \strong{Required.} A character scalar with the name of
#'   indicator to be sourced from the
#'   \href{http://statistics.gov.scot}{statistics.gov.scot}.
#'
#' @param geography \strong{Required.} A character vector corresponding to
#'   geographies for which to derive the data or URI corresponding to a
#'   specific geography.
#'
#' @param measure \strong{Required.} A character scalar pertaining to
#'   measure to be sources, usually \emph{count} or \emph{ratio}. Defaults to
#'   \emph{count}.
#'
#' @param period \strong{Optional.} A character vector corresponding to period(s)
#'  for which the data should be sourced. This has to correspond to the
#'  available periods. For instance \code{c("2001Q2", "2010Q2")} would get the
#'  data for 2nd quarters of 2001 and 2010 respectively.
#'
#' @param pre_process_results A logical, defaults to \code{FALSE}. If
#'   \code{TRUE} will apply \code{\link[SmarterScotland]{pre_process_data}}
#'   function to the results, removing from data frame columns with
#'   information on types \emph{(uri, literal, ...)} and leaving only value
#'   and dimension columns. The function can also remove initial part of URI for
#'   value and dimension columns.
#'
#' @param ... Options passed to \code{\link[SmarterScotland]{pre_process_data}},
#'   ignored if \code{pre_process_results} is \code{FALSE}.
#'
#' @details The following function attempts to source all data available for
#'   specific geographies for the provided period periods. Colloquial
#'   geography names are matched with URI values in
#'   \emph{Standard Geography Code Register}. A well known locality name like
#'   \emph{Edinburgh} may match multiple URIs. Have a look at
#'   \code{\link[SmarterScotland]{available_data_sets}} documentations to
#'   explore data available in the package and
#'   \code{\link[SmarterScotland]{get_available_datasets}} function to get most
#'   recent list.
#'
#'
#' @return A data frame.
#'
#' @export
#'
#' @examples
#' \dontrun{
#' get_geography_data(data_set = "alcohol-related-discharge",
#'                    geography = "East Lothian")
#' }
get_geography_data <-
  function(data_set,
           geography,
           period,
           measure = "count",
           pre_process_results = FALSE,
           ...) {
    # Check if all arguments were specified
    assert_string(x = data_set,
                  na.ok = FALSE,
                  null.ok = FALSE)
    assert_character(x = geography,
                     null.ok = FALSE)

    # Import basic SPARQL query
    query_orig <-
      read_query_file(query_file("qry_get_geography_data"))

    # Check if period was provided and input in the query
    if (missing(period)) {
      filter_time <- ""
    } else {
      # Expand query with the provided period data to filter per period.
      filter_time <- construct_filter(sparql_variable = "time",
                                      filter_values = period)
    }

    # Create a filter for geography
    if (missing(geography)) {
      geographies_URIs <- ""
    } else {
      # Get list of URIs corresponding to the matched geographies
      geographies_URIs <- find_geography_URI(geography = geography,
                                             database = "internal")
      filter_geographies <-
        construct_filter(
          sparql_variable = "reference_area",
          filter_values = geographies_URIs,
          use_str = TRUE
        )
    }

    # Check data set properties and construct relevant SPARQL calls for each
    dta_properties <- get_data_properties(data_set = data_set)

    # Get data set URI
    data_set_URI <-
      paste0("<http://statistics.gov.scot/data/", data_set, ">")

    # Construct variable names
    var_nms <-
      make_SPARQL_variable_names(x = dta_properties$label.value)

    # Selected properties
    selected_properties <- sapply(
      X = dta_properties$property.value,
      FUN = function(x) {
        if (grepl(pattern = "measure-properties", x = x)) {
          if (grepl(pattern = measure, x = x)) {
            x
          } else {
            NULL
          }
        } else {
          x
        }
      }
    )

    indx_measure_remove <-
      which(sapply(selected_properties, is.null))
    var_nms <- var_nms[-indx_measure_remove]

    selected_properties <-
      Filter(Negate(is.null), selected_properties)

    # Generate where clauses
    where_vars_call <- paste(
      mapply(
        FUN = function(x, y) {
          x <- create_URI(x)
          paste("?x", x, y, ".")
        },
        selected_properties,
        var_nms,
        USE.NAMES = FALSE,
        SIMPLIFY = TRUE
      ),
      collapse = " "
    )

    # Replace ?data_set with the call to URI
    where_vars_call <-
      sub(
        pattern = "?data_set",
        replacement = data_set_URI,
        x = where_vars_call,
        fixed = TRUE
      )

    # Collapse variable names for SELECT statement
    var_nms <- paste(var_nms, collapse = " ")

    # Replace query content
    query <- glue(query_orig, .open = "[", .close = "]")

    # Query Scotstat
    response <- query_scotstat(query)

    # Prepare response
    results <- parse_response(response)

    # Check whether to pre-process
    if (pre_process_results) {
      results <- pre_process_data(x = results, ...)
    }

    # Return query results
    return(results)
  }