2021-05-15 21:36:22 +02:00
|
|
|
# ==================================================================== #
|
2023-07-08 17:30:05 +02:00
|
|
|
# TITLE: #
|
2022-10-05 09:12:22 +02:00
|
|
|
# AMR: An R Package for Working with Antimicrobial Resistance Data #
|
2021-05-15 21:36:22 +02:00
|
|
|
# #
|
2023-07-08 17:30:05 +02:00
|
|
|
# SOURCE CODE: #
|
2021-05-15 21:36:22 +02:00
|
|
|
# https://github.com/msberends/AMR #
|
|
|
|
# #
|
2023-07-08 17:30:05 +02:00
|
|
|
# PLEASE CITE THIS SOFTWARE AS: #
|
2022-10-05 09:12:22 +02:00
|
|
|
# Berends MS, Luz CF, Friedrich AW, Sinha BNM, Albers CJ, Glasner C #
|
|
|
|
# (2022). AMR: An R Package for Working with Antimicrobial Resistance #
|
|
|
|
# Data. Journal of Statistical Software, 104(3), 1-31. #
|
2023-05-27 10:39:22 +02:00
|
|
|
# https://doi.org/10.18637/jss.v104.i03 #
|
2022-10-05 09:12:22 +02:00
|
|
|
# #
|
2022-12-27 15:16:15 +01:00
|
|
|
# Developed at the University of Groningen and the University Medical #
|
|
|
|
# Center Groningen in The Netherlands, in collaboration with many #
|
|
|
|
# colleagues from around the world, see our website. #
|
2021-05-15 21:36:22 +02:00
|
|
|
# #
|
|
|
|
# This R package is free software; you can freely use and distribute #
|
|
|
|
# it for both personal and commercial purposes under the terms of the #
|
|
|
|
# GNU General Public License version 2.0 (GNU GPL-2), as published by #
|
|
|
|
# the Free Software Foundation. #
|
|
|
|
# We created this package for both routine data analysis and academic #
|
|
|
|
# research and it was publicly released in the hope that it will be #
|
|
|
|
# useful, but it comes WITHOUT ANY WARRANTY OR LIABILITY. #
|
|
|
|
# #
|
|
|
|
# Visit our website for the full manual and a complete tutorial about #
|
|
|
|
# how to conduct AMR data analysis: https://msberends.github.io/AMR/ #
|
|
|
|
# ==================================================================== #
|
|
|
|
|
2023-02-18 14:56:06 +01:00
|
|
|
if (AMR:::pkg_is_available("dplyr", min_version = "1.0.0", also_load = TRUE) &&
|
|
|
|
AMR:::pkg_is_available("ggplot2", also_load = TRUE)) {
|
2021-05-15 21:36:22 +02:00
|
|
|
pdf(NULL) # prevent Rplots.pdf being created
|
2022-08-28 10:31:50 +02:00
|
|
|
|
2021-05-15 21:36:22 +02:00
|
|
|
# data should be equal
|
|
|
|
expect_equal(
|
|
|
|
(example_isolates %>%
|
2022-08-28 10:31:50 +02:00
|
|
|
select(AMC, CIP) %>%
|
2023-01-21 23:47:20 +01:00
|
|
|
ggplot_sir())$data %>%
|
2021-05-15 21:36:22 +02:00
|
|
|
summarise_all(resistance) %>%
|
|
|
|
as.double(),
|
|
|
|
example_isolates %>%
|
|
|
|
select(AMC, CIP) %>%
|
|
|
|
summarise_all(resistance) %>%
|
|
|
|
as.double()
|
|
|
|
)
|
2022-08-28 10:31:50 +02:00
|
|
|
|
2021-05-15 21:36:22 +02:00
|
|
|
expect_stdout(print(example_isolates %>%
|
2022-08-28 10:31:50 +02:00
|
|
|
select(AMC, CIP) %>%
|
2023-01-21 23:47:20 +01:00
|
|
|
ggplot_sir(x = "interpretation", facet = "antibiotic")))
|
2021-05-15 21:36:22 +02:00
|
|
|
expect_stdout(print(example_isolates %>%
|
2022-08-28 10:31:50 +02:00
|
|
|
select(AMC, CIP) %>%
|
2023-01-21 23:47:20 +01:00
|
|
|
ggplot_sir(x = "antibiotic", facet = "interpretation")))
|
2022-08-28 10:31:50 +02:00
|
|
|
|
2021-05-15 21:36:22 +02:00
|
|
|
expect_equal(
|
|
|
|
(example_isolates %>%
|
2022-08-28 10:31:50 +02:00
|
|
|
select(AMC, CIP) %>%
|
2023-01-21 23:47:20 +01:00
|
|
|
ggplot_sir(x = "interpretation", facet = "antibiotic"))$data %>%
|
2021-05-15 21:36:22 +02:00
|
|
|
summarise_all(resistance) %>%
|
|
|
|
as.double(),
|
|
|
|
example_isolates %>%
|
|
|
|
select(AMC, CIP) %>%
|
|
|
|
summarise_all(resistance) %>%
|
|
|
|
as.double()
|
|
|
|
)
|
2022-08-28 10:31:50 +02:00
|
|
|
|
2021-05-15 21:36:22 +02:00
|
|
|
expect_equal(
|
|
|
|
(example_isolates %>%
|
2022-08-28 10:31:50 +02:00
|
|
|
select(AMC, CIP) %>%
|
2023-01-21 23:47:20 +01:00
|
|
|
ggplot_sir(x = "antibiotic", facet = "interpretation"))$data %>%
|
2021-05-15 21:36:22 +02:00
|
|
|
summarise_all(resistance) %>%
|
|
|
|
as.double(),
|
|
|
|
example_isolates %>%
|
|
|
|
select(AMC, CIP) %>%
|
|
|
|
summarise_all(resistance) %>%
|
|
|
|
as.double()
|
|
|
|
)
|
2022-08-28 10:31:50 +02:00
|
|
|
|
2021-05-15 21:36:22 +02:00
|
|
|
expect_equal(
|
|
|
|
(example_isolates %>%
|
2022-08-28 10:31:50 +02:00
|
|
|
select(AMC, CIP) %>%
|
2023-01-21 23:47:20 +01:00
|
|
|
ggplot_sir(x = "antibiotic", facet = "interpretation"))$data %>%
|
2021-05-15 21:36:22 +02:00
|
|
|
summarise_all(count_resistant) %>%
|
|
|
|
as.double(),
|
|
|
|
example_isolates %>%
|
|
|
|
select(AMC, CIP) %>%
|
|
|
|
summarise_all(count_resistant) %>%
|
|
|
|
as.double()
|
|
|
|
)
|
2022-08-28 10:31:50 +02:00
|
|
|
|
2021-05-15 21:36:22 +02:00
|
|
|
# support for scale_type ab and mo
|
2022-08-28 10:31:50 +02:00
|
|
|
expect_inherits(
|
|
|
|
(data.frame(
|
|
|
|
mo = as.mo(c("e. coli", "s aureus")),
|
|
|
|
n = c(40, 100)
|
|
|
|
) %>%
|
|
|
|
ggplot(aes(x = mo, y = n)) +
|
|
|
|
geom_col())$data,
|
|
|
|
"data.frame"
|
|
|
|
)
|
|
|
|
expect_inherits(
|
|
|
|
(data.frame(
|
|
|
|
ab = as.ab(c("amx", "amc")),
|
|
|
|
n = c(40, 100)
|
|
|
|
) %>%
|
|
|
|
ggplot(aes(x = ab, y = n)) +
|
|
|
|
geom_col())$data,
|
|
|
|
"data.frame"
|
|
|
|
)
|
|
|
|
|
|
|
|
expect_inherits(
|
|
|
|
(data.frame(
|
|
|
|
ab = as.ab(c("amx", "amc")),
|
|
|
|
n = c(40, 100)
|
|
|
|
) %>%
|
|
|
|
ggplot(aes(x = ab, y = n)) +
|
|
|
|
geom_col())$data,
|
|
|
|
"data.frame"
|
|
|
|
)
|
|
|
|
|
2021-05-15 21:36:22 +02:00
|
|
|
# support for manual colours
|
2022-08-28 10:31:50 +02:00
|
|
|
expect_inherits(
|
|
|
|
(ggplot(data.frame(
|
|
|
|
x = c("Value1", "Value2", "Value3"),
|
|
|
|
y = c(1, 2, 3),
|
|
|
|
z = c("Value4", "Value5", "Value6")
|
|
|
|
)) +
|
|
|
|
geom_col(aes(x = x, y = y, fill = z)) +
|
2023-01-21 23:47:20 +01:00
|
|
|
scale_sir_colours(Value4 = "S", Value5 = "I", Value6 = "R"))$data,
|
2022-08-28 10:31:50 +02:00
|
|
|
"data.frame"
|
|
|
|
)
|
2021-05-15 21:36:22 +02:00
|
|
|
}
|