How to determine multi-drug resistance (MDR)

Matthijs S. Berends

28 May 2020

With the function mdro(), you can determine which micro-organisms are multi-drug resistant organisms (MDRO).

Type of input

The mdro() function takes a data set as input, such as a regular data.frame. It tries to automatically determine the right columns for info about your isolates, like the name of the species and all columns with results of antimicrobial agents. See the help page for more info about how to set the right settings for your data with the command ?mdro.

For WHONET data (and most other data), all settings are automatically set correctly.

Guidelines

The function support multiple guidelines. You can select a guideline with the guideline parameter. Currently supported guidelines are (case-insensitive):

Examples

The mdro() function always returns an ordered factor. For example, the output of the default guideline by Magiorakos et al. returns a factor with levels ‘Negative’, ‘MDR’, ‘XDR’ or ‘PDR’ in that order.

The next example uses the example_isolates data set. This is a data set included with this package and contains 2,000 microbial isolates with their full antibiograms. It reflects reality and can be used to practice AMR analysis. If we test the MDR/XDR/PDR guideline on this data set, we get:

library(dplyr)   # to support pipes: %>%
library(cleaner) # to create frequency tables
example_isolates %>% 
  mdro() %>% 
  freq() # show frequency table of the result
# NOTE: Using column `mo` as input for `col_mo`.
# NOTE: Auto-guessing columns suitable for analysis...OK.
# NOTE: Reliability would be improved if these antimicrobial results would be available too: ceftaroline (CPT), fusidic acid (FUS), telavancin (TLV), daptomycin (DAP), quinupristin/dalfopristin (QDA), minocycline (MNO), gentamicin-high (GEH), streptomycin-high (STH), doripenem (DOR), levofloxacin (LVX), netilmicin (NET), ticarcillin/clavulanic acid (TCC), ertapenem (ETP), cefotetan (CTT), aztreonam (ATM), ampicillin/sulbactam (SAM), polymyxin B (PLB)
# Warning in mdro(.): NA introduced for isolates where the available percentage of
# antimicrobial classes was below 50% (set with `pct_required_classes`)

Frequency table

Class: factor > ordered (numeric)
Length: 2,000
Levels: 4: Negative < Multi-drug-resistant (MDR) < Extensively drug-resistant …
Available: 1,711 (85.55%, NA: 289 = 14.45%)
Unique: 2

Item Count Percent Cum. Count Cum. Percent
1 Negative 1595 93.22% 1595 93.22%
2 Multi-drug-resistant (MDR) 116 6.78% 1711 100.00%

For another example, I will create a data set to determine multi-drug resistant TB:

# a helper function to get a random vector with values S, I and R
# with the probabilities 50% - 10% - 40%
sample_rsi <- function() {
  sample(c("S", "I", "R"),
         size = 5000,
         prob = c(0.5, 0.1, 0.4),
         replace = TRUE)
}

my_TB_data <- data.frame(rifampicin = sample_rsi(),
                         isoniazid = sample_rsi(),
                         gatifloxacin = sample_rsi(),
                         ethambutol = sample_rsi(),
                         pyrazinamide = sample_rsi(),
                         moxifloxacin = sample_rsi(),
                         kanamycin = sample_rsi())

Because all column names are automatically verified for valid drug names or codes, this would have worked exactly the same:

my_TB_data <- data.frame(RIF = sample_rsi(),
                         INH = sample_rsi(),
                         GAT = sample_rsi(),
                         ETH = sample_rsi(),
                         PZA = sample_rsi(),
                         MFX = sample_rsi(),
                         KAN = sample_rsi())

The data set now looks like this:

head(my_TB_data)
#   rifampicin isoniazid gatifloxacin ethambutol pyrazinamide moxifloxacin
# 1          S         S            R          R            S            S
# 2          S         S            S          I            S            R
# 3          R         R            R          I            R            R
# 4          S         R            S          S            S            I
# 5          I         R            S          I            R            S
# 6          S         S            R          R            S            I
#   kanamycin
# 1         S
# 2         R
# 3         S
# 4         S
# 5         R
# 6         I

We can now add the interpretation of MDR-TB to our data set. You can use:

mdro(my_TB_data, guideline = "TB")

or its shortcut mdr_tb():

my_TB_data$mdr <- mdr_tb(my_TB_data)
# NOTE: No column found as input for `col_mo`, assuming all records contain Mycobacterium tuberculosis.
# NOTE: Auto-guessing columns suitable for analysis...OK.
# NOTE: Reliability would be improved if these antimicrobial results would be available too: capreomycin (CAP), rifabutin (RIB), rifapentine (RFP)

Create a frequency table of the results:

freq(my_TB_data$mdr)

Frequency table

Class: factor > ordered (numeric)
Length: 5,000
Levels: 5: Negative < Mono-resistant < Poly-resistant < Multi-drug-resistant <…
Available: 5,000 (100%, NA: 0 = 0%)
Unique: 5

Item Count Percent Cum. Count Cum. Percent
1 Mono-resistant 3302 66.04% 3302 66.04%
2 Negative 634 12.68% 3936 78.72%
3 Multi-drug-resistant 578 11.56% 4514 90.28%
4 Poly-resistant 287 5.74% 4801 96.02%
5 Extensively drug-resistant 199 3.98% 5000 100.00%