Calculates a normalised mean for antimicrobial resistance between multiple observations, to help to identify similar isolates without comparing antibiograms by hand.
Usage
mean_amr_distance(x, ...)
# S3 method for rsi
mean_amr_distance(x, ..., combine_SI = TRUE)
# S3 method for data.frame
mean_amr_distance(x, ..., combine_SI = TRUE)
amr_distance_from_row(amr_distance, row)Arguments
- x
a vector of class rsi, mic or disk, or a data.frame containing columns of any of these classes
- ...
variables to select (supports tidyselect language such as
column1:column4andwhere(is.mic), and can thus also be antibiotic selectors- combine_SI
a logical to indicate whether all values of S and I must be merged into one, so the input only consists of S+I vs. R (susceptible vs. resistant), defaults to
TRUE- amr_distance
the outcome of
mean_amr_distance()- row
an index, such as a row number
Details
The mean AMR distance is effectively the Z-score; a normalised numeric value to compare AMR test results which can help to identify similar isolates, without comparing antibiograms by hand.
MIC values (see as.mic()) are transformed with log2() first; their distance is thus calculated as (log2(x) - mean(log2(x))) / sd(log2(x)).
R/SI values (see as.rsi()) are transformed using "S" = 1, "I" = 2, and "R" = 3. If combine_SI is TRUE (default), the "I" will be considered to be 1.
For data sets, the mean AMR distance will be calculated per column, after which the mean per row will be returned, see Examples.
Use amr_distance_from_row() to subtract distances from the distance of one row, see Examples.
Interpretation
Isolates with distances less than 0.01 difference from each other should be considered similar. Differences lower than 0.025 should be considered suspicious.
Examples
rsi <- random_rsi(10)
rsi
#> Class 'rsi'
#> [1] S I S R S I S S S S
mean_amr_distance(rsi)
#> [1] -0.3162278 -0.3162278 -0.3162278 2.8460499 -0.3162278 -0.3162278
#> [7] -0.3162278 -0.3162278 -0.3162278 -0.3162278
mic <- random_mic(10)
mic
#> Class 'mic'
#> [1] 0.005 32 >=64 0.001 32 0.005 >=64 0.001 >=64 32
mean_amr_distance(mic)
#> [1] -0.9977554 0.7020018 0.8364352 -1.3099002 0.7020018 -0.9977554
#> [7] 0.8364352 -1.3099002 0.8364352 0.7020018
# equal to the Z-score of their log2:
(log2(mic) - mean(log2(mic))) / sd(log2(mic))
#> [1] -0.9977554 0.7020018 0.8364352 -1.3099002 0.7020018 -0.9977554
#> [7] 0.8364352 -1.3099002 0.8364352 0.7020018
disk <- random_disk(10)
disk
#> Class 'disk'
#> [1] 43 13 6 50 43 19 14 23 8 44
mean_amr_distance(disk)
#> [1] 0.9878564 -0.7867359 -1.2008074 1.4019279 0.9878564 -0.4318175
#> [7] -0.7275828 -0.1952051 -1.0825013 1.0470094
y <- data.frame(
id = LETTERS[1:10],
amox = random_rsi(10, ab = "amox", mo = "Escherichia coli"),
cipr = random_disk(10, ab = "cipr", mo = "Escherichia coli"),
gent = random_mic(10, ab = "gent", mo = "Escherichia coli"),
tobr = random_mic(10, ab = "tobr", mo = "Escherichia coli")
)
y
#> id amox cipr gent tobr
#> 1 A I 22 4 2
#> 2 B I 24 0.5 2
#> 3 C R 30 2 1
#> 4 D S 29 0.5 1
#> 5 E I 30 4 2
#> 6 F R 19 0.25 4
#> 7 G I 26 4 4
#> 8 H S 22 2 2
#> 9 I S 29 2 1
#> 10 J I 28 0.5 1
mean_amr_distance(y)
#> ℹ Calculating mean AMR distance based on columns "amox", "cipr", "gent" and
#> "tobr"
#> [1] -0.03769893 -0.40858317 0.58133901 -0.40733070 0.47139759 0.01703862
#> [7] 0.53378217 -0.20375172 -0.07522511 -0.47096776
y$amr_distance <- mean_amr_distance(y, where(is.mic))
#> ℹ Calculating mean AMR distance based on columns "gent" and "tobr"
y[order(y$amr_distance), ]
#> id amox cipr gent tobr amr_distance
#> 4 D S 29 0.5 1 -0.97204037
#> 10 J I 28 0.5 1 -0.97204037
#> 2 B I 24 0.5 2 -0.33817468
#> 3 C R 30 2 1 -0.30782920
#> 9 I S 29 2 1 -0.30782920
#> 6 F R 19 0.25 4 -0.03641457
#> 8 H S 22 2 2 0.32603649
#> 1 A I 22 4 2 0.65814207
#> 5 E I 30 4 2 0.65814207
#> 7 G I 26 4 4 1.29200776
if (require("dplyr")) {
y %>%
mutate(
amr_distance = mean_amr_distance(y),
check_id_C = amr_distance_from_row(amr_distance, id == "C")
) %>%
arrange(check_id_C)
}
#> ℹ Calculating mean AMR distance based on columns "amox", "cipr", "gent" and
#> "tobr"
#> id amox cipr gent tobr amr_distance check_id_C
#> 1 C R 30 2 1 0.58133901 0.00000000
#> 2 G I 26 4 4 0.53378217 0.04755684
#> 3 E I 30 4 2 0.47139759 0.10994142
#> 4 F R 19 0.25 4 0.01703862 0.56430040
#> 5 A I 22 4 2 -0.03769893 0.61903794
#> 6 I S 29 2 1 -0.07522511 0.65656413
#> 7 H S 22 2 2 -0.20375172 0.78509073
#> 8 D S 29 0.5 1 -0.40733070 0.98866971
#> 9 B I 24 0.5 2 -0.40858317 0.98992219
#> 10 J I 28 0.5 1 -0.47096776 1.05230677
if (require("dplyr")) {
# support for groups
example_isolates %>%
filter(mo_genus() == "Enterococcus" & mo_species() != "") %>%
select(mo, TCY, carbapenems()) %>%
group_by(mo) %>%
mutate(dist = mean_amr_distance(.)) %>%
arrange(mo, dist)
}
#> ℹ Using column 'mo' as input for mo_genus()
#> ℹ Using column 'mo' as input for mo_species()
#> ℹ For carbapenems() using columns 'IPM' (imipenem) and 'MEM' (meropenem)
#> ℹ Calculating mean AMR distance based on columns "TCY", "IPM" and "MEM"
#> # A tibble: 63 × 5
#> # Groups: mo [4]
#> mo TCY IPM MEM dist
#> <mo> <rsi> <rsi> <rsi> <dbl>
#> 1 B_ENTRC_AVIM S S NA 0
#> 2 B_ENTRC_AVIM S S NA 0
#> 3 B_ENTRC_CSSL NA S NA NA
#> 4 B_ENTRC_FACM S S NA -2.66
#> 5 B_ENTRC_FACM S R R -0.423
#> 6 B_ENTRC_FACM S R R -0.423
#> 7 B_ENTRC_FACM NA R R 0.224
#> 8 B_ENTRC_FACM NA R R 0.224
#> 9 B_ENTRC_FACM NA R R 0.224
#> 10 B_ENTRC_FACM NA R R 0.224
#> # … with 53 more rows