AMR/R/zzz.R

68 lines
3.7 KiB
R
Raw Normal View History

# ==================================================================== #
# TITLE #
# Antimicrobial Resistance (AMR) Analysis #
# #
2019-01-02 23:24:07 +01:00
# SOURCE #
2020-07-08 14:48:06 +02:00
# https://github.com/msberends/AMR #
# #
# LICENCE #
# (c) 2018-2020 Berends MS, Luz CF et al. #
# #
2019-01-02 23:24:07 +01: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. #
2020-07-08 14:48:06 +02:00
# Visit our website for more info: https://msberends.github.io/AMR. #
# ==================================================================== #
2018-04-19 14:10:57 +02:00
.onLoad <- function(libname, pkgname) {
2020-05-16 13:05:47 +02:00
assign(x = "MO_lookup",
value = create_MO_lookup(),
2019-03-12 12:19:27 +01:00
envir = asNamespace("AMR"))
2020-05-16 13:05:47 +02:00
assign(x = "MO.old_lookup",
value = create_MO.old_lookup(),
2019-03-12 12:19:27 +01:00
envir = asNamespace("AMR"))
2018-04-19 14:10:57 +02:00
}
2019-02-21 23:32:30 +01:00
2019-10-11 17:21:02 +02:00
# maybe add survey later: "https://www.surveymonkey.com/r/AMR_for_R"
2020-05-16 13:05:47 +02:00
create_MO_lookup <- function() {
MO_lookup <- AMR::microorganisms
MO_lookup$kingdom_index <- 99
MO_lookup[which(MO_lookup$kingdom == "Bacteria" | MO_lookup$mo == "UNKNOWN"), "kingdom_index"] <- 1
MO_lookup[which(MO_lookup$kingdom == "Fungi"), "kingdom_index"] <- 2
MO_lookup[which(MO_lookup$kingdom == "Protozoa"), "kingdom_index"] <- 3
MO_lookup[which(MO_lookup$kingdom == "Archaea"), "kingdom_index"] <- 4
# use this paste instead of `fullname` to work with Viridans Group Streptococci, etc.
2020-05-16 13:05:47 +02:00
MO_lookup$fullname_lower <- tolower(trimws(paste(MO_lookup$genus,
MO_lookup$species,
MO_lookup$subspecies)))
2020-05-16 13:05:47 +02:00
MO_lookup[MO_lookup$genus == "" | grepl("^[(]unknown ", MO_lookup$fullname), "fullname_lower"] <- tolower(trimws(MO_lookup[MO_lookup$genus == "" | grepl("^[(]unknown ", MO_lookup$fullname),
"fullname"]))
2020-05-18 13:59:34 +02:00
MO_lookup$fullname_lower <- gsub("[^.a-z0-9/ \\-]+", "", MO_lookup$fullname_lower)
2020-05-16 13:05:47 +02:00
# add a column with only "e coli" like combinations
MO_lookup$g_species <- gsub("^([a-z])[a-z]+ ([a-z]+) ?.*", "\\1 \\2", MO_lookup$fullname_lower)
2019-11-15 16:50:46 +01:00
2019-11-15 15:25:03 +01:00
# so arrange data on prevalence first, then kingdom, then full name
2020-05-18 13:59:34 +02:00
MO_lookup[order(MO_lookup$prevalence, MO_lookup$kingdom_index, MO_lookup$fullname_lower), ]
2019-02-21 23:32:30 +01:00
}
2020-05-16 13:05:47 +02:00
create_MO.old_lookup <- function() {
MO.old_lookup <- AMR::microorganisms.old
MO.old_lookup$fullname_lower <- gsub("[^.a-z0-9/ \\-]+", "", tolower(trimws(MO.old_lookup$fullname)))
# add a column with only "e coli" like combinations
MO.old_lookup$g_species <- gsub("^([a-z])[a-z]+ ([a-z]+) ?.*", "\\1 \\2", MO.old_lookup$fullname_lower)
2019-12-21 10:56:06 +01:00
# so arrange data on prevalence first, then full name
2020-05-18 13:59:34 +02:00
MO.old_lookup[order(MO.old_lookup$prevalence, MO.old_lookup$fullname_lower), ]
2019-12-21 10:56:06 +01:00
}