REDCapCAST/R/REDCap_split.r

204 lines
6.6 KiB
R
Raw Permalink Normal View History

2018-01-19 17:50:08 +01:00
#' Split REDCap repeating instruments table into multiple tables
#'
2018-06-03 22:08:26 +02:00
#' This will take output from a REDCap export and split it into a base table
#' and child tables for each repeating instrument. Metadata
2018-01-19 17:50:08 +01:00
#' is used to determine which fields should be included in each resultant table.
#'
#' @param records Exported project records. May be a \code{data.frame},
#' \code{response}, or \code{character} vector containing JSON from an API
2018-06-03 22:08:26 +02:00
#' call.
#' @param metadata Project metadata (the data dictionary). May be a
#' \code{data.frame}, \code{response}, or \code{character} vector containing
#' JSON from an API call.
#' @param primary_table_name Name given to the list element for the primary
#' output table (as described in \emph{README.md}). Ignored if
#' \code{forms = 'all'}.
#' @param forms Indicate whether to create separate tables for repeating
#' instruments only or for all forms.
2018-01-19 17:50:08 +01:00
#' @author Paul W. Egeler, M.S., GStat
#' @examples
#' \dontrun{
#' # Using an API call -------------------------------------------------------
#'
2018-01-19 17:50:08 +01:00
#' library(RCurl)
#'
2018-01-19 17:50:08 +01:00
#' # Get the records
#' records <- postForm(
#' uri = api_url, # Supply your site-specific URI
#' token = api_token, # Supply your own API token
#' content = 'record',
#' format = 'json',
#' returnFormat = 'json'
2018-01-19 17:50:08 +01:00
#' )
#'
#' # Get the metadata
#' metadata <- postForm(
#' uri = api_url, # Supply your site-specific URI
#' token = api_token, # Supply your own API token
#' content = 'metadata',
#' format = 'json'
#' )
#'
2018-06-03 22:08:26 +02:00
#' # Convert exported JSON strings into a list of data.frames
#' REDCapRITS::REDCap_split(records, metadata)
#'
#' # Using a raw data export -------------------------------------------------
#'
#' # Get the records
#' records <- read.csv("/path/to/data/ExampleProject_DATA_2018-06-03_1700.csv")
#'
#' # Get the metadata
2023-04-13 10:57:04 +02:00
#' metadata <- read.csv(
#' "/path/to/data/ExampleProject_DataDictionary_2018-06-03.csv")
#'
#' # Split the tables
#' REDCapRITS::REDCap_split(records, metadata)
#'
#' # In conjunction with the R export script ---------------------------------
#'
2023-04-13 10:57:04 +02:00
#' # You must set the working directory first since the REDCap data export
#' # script contains relative file references.
#' old <- getwd()
#' setwd("/path/to/data/")
#'
#' # Run the data export script supplied by REDCap.
#' # This will create a data.frame of your records called 'data'
#' source("ExampleProject_R_2018-06-03_1700.r")
#'
#' # Get the metadatan
#' metadata <- read.csv("ExampleProject_DataDictionary_2018-06-03.csv")
#'
#' # Split the tables
#' REDCapRITS::REDCap_split(data, metadata)
#' setwd(old)
2018-01-19 17:50:08 +01:00
#' }
#' @return A list of \code{"data.frame"}s. The number of tables will differ
#' depending on the \code{forms} option selected.
#' \itemize{
#' \item \code{'repeating'}: one base table and one or more
#' tables for each repeating instrument.
#' \item \code{'all'}: a data.frame for each instrument, regardless of
#' whether it is a repeating instrument or not.
#' }
#' @include process_user_input.r utils.r
2018-01-19 17:50:08 +01:00
#' @export
REDCap_split <- function(records,
metadata,
primary_table_name = "",
2023-01-16 09:49:17 +01:00
forms = c("repeating", "all")) {
2023-03-06 14:36:32 +01:00
# Process user input
records <- process_user_input(records)
2023-01-16 09:49:17 +01:00
metadata <-
as.data.frame(process_user_input(metadata))
# Get the variable names in the dataset
vars_in_data <- names(records)
# Process repeat instrument names to match the redcap naming
2024-02-07 10:09:36 +01:00
if (is_repeated_longitudinal(records)){
records$redcap_repeat_instrument <- clean_redcap_name(records$redcap_repeat_instrument)
# Match arg for forms
forms <- match.arg(forms, c("repeating", "all"))
# Check to see if there were any repeating instruments
2023-01-16 09:49:17 +01:00
if (forms == "repeating" &&
!"redcap_repeat_instrument" %in% vars_in_data) {
stop("There are no repeating instruments in this dataset.")
}
2018-01-19 17:50:08 +01:00
2019-07-26 23:26:33 +02:00
# Remove NAs from `redcap_repeat_instrument` (see issue #12)
2023-01-16 09:49:17 +01:00
if (any(is.na(records$redcap_repeat_instrument))) {
2019-07-26 23:26:33 +02:00
records$redcap_repeat_instrument <- ifelse(
is.na(records$redcap_repeat_instrument),
"",
as.character(records$redcap_repeat_instrument)
)
}
}
2019-07-26 23:26:33 +02:00
# Standardize variable names for metadata
2023-03-06 14:36:32 +01:00
# names(metadata) <- metadata_names
# Make sure that no metadata columns are factors
2023-01-16 09:49:17 +01:00
metadata <-
rapply(metadata, as.character, classes = "factor", how = "replace")
# Find the fields and associated form
fields <- match_fields_to_form(metadata, vars_in_data)
2018-01-19 17:50:08 +01:00
# Variables to be present in each output table
universal_fields <- c(
vars_in_data[1],
grep(
"^redcap_(?!(repeat)).*",
vars_in_data,
value = TRUE,
perl = TRUE
)
)
if ("redcap_repeat_instrument" %in% vars_in_data) {
# Variables to be at the beginning of each repeating instrument
2023-01-16 09:49:17 +01:00
repeat_instrument_fields <- grep("^redcap_repeat.*",
vars_in_data,
value = TRUE)
# Identify the subtables in the data
subtables <- unique(records$redcap_repeat_instrument)
subtables <- subtables[subtables != ""]
# Split the table based on instrument
2023-01-16 09:49:17 +01:00
out <-
split.data.frame(records, records$redcap_repeat_instrument)
primary_table_index <- which(names(out) == "")
2018-01-19 17:50:08 +01:00
if (forms == "repeating" && primary_table_name %in% subtables) {
2023-01-16 09:49:17 +01:00
warning(
2023-04-13 10:57:04 +02:00
"The label given to the primary table is already used by a repeating
instrument. The primary table label will be left blank."
2023-01-16 09:49:17 +01:00
)
primary_table_name <- ""
} else if (primary_table_name > "") {
names(out)[[primary_table_index]] <- primary_table_name
}
# Delete the variables that are not relevant
for (i in names(out)) {
if (i == primary_table_name) {
2023-01-16 09:49:17 +01:00
out_fields <- which(vars_in_data %in% c(universal_fields,
2023-04-13 10:57:04 +02:00
fields[!fields[, 2] %in%
subtables, 1]))
2023-01-16 09:49:17 +01:00
out[[primary_table_index]] <-
out[[primary_table_index]][out_fields]
2018-01-19 17:50:08 +01:00
} else {
2023-01-16 09:49:17 +01:00
out_fields <- which(vars_in_data %in% c(universal_fields,
repeat_instrument_fields,
fields[fields[, 2] == i, 1]))
out[[i]] <- out[[i]][out_fields]
2018-01-19 17:50:08 +01:00
}
2018-01-19 17:50:08 +01:00
}
if (forms == "all") {
2023-01-16 09:49:17 +01:00
out <- c(split_non_repeating_forms(out[[primary_table_index]],
universal_fields,
fields[!fields[, 2] %in% subtables, ]),
out[-primary_table_index])
}
} else {
out <- split_non_repeating_forms(records, universal_fields, fields)
}
2018-01-19 17:50:08 +01:00
2018-06-03 22:08:26 +02:00
out
2018-01-19 17:50:08 +01:00
}