The daiquiri package generates data quality reports that enable quick visual review of temporal shifts in record-level data. Time series plots showing aggregated values are automatically created for each data field (column) depending on its contents (e.g. min/max/mean values for numeric data, no. of distinct values for categorical data), as well as overviews for missing values, non-conformant values, and duplicated rows.

Essentially, it takes input such as this:

And outputs this:

The resulting html reports are shareable and can contribute to forming a transparent record of the entire analysis process. It is designed with electronic health records in mind, but can be used for any type of record-level temporal data.

## Why should I use it?

Large routinely-collected datasets are increasingly being used in research. However, given their data are collected for operational rather than research purposes, there is a greater-than-usual need for them to be checked for data quality issues before any analyses are conducted. Events occurring at the institutional level such as software updates, new machinery or processes can cause temporal artefacts that, if not identified and taken into account, can lead to biased results and incorrect conclusions. For example, the figures below show real data from a large hospital in the UK, and how it has changed over time.

The first figure shows the percentage of missing values in the ‘Duration’ field of a dataset containing antibiotic prescriptions, and the second figure shows the mean value of all laboratory tests checking for levels of ‘creatinine’ in the blood. As you can see, there are points in time where these values shift up or down suddenly and unnaturally, indicating that something changed in the way the data was collected or processed. A careful researcher needs to take these sudden changes into account, particularly if comparing or combining the data before and after these ‘change points’.

While these checks should theoretically be conducted by the researcher at the initial data analysis stage, in practice it is unclear to what extent this is actually done, since it is rarely, if ever, reported in published papers. With the increasing drive towards greater transparency and reproducibility within the scientific community, this essential yet often-overlooked part of the analysis process will inevitably begin to come under greater scrutiny. The daiquiri package helps researchers conduct this part of the process more thoroughly, consistently and transparently, hence increasing the quality of their studies as well as trust in the scientific process.

## Installation

# install the latest release from CRAN
install.packages("daiquiri")

# or install the development version from rOpenSci
install.packages("daiquiri", repos = "https://ropensci.r-universe.dev")

# or install direct from source
# install.packages("remotes")
remotes::install_github("ropensci/daiquiri")

## Usage

library(daiquiri)

# load delimited file into a data.frame without doing any datatype conversion
path <- system.file("extdata", "example_prescriptions.csv", package = "daiquiri")
raw_data <- read_data(path, show_progress = FALSE)

head(raw_data)
## # A tibble: 6 × 8
##   PrescriptionID PrescriptionDate    Admis…¹ Drug  Dose  DoseU…² Patie…³ Locat…⁴
##   <chr>          <chr>               <chr>   <chr> <chr> <chr>   <chr>   <chr>
## 1 6000           2021-01-01 00:00:00 2020-1… Ceft… 500   mg      4993679 SITE1
## 2 6001           NULL                2020-1… Fluc… 1000  mg      819452  SITE1
## 3 6002           NULL                2020-1… Teic… 400   mg      275597  SITE1
## 4 6003           2021-01-01 01:00:00 2020-1… Fluc… 1000  NULL    819452  SITE1
## 5 6004           2021-01-01 02:00:00 2020-1… Fluc… 1000  NULL    528071  SITE1
## 6 6005           2021-01-01 03:00:00 2020-1… Co-a… 1.2   g       1001434 SITE1
## # … with abbreviated variable names ¹​AdmissionDate, ²​DoseUnit, ³​PatientID,
## #   ⁴​Location
# specify the type of data expected in each column of the data.frame
fts <- field_types(
PrescriptionID = ft_uniqueidentifier(),
PrescriptionDate = ft_timepoint(),
Drug = ft_freetext(),
Dose = ft_numeric(),
DoseUnit = ft_categorical(),
PatientID = ft_ignore(),
Location = ft_categorical(aggregate_by_each_category = TRUE)
)
# create a report in the current directory
daiq_obj <- daiquiri_report(
raw_data,
field_types = fts
)

An example report is available from the package website.

More detailed guidance can be found in the walkthrough vignette:

vignette("daiquiri", package = "daiquiri")

## How to cite this package

Please remember to update the version number to match the version you used.

Quan TP (2022). daiquiri: Data quality reporting for temporal datasets. R package version v1.0.0. Zenodo. https://doi.org/10.5281/zenodo.6334341. URL: https://github.com/ropensci/daiquiri

## Acknowledgements

This work was supported by the National Institute for Health Research Health Protection Research Unit (NIHR HPRU) in Healthcare Associated Infections and Antimicrobial Resistance at the University of Oxford in partnership with Public Health England (PHE) (NIHR200915), and by the NIHR Oxford Biomedical Research Centre.

## Contributing to this package

Please report any bugs or suggestions by opening a github issue.

Please note that this package is released with a Contributor Code of Conduct. By contributing to this project, you agree to abide by its terms.