eucast_rules(): add add_if_missing argument to control NA imputation (#259)

When `add_if_missing = FALSE`, rules are only applied to cells that already
contain an SIR value; `NA` cells are left untouched. This is useful with
`overwrite = TRUE` to update reported results without imputing values for
drugs that were not tested.

https://claude.ai/code/session_01Nucc8nXGLqNUjtuC9GrhTc

CLAUDE.md

@@ -148,24 +148,34 @@ Version format: `major.minor.patch.dev` (e.g., `3.0.1.9021`)
 
 ### Version and date bump required for every PR
 
-Before opening a pull request, always increment the four-digit dev counter by 1 in **both** of these files:
+All PRs are **squash-merged**, so each PR lands as exactly **one commit** on the default branch. Version numbers are kept in sync with the cumulative commit count since the last released tag. Therefore **exactly one version bump is allowed per PR**, regardless of how many intermediate commits are made on the branch.
 
-1. **`DESCRIPTION`** — the `Version:` field:
-   ```
-   Version: 3.0.1.9021  →  Version: 3.0.1.9022
-   ```
+#### Computing the correct version number
 
-2. **`NEWS.md`** — the top-level heading:
-   ```
-   # AMR 3.0.1.9021  →  # AMR 3.0.1.9022
-   ```
+Run the following from the repo root to determine the version string to use:
 
-Read the current version from `DESCRIPTION`, add 1 to the last numeric component, and write the new version to both files in the same commit as the rest of the PR changes.
+```bash
+currenttag=$(git describe --tags --abbrev=0 | sed 's/v//')
+currenttagfull=$(git describe --tags --abbrev=0)
+defaultbranch=$(git branch | cut -c 3- | grep -E '^master$|^main$')
+currentcommit=$(git rev-list --count ${currenttagfull}..${defaultbranch})
+currentversion="${currenttag}.$((currentcommit + 9001 + 1))"
+echo "$currentversion"
+```
 
-Also bump the date to the current date in **`DESCRIPTION`**, where it's in the `Date:` field in ISO format:
+The `+ 1` accounts for the fact that this PR's squash commit is not yet on the default branch. Set **both** of these files to the resulting version string (and only once per PR, even across multiple commits):
+
+1. **`DESCRIPTION`** — the `Version:` field
+2. **`NEWS.md`** — the top-level heading `# AMR <version>`
+
+If `git describe` fails (e.g. no tags exist in the environment), fall back to reading the current version from `DESCRIPTION` and adding 1 to the last numeric component — but only if no bump has already been made in this PR.
+
+#### Date field
+
+The `Date:` field in `DESCRIPTION` must reflect the date of the **last commit to the PR** (not the first), in ISO format. Update it with every commit so it is always current:
 
 ```
-Date: 2025-12-31
+Date: 2026-03-07
 ```
 
 ## Internal State

DESCRIPTION

@@ -1,6 +1,6 @@
 Package: AMR
-Version: 3.0.1.9029
-Date: 2026-03-06
+Version: 3.0.1.9034
+Date: 2026-03-09
 Title: Antimicrobial Resistance Data Analysis
 Description: Functions to simplify and standardise antimicrobial resistance (AMR)
   data analysis and to work with microbial and antimicrobial properties by

NAMESPACE

@@ -248,6 +248,7 @@ export(glycopeptides)
 export(guess_ab_col)
 export(inner_join_microorganisms)
 export(interpretive_rules)
+export(ionophores)
 export(is.ab)
 export(is.av)
 export(is.disk)

NEWS.md

@@ -1,4 +1,4 @@
-# AMR 3.0.1.9029
+# AMR 3.0.1.9034
 
 ### New
 * Integration with the **tidymodels** framework to allow seamless use of SIR, MIC and disk data in modelling pipelines via `recipes`
@@ -8,17 +8,19 @@
     - `all_mic()`, `all_mic_predictors()`
     - `all_disk()`, `all_disk_predictors()`
 * Data set `esbl_isolates` to practise with AMR modelling
-* AMR selectors `peptides()`, `phosphonics()` and `spiropyrimidinetriones()`
+* AMR selectors `ionophores()`, `peptides()`, `phosphonics()` and `spiropyrimidinetriones()`
 * Support for Wildtype (WT) / Non-wildtype (NWT) in `as.sir()`, all plotting functions, and all susceptibility/resistance functions.
   - `as.sir()` gained an argument `as_wt_nwt`, which defaults to `TRUE` only when `breakpoint_type = "ECOFF"` (#254)
   - This transforms the output from S/R to WT/NWT
   - Functions such as `susceptibility()` count WT as S and NWT as R
 * `interpretive_rules()`, which allows future implementation of CLSI interpretive rules (#235)
   - `eucast_rules()` has become a wrapper around that function
+* `eucast_rules()` / `interpretive_rules()` gained argument `add_if_missing` (default: `TRUE`). When set to `FALSE`, rules are only applied to cells that already contain an SIR value; `NA` cells are left untouched. This is useful with `overwrite = TRUE` to update reported results without imputing values for drugs that were not tested (#259)
 * Two new `NA` objects, `NA_ab_` and `NA_mo_`, analogous to base R's `NA_character_` and `NA_integer_`, for use in pipelines that require typed missing values
 
 ### Fixes
-* `mdro()`: when a base beta-lactam drug column is missing but a corresponding drug+inhibitor combination is present in the data and resistant (e.g., piperacillin/tazobactam = R while piperacillin is absent), the base drug is now correctly inferred as resistant. This ensures MDRO classification is not missed due to test-ordering differences in the laboratory. The reverse direction is also valid: susceptibility in a combination does not imply susceptibility in the base drug (the inhibitor may be responsible), so only resistance is propagated. Closes #209
+* Fixed a bug in `as.sir()` where values that were purely numeric (e.g., `"1"`) and matched the broad SIR-matching regex would be incorrectly stripped of all content by the Unicode letter filter
+* Fixed a bug in `as.mic()` where MIC values in scientific notation (e.g., `"1e-3"`) were incorrectly handled because the letter `e` was removed along with other Unicode letters; scientific notation `e` is now preserved
 * Fixed a bug in `as.ab()` where certain AB codes containing "PH" or "TH" (such as `ETH`, `MTH`, `PHE`, `PHN`, `STH`, `THA`, `THI1`) would incorrectly return `NA` when combined in a vector with any untranslatable value (#245)
 * Fixed a bug in `antibiogram()` for when no antimicrobials are set
 * Fixed a bug in `as.sir()` where for numeric input the arguments `S`,  `I`,  and `R` would not be considered (#244)
@@ -29,11 +31,12 @@
 * Fixed SIR and MIC coercion of combined values, e.g. `as.sir("<= 0.002; S") ` or `as.mic("S; 0.002")` (#252)
 
 ### Updates
+* `mdro()` now infers resistance for a _missing_ base drug column from an _available_ corresponding drug+inhibitor combination showing resistance (e.g., piperacillin is absent but required, while piperacillin/tazobactam available and resistant). Can be set with the new argument `infer_from_combinations`, which defaults to `TRUE` (#209). Note that this can yield a higher MDRO detection (which is a good thing as it has become more reliable).
 * `susceptibility()` and `resistance()` gained the argument `guideline`, which defaults to EUCAST, for interpreting the 'I' category correctly.
+* Added to the `antimicrobials` data set: cefepime/taniborbactam (`FTA`), ceftibuten/avibactam (`CTA`), clorobiocin (`CLB`), kasugamycin (`KAS`), ostreogrycin (`OST`), taniborbactam (`TAN`), thiostrepton (`THS`), xeruborbactam (`XER`), and zorbamycin (`ZOR`)
 * `as.mic()` and `rescale_mic()` gained the argument `round_to_next_log2`, which can be set to `TRUE` to round all values up to the nearest next log2 level (#255)
 * `antimicrobials$group` is now a `list` instead of a `character`, to contain any group the drug is in (#246)
 * `ab_group()` gained an argument `all_groups` to return all groups the antimicrobial drug is in (#246)
-* Added to the `antimicrobials` data set: cefepime/taniborbactam (`FTA`), ceftibuten/avibactam (`CTA`), kasugamycin (`KAS`), ostreogrycin (`OST`), taniborbactam (`TAN`), thiostrepton (`THS`), xeruborbactam (`XER`), and zorbamycin (`ZOR`)
 * Added explaining message to `as.sir()` when interpreting numeric values (e.g., 1 for S, 2 for I, 3 for R) (#244)
 * Updated handling of capped MIC values (`<`, `<=`, `>`, `>=`) in `as.sir()` in the argument `capped_mic_handling`: (#243)
   * Introduced four clearly defined options: `"none"`, `"conservative"` (default), `"standard"`, and `"lenient"`

R/aa_helper_functions.R

@@ -387,6 +387,10 @@ import_fn <- function(name, pkg, error_on_fail = TRUE) {
   if (isTRUE(error_on_fail)) {
     stop_ifnot_installed(pkg)
   }
+  if (pkg == "rstudioapi" && tryCatch(!rstudioapi::isAvailable(), error = function(e) TRUE)) {
+    # only allow rstudioapi to be imported if RStudio is available
+    return(NULL)
+  }
   tryCatch(
     # don't use get() to avoid fetching non-API functions
     getExportedValue(name = name, ns = asNamespace(pkg)),
@@ -1220,10 +1224,14 @@ try_colour <- function(..., before, after, collapse = " ") {
   }
 }
 is_dark <- function() {
-  AMR_env$current_theme <- tryCatch(getExportedValue("getThemeInfo", ns = asNamespace("rstudioapi"))()$editor, error = function(e) NULL)
+  AMR_env$current_theme <- NULL
+  current_theme_fn <- import_fn("getThemeInfo", "rstudioapi", error_on_fail = FALSE)
+  if (!is.null(current_theme_fn)) {
+    AMR_env$current_theme <- current_theme_fn()$editor
+  }
   if (!identical(AMR_env$current_theme, AMR_env$former_theme) || is.null(AMR_env$is_dark_theme)) {
     AMR_env$former_theme <- AMR_env$current_theme
-    AMR_env$is_dark_theme <- !has_colour() || tryCatch(isTRUE(getExportedValue("getThemeInfo", ns = asNamespace("rstudioapi"))()$dark), error = function(e) TRUE)
+    AMR_env$is_dark_theme <- !has_colour() || tryCatch(isTRUE(current_theme_fn()$dark), error = function(e) TRUE)
   }
   isTRUE(AMR_env$is_dark_theme)
 }

R/amr_selectors.R

@@ -352,6 +352,14 @@ glycopeptides <- function(only_sir_columns = FALSE, return_all = TRUE, ...) {
   amr_select_exec("glycopeptides", only_sir_columns = only_sir_columns, return_all = return_all)
 }
 
+#' @rdname antimicrobial_selectors
+#' @export
+ionophores <- function(only_sir_columns = FALSE, return_all = TRUE, ...) {
+  meet_criteria(only_sir_columns, allow_class = "logical", has_length = 1)
+  meet_criteria(return_all, allow_class = "logical", has_length = 1)
+  amr_select_exec("ionophores", only_sir_columns = only_sir_columns, return_all = return_all)
+}
+
 #' @rdname antimicrobial_selectors
 #' @export
 isoxazolylpenicillins <- function(only_sir_columns = FALSE, only_treatable = TRUE, return_all = TRUE, ...) {

R/interpretive_rules.R

@@ -79,6 +79,7 @@ format_eucast_version_nr <- function(version, markdown = TRUE) {
 #' @param only_sir_columns A [logical] to indicate whether only antimicrobial columns must be included that were transformed to class [sir][as.sir()] on beforehand. Defaults to `FALSE` if no columns of `x` have a class [sir][as.sir()].
 #' @param custom_rules Custom rules to apply, created with [custom_eucast_rules()].
 #' @param overwrite A [logical] indicating whether to overwrite existing SIR values (default: `FALSE`). When `FALSE`, only non-SIR values are modified (i.e., any value that is not already S, I or R). To ensure compliance with EUCAST guidelines, **this should remain** `FALSE`, as EUCAST notes often state that an organism "should be tested for susceptibility to individual agents or be reported resistant".
+#' @param add_if_missing A [logical] indicating whether rules should also be applied to missing (`NA`) values (default: `TRUE`). When `FALSE`, rules are only applied to cells that already contain an SIR value; cells with `NA` are left untouched. This is particularly useful when using `overwrite = TRUE` with custom rules and you want to update reported results without imputing values for untested drugs.
 #' @inheritParams first_isolate
 #' @details
 #' **Note:** This function does not translate MIC values to SIR values. Use [as.sir()] for that. \cr
@@ -175,6 +176,7 @@ interpretive_rules <- function(x,
                                only_sir_columns = any(is.sir(x)),
                                custom_rules = NULL,
                                overwrite = FALSE,
+                               add_if_missing = TRUE,
                                ...) {
   meet_criteria(x, allow_class = "data.frame")
   meet_criteria(col_mo, allow_class = "character", has_length = 1, is_in = colnames(x), allow_NULL = TRUE)
@@ -189,6 +191,7 @@ interpretive_rules <- function(x,
   meet_criteria(only_sir_columns, allow_class = "logical", has_length = 1)
   meet_criteria(custom_rules, allow_class = "custom_eucast_rules", allow_NULL = TRUE)
   meet_criteria(overwrite, allow_class = "logical", has_length = 1)
+  meet_criteria(add_if_missing, allow_class = "logical", has_length = 1)
 
   stop_if(
     guideline == "CLSI",
@@ -538,7 +541,8 @@ interpretive_rules <- function(x,
           warned = warned,
           info = info,
           verbose = verbose,
-          overwrite = overwrite
+          overwrite = overwrite,
+          add_if_missing = add_if_missing
         )
         n_added <- n_added + run_changes$added
         n_changed <- n_changed + run_changes$changed
@@ -580,7 +584,8 @@ interpretive_rules <- function(x,
           warned = warned,
           info = info,
           verbose = verbose,
-          overwrite = overwrite
+          overwrite = overwrite,
+          add_if_missing = add_if_missing
         )
         n_added <- n_added + run_changes$added
         n_changed <- n_changed + run_changes$changed
@@ -877,7 +882,8 @@ interpretive_rules <- function(x,
       warned = warned,
       info = info,
       verbose = verbose,
-      overwrite = overwrite
+      overwrite = overwrite,
+      add_if_missing = add_if_missing
     )
     n_added <- n_added + run_changes$added
     n_changed <- n_changed + run_changes$changed
@@ -947,7 +953,8 @@ interpretive_rules <- function(x,
         warned = warned,
         info = info,
         verbose = verbose,
-        overwrite = overwrite
+        overwrite = overwrite,
+        add_if_missing = add_if_missing
       )
       n_added <- n_added + run_changes$added
       n_changed <- n_changed + run_changes$changed
@@ -1139,7 +1146,8 @@ edit_sir <- function(x,
                      warned,
                      info,
                      verbose,
-                     overwrite) {
+                     overwrite,
+                     add_if_missing) {
   cols <- unique(cols[!is.na(cols) & !is.null(cols)])
 
   # for Verbose Mode, keep track of all changes and return them
@@ -1172,13 +1180,15 @@ edit_sir <- function(x,
     if (isFALSE(overwrite) && any(isSIR) && message_not_thrown_before("edit_sir.warning_overwrite")) {
       warning_("Some values had SIR values and were not overwritten, since `overwrite = FALSE`.")
     }
+    # determine which cells to modify based on overwrite and add_if_missing
+    apply_mask <- if (isTRUE(overwrite)) {
+      if (isFALSE(add_if_missing)) !isNA else rep(TRUE, length(isNA))
+    } else {
+      if (isFALSE(add_if_missing)) isSIR else non_SIR
+    }
     tryCatch(
       # insert into original table
-      if (isTRUE(overwrite)) {
-        new_edits[rows, cols] <- to
-      } else {
-        new_edits[rows, cols][non_SIR] <- to
-      },
+      new_edits[rows, cols][apply_mask] <- to,
       warning = function(w) {
         if (w$message %like% "invalid factor level") {
           xyz <- vapply(FUN.VALUE = logical(1), cols, function(col) {
@@ -1188,11 +1198,7 @@ edit_sir <- function(x,
             )
             TRUE
           })
-          if (isTRUE(overwrite)) {
-            suppressWarnings(new_edits[rows, cols] <<- to)
-          } else {
-            suppressWarnings(new_edits[rows, cols][non_SIR] <<- to)
-          }
+          suppressWarnings(new_edits[rows, cols][apply_mask] <<- to)
           warning_(
             "in `eucast_rules()`: value \"", to, "\" added to the factor levels of column",
             ifelse(length(cols) == 1, "", "s"),

R/mdro.R

@@ -31,7 +31,7 @@
 #'
 #' Determine which isolates are multidrug-resistant organisms (MDRO) according to international, national, or custom guidelines.
 #' @param x A [data.frame] with antimicrobials columns, like `AMX` or `amox`. Can be left blank for automatic determination.
-#' @param guideline A specific guideline to follow, see sections *Supported international / national guidelines* and *Using Custom Guidelines* below. When left empty, the publication by Magiorakos *et al.* (see below) will be followed.
+#' @param guideline A specific guideline to follow, see sections *Supported International / National Guidelines* and *Using Custom Guidelines* below. When left empty, the publication by Magiorakos *et al.* (see below) will be followed.
 #' @param esbl [logical] values, or a column name containing logical values, indicating the presence of an ESBL gene (or production of its proteins).
 #' @param carbapenemase [logical] values, or a column name containing logical values, indicating the presence of a carbapenemase gene (or production of its proteins).
 #' @param mecA [logical] values, or a column name containing logical values, indicating the presence of a *mecA* gene (or production of its proteins).
@@ -42,6 +42,7 @@
 #' @param pct_required_classes Minimal required percentage of antimicrobial classes that must be available per isolate, rounded down. For example, with the default guideline, 17 antimicrobial classes must be available for *S. aureus*. Setting this `pct_required_classes` argument to `0.5` (default) means that for every *S. aureus* isolate at least 8 different classes must be available. Any lower number of available classes will return `NA` for that isolate.
 #' @param combine_SI A [logical] to indicate whether all values of S and I must be merged into one, so resistance is only considered when isolates are R, not I. As this is the default behaviour of the [mdro()] function, it follows the redefinition by EUCAST about the interpretation of I (increased exposure) in 2019, see section 'Interpretation of S, I and R' below. When using `combine_SI = FALSE`, resistance is considered when isolates are R or I.
 #' @param verbose A [logical] to turn Verbose mode on and off (default is off). In Verbose mode, the function returns a data set with the MDRO results in logbook form with extensive info about which isolates would be MDRO-positive, or why they are not.
+#' @param infer_from_combinations A [logical] to indicate whether resistance for a missing base beta-lactam drug should be inferred from an available drug+inhibitor combination (e.g., piperacillin from piperacillin/tazobactam). The clinical basis is that resistance in a combination always implies resistance in the base drug, since the enzyme inhibitor provides no benefit when the organism is truly resistant. Only resistance is inferred; susceptibility in a combination does **not** imply susceptibility in the base drug (the inhibitor may be responsible). Defaults to `TRUE`.
 #' @details
 #' These functions are context-aware. This means that the `x` argument can be left blank if used inside a [data.frame] call, see *Examples*.
 #'
@@ -143,6 +144,7 @@ mdro <- function(x = NULL,
                  combine_SI = TRUE,
                  verbose = FALSE,
                  only_sir_columns = any(is.sir(x)),
+                 infer_from_combinations = TRUE,
                  ...) {
   if (is_null_or_grouped_tbl(x)) {
     # when `x` is left blank, auto determine it (get_current_data() searches underlying data within call)
@@ -165,7 +167,7 @@ mdro <- function(x = NULL,
   meet_criteria(combine_SI, allow_class = "logical", has_length = 1)
   meet_criteria(verbose, allow_class = "logical", has_length = 1)
   meet_criteria(only_sir_columns, allow_class = "logical", has_length = 1)
-
+  meet_criteria(infer_from_combinations, allow_class = "logical", has_length = 1)
 
   if (isTRUE(only_sir_columns) && !any(is.sir(x))) {
     stop_("There were no SIR columns found in the data set, despite `only_sir_columns` being `TRUE`. Transform columns with `as.sir()` for valid antimicrobial interpretations.")
@@ -480,12 +482,8 @@ mdro <- function(x = NULL,
   }
   cols_ab <- cols_ab[!duplicated(cols_ab)]
 
-  # Infer resistance for missing base drugs from available drug+inhibitor combination columns.
-  # Clinical principle: resistance in drug+inhibitor (e.g., piperacillin/tazobactam = R)
-  # always implies resistance in the base drug (e.g., piperacillin = R), because the
-  # enzyme inhibitor adds nothing when the organism is truly resistant to the base drug.
-  # NOTE: susceptibility in a combination does NOT imply susceptibility in the base drug
-  # (the inhibitor may be responsible), so synthetic proxy columns only propagate R, not S/I.
+  # Infer resistance for missing base drugs ----
+  if (isTRUE(infer_from_combinations)) {
     .combos_in_data <- AB_BETALACTAMS_WITH_INHIBITOR[AB_BETALACTAMS_WITH_INHIBITOR %in% names(cols_ab)]
     if (length(.combos_in_data) > 0) {
       .base_drugs <- suppressMessages(
@@ -505,16 +503,22 @@ mdro <- function(x = NULL,
               lapply(x[, .combo_cols, drop = FALSE], function(col) as.character(as.sir(col))),
               stringsAsFactors = FALSE
             )
-          .new_col <- paste0(".sir_proxy_", .base_code)
+            .new_col <- paste0(.base_code, ".inferred_sir_proxy_from#", paste0(.combos, collapse = "/"), "#")
             x[[.new_col]] <- ifelse(rowSums(.sir_chars == "R", na.rm = TRUE) > 0L, "R", NA_character_)
-          cols_ab <- c(cols_ab, setNames(.new_col, .base_code))
-          if (isTRUE(verbose)) {
+            cols_ab <- c(cols_ab, stats::setNames(.new_col, .base_code))
+            if (isTRUE(info.bak)) {
               message_(
-              "Inferring resistance for ", ab_name(.base_code, language = NULL),
-              " from available drug+inhibitor combination(s): ",
-              paste(ab_name(as.character(.combos), language = NULL), collapse = ", "),
-              " (resistance in a combination always implies resistance in the base drug)",
-              add_fn = font_blue
+                "Inferring resistance for ",
+                ab_name(.base_code, language = NULL, tolower = TRUE),
+                " (", font_bold(.base_code, collapse = NULL), ", ", font_italic("missing"), ") from ",
+                vector_or(
+                  quotes = FALSE,
+                  last_sep = " and/or ",
+                  paste0(
+                    ab_name(.combos, language = NULL, tolower = TRUE),
+                    " (", font_bold(.combos, collapse = NULL), ", ", font_italic("available"), ")"
+                  )
+                )
               )
             }
           }
@@ -522,7 +526,7 @@ mdro <- function(x = NULL,
       }
       cols_ab <- cols_ab[!duplicated(names(cols_ab))]
     }
-  rm(list = intersect(ls(), c(".combos_in_data", ".base_drugs", ".unique_bases", ".base", ".base_code", ".combos", ".combo_cols", ".sir_chars", ".new_col")))
+  }
 
   # nolint start
   AMC <- cols_ab["AMC"]
@@ -1937,7 +1941,8 @@ mdro <- function(x = NULL,
     # format data set
     colnames(x)[colnames(x) == col_mo] <- "microorganism"
     x$microorganism <- mo_name(x$microorganism, language = NULL)
-    x$guideline <- paste0(guideline$author, " - ", guideline$name, ", ", guideline$version, ")")
+    x$guideline <- paste0(guideline$author, " - ", guideline$name, ifelse(is.na(guideline$version), "", paste0(" (", guideline$version, ")")))
+    x$all_nonsusceptible_columns <- gsub(".inferred_sir_proxy_from#(.*?)#", " (inferred from \\1)", x$all_nonsusceptible_columns, perl = TRUE)
     x[, c(
       "row_number",
       "microorganism",

R/mic.R

@@ -217,9 +217,9 @@ as.mic <- function(x, na.rm = FALSE, keep_operators = "all", round_to_next_log2
     warning_("Some MICs were combined values, only the first values are kept")
     x[x %like% "[0-9]/.*[0-9]"] <- gsub("/.*", "", x[x %like% "[0-9]/.*[0-9]"])
   }
-  x <- trimws2(gsub("[\\p{L}]", "", x, perl = TRUE)) # \p{L} is the Unicode category for all letters, including those with diacritics
+  x <- trimws2(gsub("[^e\\P{L}]", "", x, perl = TRUE)) # \p{L} is the Unicode category for all letters, including those with diacritics
   # remove other invalid characters
-  x <- gsub("[^0-9.><= -]+", "", x, perl = TRUE)
+  x <- gsub("[^0-9e.><= -]+", "", x, perl = TRUE)
   # transform => to >= and =< to <=
   x <- gsub("=<", "<=", x, fixed = TRUE)
   x <- gsub("=>", ">=", x, fixed = TRUE)

R/sir.R

@@ -568,7 +568,7 @@ as.sir.default <- function(x,
     x[x %like% "dose"] <- "SDD"
     mtch <- grepl(paste0("(", S, "|", I, "|", R, "|", NI, "|", SDD, "|", WT, "|", NWT, "|", NS, "|[A-Z]+)"), x, perl = TRUE)
     x[!mtch] <- ""
-    x[mtch] <- trimws2(gsub("[^\\p{L}]", "", x[mtch], perl = TRUE)) # \p{L} is the Unicode category for all letters, including those with diacritics
+    x[mtch & x %unlike% "^[0-9+]$"] <- trimws2(gsub("[^\\p{L}]", "", x[mtch & x %unlike% "^[0-9+]$"], perl = TRUE)) # \p{L} is the Unicode category for all letters, including those with diacritics
     # apply regexes set by user
     x[x %like% S] <- "S"
     x[x %like% I] <- "I"

data-raw/_pre_commit_checks.R

@@ -369,6 +369,9 @@ pre_commit_lst$AB_AMINOGLYCOSIDES <- antimicrobials %>%
   filter(group %like% "aminoglycoside|paromomycin|spectinomycin") %>%
   pull(ab)
 pre_commit_lst$AB_AMINOPENICILLINS <- as.ab(c("AMP", "AMX", "AMC"))
+pre_commit_lst$AB_AMINOCOUMARINS <- antimicrobials %>%
+  filter(name %like% "novobiocin|clorobiocin") %>%
+  pull(ab)
 pre_commit_lst$AB_ANTIFUNGALS <- antimicrobials %>%
   filter(group %like% "antifungal") %>%
   pull(ab)
@@ -486,6 +489,18 @@ pre_commit_lst$AB_BETALACTAMS_WITH_INHIBITOR <- antimicrobials %>%
 # this will be used for documentation:
 pre_commit_lst$DEFINED_AB_GROUPS <- sort(names(pre_commit_lst)[names(pre_commit_lst) %like% "^AB_" & names(pre_commit_lst) != "AB_LOOKUP"])
 
+# Check that all AB_* groups with >= 4 members have a corresponding function
+for (grp in pre_commit_lst$DEFINED_AB_GROUPS[pre_commit_lst$DEFINED_AB_GROUPS %unlike% "BETALACTAMASE_INHIBITORS|EXCEPT"]) {
+  if (length(pre_commit_lst[[grp]]) >= 4) {
+    fn_name <- tolower(gsub("^AB_", "", grp))
+    if (!fn_name %in% ls(envir = asNamespace("AMR"))) {
+      stop("Group '", grp, "' has ", length(pre_commit_lst[[grp]]),
+           " members (", toString(ab_name(pre_commit_lst[[grp]], tolower = T)), ") but no corresponding function '", fn_name, "()' exists in the AMR namespace.",
+           call. = FALSE)
+    }
+  }
+}
+
 # Update the antimicrobials$group column
 usethis::ui_info("Updating 'group' column in antimicrobials data set from AB_* vectors")
 prettify_group_name <- function(name) {
@@ -557,6 +572,7 @@ pre_commit_lst$ABX_PRIORITY_LIST <- c("Aminopenicillins",
                                       "Beta-lactams",
                                       "Beta-lactamase inhibitors",
                                       "Pleuromutilins",
+                                      "Aminocoumarins",
                                       "Other")
 if (!all(unlist(antimicrobials$group) %in% pre_commit_lst$ABX_PRIORITY_LIST)) {
   stop("Missing group(s) in priority list: ", paste(setdiff(unlist(antimicrobials$group), pre_commit_lst$ABX_PRIORITY_LIST), collapse = ", "))

data-raw/_reproduction_scripts/reproduction_of_antimicrobials.R

@@ -972,6 +972,17 @@ antimicrobials <- antimicrobials |>
       select(1:4),
   )
 
+antimicrobials <- antimicrobials |>
+  mutate(ab = as.character(ab)) |>
+  bind_rows(
+    antimicrobials |>
+      filter(ab == "NOV") |>
+      mutate(ab = "CLB",
+             cid = 54706138,
+             name = "Clorobiocin") |>
+      select(1:4),
+  )
+
 # update ATC codes from WHOCC website -------------------------------------
 
 library(rvest)
@@ -1171,6 +1182,11 @@ for (i in 1:nrow(antimicrobials)) {
     antimicrobials[i, "loinc"][[1]] <- ifelse(length(loinc) == 0, list(NA_character_), list(loinc))
   }
 }
+antimicrobials$group <- unname(antimicrobials$group)
+antimicrobials$atc <- unname(antimicrobials$atc)
+antimicrobials$abbreviations <- unname(antimicrobials$abbreviations)
+antimicrobials$synonyms <- unname(antimicrobials$synonyms)
+antimicrobials$loinc <- unname(antimicrobials$loinc)
 
 
 usethis::use_data(antimicrobials, overwrite = TRUE, version = 2, compress = "xz")

data-raw/ab.md5

@@ -1 +1 @@
-147709d2fb1b31b013c6ffb387f4d3ba
+11aade8a39bfdff02d01fb52b04eacdc

data-raw/datasets/antimicrobials.txt

@@ -157,6 +157,7 @@
 "CLF1"	2799	"Clofoctol"	"Other"	"J01XX03,QJ01XX03"	"Other antibacterials"	"Other antibacterials"	"NA"	"clofoctolo,clofoctolum,gramplus,octofene,phenol"					"NA"
 "CLM"	71807	"Clometocillin"	"Penicillins,Beta-lactams"	"J01CE07,QJ01CE07"	"Beta-lactam antibacterials, penicillins"	"Beta-lactamase sensitive penicillins"	"NA"	"chlomethocillin,clometacillin,clomethacillin,clomethocillin,clometocilina,clometocilline,clometocillinsalt,clometocillinum,penicilline,rixapen"	1	"g"			"NA"
 "CLM1"	54680675	"Clomocycline"	"Tetracyclines"	"J01AA11,QJ01AA11"	"Tetracyclines"	"Tetracyclines"	"NA"	"clomociclina,clomocyclinum,megaclor"	1	"g"			"NA"
+"CLB"	54706138	"Clorobiocin"	"Aminocoumarins"	"NA"			"NA"	"chlorobiocin"					"NA"
 "CTR"	2812	"Clotrimazole"	"Antifungals"	"A01AB18,D01AC01,G01AF02,QA01AB18,QD01AC01,QG01AF02,QJ02AB90"			"clot"	"alevazol,bisphenyl,canesten,canestene,canestine,canifug,chlotrimazole,clomatin,clotrimaderm,clotrimazol,clotrimazolum,coltrimazole,cutistad,diphenylmethane,empecid,esparol,femmesil,footlogix,fortinia,gynix,imidazole,jidesheng,klotrimazole,lakesia,lombazol,lombazole,lombazolum,lotrimax,lotrimin,monobaycuten,mycelax,mycelex,mycofug,mycosporin,mykosporin,nalbix,otomax,pedesil,pedisafe,ringworm,stiemazol,tibatin,trimysten,trivagizole"					"10653-4,10654-2,18909-2,54177-1,55663-9"
 "CLO"	6098	"Cloxacillin"	"Isoxazolylpenicillins,Penicillins,Beta-lactams"	"J01CF02,QJ01CF02,QJ51CF02,QS01AA90"	"Beta-lactam antibacterials, penicillins"	"Beta-lactamase resistant penicillins"	"clox,cloxac"	"ankerbin,austrastaph,biocloxin,brispen,chloroxacillin,ciclex,clocil,clossacillina,cloxacilina,cloxacillinanhydrous,cloxacilline,cloxacillinsalt,cloxacillinum,cloxapen,constaphyl,dariclox,dichlorstapenor,diclocil,dicloxacillinhydrate,diflor,digloxilline,dynapen,ekvacillin,gelstaph,novapen,noxaben,orbenin,pathocil,stampen,staphybiotic,syntarpen,syntarpensalt,tegopen"	2	"g"	2	"g"	"16628-0,18910-0,196-6,197-4,198-2,199-0,25250-2,55664-7"
 "COL"	5311054	"Colistin"	"Polymyxins"	"A07AA10,J01XB01,QA07AA10,QJ01XB01,QJ51XB01"	"Other antibacterials"	"Polymyxins"	"cl,coli,colist,cs,cst,ct"	"colimycin,colisticin,colisticina,colistina,colistine,colistinum,colobreathe,colomycin,kangdisu,kolimitsin,kolimycin,promixin,sogecoli,totazina"	9	"MU"	9	"MU"	"16645-4,18912-6,204-8,205-5,206-3,207-1,29493-4,33333-6"
@@ -307,7 +308,7 @@
 "NME"		"Norfloxacin/metronidazole"	"Fluoroquinolones,Quinolones"	"J01RA14,QJ01RA14"	"Combinations of antibacterials"	"Combinations of antibacterials"	"NA"	"NA"					"NA"
 "NTI"		"Norfloxacin/tinidazole"	"Fluoroquinolones,Quinolones"	"J01RA13,QJ01RA13"	"Combinations of antibacterials"	"Combinations of antibacterials"	"NA"	"NA"					"NA"
 "NVA"	10419027	"Norvancomycin"	"Glycopeptides,Peptides"	"NA"			"NA"	"NA"					"NA"
-"NOV"	54675769	"Novobiocin"	"Other"	"QJ01XX95"			"novo,novobi"	"albadry,albamix,albamycin,biotexin,cardelmycin,cardelmycinsalt,cathocin,cathomycin,inabiocin,novobiocina,novobiocine,novobiocinsalt,novobiocinum,robiocina,sirbiocina,spheromycin,stilbiocina,streptonivicin,streptonivicinsalt,vulcamicina,vulcamycin,vulkamycin"					"17378-1,18957-1,370-7,371-5,372-3,373-1,41706-3"
+"NOV"	54675769	"Novobiocin"	"Aminocoumarins"	"QJ01XX95"			"novo,novobi"	"albadry,albamix,albamycin,biotexin,cardelmycin,cardelmycinsalt,cathocin,cathomycin,inabiocin,novobiocina,novobiocine,novobiocinsalt,novobiocinum,robiocina,sirbiocina,spheromycin,stilbiocina,streptonivicin,streptonivicinsalt,vulcamicina,vulcamycin,vulkamycin"					"17378-1,18957-1,370-7,371-5,372-3,373-1,41706-3"
 "NYS"	6433272	"Nystatin"	"Ionophores,Antifungals"	"A07AA02,D01AA01,G01AA01,QA07AA02,QD01AA01,QG01AA01"			"nyst,nystan"	"biofanal,diastatin,herniocid,moronal,myconystatin,mycostatin,mykostatyna,nilstat,nistatin,nistatina,nyotran,nystan,nystatyna,nystavescent,nystex"	1.5	"MU"			"10697-1,10698-9,18958-9,35824-2,55689-4"
 "OFX"	4583	"Ofloxacin"	"Fluoroquinolones,Quinolones"	"J01MA01,QJ01MA01,QS01AE01,QS02AA16,S01AE01,S02AA16"	"Quinolone antibacterials"	"Fluoroquinolones"	"of,ofl,oflo,ofloxa,ofx"	"exocin,exocine,flobacin,floxil,floxin,monoflocet,oflocet,ofloxacina,ofloxacine,ofloxacino,ofloxacinum,ofloxaxin,oxaldin,tarivid,visiren,zanocin"	0.4	"g"	0.4	"g"	"18959-7,20384-4,23948-3,25264-3,374-9,375-6,376-4,377-2,3877-8,41408-6,41409-4,41410-2,42653-6,7038-3,72168-8"
 "OOR"		"Ofloxacin/ornidazole"	"Fluoroquinolones,Quinolones"	"J01RA09,QJ01RA09"	"Combinations of antibacterials"	"Combinations of antibacterials"	"NA"	"NA"					"NA"
@@ -332,7 +333,7 @@
 "PEF"	51081	"Pefloxacin"	"Fluoroquinolones,Quinolones"	"J01MA03,QJ01MA03"	"Quinolone antibacterials"	"Fluoroquinolones"	"pefl,perflo"	"labocton,pefbid,pefloxacine,pefloxacinium,pefloxacino,pefloxacinum,pefocin,pefran,pelox"	0.8	"g"	0.8	"g"	"18963-9,35828-3,390-5,3906-5,7040-9"
 "PEF-S"		"Pefloxacin screening test"	"Fluoroquinolones,Quinolones"	"NA"			"pef screen"	"NA"					"NA"
 "PNM"	10250769	"Penamecillin"	"Penicillins,Beta-lactams"	"J01CE06,QJ01CE06"	"Beta-lactam antibacterials, penicillins"	"Beta-lactamase sensitive penicillins"	"NA"	"havapen,hydroxymethyl,penamecilina,penamecillina,penamecilline,penamecillinum"	1.05	"g"			"NA"
-"PNO"		"Penicillin/novobiocin"	"Penicillins,Beta-lactams"	"NA"			"pennov"	"NA"					"35872-1,35873-9,35874-7"
+"PNO"		"Penicillin/novobiocin"	"Penicillins,Beta-lactams,Aminocoumarins"	"NA"			"pennov"	"NA"					"35872-1,35873-9,35874-7"
 "PSU"		"Penicillin/sulbactam"	"Penicillins,Beta-lactams,Beta-lactamase inhibitors"	"NA"			"NA"	"NA"					"NA"
 "PNM1"	54686187	"Penimepicycline"	"Tetracyclines"	"J01AA10,QJ01AA10"	"Tetracyclines"	"Tetracyclines"	"NA"	"criseocil,duamine,geotricyn,hydrocycline,penetracyne,penimepiciclina,penimepicyclinum"					"NA"
 "PIM"	65453	"Pentisomicin"	"Aminoglycosides"	"NA"			"NA"	"mutamicin,mutamycin,pentisomicina,pentisomicine,pentisomicinum"					"NA"

man/antimicrobial_selectors.Rd

@@ -17,6 +17,7 @@
 \alias{cephalosporins_5th}
 \alias{fluoroquinolones}
 \alias{glycopeptides}
+\alias{ionophores}
 \alias{isoxazolylpenicillins}
 \alias{lincosamides}
 \alias{lipoglycopeptides}
@@ -81,6 +82,8 @@ fluoroquinolones(only_sir_columns = FALSE, only_treatable = TRUE,
 
 glycopeptides(only_sir_columns = FALSE, return_all = TRUE, ...)
 
+ionophores(only_sir_columns = FALSE, return_all = TRUE, ...)
+
 isoxazolylpenicillins(only_sir_columns = FALSE, only_treatable = TRUE,
   return_all = TRUE, ...)
 
@@ -202,6 +205,7 @@ The \code{\link[=not_intrinsic_resistant]{not_intrinsic_resistant()}} function c
 \item \code{\link[=cephalosporins_5th]{cephalosporins_5th()}} can select: \cr ceftaroline (CPT), ceftaroline/avibactam (CPA), ceftobiprole (BPR), ceftobiprole medocaril (CFM1), and ceftolozane/tazobactam (CZT)
 \item \code{\link[=fluoroquinolones]{fluoroquinolones()}} can select: \cr besifloxacin (BES), ciprofloxacin (CIP), ciprofloxacin/metronidazole (CIM), ciprofloxacin/ornidazole (CIO), ciprofloxacin/tinidazole (CIT), clinafloxacin (CLX), danofloxacin (DAN), delafloxacin (DFX), difloxacin (DIF), enoxacin (ENX), enrofloxacin (ENR), finafloxacin (FIN), fleroxacin (FLE), garenoxacin (GRN), gatifloxacin (GAT), gemifloxacin (GEM), grepafloxacin (GRX), lascufloxacin (LSC), levofloxacin (LVX), levofloxacin/ornidazole (LEO), levonadifloxacin (LND), lomefloxacin (LOM), marbofloxacin (MAR), metioxate (MXT), miloxacin (MIL), moxifloxacin (MFX), nadifloxacin (NAD), nemonoxacin (NEM), nifuroquine (NIF), nitroxoline (NTR), norfloxacin (NOR), norfloxacin screening test (NOR-S), norfloxacin/metronidazole (NME), norfloxacin/tinidazole (NTI), ofloxacin (OFX), ofloxacin/ornidazole (OOR), orbifloxacin (ORB), pazufloxacin (PAZ), pefloxacin (PEF), pefloxacin screening test (PEF-S), pradofloxacin (PRA), premafloxacin (PRX), prulifloxacin (PRU), rufloxacin (RFL), sarafloxacin (SAR), sitafloxacin (SIT), sparfloxacin (SPX), temafloxacin (TMX), tilbroquinol (TBQ), tioxacin (TXC), tosufloxacin (TFX), and trovafloxacin (TVA)
 \item \code{\link[=glycopeptides]{glycopeptides()}} can select: \cr avoparcin (AVO), bleomycin (BLM), dalbavancin (DAL), norvancomycin (NVA), oritavancin (ORI), ramoplanin (RAM), teicoplanin (TEC), teicoplanin-macromethod (TCM), telavancin (TLV), vancomycin (VAN), vancomycin-macromethod (VAM), and zorbamycin (ZOR)
+\item \code{\link[=ionophores]{ionophores()}} can select: \cr lasalocid (LAS), monensin sodium (MON), narasin (NAR), nystatin (NYS), and salinomycin (SAL)
 \item \code{\link[=isoxazolylpenicillins]{isoxazolylpenicillins()}} can select: \cr cloxacillin (CLO), dicloxacillin (DIC), flucloxacillin (FLC), meticillin (MET), oxacillin (OXA), and oxacillin screening test (OXA-S)
 \item \code{\link[=lincosamides]{lincosamides()}} can select: \cr clindamycin (CLI), clindamycin inducible screening test (CLI-S), lincomycin (LIN), and pirlimycin (PRL)
 \item \code{\link[=lipoglycopeptides]{lipoglycopeptides()}} can select: \cr dalbavancin (DAL), oritavancin (ORI), and telavancin (TLV)

man/antimicrobials.Rd

@@ -5,9 +5,9 @@
 \alias{antimicrobials}
 \alias{antibiotics}
 \alias{antivirals}
-\title{Data Sets with 624 Antimicrobial Drugs}
+\title{Data Sets with 625 Antimicrobial Drugs}
 \format{
-\subsection{For the \link{antimicrobials} data set: a \link[tibble:tibble]{tibble} with 504 observations and 14 variables:}{
+\subsection{For the \link{antimicrobials} data set: a \link[tibble:tibble]{tibble} with 505 observations and 14 variables:}{
 \itemize{
 \item \code{ab}\cr antimicrobial ID as used in this package (such as \code{AMC}), using the official EARS-Net (European Antimicrobial Resistance Surveillance Network) codes where available. \emph{\strong{This is a unique identifier.}}
 \item \code{cid}\cr Compound ID as found in PubChem. \emph{\strong{This is a unique identifier.}}
@@ -50,7 +50,7 @@ LOINC:
 }
 }
 
-An object of class \code{deprecated_amr_dataset} (inherits from \code{tbl_df}, \code{tbl}, \code{data.frame}) with 504 rows and 14 columns.
+An object of class \code{deprecated_amr_dataset} (inherits from \code{tbl_df}, \code{tbl}, \code{data.frame}) with 505 rows and 14 columns.
 
 An object of class \code{tbl_df} (inherits from \code{tbl}, \code{data.frame}) with 120 rows and 11 columns.
 }

man/custom_eucast_rules.Rd

@@ -98,8 +98,9 @@ x
 #>      amoxicillin (AMX), ampicillin (AMP), azlocillin (AZL), mezlocillin (MEZ), piperacillin (PIP), piperacillin/tazobactam (TZP)
 }\if{html}{\out{</div>}}
 
-These 42 antimicrobial groups are allowed in the rules (case-insensitive) and can be used in any combination:
+These 43 antimicrobial groups are allowed in the rules (case-insensitive) and can be used in any combination:
 \itemize{
+\item aminocoumarins\cr(clorobiocin, novobiocin, and penicillin/novobiocin)
 \item aminoglycosides\cr(amikacin, amikacin/fosfomycin, apramycin, arbekacin, astromicin, bekanamycin, dibekacin, framycetin, gentamicin, gentamicin-high, habekacin, hygromycin, isepamicin, kanamycin, kanamycin-high, kanamycin/cephalexin, kasugamycin, micronomicin, neomycin, netilmicin, pentisomicin, plazomicin, propikacin, ribostamycin, sisomicin, streptoduocin, streptomycin, streptomycin-high, tobramycin, and tobramycin-high)
 \item aminopenicillins\cr(amoxicillin, amoxicillin/clavulanic acid, and ampicillin)
 \item antifungals\cr(amorolfine, amphotericin B, amphotericin B-high, anidulafungin, butoconazole, caspofungin, ciclopirox, clotrimazole, econazole, fluconazole, flucytosine, fosfluconazole, griseofulvin, hachimycin, ibrexafungerp, isavuconazole, isoconazole, itraconazole, ketoconazole, manogepix, micafungin, miconazole, nystatin, oteseconazole, pimaricin, posaconazole, rezafungin, ribociclib, sulconazole, terbinafine, terconazole, and voriconazole)

man/custom_mdro_guideline.Rd

@@ -94,7 +94,7 @@ You can define antimicrobial groups instead of single antimicrobials for the rul
 )
 }\if{html}{\out{</div>}}
 
-All 42 antimicrobial selectors are supported for use in the rules:
+All 43 antimicrobial selectors are supported for use in the rules:
 \itemize{
 \item \code{\link[=aminoglycosides]{aminoglycosides()}} can select: \cr amikacin, amikacin/fosfomycin, apramycin, arbekacin, astromicin, bekanamycin, dibekacin, framycetin, gentamicin, gentamicin-high, habekacin, hygromycin, isepamicin, kanamycin, kanamycin-high, kanamycin/cephalexin, kasugamycin, micronomicin, neomycin, netilmicin, pentisomicin, plazomicin, propikacin, ribostamycin, sisomicin, streptoduocin, streptomycin, streptomycin-high, tobramycin, and tobramycin-high
 \item \code{\link[=aminopenicillins]{aminopenicillins()}} can select: \cr amoxicillin, amoxicillin/clavulanic acid, and ampicillin
@@ -111,6 +111,7 @@ All 42 antimicrobial selectors are supported for use in the rules:
 \item \code{\link[=cephalosporins_5th]{cephalosporins_5th()}} can select: \cr ceftaroline, ceftaroline/avibactam, ceftobiprole, ceftobiprole medocaril, and ceftolozane/tazobactam
 \item \code{\link[=fluoroquinolones]{fluoroquinolones()}} can select: \cr besifloxacin, ciprofloxacin, ciprofloxacin/metronidazole, ciprofloxacin/ornidazole, ciprofloxacin/tinidazole, clinafloxacin, danofloxacin, delafloxacin, difloxacin, enoxacin, enrofloxacin, finafloxacin, fleroxacin, garenoxacin, gatifloxacin, gemifloxacin, grepafloxacin, lascufloxacin, levofloxacin, levofloxacin/ornidazole, levonadifloxacin, lomefloxacin, marbofloxacin, metioxate, miloxacin, moxifloxacin, nadifloxacin, nemonoxacin, nifuroquine, nitroxoline, norfloxacin, norfloxacin screening test, norfloxacin/metronidazole, norfloxacin/tinidazole, ofloxacin, ofloxacin/ornidazole, orbifloxacin, pazufloxacin, pefloxacin, pefloxacin screening test, pradofloxacin, premafloxacin, prulifloxacin, rufloxacin, sarafloxacin, sitafloxacin, sparfloxacin, temafloxacin, tilbroquinol, tioxacin, tosufloxacin, and trovafloxacin
 \item \code{\link[=glycopeptides]{glycopeptides()}} can select: \cr avoparcin, bleomycin, dalbavancin, norvancomycin, oritavancin, ramoplanin, teicoplanin, teicoplanin-macromethod, telavancin, vancomycin, vancomycin-macromethod, and zorbamycin
+\item \code{\link[=ionophores]{ionophores()}} can select: \cr lasalocid, monensin sodium, narasin, nystatin, and salinomycin
 \item \code{\link[=isoxazolylpenicillins]{isoxazolylpenicillins()}} can select: \cr cloxacillin, dicloxacillin, flucloxacillin, meticillin, oxacillin, and oxacillin screening test
 \item \code{\link[=lincosamides]{lincosamides()}} can select: \cr clindamycin, clindamycin inducible screening test, lincomycin, and pirlimycin
 \item \code{\link[=lipoglycopeptides]{lipoglycopeptides()}} can select: \cr dalbavancin, oritavancin, and telavancin

man/mdro.Rd

@@ -18,7 +18,8 @@
 mdro(x = NULL, guideline = "CMI 2012", col_mo = NULL, esbl = NA,
   carbapenemase = NA, mecA = NA, mecC = NA, vanA = NA, vanB = NA,
   info = interactive(), pct_required_classes = 0.5, combine_SI = TRUE,
-  verbose = FALSE, only_sir_columns = any(is.sir(x)), ...)
+  verbose = FALSE, only_sir_columns = any(is.sir(x)),
+  infer_from_combinations = TRUE, ...)
 
 brmo(x = NULL, only_sir_columns = any(is.sir(x)), ...)
 
@@ -35,7 +36,7 @@ eucast_exceptional_phenotypes(x = NULL, only_sir_columns = any(is.sir(x)),
 \arguments{
 \item{x}{A \link{data.frame} with antimicrobials columns, like \code{AMX} or \code{amox}. Can be left blank for automatic determination.}
 
-\item{guideline}{A specific guideline to follow, see sections \emph{Supported international / national guidelines} and \emph{Using Custom Guidelines} below. When left empty, the publication by Magiorakos \emph{et al.} (see below) will be followed.}
+\item{guideline}{A specific guideline to follow, see sections \emph{Supported International / National Guidelines} and \emph{Using Custom Guidelines} below. When left empty, the publication by Magiorakos \emph{et al.} (see below) will be followed.}
 
 \item{col_mo}{Column name of the names or codes of the microorganisms (see \code{\link[=as.mo]{as.mo()}}) - the default is the first column of class \code{\link{mo}}. Values will be coerced using \code{\link[=as.mo]{as.mo()}}.}
 
@@ -61,6 +62,8 @@ eucast_exceptional_phenotypes(x = NULL, only_sir_columns = any(is.sir(x)),
 
 \item{only_sir_columns}{A \link{logical} to indicate whether only antimicrobial columns must be included that were transformed to class \link[=as.sir]{sir} on beforehand. Defaults to \code{FALSE} if no columns of \code{x} have a class \link[=as.sir]{sir}.}
 
+\item{infer_from_combinations}{A \link{logical} to indicate whether resistance for a missing base beta-lactam drug should be inferred from an available drug+inhibitor combination (e.g., piperacillin from piperacillin/tazobactam). The clinical basis is that resistance in a combination always implies resistance in the base drug, since the enzyme inhibitor provides no benefit when the organism is truly resistant. Only resistance is inferred; susceptibility in a combination does \strong{not} imply susceptibility in the base drug (the inhibitor may be responsible). Defaults to \code{TRUE}.}
+
 \item{...}{Column names of antimicrobials. To automatically detect antimicrobial column names, do not provide any named arguments; \code{\link[=guess_ab_col]{guess_ab_col()}} will then be used for detection. To manually specify a column, provide its name (case-insensitive) as an argument, e.g. \code{AMX = "amoxicillin"}. To skip a specific antimicrobial, set it to \code{NULL}, e.g. \code{TIC = NULL} to exclude ticarcillin. If a manually defined column does not exist in the data, it will be skipped with a warning.}
 }
 \value{

tests/testthat/test-data.R

@@ -35,7 +35,10 @@ test_that("test-data.R", {
   expect_identical(class(microorganisms$mo), c("mo", "character"))
   expect_identical(nrow(antimicrobials), length(unique(AMR::antimicrobials$ab)))
   expect_identical(class(AMR::antimicrobials$ab), c("ab", "character"))
-  expect_identical(nrow(antimicrobials[!is.na(antimicrobials$cid), ]), length(unique(AMR::antimicrobials$cid[!is.na(antimicrobials$cid)])), label = "nr of rows with CIDs", expected.label = "unique nr of CIDs")
+  expect_identical(
+    nrow(antimicrobials[!is.na(antimicrobials$cid), ]),
+    length(unique(AMR::antimicrobials$cid[!is.na(antimicrobials$cid)]))
+  )
 
   # check cross table reference
   expect_true(all(microorganisms.codes$mo %in% microorganisms$mo))

tests/testthat/test-mdro.R

@@ -296,4 +296,55 @@ test_that("test-mdro.R", {
     expect_output(x <- mdro(example_isolates %>% group_by(ward), info = TRUE, pct_required_classes = 0))
     expect_output(x <- mdro(example_isolates %>% group_by(ward), guideline = custom, info = TRUE))
   }
+
+  # drug+inhibitor inference for missing base drug columns (issue #209) -------
+  # Resistance in drug+inhibitor implies resistance in the base drug.
+  # MRGN guideline is used because it explicitly requires PIP=R (not PIP OR TZP)
+  # for Pseudomonas aeruginosa 4MRGN, making the proxy effect directly testable.
+  pseud_no_pip <- data.frame(
+    mo = as.mo("Pseudomonas aeruginosa"),
+    TZP = as.sir("R"), # piperacillin/tazobactam; no PIP column
+    CAZ = as.sir("R"),
+    IPM = as.sir("R"),
+    MEM = as.sir("R"),
+    CIP = as.sir("R"),
+    stringsAsFactors = FALSE
+  )
+  # Inference message goes to message() / stderr, not stdout
+  # -> must use expect_message(), NOT expect_output()
+  expect_message(
+    suppressWarnings(mdro(pseud_no_pip, guideline = "mrgn", info = TRUE)),
+    "Inferring resistance"
+  )
+  inferred <- suppressWarnings(mdro(pseud_no_pip, guideline = "mrgn", info = FALSE))
+  not_inferred <- suppressWarnings(mdro(pseud_no_pip, guideline = "mrgn", info = FALSE, infer_from_combinations = FALSE))
+  expect_equal(as.character(inferred), "4MRGN")
+  expect_equal(as.character(not_inferred), "Negative")
+
+  # Susceptibility in combo does NOT propagate: proxy = NA, not S
+  # -> 4MRGN criteria no longer met -> lower level than when TZP=R
+  pseud_tzp_s <- pseud_no_pip
+  pseud_tzp_s$TZP <- as.sir("S")
+  result_tzp_s <- suppressMessages(suppressWarnings(
+    mdro(pseud_tzp_s, guideline = "mrgn", info = FALSE)
+  ))
+  expect_true(as.integer(result_tzp_s) < as.integer(inferred))
+
+  # Multiple combos for the same base drug: AMX can come from AMC (amoxi/clavulanic acid) and AXS (amoxi/sulbactam)
+  ente_no_amx <- data.frame(
+    mo = as.mo("Enterococcus faecium"),
+    AMC = as.sir("R"), # amoxicillin/clavulanic acid
+    AXS = as.sir("R"), # amoxicillin/sulbactam
+    VAN = as.sir("R"),
+    TEC = as.sir("R"),
+    LNZ = as.sir("R"),
+    DAP = as.sir("R"),
+    stringsAsFactors = FALSE
+  )
+  # Should have multiple columns in the verbose explanation
+  out <- mdro(ente_no_amx, guideline = "EUCAST 3.3", info = FALSE, verbose = TRUE)
+  expect_identical(
+    out$all_nonsusceptible_columns,
+    "AMC, AMX (inferred from AMC/AXS), AXS, DAP, LNZ, TEC, VAN"
+  )
 })

tests/testthat/test-zzz.R

@@ -178,9 +178,15 @@ test_that("test-zzz.R", {
       also_load = FALSE,
       min_version = if (pkg == "dplyr") "1.0.0" else NULL
     )) {
+      if (pkg == "rstudioapi") {
+        expect_true(is.function(tryCatch(get(fn, envir = asNamespace(pkg)), error = function(e) NULL)),
+          info = paste0("Function does not exist (anymore): function `", pkg, "::", fn, "()`")
+        )
+      } else {
         expect_true(!is.null(AMR:::import_fn(name = fn, pkg = pkg, error_on_fail = FALSE)),
           info = paste0("Function does not exist (anymore): function `", pkg, "::", fn, "()`")
         )
+      }
     } else if (pkg != "rstudioapi") {
       warning("Package '", pkg, "' not available")
     }