# Install the development version from GitHub:
# install.packages("pak")
pak::pak("mikejohnpage/finddxgap")

library(finddxgap)

# Compute a diagnostic gap
cases <- data.frame(
  estimated = c(100, 80, 100),
  notified = c(10, 60, 70)
)

dxgap(cases, estimated, notified)
#>   estimated notified dxgap
#> 1       100       10    90
#> 2        80       60    25
#> 3       100       70    30