@@ -401,19 +401,19 @@ names or codes, this would have worked exactly the same way:
head ( my_TB_data )
#> rifampicin isoniazid gatifloxacin ethambutol pyrazinamide moxifloxacin
-#> 1 I I S I I I
-#> 2 I I I I I I
-#> 3 I S R R S S
-#> 4 I I R I R I
-#> 5 I S I I R I
-#> 6 S I R R R I
+#> 1 S I R I S I
+#> 2 S R I I R I
+#> 3 S I R R S I
+#> 4 S S I R S I
+#> 5 S R S I S S
+#> 6 S I R I R R
#> kanamycin
-#> 1 I
-#> 2 I
-#> 3 S
-#> 4 R
+#> 1 S
+#> 2 S
+#> 3 I
+#> 4 I
#> 5 S
-#> 6 S
+#> 6 R
We can now add the interpretation of MDR-TB to our data set. You can
use:
@@ -454,32 +454,32 @@ Unique: 5
1
Mono-resistant
-3223
-64.46%
-3223
-64.46%
+3233
+64.66%
+3233
+64.66%
2
Negative
-1033
-20.66%
-4256
-85.12%
+945
+18.90%
+4178
+83.56%
3
Multi-drug-resistant
-402
-8.04%
-4658
-93.16%
+446
+8.92%
+4624
+92.48%
4
Poly-resistant
-241
-4.82%
+275
+5.50%
4899
97.98%
diff --git a/articles/PCA.html b/articles/PCA.html
index aaa0a342..b94fd2b6 100644
--- a/articles/PCA.html
+++ b/articles/PCA.html
@@ -14,8 +14,8 @@
-
-
+
+
@@ -38,7 +38,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -307,7 +307,7 @@ function:
#> [1] "Caryophanales" "Enterobacterales" "Lactobacillales" "Pseudomonadales"
#> Importance of components:
#> PC1 PC2 PC3 PC4 PC5 PC6 PC7
-#> Standard deviation 2.1539 1.6807 0.6138 0.33879 0.20808 0.03140 9.577e-17
+#> Standard deviation 2.1539 1.6807 0.6138 0.33879 0.20808 0.03140 1.232e-16
#> Proportion of Variance 0.5799 0.3531 0.0471 0.01435 0.00541 0.00012 0.000e+00
#> Cumulative Proportion 0.5799 0.9330 0.9801 0.99446 0.99988 1.00000 1.000e+00
#> Groups (n=4, named as 'order'):
diff --git a/articles/WHONET.html b/articles/WHONET.html
index f1ae169f..ef94a4dd 100644
--- a/articles/WHONET.html
+++ b/articles/WHONET.html
@@ -14,8 +14,8 @@
-
-
+
+
@@ -38,7 +38,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/articles/datasets.html b/articles/datasets.html
index 0bf7ce87..6d615cdd 100644
--- a/articles/datasets.html
+++ b/articles/datasets.html
@@ -14,8 +14,8 @@
-
-
+
+
@@ -38,7 +38,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -193,7 +193,7 @@
Automatically, using the rOpenSci R-universe platform , by adding our R-universe address to your list of repositories (‘repos’):
diff --git a/news/index.html b/news/index.html
index e4756377..f9214a5b 100644
--- a/news/index.html
+++ b/news/index.html
@@ -1,5 +1,5 @@
-Changelog • AMR (for R) Changelog • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -158,6 +158,20 @@
Source: NEWS.md
+
+
AMR 2.1.1.9001
+
+
New
+
Function scale_x_mic()
, an advanced function to use in ggplot, to allow plotting of MIC values on the x axis. It allow for manual range definition and plotting missing intermediate log2 levels.
+
+
Changed
+
For MICs:
+Added 4096 and 5 powers of 192 as valid levels (192, 384, 576, 768, 960)
+Added new argument keep_operators
to as.mic()
. This can be "all"
(default), "none"
, or "edges"
.
+
+
+
+
AMR 2.1.1 CRAN release: 2023-10-21
Fix for selecting first isolates using the phenotype-based method
diff --git a/pkgdown.yml b/pkgdown.yml
index f357feeb..2cb81160 100644
--- a/pkgdown.yml
+++ b/pkgdown.yml
@@ -11,7 +11,7 @@ articles:
other_pkg: other_pkg.html
resistance_predict: resistance_predict.html
welcome_to_AMR: welcome_to_AMR.html
-last_built: 2023-10-22T13:55Z
+last_built: 2023-12-03T00:11Z
urls:
reference: https://msberends.github.io/AMR/reference
article: https://msberends.github.io/AMR/articles
diff --git a/reference/AMR-deprecated.html b/reference/AMR-deprecated.html
index b5afd08b..fdb18a9b 100644
--- a/reference/AMR-deprecated.html
+++ b/reference/AMR-deprecated.html
@@ -1,5 +1,5 @@
-Deprecated Functions — AMR-deprecated • AMR (for R) Deprecated Functions — AMR-deprecated • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/AMR-options.html b/reference/AMR-options.html
index 0e8f240a..8d3f3a07 100644
--- a/reference/AMR-options.html
+++ b/reference/AMR-options.html
@@ -1,5 +1,5 @@
-Options for the AMR package — AMR-options • AMR (for R) Options for the AMR package — AMR-options • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/AMR.html b/reference/AMR.html
index 35775906..b68ccf39 100644
--- a/reference/AMR.html
+++ b/reference/AMR.html
@@ -6,7 +6,7 @@ This work was published in the Journal of Statistical Software (Volume 104(3); d
and doi:10.33612/diss.192486375
).
After installing this package, R knows ~52 000 microorganisms (updated December 2022) and all ~600 antibiotic, antimycotic and antiviral drugs by name and code (including ATC, EARS-Net, ASIARS-Net, PubChem, LOINC and SNOMED CT), and knows all about valid SIR and MIC values. The integral clinical breakpoint guidelines from CLSI and EUCAST are included, even with epidemiological cut-off (ECOFF) values. It supports and can read any data format, including WHONET data. This package works on Windows, macOS and Linux with all versions of R since R-3.0 (April 2013). It was designed to work in any setting, including those with very limited resources. It was created for both routine data analysis and academic research at the Faculty of Medical Sciences of the public University of Groningen, in collaboration with non-profit organisations Certe Medical Diagnostics and Advice Foundation and University Medical Center Groningen.
-The AMR package is available in English, Chinese, Czech, Danish, Dutch, Finnish, French, German, Greek, Italian, Japanese, Norwegian, Polish, Portuguese, Romanian, Russian, Spanish, Swedish, Turkish, and Ukrainian. Antimicrobial drug (group) names and colloquial microorganism names are provided in these languages.">The AMR Package — AMR • AMR (for R) The AMR Package — AMR • AMR (for R) AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/Rplot002.png b/reference/Rplot002.png
index 8c752534..971d6676 100644
Binary files a/reference/Rplot002.png and b/reference/Rplot002.png differ
diff --git a/reference/Rplot003.png b/reference/Rplot003.png
index 73fc7bc8..ad853f36 100644
Binary files a/reference/Rplot003.png and b/reference/Rplot003.png differ
diff --git a/reference/Rplot004.png b/reference/Rplot004.png
index 1414bc34..4e8840a7 100644
Binary files a/reference/Rplot004.png and b/reference/Rplot004.png differ
diff --git a/reference/Rplot005.png b/reference/Rplot005.png
index 2e4f99f9..62fdfee5 100644
Binary files a/reference/Rplot005.png and b/reference/Rplot005.png differ
diff --git a/reference/Rplot006.png b/reference/Rplot006.png
index af5a59f3..3859203c 100644
Binary files a/reference/Rplot006.png and b/reference/Rplot006.png differ
diff --git a/reference/Rplot007.png b/reference/Rplot007.png
index 5faaddb2..d0d49549 100644
Binary files a/reference/Rplot007.png and b/reference/Rplot007.png differ
diff --git a/reference/Rplot008.png b/reference/Rplot008.png
index df21079c..2c3980f9 100644
Binary files a/reference/Rplot008.png and b/reference/Rplot008.png differ
diff --git a/reference/Rplot009.png b/reference/Rplot009.png
index 73fc7bc8..b88f3ad8 100644
Binary files a/reference/Rplot009.png and b/reference/Rplot009.png differ
diff --git a/reference/Rplot010.png b/reference/Rplot010.png
index 46064676..08ff673d 100644
Binary files a/reference/Rplot010.png and b/reference/Rplot010.png differ
diff --git a/reference/Rplot011.png b/reference/Rplot011.png
new file mode 100644
index 00000000..7ac9eb77
Binary files /dev/null and b/reference/Rplot011.png differ
diff --git a/reference/Rplot012.png b/reference/Rplot012.png
new file mode 100644
index 00000000..21c3a434
Binary files /dev/null and b/reference/Rplot012.png differ
diff --git a/reference/Rplot013.png b/reference/Rplot013.png
new file mode 100644
index 00000000..b3e0b9b5
Binary files /dev/null and b/reference/Rplot013.png differ
diff --git a/reference/WHOCC.html b/reference/WHOCC.html
index 417a80d2..aa57037a 100644
--- a/reference/WHOCC.html
+++ b/reference/WHOCC.html
@@ -1,5 +1,5 @@
-WHOCC: WHO Collaborating Centre for Drug Statistics Methodology — WHOCC • AMR (for R) WHOCC: WHO Collaborating Centre for Drug Statistics Methodology — WHOCC • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/WHONET.html b/reference/WHONET.html
index e801aecc..a0017d78 100644
--- a/reference/WHONET.html
+++ b/reference/WHONET.html
@@ -1,5 +1,5 @@
-Data Set with 500 Isolates - WHONET Example — WHONET • AMR (for R) Data Set with 500 Isolates - WHONET Example — WHONET • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/ab_from_text.html b/reference/ab_from_text.html
index f1c1f2c8..d4381264 100644
--- a/reference/ab_from_text.html
+++ b/reference/ab_from_text.html
@@ -1,5 +1,5 @@
-Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text • AMR (for R) Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/ab_property.html b/reference/ab_property.html
index f0c93ac2..73a77ecc 100644
--- a/reference/ab_property.html
+++ b/reference/ab_property.html
@@ -1,5 +1,5 @@
-Get Properties of an Antibiotic — ab_property • AMR (for R) Get Properties of an Antibiotic — ab_property • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/add_custom_antimicrobials.html b/reference/add_custom_antimicrobials.html
index 21084e61..63daee96 100644
--- a/reference/add_custom_antimicrobials.html
+++ b/reference/add_custom_antimicrobials.html
@@ -1,5 +1,5 @@
-Add Custom Antimicrobials — add_custom_antimicrobials • AMR (for R) Add Custom Antimicrobials — add_custom_antimicrobials • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/add_custom_microorganisms.html b/reference/add_custom_microorganisms.html
index dfd47a1a..e45d1d18 100644
--- a/reference/add_custom_microorganisms.html
+++ b/reference/add_custom_microorganisms.html
@@ -1,5 +1,5 @@
-Add Custom Microorganisms — add_custom_microorganisms • AMR (for R) Add Custom Microorganisms — add_custom_microorganisms • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/age.html b/reference/age.html
index d58a5b5d..f918976a 100644
--- a/reference/age.html
+++ b/reference/age.html
@@ -1,5 +1,5 @@
-Age in Years of Individuals — age • AMR (for R) Age in Years of Individuals — age • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -222,16 +222,16 @@
df
#> birth_date age age_exact age_at_y2k
-#> 1 1950-06-04 73 73.38356 49
-#> 2 1953-03-27 70 70.57260 46
-#> 3 1968-11-29 54 54.89589 31
-#> 4 1997-09-02 26 26.13699 2
-#> 5 1981-09-03 42 42.13425 18
-#> 6 1974-02-03 49 49.71507 25
-#> 7 1947-03-27 76 76.57260 52
-#> 8 1972-02-11 51 51.69315 27
-#> 9 1959-12-06 63 63.87671 40
-#> 10 1935-10-25 87 87.99178 64
+#> 1 1968-10-06 55 55.15890 31
+#> 2 1954-09-27 69 69.18356 45
+#> 3 1984-09-22 39 39.19726 15
+#> 4 1939-07-03 84 84.41918 60
+#> 5 1984-07-23 39 39.36438 15
+#> 6 1998-04-21 25 25.61918 1
+#> 7 1940-12-19 82 82.95616 59
+#> 8 1993-02-19 30 30.78630 6
+#> 9 1968-03-29 55 55.68219 31
+#> 10 1969-09-11 54 54.22740 30
On this page
diff --git a/reference/age_groups.html b/reference/age_groups.html
index 050adfcb..223e4603 100644
--- a/reference/age_groups.html
+++ b/reference/age_groups.html
@@ -1,5 +1,5 @@
-Split Ages into Age Groups — age_groups • AMR (for R) Split Ages into Age Groups — age_groups • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/antibiogram.html b/reference/antibiogram.html
index 2e68d7a5..5487c965 100644
--- a/reference/antibiogram.html
+++ b/reference/antibiogram.html
@@ -1,5 +1,5 @@
-Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram • AMR (for R) Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/antibiotic_class_selectors.html b/reference/antibiotic_class_selectors.html
index 294d2377..3bda475d 100644
--- a/reference/antibiotic_class_selectors.html
+++ b/reference/antibiotic_class_selectors.html
@@ -1,6 +1,6 @@
Antibiotic Selectors — antibiotic_class_selectors • AMR (for R) Antibiotic Selectors — antibiotic_class_selectors • AMR (for R) Data Sets with 603 Antimicrobial Drugs — antibiotics • AMR (for R) Data Sets with 603 Antimicrobial Drugs — antibiotics • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/as.ab.html b/reference/as.ab.html
index a68f3046..5df52c2a 100644
--- a/reference/as.ab.html
+++ b/reference/as.ab.html
@@ -1,5 +1,5 @@
-Transform Input to an Antibiotic ID — as.ab • AMR (for R) Transform Input to an Antibiotic ID — as.ab • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/as.av.html b/reference/as.av.html
index 1160b583..7daa7d57 100644
--- a/reference/as.av.html
+++ b/reference/as.av.html
@@ -1,5 +1,5 @@
-Transform Input to an Antiviral Drug ID — as.av • AMR (for R) Transform Input to an Antiviral Drug ID — as.av • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/as.disk.html b/reference/as.disk.html
index 436c6a4f..725daf5d 100644
--- a/reference/as.disk.html
+++ b/reference/as.disk.html
@@ -1,5 +1,5 @@
-Transform Input to Disk Diffusion Diameters — as.disk • AMR (for R) Transform Input to Disk Diffusion Diameters — as.disk • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/as.mic.html b/reference/as.mic.html
index 37581e9a..07938b15 100644
--- a/reference/as.mic.html
+++ b/reference/as.mic.html
@@ -1,5 +1,5 @@
-Transform Input to Minimum Inhibitory Concentrations (MIC) — as.mic • AMR (for R) Transform Input to Minimum Inhibitory Concentrations (MIC) — as.mic • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -165,7 +165,7 @@
Usage
-
as.mic ( x , na.rm = FALSE )
+ as.mic ( x , na.rm = FALSE , keep_operators = "all" )
NA_mic_
@@ -185,6 +185,10 @@
a logical indicating whether missing values should be removed
+keep_operators
+a character specifying how to handle operators (such as >
and <=
) in the input. Accepts one of three values: "all"
(or TRUE
) to keep all operators, "none"
(or FALSE
) to remove all operators, or "edges"
to keep operators only at both ends of the range.
+
+
as.mic
a logical to indicate whether the mic
class should be kept - the default is FALSE
diff --git a/reference/as.mo.html b/reference/as.mo.html
index afd81595..e0f19944 100644
--- a/reference/as.mo.html
+++ b/reference/as.mo.html
@@ -1,5 +1,5 @@
-Transform Arbitrary Input to Valid Microbial Taxonomy — as.mo • AMR (for R) Transform Arbitrary Input to Valid Microbial Taxonomy — as.mo • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/as.sir.html b/reference/as.sir.html
index 04cdfbfd..b66e653e 100644
--- a/reference/as.sir.html
+++ b/reference/as.sir.html
@@ -1,7 +1,7 @@
Translate MIC and Disk Diffusion to SIR, or Clean Existing SIR Data — as.sir • AMR (for R) Translate MIC and Disk Diffusion to SIR, or Clean Existing SIR Data — as.sir • AMR (for R) Get ATC Properties from WHOCC Website — atc_online_property • AMR (for R) Get ATC Properties from WHOCC Website — atc_online_property • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/av_from_text.html b/reference/av_from_text.html
index 5a427a30..16aae058 100644
--- a/reference/av_from_text.html
+++ b/reference/av_from_text.html
@@ -1,5 +1,5 @@
-Retrieve Antiviral Drug Names and Doses from Clinical Text — av_from_text • AMR (for R) Retrieve Antiviral Drug Names and Doses from Clinical Text — av_from_text • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/av_property.html b/reference/av_property.html
index 5ba5c418..fa8a2fb5 100644
--- a/reference/av_property.html
+++ b/reference/av_property.html
@@ -1,5 +1,5 @@
-Get Properties of an Antiviral Drug — av_property • AMR (for R) Get Properties of an Antiviral Drug — av_property • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/availability.html b/reference/availability.html
index 783abe41..a8045977 100644
--- a/reference/availability.html
+++ b/reference/availability.html
@@ -1,5 +1,5 @@
-Check Availability of Columns — availability • AMR (for R) Check Availability of Columns — availability • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/bug_drug_combinations.html b/reference/bug_drug_combinations.html
index 794022c6..8b6e4ca5 100644
--- a/reference/bug_drug_combinations.html
+++ b/reference/bug_drug_combinations.html
@@ -1,5 +1,5 @@
-Determine Bug-Drug Combinations — bug_drug_combinations • AMR (for R) Determine Bug-Drug Combinations — bug_drug_combinations • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/clinical_breakpoints.html b/reference/clinical_breakpoints.html
index aa50ed1f..19afa6f8 100644
--- a/reference/clinical_breakpoints.html
+++ b/reference/clinical_breakpoints.html
@@ -1,5 +1,5 @@
-Data Set with Clinical Breakpoints for SIR Interpretation — clinical_breakpoints • AMR (for R) Data Set with Clinical Breakpoints for SIR Interpretation — clinical_breakpoints • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/count.html b/reference/count.html
index e7107bdf..e9ab79f4 100644
--- a/reference/count.html
+++ b/reference/count.html
@@ -1,6 +1,6 @@
Count Available Isolates — count • AMR (for R) Count Available Isolates — count • AMR (for R) Define Custom EUCAST Rules — custom_eucast_rules • AMR (for R) Define Custom EUCAST Rules — custom_eucast_rules • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/dosage.html b/reference/dosage.html
index aed85f6a..84a16cd9 100644
--- a/reference/dosage.html
+++ b/reference/dosage.html
@@ -1,5 +1,5 @@
-Data Set with Treatment Dosages as Defined by EUCAST — dosage • AMR (for R) Data Set with Treatment Dosages as Defined by EUCAST — dosage • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/eucast_rules.html b/reference/eucast_rules.html
index 6bd1d2f4..0363c376 100644
--- a/reference/eucast_rules.html
+++ b/reference/eucast_rules.html
@@ -1,6 +1,6 @@
Apply EUCAST Rules — eucast_rules • AMR (for R) Apply EUCAST Rules — eucast_rules • AMR (for R) Data Set with 2 000 Example Isolates — example_isolates • AMR (for R) Data Set with 2 000 Example Isolates — example_isolates • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/example_isolates_unclean.html b/reference/example_isolates_unclean.html
index 93ff9671..fc337e91 100644
--- a/reference/example_isolates_unclean.html
+++ b/reference/example_isolates_unclean.html
@@ -1,5 +1,5 @@
-Data Set with Unclean Data — example_isolates_unclean • AMR (for R) Data Set with Unclean Data — example_isolates_unclean • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/first_isolate.html b/reference/first_isolate.html
index cfa9b4db..3cb3c794 100644
--- a/reference/first_isolate.html
+++ b/reference/first_isolate.html
@@ -1,6 +1,6 @@
Determine First Isolates — first_isolate • AMR (for R) Determine First Isolates — first_isolate • AMR (for R) G-test for Count Data — g.test • AMR (for R) G-test for Count Data — g.test • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/get_episode.html b/reference/get_episode.html
index 06a421a9..b5e6f7c6 100644
--- a/reference/get_episode.html
+++ b/reference/get_episode.html
@@ -1,5 +1,5 @@
-Determine Clinical or Epidemic Episodes — get_episode • AMR (for R) Determine Clinical or Epidemic Episodes — get_episode • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -263,27 +263,28 @@
df <- example_isolates [ sample ( seq_len ( 2000 ) , size = 100 ) , ]
get_episode ( df $ date , episode_days = 60 ) # indices
-#> [1] 39 42 25 32 10 25 17 9 4 9 41 18 13 36 30 24 41 14 38 12 23 12 24 25 42
-#> [26] 20 26 19 16 14 4 5 31 9 20 26 22 17 5 6 30 18 34 27 43 28 37 29 38 31
-#> [51] 11 3 41 39 2 42 33 11 34 40 16 27 13 23 4 11 9 10 1 41 30 24 34 8 11
-#> [76] 30 35 31 2 32 9 40 11 17 35 21 7 38 15 42 21 38 39 44 1 15 22 11 22 12
+#> [1] 8 23 24 26 5 8 42 39 25 48 36 19 7 26 4 37 29 15 40 5 16 29 3 34 44
+#> [26] 16 3 48 1 9 10 26 38 13 20 30 31 7 18 5 16 22 2 2 41 2 4 18 37 43
+#> [51] 33 28 10 16 41 15 43 28 12 12 9 47 21 8 46 7 10 8 9 27 30 30 28 8 26
+#> [76] 13 44 10 16 35 38 42 17 32 10 39 34 35 25 13 12 19 26 14 13 6 11 45 7 27
is_new_episode ( df $ date , episode_days = 60 ) # TRUE/FALSE
-#> [1] TRUE TRUE TRUE TRUE TRUE FALSE TRUE TRUE TRUE FALSE TRUE TRUE
-#> [13] TRUE TRUE TRUE TRUE FALSE TRUE TRUE TRUE TRUE FALSE FALSE FALSE
-#> [25] FALSE TRUE TRUE TRUE TRUE FALSE FALSE TRUE TRUE FALSE FALSE FALSE
-#> [37] TRUE FALSE FALSE TRUE FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE
-#> [49] FALSE FALSE TRUE TRUE FALSE FALSE TRUE FALSE TRUE FALSE FALSE TRUE
-#> [61] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE FALSE FALSE FALSE
-#> [73] FALSE TRUE FALSE FALSE TRUE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
-#> [85] FALSE TRUE TRUE FALSE TRUE FALSE FALSE FALSE FALSE TRUE FALSE FALSE
-#> [97] FALSE FALSE FALSE FALSE
+#> [1] TRUE TRUE TRUE TRUE TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE
+#> [13] TRUE FALSE TRUE TRUE TRUE TRUE TRUE FALSE TRUE FALSE TRUE TRUE
+#> [25] TRUE FALSE FALSE FALSE TRUE TRUE TRUE FALSE TRUE TRUE TRUE TRUE
+#> [37] TRUE FALSE TRUE FALSE FALSE TRUE TRUE FALSE TRUE FALSE FALSE FALSE
+#> [49] FALSE TRUE TRUE TRUE FALSE FALSE FALSE FALSE FALSE FALSE TRUE FALSE
+#> [61] FALSE TRUE TRUE FALSE TRUE FALSE FALSE FALSE FALSE TRUE FALSE FALSE
+#> [73] FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE FALSE FALSE TRUE TRUE
+#> [85] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE FALSE TRUE
+#> [97] TRUE TRUE FALSE FALSE
# filter on results from the third 60-day episode only, using base R
df [ which ( get_episode ( df $ date , 60 ) == 3 ) , ]
-#> # A tibble: 1 × 46
+#> # A tibble: 2 × 46
#> date patient age gender ward mo PEN OXA FLC AMX
#> <date> <chr> <dbl> <chr> <chr> <mo> <sir> <sir> <sir> <sir>
-#> 1 2003-01-14 553529 89 F Clinical B_ESCHR_COLI R NA NA NA
+#> 1 2002-10-14 FCC668 54 F ICU B_ACNTB R NA NA NA
+#> 2 2002-11-18 956065 89 F Clinical B_ESCHR_COLI R NA NA NA
#> # ℹ 36 more variables: AMC <sir>, AMP <sir>, TZP <sir>, CZO <sir>, FEP <sir>,
#> # CXM <sir>, FOX <sir>, CTX <sir>, CAZ <sir>, CRO <sir>, GEN <sir>,
#> # TOB <sir>, AMK <sir>, KAN <sir>, TMP <sir>, SXT <sir>, NIT <sir>,
@@ -317,19 +318,19 @@
arrange ( patient , condition , date )
}
#> # A tibble: 100 × 4
-#> # Groups: patient, condition [97]
+#> # Groups: patient, condition [98]
#> patient date condition new_episode
#> <chr> <date> <chr> <lgl>
-#> 1 04C169 2012-07-18 A TRUE
-#> 2 051150 2007-01-02 C TRUE
-#> 3 067927 2002-02-05 A TRUE
-#> 4 097186 2015-10-28 C TRUE
-#> 5 0DBB93 2009-05-08 B TRUE
-#> 6 0DBF93 2015-10-12 C TRUE
-#> 7 0E2483 2007-11-10 C TRUE
-#> 8 122506 2007-08-10 B TRUE
-#> 9 126334 2009-11-26 A TRUE
-#> 10 187841 2012-09-30 B TRUE
+#> 1 077552 2002-05-14 B TRUE
+#> 2 080086 2010-08-08 A TRUE
+#> 3 0DBB93 2009-05-08 C TRUE
+#> 4 0DBB93 2009-05-08 C FALSE
+#> 5 101305 2006-12-13 A TRUE
+#> 6 114570 2003-04-08 A TRUE
+#> 7 124128 2006-10-02 A TRUE
+#> 8 136315 2004-02-02 B TRUE
+#> 9 138869 2009-05-20 A TRUE
+#> 10 1435C8 2004-03-03 B TRUE
#> # ℹ 90 more rows
if ( require ( "dplyr" ) ) {
@@ -343,19 +344,19 @@
arrange ( patient , ward , date )
}
#> # A tibble: 100 × 5
-#> # Groups: ward, patient [95]
+#> # Groups: ward, patient [97]
#> ward date patient new_index new_logical
#> <chr> <date> <chr> <int> <lgl>
-#> 1 Clinical 2012-07-18 04C169 1 TRUE
-#> 2 ICU 2007-01-02 051150 1 TRUE
-#> 3 ICU 2002-02-05 067927 1 TRUE
-#> 4 Clinical 2015-10-28 097186 1 TRUE
-#> 5 Clinical 2009-05-08 0DBB93 1 TRUE
-#> 6 Clinical 2015-10-12 0DBF93 1 TRUE
-#> 7 Clinical 2007-11-10 0E2483 1 TRUE
-#> 8 Clinical 2007-08-10 122506 1 TRUE
-#> 9 Outpatient 2009-11-26 126334 1 TRUE
-#> 10 Clinical 2012-09-30 187841 1 TRUE
+#> 1 Clinical 2002-05-14 077552 1 TRUE
+#> 2 Clinical 2010-08-08 080086 1 TRUE
+#> 3 Clinical 2009-05-08 0DBB93 1 TRUE
+#> 4 ICU 2009-05-08 0DBB93 1 TRUE
+#> 5 Clinical 2006-12-13 101305 1 TRUE
+#> 6 ICU 2003-04-08 114570 1 TRUE
+#> 7 Clinical 2006-10-02 124128 1 TRUE
+#> 8 Clinical 2004-02-02 136315 1 TRUE
+#> 9 Outpatient 2009-05-20 138869 1 TRUE
+#> 10 Clinical 2004-03-03 1435C8 1 TRUE
#> # ℹ 90 more rows
if ( require ( "dplyr" ) ) {
@@ -371,9 +372,9 @@
#> # A tibble: 3 × 5
#> ward n_patients n_episodes_365 n_episodes_60 n_episodes_30
#> <chr> <int> <int> <int> <int>
-#> 1 Clinical 60 13 34 45
-#> 2 ICU 26 9 19 22
-#> 3 Outpatient 9 6 8 8
+#> 1 Clinical 57 13 36 46
+#> 2 ICU 33 8 20 22
+#> 3 Outpatient 7 6 7 7
# grouping on patients and microorganisms leads to the same
# results as first_isolate() when using 'episode-based':
@@ -402,19 +403,19 @@
select ( group_vars ( . ) , flag_episode )
}
#> # A tibble: 100 × 4
-#> # Groups: patient, mo, ward [97]
-#> patient mo ward flag_episode
-#> <chr> <mo> <chr> <lgl>
-#> 1 474360 B_STPHY_HMNS ICU TRUE
-#> 2 468282 B_STPHY_AURS ICU TRUE
-#> 3 AB0003 B_STRPT_ANGN Clinical TRUE
-#> 4 A95245 B_STPHY_HMNS Outpatient TRUE
-#> 5 219024 B_ESCHR_COLI Clinical TRUE
-#> 6 192353 B_STRPT_PNMN Clinical TRUE
-#> 7 8B8705 B_ESCHR_COLI Clinical TRUE
-#> 8 BB8157 B_STPHY_CONS Clinical TRUE
-#> 9 E44854 B_STRPT_PNMN ICU TRUE
-#> 10 400169 B_SERRT_MRCS ICU TRUE
+#> # Groups: patient, mo, ward [98]
+#> patient mo ward flag_episode
+#> <chr> <mo> <chr> <lgl>
+#> 1 1435C8 B_ESCHR_COLI Clinical TRUE
+#> 2 493589 B_STRPT_PNMN Clinical TRUE
+#> 3 B93C54 B_KLBSL_PNMN Clinical TRUE
+#> 4 889500 B_ESCHR_COLI Clinical TRUE
+#> 5 501429 B_HMPHL_PRNF ICU TRUE
+#> 6 E72482 B_STPHY_CONS ICU TRUE
+#> 7 917895 B_STPHY_CPTS ICU TRUE
+#> 8 534091 B_ESCHR_COLI Clinical TRUE
+#> 9 F79032 B_ESCHR_COLI Clinical TRUE
+#> 10 527306 B_STPHY_EPDR ICU TRUE
#> # ℹ 90 more rows
# }
diff --git a/reference/ggplot_pca.html b/reference/ggplot_pca.html
index f3c408a9..be37327a 100644
--- a/reference/ggplot_pca.html
+++ b/reference/ggplot_pca.html
@@ -1,5 +1,5 @@
-PCA Biplot with ggplot2 — ggplot_pca • AMR (for R) PCA Biplot with ggplot2 — ggplot_pca • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/ggplot_sir.html b/reference/ggplot_sir.html
index 0f9f13f9..14200469 100644
--- a/reference/ggplot_sir.html
+++ b/reference/ggplot_sir.html
@@ -1,5 +1,5 @@
-AMR Plots with ggplot2 — ggplot_sir • AMR (for R) AMR Plots with ggplot2 — ggplot_sir • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/guess_ab_col.html b/reference/guess_ab_col.html
index 50b4d529..9ecaae5a 100644
--- a/reference/guess_ab_col.html
+++ b/reference/guess_ab_col.html
@@ -1,5 +1,5 @@
-Guess Antibiotic Column — guess_ab_col • AMR (for R) Guess Antibiotic Column — guess_ab_col • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/index.html b/reference/index.html
index 9b8a7952..f2e9af0f 100644
--- a/reference/index.html
+++ b/reference/index.html
@@ -1,5 +1,5 @@
-Function reference • AMR (for R) Function reference • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -324,7 +324,7 @@
Determine Multidrug-Resistant Organisms (MDRO)
- plot(<mic> )
autoplot(<mic> )
fortify(<mic> )
plot(<disk> )
autoplot(<disk> )
fortify(<disk> )
plot(<sir> )
autoplot(<sir> )
fortify(<sir> )
+ plot(<mic> )
autoplot(<mic> )
fortify(<mic> )
scale_x_mic()
plot(<disk> )
autoplot(<disk> )
fortify(<disk> )
plot(<sir> )
autoplot(<sir> )
fortify(<sir> )
Plotting for Classes sir
, mic
and disk
diff --git a/reference/intrinsic_resistant.html b/reference/intrinsic_resistant.html
index f126d638..6b2d9fa4 100644
--- a/reference/intrinsic_resistant.html
+++ b/reference/intrinsic_resistant.html
@@ -1,5 +1,5 @@
-Data Set with Bacterial Intrinsic Resistance — intrinsic_resistant • AMR (for R) Data Set with Bacterial Intrinsic Resistance — intrinsic_resistant • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/italicise_taxonomy.html b/reference/italicise_taxonomy.html
index 035d4f0d..05d78f6b 100644
--- a/reference/italicise_taxonomy.html
+++ b/reference/italicise_taxonomy.html
@@ -1,5 +1,5 @@
-Italicise Taxonomic Families, Genera, Species, Subspecies — italicise_taxonomy • AMR (for R) Italicise Taxonomic Families, Genera, Species, Subspecies — italicise_taxonomy • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/join.html b/reference/join.html
index 8dec24f7..d1e6fdbe 100644
--- a/reference/join.html
+++ b/reference/join.html
@@ -1,5 +1,5 @@
-Join microorganisms to a Data Set — join • AMR (for R) Join microorganisms to a Data Set — join • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/key_antimicrobials.html b/reference/key_antimicrobials.html
index 8a4d39f5..d01657a7 100644
--- a/reference/key_antimicrobials.html
+++ b/reference/key_antimicrobials.html
@@ -1,5 +1,5 @@
-(Key) Antimicrobials for First Weighted Isolates — key_antimicrobials • AMR (for R) (Key) Antimicrobials for First Weighted Isolates — key_antimicrobials • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
diff --git a/reference/kurtosis.html b/reference/kurtosis.html
index 861693c0..7b5065e1 100644
--- a/reference/kurtosis.html
+++ b/reference/kurtosis.html
@@ -1,5 +1,5 @@
-Kurtosis of the Sample — kurtosis • AMR (for R) Kurtosis of the Sample — kurtosis • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -199,9 +199,9 @@
Examples
kurtosis ( rnorm ( 10000 ) )
-#> [1] 3.065346
+#> [1] 2.973468
kurtosis ( rnorm ( 10000 ) , excess = TRUE )
-#> [1] -0.03762898
+#> [1] -0.03270321
On this page
diff --git a/reference/like.html b/reference/like.html
index ba54ee03..9c681adb 100644
--- a/reference/like.html
+++ b/reference/like.html
@@ -1,5 +1,5 @@
-Vectorised Pattern Matching with Keyboard Shortcut — like • AMR (for R) Vectorised Pattern Matching with Keyboard Shortcut — like • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -241,17 +241,53 @@
# \donttest{
# get isolates whose name start with 'Entero' (case-insensitive)
example_isolates [ which ( mo_name ( ) %like% "^entero" ) , ]
-#> Error: in mo_name(): argument x is missing and no column with info about
-#> microorganisms could be found.
+#> ℹ Using column 'mo' as input for mo_name()
+#> # A tibble: 106 × 46
+#> date patient age gender ward mo PEN OXA FLC AMX
+#> <date> <chr> <dbl> <chr> <chr> <mo> <sir> <sir> <sir> <sir>
+#> 1 2002-02-21 4FC193 69 M Clinic… B_ENTRC_FACM NA NA NA NA
+#> 2 2002-04-08 130252 78 M ICU B_ENTRC_FCLS NA NA NA NA
+#> 3 2002-06-23 798871 82 M Clinic… B_ENTRC_FCLS NA NA NA NA
+#> 4 2002-06-23 798871 82 M Clinic… B_ENTRC_FCLS NA NA NA NA
+#> 5 2003-04-20 6BC362 62 M ICU B_ENTRC NA NA NA NA
+#> 6 2003-04-21 6BC362 62 M ICU B_ENTRC NA NA NA NA
+#> 7 2003-08-13 F35553 52 M ICU B_ENTRBC_CLOC R NA NA R
+#> 8 2003-08-13 F35553 52 M ICU B_ENTRC_FCLS NA NA NA NA
+#> 9 2003-09-05 F35553 52 M ICU B_ENTRC NA NA NA NA
+#> 10 2003-09-05 F35553 52 M ICU B_ENTRBC_CLOC R NA NA R
+#> # ℹ 96 more rows
+#> # ℹ 36 more variables: AMC <sir>, AMP <sir>, TZP <sir>, CZO <sir>, FEP <sir>,
+#> # CXM <sir>, FOX <sir>, CTX <sir>, CAZ <sir>, CRO <sir>, GEN <sir>,
+#> # TOB <sir>, AMK <sir>, KAN <sir>, TMP <sir>, SXT <sir>, NIT <sir>,
+#> # FOS <sir>, LNZ <sir>, CIP <sir>, MFX <sir>, VAN <sir>, TEC <sir>,
+#> # TCY <sir>, TGC <sir>, DOX <sir>, ERY <sir>, CLI <sir>, AZM <sir>,
+#> # IPM <sir>, MEM <sir>, MTR <sir>, CHL <sir>, COL <sir>, MUP <sir>, …
if ( require ( "dplyr" ) ) {
example_isolates %>%
filter ( mo_name ( ) %like% "^ent" )
}
-#> Error in filter(., mo_name() %like% "^ent"): ℹ In argument: `mo_name() %like% "^ent"`.
-#> Caused by error:
-#> ! in mo_name(): argument x is missing and no column with info about
-#> microorganisms could be found.
+#> ℹ Using column 'mo' as input for mo_name()
+#> # A tibble: 106 × 46
+#> date patient age gender ward mo PEN OXA FLC AMX
+#> <date> <chr> <dbl> <chr> <chr> <mo> <sir> <sir> <sir> <sir>
+#> 1 2002-02-21 4FC193 69 M Clinic… B_ENTRC_FACM NA NA NA NA
+#> 2 2002-04-08 130252 78 M ICU B_ENTRC_FCLS NA NA NA NA
+#> 3 2002-06-23 798871 82 M Clinic… B_ENTRC_FCLS NA NA NA NA
+#> 4 2002-06-23 798871 82 M Clinic… B_ENTRC_FCLS NA NA NA NA
+#> 5 2003-04-20 6BC362 62 M ICU B_ENTRC NA NA NA NA
+#> 6 2003-04-21 6BC362 62 M ICU B_ENTRC NA NA NA NA
+#> 7 2003-08-13 F35553 52 M ICU B_ENTRBC_CLOC R NA NA R
+#> 8 2003-08-13 F35553 52 M ICU B_ENTRC_FCLS NA NA NA NA
+#> 9 2003-09-05 F35553 52 M ICU B_ENTRC NA NA NA NA
+#> 10 2003-09-05 F35553 52 M ICU B_ENTRBC_CLOC R NA NA R
+#> # ℹ 96 more rows
+#> # ℹ 36 more variables: AMC <sir>, AMP <sir>, TZP <sir>, CZO <sir>, FEP <sir>,
+#> # CXM <sir>, FOX <sir>, CTX <sir>, CAZ <sir>, CRO <sir>, GEN <sir>,
+#> # TOB <sir>, AMK <sir>, KAN <sir>, TMP <sir>, SXT <sir>, NIT <sir>,
+#> # FOS <sir>, LNZ <sir>, CIP <sir>, MFX <sir>, VAN <sir>, TEC <sir>,
+#> # TCY <sir>, TGC <sir>, DOX <sir>, ERY <sir>, CLI <sir>, AZM <sir>,
+#> # IPM <sir>, MEM <sir>, MTR <sir>, CHL <sir>, COL <sir>, MUP <sir>, …
# }
diff --git a/reference/mdro.html b/reference/mdro.html
index 6d2e17b2..8756b34c 100644
--- a/reference/mdro.html
+++ b/reference/mdro.html
@@ -1,5 +1,5 @@
-Determine Multidrug-Resistant Organisms (MDRO) — mdro • AMR (for R) Determine Multidrug-Resistant Organisms (MDRO) — mdro • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -335,11 +335,12 @@ A microorganism is categorised as "Resistant" when there is a high likelihood of
Examples
out <- mdro ( example_isolates , guideline = "EUCAST" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
str ( out )
-#> Error in eval(expr, envir, enclos): object 'out' not found
+#> Ord.factor w/ 3 levels "Negative"<"Positive, unconfirmed"<..: 1 1 1 1 1 1 1 1 1 1 ...
table ( out )
-#> Error in eval(expr, envir, enclos): object 'out' not found
+#> out
+#> Negative Positive, unconfirmed Positive
+#> 1978 0 6
out <- mdro ( example_isolates ,
guideline = custom_mdro_guideline (
@@ -347,9 +348,10 @@ A microorganism is categorised as "Resistant" when there is a high likelihood of
VAN == "R" ~ "Custom MDRO 2"
)
)
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
table ( out )
-#> Error in eval(expr, envir, enclos): object 'out' not found
+#> out
+#> Negative Custom MDRO 1 Custom MDRO 2
+#> 870 804 326
# \donttest{
if ( require ( "dplyr" ) ) {
@@ -363,7 +365,18 @@ A microorganism is categorised as "Resistant" when there is a high likelihood of
pull ( MDRO ) %>%
table ( )
}
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Warning: in mdro(): NA introduced for isolates where the available percentage of
+#> antimicrobial classes was below 50% (set with pct_required_classes)
+#> Warning: There was 1 warning in `mutate()`.
+#> ℹ In argument: `MDRO = mdro()`.
+#> Caused by warning:
+#> ! in mdro(): NA introduced for isolates where the available percentage of
+#> antimicrobial classes was below 50% (set with pct_required_classes)
+#> .
+#> Negative Multi-drug-resistant (MDR)
+#> 1601 128
+#> Extensively drug-resistant (XDR) Pandrug-resistant (PDR)
+#> 0 0
# }
diff --git a/reference/mean_amr_distance.html b/reference/mean_amr_distance.html
index fc5f1564..729d1e28 100644
--- a/reference/mean_amr_distance.html
+++ b/reference/mean_amr_distance.html
@@ -1,5 +1,5 @@
-Calculate the Mean AMR Distance — mean_amr_distance • AMR (for R) Calculate the Mean AMR Distance — mean_amr_distance • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -216,28 +216,33 @@
Examples
sir <- random_sir ( 10 )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
sir
-#> Error in eval(expr, envir, enclos): object 'sir' not found
+#> Class 'sir'
+#> [1] S R R S I I R S S R
mean_amr_distance ( sir )
-#> Error in eval(expr, envir, enclos): object 'sir' not found
+#> [1] -0.7745967 1.1618950 1.1618950 -0.7745967 -0.7745967 -0.7745967
+#> [7] 1.1618950 -0.7745967 -0.7745967 1.1618950
mic <- random_mic ( 10 )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
mic
-#> Error in eval(expr, envir, enclos): object 'mic' not found
+#> Class 'mic'
+#> [1] 0.002 >=64 0.125 >=64 2 16 0.005 <=0.001 16
+#> [10] 0.25
mean_amr_distance ( mic )
-#> Error in eval(expr, envir, enclos): object 'mic' not found
+#> [1] -1.2773949 1.1399840 -0.3137594 1.1399840 0.3323488 0.8169299
+#> [7] -1.0638678 -1.4389220 0.8169299 -0.1522324
# equal to the Z-score of their log2:
( log2 ( mic ) - mean ( log2 ( mic ) ) ) / sd ( log2 ( mic ) )
-#> Error in eval(expr, envir, enclos): object 'mic' not found
+#> [1] -1.2773949 1.1399840 -0.3137594 1.1399840 0.3323488 0.8169299
+#> [7] -1.0638678 -1.4389220 0.8169299 -0.1522324
disk <- random_disk ( 10 )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
disk
-#> Error in eval(expr, envir, enclos): object 'disk' not found
+#> Class 'disk'
+#> [1] 34 26 14 50 7 45 18 40 19 24
mean_amr_distance ( disk )
-#> Error in eval(expr, envir, enclos): object 'disk' not found
+#> [1] 0.4467074 -0.1205401 -0.9714114 1.5812024 -1.4677529 1.2266727
+#> [7] -0.6877876 0.8721430 -0.6168817 -0.2623520
y <- data.frame (
id = LETTERS [ 1 : 10 ] ,
@@ -246,15 +251,27 @@
gent = random_mic ( 10 , ab = "gent" , mo = "Escherichia coli" ) ,
tobr = random_mic ( 10 , ab = "tobr" , mo = "Escherichia coli" )
)
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
y
-#> Error in eval(expr, envir, enclos): object 'y' not found
+#> id amox cipr gent tobr
+#> 1 A R 30 8 2
+#> 2 B I 24 <=0.5 4
+#> 3 C S 19 2 <=0.5
+#> 4 D I 21 <=0.5 1
+#> 5 E I 22 >=16 1
+#> 6 F I 28 8 2
+#> 7 G I 19 4 <=0.5
+#> 8 H I 29 <=0.5 <=0.5
+#> 9 I I 26 >=16 <=0.5
+#> 10 J R 21 8 1
mean_amr_distance ( y )
-#> Error in eval(expr, envir, enclos): object 'y' not found
+#> ℹ Calculating mean AMR distance based on columns "amox", "cipr", "gent",
+#> and "tobr"
+#> [1] 1.23945841 0.03375462 -0.73801392 -0.62254823 0.04563698 0.52522389
+#> [7] -0.61650763 -0.37448920 0.05108107 0.45640401
y $ amr_distance <- mean_amr_distance ( y , where ( is.mic ) )
-#> Error in eval(expr, envir, enclos): object 'y' not found
+#> Error in .subset(x, j): invalid subscript type 'list'
y [ order ( y $ amr_distance ) , ]
-#> Error in eval(expr, envir, enclos): object 'y' not found
+#> Error in order(y$amr_distance): argument 1 is not a vector
if ( require ( "dplyr" ) ) {
y %>%
@@ -264,7 +281,19 @@
) %>%
arrange ( check_id_C )
}
-#> Error in eval(expr, envir, enclos): object 'y' not found
+#> ℹ Calculating mean AMR distance based on columns "amox", "cipr", "gent",
+#> and "tobr"
+#> id amox cipr gent tobr amr_distance check_id_C
+#> 1 C S 19 2 <=0.5 -0.73801392 0.0000000
+#> 2 D I 21 <=0.5 1 -0.62254823 0.1154657
+#> 3 G I 19 4 <=0.5 -0.61650763 0.1215063
+#> 4 H I 29 <=0.5 <=0.5 -0.37448920 0.3635247
+#> 5 B I 24 <=0.5 4 0.03375462 0.7717685
+#> 6 E I 22 >=16 1 0.04563698 0.7836509
+#> 7 I I 26 >=16 <=0.5 0.05108107 0.7890950
+#> 8 J R 21 8 1 0.45640401 1.1944179
+#> 9 F I 28 8 2 0.52522389 1.2632378
+#> 10 A R 30 8 2 1.23945841 1.9774723
if ( require ( "dplyr" ) ) {
# support for groups
example_isolates %>%
@@ -274,10 +303,25 @@
mutate ( dist = mean_amr_distance ( . ) ) %>%
arrange ( mo , dist )
}
-#> Error in filter(., mo_genus() == "Enterococcus" & mo_species() != ""): ℹ In argument: `mo_genus() == "Enterococcus" & mo_species() != ""`.
-#> Caused by error:
-#> ! in mo_genus(): argument x is missing and no column with info about
-#> microorganisms could be found.
+#> ℹ 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> <sir> <sir> <sir> <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
+#> # ℹ 53 more rows
On this page
diff --git a/reference/microorganisms.codes.html b/reference/microorganisms.codes.html
index 5c1f76a0..e4e70e08 100644
--- a/reference/microorganisms.codes.html
+++ b/reference/microorganisms.codes.html
@@ -1,5 +1,5 @@
-Data Set with 4 957 Common Microorganism Codes — microorganisms.codes • AMR (for R) Data Set with 4 957 Common Microorganism Codes — microorganisms.codes • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -185,19 +185,95 @@
Examples
microorganisms.codes
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> # A tibble: 4,957 × 2
+#> code mo
+#> <chr> <mo>
+#> 1 1011 B_GRAMP
+#> 2 1012 B_GRAMP
+#> 3 1013 B_GRAMN
+#> 4 1014 B_GRAMN
+#> 5 1015 F_YEAST
+#> 6 103 B_ESCHR_COLI
+#> 7 104 B_SLMNL_ENTR_ENTR
+#> 8 1100 B_STRPT
+#> 9 1101 B_STRPT_VIRI
+#> 10 1102 B_STRPT_HAEM
+#> # ℹ 4,947 more rows
# 'ECO' or 'eco' is the WHONET code for E. coli:
microorganisms.codes [ microorganisms.codes $ code == "ECO" , ]
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> # A tibble: 1 × 2
+#> code mo
+#> <chr> <mo>
+#> 1 ECO B_ESCHR_COLI
# and therefore, 'eco' will be understood as E. coli in this package:
mo_info ( "eco" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> $mo
+#> [1] "B_ESCHR_COLI"
+#>
+#> $kingdom
+#> [1] "Bacteria"
+#>
+#> $phylum
+#> [1] "Pseudomonadota"
+#>
+#> $class
+#> [1] "Gammaproteobacteria"
+#>
+#> $order
+#> [1] "Enterobacterales"
+#>
+#> $family
+#> [1] "Enterobacteriaceae"
+#>
+#> $genus
+#> [1] "Escherichia"
+#>
+#> $species
+#> [1] "coli"
+#>
+#> $subspecies
+#> [1] ""
+#>
+#> $status
+#> [1] "accepted"
+#>
+#> $synonyms
+#> NULL
+#>
+#> $gramstain
+#> [1] "Gram-negative"
+#>
+#> $oxygen_tolerance
+#> [1] "facultative anaerobe"
+#>
+#> $url
+#> [1] "https://lpsn.dsmz.de/species/escherichia-coli"
+#>
+#> $ref
+#> [1] "Castellani et al., 1919"
+#>
+#> $snomed
+#> [1] "1095001000112106" "715307006" "737528008" "416989002"
+#> [5] "116397003" "414097009" "414098004" "414099007"
+#> [9] "414100004" "116395006" "735270003" "116396007"
+#> [13] "83285000" "116394005" "112283007" "710886005"
+#> [17] "710887001" "710888006" "710889003" "414132004"
+#> [21] "721892009" "416812001" "416740004" "417216001"
+#> [25] "457541006" "710253004" "416530004" "417189006"
+#> [29] "409800005" "713925008" "444771000124108" "838549008"
+#>
+#> $lpsn
+#> [1] "776057"
+#>
+#> $gbif
+#> [1] "11286021"
+#>
# works for all AMR functions:
mo_is_intrinsic_resistant ( "eco" , ab = "vancomycin" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] TRUE
On this page
diff --git a/reference/microorganisms.groups.html b/reference/microorganisms.groups.html
index 342179c7..68a886fc 100644
--- a/reference/microorganisms.groups.html
+++ b/reference/microorganisms.groups.html
@@ -1,5 +1,5 @@
-Data Set with 521 Microorganisms In Species Groups — microorganisms.groups • AMR (for R) Data Set with 521 Microorganisms In Species Groups — microorganisms.groups • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -187,11 +187,31 @@
Examples
microorganisms.groups
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> # A tibble: 521 × 4
+#> mo_group mo mo_group_name mo_name
+#> <mo> <mo> <chr> <chr>
+#> 1 B_ACNTB_BMNN-C B_ACNTB_BMNN Acinetobacter baumannii complex Acinetobacter ba…
+#> 2 B_ACNTB_BMNN-C B_ACNTB_CLCC Acinetobacter baumannii complex Acinetobacter ca…
+#> 3 B_ACNTB_BMNN-C B_ACNTB_DJKS Acinetobacter baumannii complex Acinetobacter di…
+#> 4 B_ACNTB_BMNN-C B_ACNTB_NSCM Acinetobacter baumannii complex Acinetobacter no…
+#> 5 B_ACNTB_BMNN-C B_ACNTB_PITT Acinetobacter baumannii complex Acinetobacter pi…
+#> 6 B_ACNTB_BMNN-C B_ACNTB_SFRT Acinetobacter baumannii complex Acinetobacter se…
+#> 7 B_BCTRD_FRGL-C B_BCTRD_FRGL Bacteroides fragilis complex Bacteroides frag…
+#> 8 B_BCTRD_FRGL-C B_BCTRD_OVTS Bacteroides fragilis complex Bacteroides ovat…
+#> 9 B_BCTRD_FRGL-C B_BCTRD_THTT Bacteroides fragilis complex Bacteroides thet…
+#> 10 B_BCTRD_FRGL-C B_BCTRD_VLGT Bacteroides fragilis complex Bacteroides vulg…
+#> # ℹ 511 more rows
# these are all species in the Bacteroides fragilis group, as per WHONET:
microorganisms.groups [ microorganisms.groups $ mo_group == "B_BCTRD_FRGL-C" , ]
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> # A tibble: 5 × 4
+#> mo_group mo mo_group_name mo_name
+#> <mo> <mo> <chr> <chr>
+#> 1 B_BCTRD_FRGL-C B_BCTRD_FRGL Bacteroides fragilis complex Bacteroides fragilis
+#> 2 B_BCTRD_FRGL-C B_BCTRD_OVTS Bacteroides fragilis complex Bacteroides ovatus
+#> 3 B_BCTRD_FRGL-C B_BCTRD_THTT Bacteroides fragilis complex Bacteroides thetaiot…
+#> 4 B_BCTRD_FRGL-C B_BCTRD_VLGT Bacteroides fragilis complex Bacteroides vulgatus
+#> 5 B_BCTRD_FRGL-C B_PRBCT_DSTS Bacteroides fragilis complex Parabacteroides dist…
On this page
diff --git a/reference/microorganisms.html b/reference/microorganisms.html
index e2b0a0c3..6fe1ad3a 100644
--- a/reference/microorganisms.html
+++ b/reference/microorganisms.html
@@ -1,5 +1,5 @@
-Data Set with 52 171 Microorganisms — microorganisms • AMR (for R) Data Set with 52 171 Microorganisms — microorganisms • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -255,7 +255,24 @@
Examples
microorganisms
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> # A tibble: 52,171 × 23
+#> mo fullname status kingdom phylum class order family genus
+#> <mo> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr>
+#> 1 B_GRAMN (unknown Gra… accep… Bacter… (unkn… (unk… (unk… " (unk… " (un…
+#> 2 B_GRAMP (unknown Gra… accep… Bacter… (unkn… (unk… (unk… " (unk… " (un…
+#> 3 B_ANAER-NEG (unknown ana… accep… Bacter… (unkn… (unk… (unk… " (unk… " (un…
+#> 4 B_ANAER-POS (unknown ana… accep… Bacter… (unkn… (unk… (unk… " (unk… " (un…
+#> 5 B_ANAER (unknown ana… accep… Bacter… (unkn… (unk… (unk… " (unk… " (un…
+#> 6 F_FUNGUS (unknown fun… accep… Fungi (unkn… (unk… (unk… " (unk… " (un…
+#> 7 UNKNOWN (unknown nam… accep… (unkno… (unkn… (unk… (unk… " (unk… " (un…
+#> 8 F_YEAST (unknown yea… accep… Fungi (unkn… (unk… (unk… " (unk… " (un…
+#> 9 B_[FAM]_ABDTBCTR Abditibacter… accep… Bacter… Abdit… Abdi… Abdi… " Abdi… ""
+#> 10 B_[ORD]_ABDTBCTR Abditibacter… accep… Bacter… Abdit… Abdi… Abdi… "" ""
+#> # ℹ 52,161 more rows
+#> # ℹ 14 more variables: species <chr>, subspecies <chr>, rank <chr>, ref <chr>,
+#> # oxygen_tolerance <chr>, source <chr>, lpsn <chr>, lpsn_parent <chr>,
+#> # lpsn_renamed_to <chr>, gbif <chr>, gbif_parent <chr>,
+#> # gbif_renamed_to <chr>, prevalence <dbl>, snomed <list>
On this page
diff --git a/reference/mo_matching_score.html b/reference/mo_matching_score.html
index 3de7fced..efee71ca 100644
--- a/reference/mo_matching_score.html
+++ b/reference/mo_matching_score.html
@@ -1,5 +1,5 @@
-Calculate the Matching Score for Microorganisms — mo_matching_score • AMR (for R) Calculate the Matching Score for Microorganisms — mo_matching_score • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -225,16 +225,30 @@
#> ℹ Reset 17 previously matched input values.
as.mo ( "E. coli" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'mo'
+#> [1] B_ESCHR_COLI
mo_uncertainties ( )
-#> No uncertainties to show. Only uncertainties of the last call to as.mo()
-#> or any mo_*() function are stored.
+#> Matching scores are based on the resemblance between the input and the full
+#> taxonomic name, and the pathogenicity in humans. See ?mo_matching_score.
+#>
+#> --------------------------------------------------------------------------------
+#> "E. coli" -> Escherichia coli (B_ESCHR_COLI, 0.688)
+#> Also matched: Enterobacter cowanii (0.600), Eubacterium combesii
+#> (0.600), Eggerthia catenaformis (0.591), Eubacterium callanderi
+#> (0.591), Enterocloster citroniae (0.587), Eubacterium cylindroides
+#> (0.583), Enterococcus casseliflavus (0.577), Enterobacter cloacae
+#> cloacae (0.571), Enterobacter cloacae complex (0.571), and Ehrlichia
+#> canis (0.567)
+#>
+#> Only the first 10 other matches of each record are shown. Run
+#> print(mo_uncertainties(), n = ...) to view more entries, or save
+#> mo_uncertainties() to an object.
mo_matching_score (
x = "E. coli" ,
n = c ( "Escherichia coli" , "Entamoeba coli" )
)
-#> Error in .Call(Cchmatch, x, table, as.integer(nomatch[1L])): NULL value passed as symbol address
+#> [1] 0.6875000 0.3809524
On this page
diff --git a/reference/mo_property.html b/reference/mo_property.html
index a8f752c9..957ccf29 100644
--- a/reference/mo_property.html
+++ b/reference/mo_property.html
@@ -1,5 +1,5 @@
-Get Properties of a Microorganism — mo_property • AMR (for R) Get Properties of a Microorganism — mo_property • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -509,145 +509,151 @@
# taxonomic tree -----------------------------------------------------------
mo_kingdom ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Bacteria"
mo_phylum ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Pseudomonadota"
mo_class ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Gammaproteobacteria"
mo_order ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Enterobacterales"
mo_family ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Enterobacteriaceae"
mo_genus ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Klebsiella"
mo_species ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "pneumoniae"
mo_subspecies ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] ""
# full names and short names -----------------------------------------------
mo_name ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Klebsiella pneumoniae"
mo_fullname ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Klebsiella pneumoniae"
mo_shortname ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "K. pneumoniae"
# other properties ---------------------------------------------------------
mo_pathogenicity ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] Pathogenic
+#> Levels: Pathogenic < Potentially pathogenic < Non-pathogenic < Unknown
mo_gramstain ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Gram-negative"
mo_snomed ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [[1]]
+#> [1] "1098101000112102" "446870005" "1098201000112108" "409801009"
+#> [5] "56415008" "714315002" "713926009"
+#>
mo_type ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Bacteria"
mo_rank ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "species"
mo_url ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Klebsiella pneumoniae
+#> "https://lpsn.dsmz.de/species/klebsiella-pneumoniae"
mo_is_yeast ( c ( "Candida" , "Trichophyton" , "Klebsiella" ) )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] TRUE FALSE FALSE
# scientific reference -----------------------------------------------------
mo_ref ( "Klebsiella aerogenes" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Tindall et al., 2017"
mo_authors ( "Klebsiella aerogenes" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Tindall et al."
mo_year ( "Klebsiella aerogenes" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] 2017
mo_lpsn ( "Klebsiella aerogenes" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "777146"
mo_gbif ( "Klebsiella aerogenes" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "9281703"
mo_synonyms ( "Klebsiella aerogenes" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Hormaeche et al., 1960 Bascomb et al., 1971
+#> "Enterobacter aerogenes" "Klebsiella mobilis"
# abbreviations known in the field -----------------------------------------
mo_genus ( "MRSA" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Staphylococcus"
mo_species ( "MRSA" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "aureus"
mo_shortname ( "VISA" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "S. aureus"
mo_gramstain ( "VISA" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Gram-positive"
mo_genus ( "EHEC" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Escherichia"
mo_species ( "EIEC" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "coli"
mo_name ( "UPEC" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Escherichia coli"
# known subspecies ---------------------------------------------------------
mo_fullname ( "K. pneu rh" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Klebsiella pneumoniae rhinoscleromatis"
mo_shortname ( "K. pneu rh" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "K. pneumoniae"
# \donttest{
# Becker classification, see ?as.mo ----------------------------------------
mo_fullname ( "Staph epidermidis" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Staphylococcus epidermidis"
mo_fullname ( "Staph epidermidis" , Becker = TRUE )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Coagulase-negative Staphylococcus (CoNS)"
mo_shortname ( "Staph epidermidis" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "S. epidermidis"
mo_shortname ( "Staph epidermidis" , Becker = TRUE )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "CoNS"
# Lancefield classification, see ?as.mo ------------------------------------
mo_fullname ( "Strep agalactiae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Streptococcus agalactiae"
mo_fullname ( "Strep agalactiae" , Lancefield = TRUE )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Streptococcus Group B"
mo_shortname ( "Strep agalactiae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "S. agalactiae"
mo_shortname ( "Strep agalactiae" , Lancefield = TRUE )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "GBS"
# language support --------------------------------------------------------
mo_gramstain ( "Klebsiella pneumoniae" , language = "de" ) # German
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Gramnegativ"
mo_gramstain ( "Klebsiella pneumoniae" , language = "nl" ) # Dutch
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Gram-negatief"
mo_gramstain ( "Klebsiella pneumoniae" , language = "es" ) # Spanish
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Gram negativo"
mo_gramstain ( "Klebsiella pneumoniae" , language = "el" ) # Greek
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Αρνητικό κατά Gram"
mo_gramstain ( "Klebsiella pneumoniae" , language = "uk" ) # Ukrainian
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Грамнегативні"
# mo_type is equal to mo_kingdom, but mo_kingdom will remain untranslated
mo_kingdom ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Bacteria"
mo_type ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Bacteria"
mo_kingdom ( "Klebsiella pneumoniae" , language = "zh" ) # Chinese, no effect
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Bacteria"
mo_type ( "Klebsiella pneumoniae" , language = "zh" ) # Chinese, translated
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "细菌"
mo_fullname ( "S. pyogenes" , Lancefield = TRUE , language = "de" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Streptococcus Gruppe A"
mo_fullname ( "S. pyogenes" , Lancefield = TRUE , language = "uk" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Streptococcus Група A"
# other --------------------------------------------------------------------
@@ -658,28 +664,145 @@
filter ( mo_is_gram_positive ( ) ) %>%
count ( mo_genus ( ) , sort = TRUE )
}
-#> Error in filter(., mo_is_gram_positive()): ℹ In argument: `mo_is_gram_positive()`.
-#> Caused by error:
-#> ! in mo_is_gram_positive(): argument x is missing and no column with info
-#> about microorganisms could be found.
+#> ℹ Using column 'mo' as input for mo_is_gram_positive()
+#> ℹ Using column 'mo' as input for mo_genus()
+#> # A tibble: 18 × 2
+#> `mo_genus()` n
+#> <chr> <int>
+#> 1 Staphylococcus 840
+#> 2 Streptococcus 275
+#> 3 Enterococcus 83
+#> 4 Corynebacterium 17
+#> 5 Micrococcus 6
+#> 6 Gemella 3
+#> 7 Aerococcus 2
+#> 8 Cutibacterium 1
+#> 9 Dermabacter 1
+#> 10 Fusibacter 1
+#> 11 Globicatella 1
+#> 12 Granulicatella 1
+#> 13 Lactobacillus 1
+#> 14 Leuconostoc 1
+#> 15 Listeria 1
+#> 16 Paenibacillus 1
+#> 17 Rothia 1
+#> 18 Schaalia 1
if ( require ( "dplyr" ) ) {
example_isolates %>%
filter ( mo_is_intrinsic_resistant ( ab = "vanco" ) ) %>%
count ( mo_genus ( ) , sort = TRUE )
}
-#> Error in filter(., mo_is_intrinsic_resistant(ab = "vanco")): ℹ In argument: `mo_is_intrinsic_resistant(ab = "vanco")`.
-#> Caused by error:
-#> ! in mo_is_intrinsic_resistant(): argument x is missing and no column with
-#> info about microorganisms could be found.
+#> ℹ Using column 'mo' as input for mo_is_intrinsic_resistant()
+#> ℹ Using column 'mo' as input for mo_genus()
+#> # A tibble: 20 × 2
+#> `mo_genus()` n
+#> <chr> <int>
+#> 1 Escherichia 467
+#> 2 Klebsiella 77
+#> 3 Proteus 39
+#> 4 Pseudomonas 30
+#> 5 Serratia 25
+#> 6 Enterobacter 23
+#> 7 Citrobacter 11
+#> 8 Haemophilus 8
+#> 9 Acinetobacter 6
+#> 10 Morganella 6
+#> 11 Pantoea 4
+#> 12 Salmonella 3
+#> 13 Neisseria 2
+#> 14 Stenotrophomonas 2
+#> 15 Campylobacter 1
+#> 16 Enterococcus 1
+#> 17 Hafnia 1
+#> 18 Lactobacillus 1
+#> 19 Leuconostoc 1
+#> 20 Pseudescherichia 1
# get a list with the complete taxonomy (from kingdom to subspecies)
mo_taxonomy ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> $kingdom
+#> [1] "Bacteria"
+#>
+#> $phylum
+#> [1] "Pseudomonadota"
+#>
+#> $class
+#> [1] "Gammaproteobacteria"
+#>
+#> $order
+#> [1] "Enterobacterales"
+#>
+#> $family
+#> [1] "Enterobacteriaceae"
+#>
+#> $genus
+#> [1] "Klebsiella"
+#>
+#> $species
+#> [1] "pneumoniae"
+#>
+#> $subspecies
+#> [1] ""
+#>
# get a list with the taxonomy, the authors, Gram-stain,
# SNOMED codes, and URL to the online database
mo_info ( "Klebsiella pneumoniae" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> $mo
+#> [1] "B_KLBSL_PNMN"
+#>
+#> $kingdom
+#> [1] "Bacteria"
+#>
+#> $phylum
+#> [1] "Pseudomonadota"
+#>
+#> $class
+#> [1] "Gammaproteobacteria"
+#>
+#> $order
+#> [1] "Enterobacterales"
+#>
+#> $family
+#> [1] "Enterobacteriaceae"
+#>
+#> $genus
+#> [1] "Klebsiella"
+#>
+#> $species
+#> [1] "pneumoniae"
+#>
+#> $subspecies
+#> [1] ""
+#>
+#> $status
+#> [1] "accepted"
+#>
+#> $synonyms
+#> NULL
+#>
+#> $gramstain
+#> [1] "Gram-negative"
+#>
+#> $oxygen_tolerance
+#> [1] "aerobe"
+#>
+#> $url
+#> [1] "https://lpsn.dsmz.de/species/klebsiella-pneumoniae"
+#>
+#> $ref
+#> [1] "Trevisan, 1887"
+#>
+#> $snomed
+#> [1] "1098101000112102" "446870005" "1098201000112108" "409801009"
+#> [5] "56415008" "714315002" "713926009"
+#>
+#> $lpsn
+#> [1] "777151"
+#>
+#> $gbif
+#> [1] "3221874"
+#>
# }
diff --git a/reference/mo_source.html b/reference/mo_source.html
index 20d710d7..fac9fead 100644
--- a/reference/mo_source.html
+++ b/reference/mo_source.html
@@ -1,6 +1,6 @@
User-Defined Reference Data Set for Microorganisms — mo_source • AMR (for R) User-Defined Reference Data Set for Microorganisms — mo_source • AMR (for R) Principal Component Analysis (for AMR) — pca • AMR (for R) Principal Component Analysis (for AMR) — pca • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -273,9 +273,22 @@
labs ( title = "Title here" )
}
}
-#> Error in group_by(., order = mo_order(mo), genus = mo_genus(mo)): ℹ In argument: `order = mo_order(mo)`.
-#> Caused by error in `.Call()`:
-#> ! NULL value passed as symbol address
+#> Warning: There were 73 warnings in `summarise()`.
+#> The first warning was:
+#> ℹ In argument: `PEN = (function (..., minimum = 30, as_percent = FALSE,
+#> only_all_tested = FALSE) ...`.
+#> ℹ In group 5: `order = "Lactobacillales"`, `genus = "Enterococcus"`.
+#> Caused by warning:
+#> ! Introducing NA: only 14 results available for PEN in group: order =
+#> "Lactobacillales", genus = "Enterococcus" (minimum = 30).
+#> ℹ Run `dplyr::last_dplyr_warnings()` to see the 72 remaining warnings.
+#> ℹ Columns selected for PCA: "AMC", "CAZ", "CTX", "CXM", "GEN", "SXT",
+#> "TMP", and "TOB". Total observations available: 7.
+#> Groups (n=4, named as 'order'):
+#> [1] "Caryophanales" "Enterobacterales" "Lactobacillales" "Pseudomonadales"
+#>
+
+
# }
diff --git a/reference/plot-1.png b/reference/plot-1.png
new file mode 100644
index 00000000..6fd86aa3
Binary files /dev/null and b/reference/plot-1.png differ
diff --git a/reference/plot-10.png b/reference/plot-10.png
new file mode 100644
index 00000000..bf982126
Binary files /dev/null and b/reference/plot-10.png differ
diff --git a/reference/plot-11.png b/reference/plot-11.png
new file mode 100644
index 00000000..e59cc985
Binary files /dev/null and b/reference/plot-11.png differ
diff --git a/reference/plot-12.png b/reference/plot-12.png
new file mode 100644
index 00000000..5c90f2ab
Binary files /dev/null and b/reference/plot-12.png differ
diff --git a/reference/plot-13.png b/reference/plot-13.png
new file mode 100644
index 00000000..9c553b7d
Binary files /dev/null and b/reference/plot-13.png differ
diff --git a/reference/plot-2.png b/reference/plot-2.png
new file mode 100644
index 00000000..7bfd00ed
Binary files /dev/null and b/reference/plot-2.png differ
diff --git a/reference/plot-3.png b/reference/plot-3.png
new file mode 100644
index 00000000..6ddb0ca5
Binary files /dev/null and b/reference/plot-3.png differ
diff --git a/reference/plot-4.png b/reference/plot-4.png
new file mode 100644
index 00000000..2267ea9e
Binary files /dev/null and b/reference/plot-4.png differ
diff --git a/reference/plot-5.png b/reference/plot-5.png
new file mode 100644
index 00000000..50375370
Binary files /dev/null and b/reference/plot-5.png differ
diff --git a/reference/plot-6.png b/reference/plot-6.png
new file mode 100644
index 00000000..2e10087b
Binary files /dev/null and b/reference/plot-6.png differ
diff --git a/reference/plot-7.png b/reference/plot-7.png
new file mode 100644
index 00000000..6dab7052
Binary files /dev/null and b/reference/plot-7.png differ
diff --git a/reference/plot-8.png b/reference/plot-8.png
new file mode 100644
index 00000000..d16d7e5a
Binary files /dev/null and b/reference/plot-8.png differ
diff --git a/reference/plot-9.png b/reference/plot-9.png
new file mode 100644
index 00000000..71b5ca18
Binary files /dev/null and b/reference/plot-9.png differ
diff --git a/reference/plot.html b/reference/plot.html
index fd846b6d..cd4fc1ec 100644
--- a/reference/plot.html
+++ b/reference/plot.html
@@ -1,5 +1,7 @@
-Plotting for Classes sir, mic and disk — plot • AMR (for R) Plotting for Classes sir, mic and disk — plot • AMR (for R)
@@ -10,7 +12,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -161,6 +163,7 @@
Functions to plot classes sir
, mic
and disk
, with support for base R and ggplot2
.
+
Especially scale_x_mic()
is a relevant wrapper to plot MIC values for ggplot2
. It allows custom MIC ranges and to plot intermediate log2 levels for missing MIC values.
@@ -202,6 +205,16 @@
# S3 method for mic
fortify ( object , ... )
+
scale_x_mic (
+
keep_operators = "edges" ,
+
mic_range = NULL ,
+
... ,
+
drop = FALSE ,
+
guide = waiver ( ) ,
+
position = "bottom" ,
+
na.translate = TRUE
+
)
+
# S3 method for disk
plot (
x ,
@@ -313,6 +326,18 @@
...
arguments passed on to methods
+
+
keep_operators
+
a character specifying how to handle operators (such as >
and <=
) in the input. Accepts one of three values: "all"
(or TRUE
) to keep all operators, "none"
(or FALSE
) to remove all operators, or "edges"
to keep operators only at both ends of the range.
+
+
+
mic_range
+
a manual range to plot the MIC values, e.g., mic_range = c(0.001, 32)
. Use NA
to set no limit on one side, e.g., mic_range = c(NA, 32)
.
+
+
+
drop, guide, position, na.translate
+
arguments passed on to ggplot2::scale_x_discrete()
+
Value
@@ -333,40 +358,68 @@
Examples
some_mic_values <- random_mic ( size = 100 )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
some_disk_values <- random_disk ( size = 100 , mo = "Escherichia coli" , ab = "cipro" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
some_sir_values <- random_sir ( 50 , prob_SIR = c ( 0.55 , 0.05 , 0.30 ) )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
plot ( some_mic_values )
-#> Error in eval(expr, envir, enclos): object 'some_mic_values' not found
+
plot ( some_disk_values )
-#> Error in eval(expr, envir, enclos): object 'some_disk_values' not found
+
plot ( some_sir_values )
-#> Error in eval(expr, envir, enclos): object 'some_sir_values' not found
+
# when providing the microorganism and antibiotic, colours will show interpretations:
plot ( some_mic_values , mo = "S. aureus" , ab = "ampicillin" )
-#> Error in eval(expr, envir, enclos): object 'some_mic_values' not found
+
plot ( some_disk_values , mo = "Escherichia coli" , ab = "cipro" )
-#> Error in eval(expr, envir, enclos): object 'some_disk_values' not found
+
plot ( some_disk_values , mo = "Escherichia coli" , ab = "cipro" , language = "nl" )
-#> Error in eval(expr, envir, enclos): object 'some_disk_values' not found
+
# \donttest{
if ( require ( "ggplot2" ) ) {
autoplot ( some_mic_values )
}
-#> Error in eval(expr, envir, enclos): object 'some_mic_values' not found
+
if ( require ( "ggplot2" ) ) {
autoplot ( some_disk_values , mo = "Escherichia coli" , ab = "cipro" )
}
-#> Error in eval(expr, envir, enclos): object 'some_disk_values' not found
+
if ( require ( "ggplot2" ) ) {
autoplot ( some_sir_values )
}
-#> Error in eval(expr, envir, enclos): object 'some_sir_values' not found
+
+# }
+
+# Plotting using scale_x_mic()
+# \donttest{
+if ( require ( "ggplot2" ) ) {
+ mic_plot <- ggplot ( data.frame ( mics = as.mic ( c ( 0.125 , "<=4" , 4 , 8 , 32 , ">=32" ) ) ,
+ counts = c ( 1 , 1 , 2 , 2 , 3 , 3 ) ) ,
+ aes ( mics , counts ) ) +
+ geom_col ( )
+ mic_plot +
+ labs ( title = "without scale_x_mic()" )
+}
+
+if ( require ( "ggplot2" ) ) {
+ mic_plot +
+ scale_x_mic ( ) +
+ labs ( title = "with scale_x_mic()" )
+}
+
+if ( require ( "ggplot2" ) ) {
+ mic_plot +
+ scale_x_mic ( keep_operators = "all" ) +
+ labs ( title = "with scale_x_mic() keeping all operators" )
+}
+
+if ( require ( "ggplot2" ) ) {
+ mic_plot +
+ scale_x_mic ( mic_range = c ( 1 , 128 ) ) +
+ labs ( title = "with scale_x_mic() using a manual range" )
+}
+
# }
diff --git a/reference/proportion.html b/reference/proportion.html
index ba3b0177..c2f1f7f7 100644
--- a/reference/proportion.html
+++ b/reference/proportion.html
@@ -1,6 +1,6 @@
Calculate Antimicrobial Resistance — proportion • AMR (for R) Calculate Antimicrobial Resistance — proportion • AMR (for R) Random MIC Values/Disk Zones/SIR Generation — random • AMR (for R) Random MIC Values/Disk Zones/SIR Generation — random • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
-
2.1.1
+
2.1.1.9001
@@ -209,27 +209,43 @@
Examples
random_mic ( 25 )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'mic'
+#> [1] 128 2 0.25 256 8 0.025 2 0.5 0.01 0.001
+#> [11] 0.0625 0.25 1 0.5 8 0.005 0.01 16 64 256
+#> [21] 128 256 2 0.125 2
random_disk ( 25 )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'disk'
+#> [1] 28 40 26 12 43 37 35 14 33 43 16 9 34 33 41 42 46 41 41 40 23 50 44 41 27
random_sir ( 25 )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'sir'
+#> [1] I I R I S R R R R S S S R R S R I S I I S I I I I
# \donttest{
# make the random generation more realistic by setting a bug and/or drug:
random_mic ( 25 , "Klebsiella pneumoniae" ) # range 0.0625-64
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'mic'
+#> [1] 64 0.025 0.025 0.025 1 16 8 32 <=0.002
+#> [10] 64 8 0.005 2 8 0.5 32 64 0.0625
+#> [19] 0.025 0.025 0.125 0.125 1 16 32
random_mic ( 25 , "Klebsiella pneumoniae" , "meropenem" ) # range 0.0625-16
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'mic'
+#> [1] 0.5 2 2 4 16 0.5 16 16 8 8 4 8 8 0.5 16 0.5 1 1 16
+#> [20] 4 1 8 4 4 2
random_mic ( 25 , "Streptococcus pneumoniae" , "meropenem" ) # range 0.0625-4
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'mic'
+#> [1] 0.25 1 2 0.125 0.25 >=4 1 2 0.5 2 0.125 0.5
+#> [13] 0.5 0.125 >=4 2 1 2 2 0.5 0.125 >=4 1 1
+#> [25] >=4
random_disk ( 25 , "Klebsiella pneumoniae" ) # range 8-50
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'disk'
+#> [1] 36 42 24 10 21 20 40 18 45 26 22 35 22 11 43 30 48 8 43 16 36 32 41 34 32
random_disk ( 25 , "Klebsiella pneumoniae" , "ampicillin" ) # range 11-17
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'disk'
+#> [1] 15 12 14 16 12 12 12 13 14 16 12 11 13 12 12 15 12 12 17 17 14 17 15 13 15
random_disk ( 25 , "Streptococcus pneumoniae" , "ampicillin" ) # range 12-27
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> Class 'disk'
+#> [1] 27 17 27 19 22 19 27 19 21 21 19 24 20 20 27 27 22 20 18 24 19 24 17 19 18
# }
diff --git a/reference/resistance_predict-1.png b/reference/resistance_predict-1.png
new file mode 100644
index 00000000..84f1746f
Binary files /dev/null and b/reference/resistance_predict-1.png differ
diff --git a/reference/resistance_predict-2.png b/reference/resistance_predict-2.png
new file mode 100644
index 00000000..5996724e
Binary files /dev/null and b/reference/resistance_predict-2.png differ
diff --git a/reference/resistance_predict-3.png b/reference/resistance_predict-3.png
new file mode 100644
index 00000000..ed0dab97
Binary files /dev/null and b/reference/resistance_predict-3.png differ
diff --git a/reference/resistance_predict-4.png b/reference/resistance_predict-4.png
new file mode 100644
index 00000000..54f622cf
Binary files /dev/null and b/reference/resistance_predict-4.png differ
diff --git a/reference/resistance_predict.html b/reference/resistance_predict.html
index 4f4c2a71..f3bf8ab4 100644
--- a/reference/resistance_predict.html
+++ b/reference/resistance_predict.html
@@ -1,5 +1,5 @@
-Predict Antimicrobial Resistance — resistance_predict • AMR (for R) Predict Antimicrobial Resistance — resistance_predict • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -321,14 +321,13 @@ A microorganism is categorised as "Resistant" when there is a high likelihood of
year_min = 2010 ,
model = "binomial"
)
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
plot ( x )
-#> Error in eval(expr, envir, enclos): object 'x' not found
+
# \donttest{
if ( require ( "ggplot2" ) ) {
ggplot_sir_predict ( x )
}
-#> Error: object 'x' not found
+
# using dplyr:
if ( require ( "dplyr" ) ) {
@@ -342,7 +341,25 @@ A microorganism is categorised as "Resistant" when there is a high likelihood of
mymodel <- attributes ( x ) $ model
summary ( mymodel )
}
-#> Error: NULL value passed as symbol address
+
+#> NULL
+#>
+#> Call:
+#> glm(formula = df_matrix ~ year, family = binomial)
+#>
+#> Coefficients:
+#> Estimate Std. Error z value Pr(>|z|)
+#> (Intercept) 27.05312 72.24652 0.374 0.708
+#> year -0.01286 0.03601 -0.357 0.721
+#>
+#> (Dispersion parameter for binomial family taken to be 1)
+#>
+#> Null deviance: 5.7883 on 11 degrees of freedom
+#> Residual deviance: 5.6608 on 10 degrees of freedom
+#> AIC: 50.745
+#>
+#> Number of Fisher Scoring iterations: 4
+#>
# create nice plots with ggplot2 yourself
if ( require ( "dplyr" ) && require ( "ggplot2" ) ) {
@@ -358,7 +375,7 @@ A microorganism is categorised as "Resistant" when there is a high likelihood of
head ( data )
autoplot ( data )
}
-#> Error: In argument: `mo == as.mo("E. coli")`.
+
# }
diff --git a/reference/skewness.html b/reference/skewness.html
index 9a192193..91c3f0bc 100644
--- a/reference/skewness.html
+++ b/reference/skewness.html
@@ -1,6 +1,6 @@
Skewness of the Sample — skewness • AMR (for R) Skewness of the Sample — skewness • AMR (for R) Translate Strings from the AMR Package — translate • AMR (for R) Translate Strings from the AMR Package — translate • AMR (for R)
@@ -10,7 +10,7 @@
AMR (for R)
- 2.1.1
+ 2.1.1.9001
@@ -210,42 +210,43 @@
Examples
# Current settings (based on system language)
ab_name ( "Ciprofloxacin" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Ciprofloxacin"
mo_name ( "Coagulase-negative Staphylococcus (CoNS)" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Coagulase-negative Staphylococcus (CoNS)"
# setting another language
set_AMR_locale ( "Dutch" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> ℹ Using Dutch (Nederlands) for the AMR package for this session.
ab_name ( "Ciprofloxacin" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Ciprofloxacine"
mo_name ( "Coagulase-negative Staphylococcus (CoNS)" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Coagulase-negatieve Staphylococcus (CNS)"
# setting yet another language
set_AMR_locale ( "German" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> ℹ Using German (Deutsch) for the AMR package for this session.
ab_name ( "Ciprofloxacin" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Ciprofloxacin"
mo_name ( "Coagulase-negative Staphylococcus (CoNS)" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Koagulase-negative Staphylococcus (KNS)"
# set_AMR_locale() understands endonyms, English exonyms, and ISO-639-1:
set_AMR_locale ( "Deutsch" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> ℹ Using German (Deutsch) for the AMR package for this session.
set_AMR_locale ( "German" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> ℹ Using German (Deutsch) for the AMR package for this session.
set_AMR_locale ( "de" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> ℹ Using German (Deutsch) for the AMR package for this session.
ab_name ( "amox/clav" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> ℹ More than one result was found for item 1: amoxicillin and clavulansäure
+#> [1] "Amoxicillin/Clavulansäure"
# reset to system default
reset_AMR_locale ( )
#> ℹ Using the English language (English) for the AMR package for this
#> session.
ab_name ( "amox/clav" )
-#> Error in .Call(Cchin, x, table): NULL value passed as symbol address
+#> [1] "Amoxicillin/clavulanic acid"
On this page
diff --git a/search.json b/search.json
index f69b6cf1..a55a5bca 100644
--- a/search.json
+++ b/search.json
@@ -1 +1 @@
-[{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"introduction","dir":"Articles","previous_headings":"","what":"Introduction","title":"How to conduct AMR data analysis","text":"Conducting AMR data analysis unfortunately requires -depth knowledge different scientific fields, makes hard right. least, requires: Good questions (always start !) reliable data thorough understanding (clinical) epidemiology, understand clinical epidemiological relevance possible bias results thorough understanding (clinical) microbiology/infectious diseases, understand microorganisms causal infections implications pharmaceutical treatment, well understanding intrinsic acquired microbial resistance Experience data analysis microbiological tests results, understand determination limitations MIC values interpretations SIR values Availability biological taxonomy microorganisms probably normalisation factors pharmaceuticals, defined daily doses (DDD) Available (inter-)national guidelines, profound methods apply course, instantly provide knowledge experience. AMR package, aimed providing (1) tools simplify antimicrobial resistance data cleaning, transformation analysis, (2) methods easily incorporate international guidelines (3) scientifically reliable reference data, including requirements mentioned . AMR package enables standardised reproducible AMR data analysis, application evidence-based rules, determination first isolates, translation various codes microorganisms antimicrobial agents, determination (multi-drug) resistant microorganisms, calculation antimicrobial resistance, prevalence future trends.","code":""},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"preparation","dir":"Articles","previous_headings":"","what":"Preparation","title":"How to conduct AMR data analysis","text":"tutorial, create fake demonstration data work . can skip Cleaning data already data ready. start analysis, try make structure data generally look like :","code":""},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"needed-r-packages","dir":"Articles","previous_headings":"Preparation","what":"Needed R packages","title":"How to conduct AMR data analysis","text":"many uses R, need additional packages AMR data analysis. package works closely together tidyverse packages dplyr ggplot2 RStudio. tidyverse tremendously improves way conduct data science - allows natural way writing syntaxes creating beautiful plots R. also use cleaner package, can used cleaning data creating frequency tables. AMR package contains data set example_isolates_unclean, might look data users extracted laboratory systems: AMR data analysis, like microorganism column contain valid, --date taxonomy, antibiotic columns cleaned SIR values well.","code":"library(dplyr) library(ggplot2) library(AMR) # (if not yet installed, install with:) # install.packages(c(\"dplyr\", \"ggplot2\", \"AMR\")) example_isolates_unclean #> # A tibble: 3,000 × 8 #> patient_id hospital date bacteria AMX AMC CIP GEN #> #> 1 J3 A 2012-11-21 E. coli R I S S #> 2 R7 A 2018-04-03 K. pneumoniae R I S S #> 3 P3 A 2014-09-19 E. coli R S S S #> 4 P10 A 2015-12-10 E. coli S I S S #> 5 B7 A 2015-03-02 E. coli S S S S #> 6 W3 A 2018-03-31 S. aureus R S R S #> 7 J8 A 2016-06-14 E. coli R S S S #> 8 M3 A 2015-10-25 E. coli R S S S #> 9 J3 A 2019-06-19 E. coli S S S S #> 10 G6 A 2015-04-27 S. aureus S S S S #> # ℹ 2,990 more rows # we will use 'our_data' as the data set name for this tutorial our_data <- example_isolates_unclean"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"taxonomy-of-microorganisms","dir":"Articles","previous_headings":"Preparation","what":"Taxonomy of microorganisms","title":"How to conduct AMR data analysis","text":".mo(), users can transform arbitrary microorganism names codes current taxonomy. AMR package contains --date taxonomic data. specific, currently included data retrieved 11 Dec 2022. codes AMR packages come .mo() short, still human readable. importantly, .mo() supports kinds input: first character codes denote taxonomic kingdom, Bacteria (B), Fungi (F), Protozoa (P). AMR package also contain functions directly retrieve taxonomic properties, name, genus, species, family, order, even Gram-stain. start mo_ use .mo() internally, still arbitrary user input can used: Now can thus clean data: Apparently, uncertainty translation taxonomic codes. Let’s check : ’s good.","code":"as.mo(\"Klebsiella pneumoniae\") #> Class 'mo' #> [1] B_KLBSL_PNMN as.mo(\"K. pneumoniae\") #> Class 'mo' #> [1] B_KLBSL_PNMN as.mo(\"KLEPNE\") #> Class 'mo' #> [1] B_KLBSL_PNMN as.mo(\"KLPN\") #> Class 'mo' #> [1] B_KLBSL_PNMN mo_family(\"K. pneumoniae\") #> [1] \"Enterobacteriaceae\" mo_genus(\"K. pneumoniae\") #> [1] \"Klebsiella\" mo_species(\"K. pneumoniae\") #> [1] \"pneumoniae\" mo_gramstain(\"Klebsiella pneumoniae\") #> [1] \"Gram-negative\" mo_ref(\"K. pneumoniae\") #> [1] \"Trevisan, 1887\" mo_snomed(\"K. pneumoniae\") #> [[1]] #> [1] \"1098101000112102\" \"446870005\" \"1098201000112108\" \"409801009\" #> [5] \"56415008\" \"714315002\" \"713926009\" our_data$bacteria <- as.mo(our_data$bacteria, info = TRUE) #> ℹ Microorganism translation was uncertain for four microorganisms. Run #> mo_uncertainties() to review these uncertainties, or use #> add_custom_microorganisms() to add custom entries. mo_uncertainties() #> Matching scores are based on the resemblance between the input and the full #> taxonomic name, and the pathogenicity in humans. See ?mo_matching_score. #> #> -------------------------------------------------------------------------------- #> \"E. coli\" -> Escherichia coli (B_ESCHR_COLI, 0.688) #> Also matched: Enterobacter cowanii (0.600), Eubacterium combesii #> (0.600), Eggerthia catenaformis (0.591), Eubacterium callanderi #> (0.591), Enterocloster citroniae (0.587), Eubacterium cylindroides #> (0.583), Enterococcus casseliflavus (0.577), Enterobacter cloacae #> cloacae (0.571), Enterobacter cloacae complex (0.571), and Ehrlichia #> canis (0.567) #> -------------------------------------------------------------------------------- #> \"K. pneumoniae\" -> Klebsiella pneumoniae (B_KLBSL_PNMN, 0.786) #> Also matched: Klebsiella pneumoniae ozaenae (0.707), Klebsiella #> pneumoniae pneumoniae (0.688), Klebsiella pneumoniae rhinoscleromatis #> (0.658), Klebsiella pasteurii (0.500), Klebsiella planticola (0.500), #> Kingella potus (0.400), Kosakonia pseudosacchari (0.361), Kaistella #> palustris (0.333), Kocuria palustris (0.333), and Kocuria pelophila #> (0.333) #> -------------------------------------------------------------------------------- #> \"S. aureus\" -> Staphylococcus aureus (B_STPHY_AURS, 0.690) #> Also matched: Staphylococcus aureus aureus (0.643), Staphylococcus #> argenteus (0.625), Staphylococcus aureus anaerobius (0.625), Salmonella #> Aurelianis (0.595), Salmonella Aarhus (0.588), Salmonella Amounderness #> (0.587), Selenomonas artemidis (0.571), Salmonella choleraesuis #> arizonae (0.562), Streptococcus anginosus anginosus (0.561), and #> Salmonella Abaetetuba (0.548) #> -------------------------------------------------------------------------------- #> \"S. pneumoniae\" -> Streptococcus pneumoniae (B_STRPT_PNMN, 0.750) #> Also matched: Streptococcus pseudopneumoniae (0.700), Serratia #> proteamaculans quinovora (0.545), Streptococcus pseudoporcinus (0.536), #> Staphylococcus pseudintermedius (0.532), Serratia proteamaculans #> proteamaculans (0.526), Salmonella Portanigra (0.524), Sphingomonas #> paucimobilis (0.520), Streptococcus pluranimalium (0.519), #> Streptococcus constellatus pharyngis (0.514), and Salmonella Pakistan #> (0.500) #> #> Only the first 10 other matches of each record are shown. Run #> print(mo_uncertainties(), n = ...) to view more entries, or save #> mo_uncertainties() to an object."},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"antibiotic-results","dir":"Articles","previous_headings":"Preparation","what":"Antibiotic results","title":"How to conduct AMR data analysis","text":"column antibiotic test results must also cleaned. AMR package comes three new data types work test results: mic minimal inhibitory concentrations (MIC), disk disk diffusion diameters, sir SIR data interpreted already. package can also determine SIR values based MIC disk diffusion values, read .sir() page. now, just clean SIR columns data using dplyr: basically cleaning, time start data inclusion.","code":"# method 1, be explicit about the columns: our_data <- our_data %>% mutate_at(vars(AMX:GEN), as.sir) # method 2, let the AMR package determine the eligible columns our_data <- our_data %>% mutate_if(is_sir_eligible, as.sir) # result: our_data #> # A tibble: 3,000 × 8 #> patient_id hospital date bacteria AMX AMC CIP GEN #> #> 1 J3 A 2012-11-21 B_ESCHR_COLI R I S S #> 2 R7 A 2018-04-03 B_KLBSL_PNMN R I S S #> 3 P3 A 2014-09-19 B_ESCHR_COLI R S S S #> 4 P10 A 2015-12-10 B_ESCHR_COLI S I S S #> 5 B7 A 2015-03-02 B_ESCHR_COLI S S S S #> 6 W3 A 2018-03-31 B_STPHY_AURS R S R S #> 7 J8 A 2016-06-14 B_ESCHR_COLI R S S S #> 8 M3 A 2015-10-25 B_ESCHR_COLI R S S S #> 9 J3 A 2019-06-19 B_ESCHR_COLI S S S S #> 10 G6 A 2015-04-27 B_STPHY_AURS S S S S #> # ℹ 2,990 more rows"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"first-isolates","dir":"Articles","previous_headings":"Preparation","what":"First isolates","title":"How to conduct AMR data analysis","text":"need know isolates can actually use analysis without repetition bias. conduct analysis antimicrobial resistance, must include first isolate every patient per episode (Hindler et al., Clin Infect Dis. 2007). , easily get overestimate underestimate resistance antibiotic. Imagine patient admitted MRSA found 5 different blood cultures following weeks (yes, countries like Netherlands blood drawing policies). resistance percentage oxacillin isolates overestimated, included MRSA . clearly selection bias. Clinical Laboratory Standards Institute (CLSI) appoints follows: (…) preparing cumulative antibiogram guide clinical decisions empirical antimicrobial therapy initial infections, first isolate given species per patient, per analysis period (eg, one year) included, irrespective body site, antimicrobial susceptibility profile, phenotypical characteristics (eg, biotype). first isolate easily identified, cumulative antimicrobial susceptibility test data prepared using first isolate generally comparable cumulative antimicrobial susceptibility test data calculated methods, providing duplicate isolates excluded. M39-A4 Analysis Presentation Cumulative Antimicrobial Susceptibility Test Data, 4th Edition. CLSI, 2014. Chapter 6.4 AMR package includes methodology first_isolate() function able apply four different methods defined Hindler et al. 2007: phenotype-based, episode-based, patient-based, isolate-based. right method depends goals analysis, default phenotype-based method case method properly correct duplicate isolates. Read methods first_isolate() page. outcome function can easily added data: 90% suitable resistance analysis! can now filter filter() function, also dplyr package: future use, two syntaxes can shortened: end 2 712 isolates analysis. Now data looks like: Time analysis.","code":"our_data <- our_data %>% mutate(first = first_isolate(info = TRUE)) #> ℹ Determining first isolates using an episode length of 365 days #> ℹ Using column 'bacteria' as input for col_mo. #> ℹ Using column 'date' as input for col_date. #> ℹ Using column 'patient_id' as input for col_patient_id. #> ℹ Basing inclusion on all antimicrobial results, using a points threshold #> of 2 #> => Found 2,712 'phenotype-based' first isolates (90.4% of total where a #> microbial ID was available) our_data_1st <- our_data %>% filter(first == TRUE) our_data_1st <- our_data %>% filter_first_isolate() our_data_1st #> # A tibble: 2,712 × 9 #> patient_id hospital date bacteria AMX AMC CIP GEN first #> #> 1 J3 A 2012-11-21 B_ESCHR_COLI R I S S TRUE #> 2 R7 A 2018-04-03 B_KLBSL_PNMN R I S S TRUE #> 3 P10 A 2015-12-10 B_ESCHR_COLI S I S S TRUE #> 4 B7 A 2015-03-02 B_ESCHR_COLI S S S S TRUE #> 5 W3 A 2018-03-31 B_STPHY_AURS R S R S TRUE #> 6 M3 A 2015-10-25 B_ESCHR_COLI R S S S TRUE #> 7 J3 A 2019-06-19 B_ESCHR_COLI S S S S TRUE #> 8 G6 A 2015-04-27 B_STPHY_AURS S S S S TRUE #> 9 P4 A 2011-06-21 B_ESCHR_COLI S S S S TRUE #> 10 Z1 A 2014-09-05 B_ESCHR_COLI S S S S TRUE #> # ℹ 2,702 more rows"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"analysing-the-data","dir":"Articles","previous_headings":"","what":"Analysing the data","title":"How to conduct AMR data analysis","text":"base R summary() function gives good first impression, comes support new mo sir classes now data set:","code":"summary(our_data_1st) #> patient_id hospital date #> Length:2712 Length:2712 Min. :2011-01-01 #> Class :character Class :character 1st Qu.:2013-05-03 #> Mode :character Mode :character Median :2015-06-16 #> Mean :2015-06-21 #> 3rd Qu.:2017-08-24 #> Max. :2019-12-27 #> bacteria AMX AMC #> Class :mo Class:sir Class:sir #> :0 %R :42.9% (n=1163) %R :36.1% (n=978) #> Unique:4 %SI :57.1% (n=1549) %SI :63.9% (n=1734) #> #1 :B_ESCHR_COLI - %S :41.0% (n=1112) - %S :52.0% (n=1409) #> #2 :B_STPHY_AURS - %I :16.1% (n=437) - %I :12.0% (n=325) #> #3 :B_STRPT_PNMN #> CIP GEN first #> Class:sir Class:sir Mode:logical #> %R :42.0% (n=1138) %R :37.3% (n=1011) TRUE:2712 #> %SI :58.0% (n=1574) %SI :62.7% (n=1701) #> - %S :51.5% (n=1396) - %S :59.6% (n=1616) #> - %I : 6.6% (n=178) - %I : 3.1% (n=85) #> glimpse(our_data_1st) #> Rows: 2,712 #> Columns: 9 #> $ patient_id \"J3\", \"R7\", \"P10\", \"B7\", \"W3\", \"M3\", \"J3\", \"G6\", \"P4\", \"Z1\"… #> $ hospital \"A\", \"A\", \"A\", \"A\", \"A\", \"A\", \"A\", \"A\", \"A\", \"A\", \"A\", \"A\",… #> $ date 2012-11-21, 2018-04-03, 2015-12-10, 2015-03-02, 2018-03-31… #> $ bacteria \"B_ESCHR_COLI\", \"B_KLBSL_PNMN\", \"B_ESCHR_COLI\", \"B_ESCHR_COL… #> $ AMX R, R, S, S, R, R, S, S, S, S, R, S, S, S, R, R, R, R, S, R,… #> $ AMC I, I, I, S, S, S, S, S, S, S, S, S, I, S, S, S, S, R, S, S,… #> $ CIP S, S, S, S, R, S, S, S, S, S, S, S, S, S, S, S, S, S, S, S,… #> $ GEN S, S, S, S, S, S, S, S, S, S, R, S, S, S, S, S, S, S, S, S,… #> $ first TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE,… # number of unique values per column: sapply(our_data_1st, n_distinct) #> patient_id hospital date bacteria AMX AMC CIP #> 260 3 1852 4 3 3 3 #> GEN first #> 3 1"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"availability-of-species","dir":"Articles","previous_headings":"Analysing the data","what":"Availability of species","title":"How to conduct AMR data analysis","text":"just get idea species distributed, create frequency table count() based name microorganisms:","code":"our_data %>% count(mo_name(bacteria), sort = TRUE) #> # A tibble: 4 × 2 #> `mo_name(bacteria)` n #> #> 1 Escherichia coli 1518 #> 2 Staphylococcus aureus 730 #> 3 Streptococcus pneumoniae 426 #> 4 Klebsiella pneumoniae 326 our_data_1st %>% count(mo_name(bacteria), sort = TRUE) #> # A tibble: 4 × 2 #> `mo_name(bacteria)` n #> #> 1 Escherichia coli 1319 #> 2 Staphylococcus aureus 676 #> 3 Streptococcus pneumoniae 400 #> 4 Klebsiella pneumoniae 317"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"select-and-filter-with-antibiotic-selectors","dir":"Articles","previous_headings":"Analysing the data","what":"Select and filter with antibiotic selectors","title":"How to conduct AMR data analysis","text":"Using -called antibiotic class selectors, can select filter columns based antibiotic class antibiotic results :","code":"our_data_1st %>% select(date, aminoglycosides()) #> ℹ For aminoglycosides() using column 'GEN' (gentamicin) #> # A tibble: 2,712 × 2 #> date GEN #> #> 1 2012-11-21 S #> 2 2018-04-03 S #> 3 2015-12-10 S #> 4 2015-03-02 S #> 5 2018-03-31 S #> 6 2015-10-25 S #> 7 2019-06-19 S #> 8 2015-04-27 S #> 9 2011-06-21 S #> 10 2014-09-05 S #> # ℹ 2,702 more rows our_data_1st %>% select(bacteria, betalactams()) #> ℹ For betalactams() using columns 'AMX' (amoxicillin) and 'AMC' #> (amoxicillin/clavulanic acid) #> # A tibble: 2,712 × 3 #> bacteria AMX AMC #> #> 1 B_ESCHR_COLI R I #> 2 B_KLBSL_PNMN R I #> 3 B_ESCHR_COLI S I #> 4 B_ESCHR_COLI S S #> 5 B_STPHY_AURS R S #> 6 B_ESCHR_COLI R S #> 7 B_ESCHR_COLI S S #> 8 B_STPHY_AURS S S #> 9 B_ESCHR_COLI S S #> 10 B_ESCHR_COLI S S #> # ℹ 2,702 more rows our_data_1st %>% select(bacteria, where(is.sir)) #> # A tibble: 2,712 × 5 #> bacteria AMX AMC CIP GEN #> #> 1 B_ESCHR_COLI R I S S #> 2 B_KLBSL_PNMN R I S S #> 3 B_ESCHR_COLI S I S S #> 4 B_ESCHR_COLI S S S S #> 5 B_STPHY_AURS R S R S #> 6 B_ESCHR_COLI R S S S #> 7 B_ESCHR_COLI S S S S #> 8 B_STPHY_AURS S S S S #> 9 B_ESCHR_COLI S S S S #> 10 B_ESCHR_COLI S S S S #> # ℹ 2,702 more rows # filtering using AB selectors is also possible: our_data_1st %>% filter(any(aminoglycosides() == \"R\")) #> ℹ For aminoglycosides() using column 'GEN' (gentamicin) #> # A tibble: 1,011 × 9 #> patient_id hospital date bacteria AMX AMC CIP GEN first #> #> 1 J5 A 2017-12-25 B_STRPT_PNMN R S S R TRUE #> 2 X1 A 2017-07-04 B_STPHY_AURS R S S R TRUE #> 3 B3 A 2016-07-24 B_ESCHR_COLI S S S R TRUE #> 4 V7 A 2012-04-03 B_ESCHR_COLI S S S R TRUE #> 5 C9 A 2017-03-23 B_ESCHR_COLI S S S R TRUE #> 6 R1 A 2018-06-10 B_STPHY_AURS S S S R TRUE #> 7 S2 A 2013-07-19 B_STRPT_PNMN S S S R TRUE #> 8 P5 A 2019-03-09 B_STPHY_AURS S S S R TRUE #> 9 Q8 A 2019-08-10 B_STPHY_AURS S S S R TRUE #> 10 K5 A 2013-03-15 B_STRPT_PNMN S S S R TRUE #> # ℹ 1,001 more rows our_data_1st %>% filter(all(betalactams() == \"R\")) #> ℹ For betalactams() using columns 'AMX' (amoxicillin) and 'AMC' #> (amoxicillin/clavulanic acid) #> # A tibble: 483 × 9 #> patient_id hospital date bacteria AMX AMC CIP GEN first #> #> 1 M7 A 2013-07-22 B_STRPT_PNMN R R S S TRUE #> 2 R10 A 2013-12-20 B_STPHY_AURS R R S S TRUE #> 3 R7 A 2015-10-25 B_STPHY_AURS R R S S TRUE #> 4 R8 A 2019-10-25 B_STPHY_AURS R R S S TRUE #> 5 B6 A 2016-11-20 B_ESCHR_COLI R R R R TRUE #> 6 I7 A 2015-08-19 B_ESCHR_COLI R R S S TRUE #> 7 N3 A 2014-12-29 B_STRPT_PNMN R R R S TRUE #> 8 Q2 A 2019-09-22 B_ESCHR_COLI R R S S TRUE #> 9 X7 A 2011-03-20 B_ESCHR_COLI R R S R TRUE #> 10 V1 A 2018-08-07 B_STPHY_AURS R R S S TRUE #> # ℹ 473 more rows # even works in base R (since R 3.0): our_data_1st[all(betalactams() == \"R\"), ] #> ℹ For betalactams() using columns 'AMX' (amoxicillin) and 'AMC' #> (amoxicillin/clavulanic acid) #> # A tibble: 483 × 9 #> patient_id hospital date bacteria AMX AMC CIP GEN first #> #> 1 M7 A 2013-07-22 B_STRPT_PNMN R R S S TRUE #> 2 R10 A 2013-12-20 B_STPHY_AURS R R S S TRUE #> 3 R7 A 2015-10-25 B_STPHY_AURS R R S S TRUE #> 4 R8 A 2019-10-25 B_STPHY_AURS R R S S TRUE #> 5 B6 A 2016-11-20 B_ESCHR_COLI R R R R TRUE #> 6 I7 A 2015-08-19 B_ESCHR_COLI R R S S TRUE #> 7 N3 A 2014-12-29 B_STRPT_PNMN R R R S TRUE #> 8 Q2 A 2019-09-22 B_ESCHR_COLI R R S S TRUE #> 9 X7 A 2011-03-20 B_ESCHR_COLI R R S R TRUE #> 10 V1 A 2018-08-07 B_STPHY_AURS R R S S TRUE #> # ℹ 473 more rows"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"generate-antibiograms","dir":"Articles","previous_headings":"Analysing the data","what":"Generate antibiograms","title":"How to conduct AMR data analysis","text":"Since AMR v2.0 (March 2023), easy create different types antibiograms, support 20 different languages. four antibiogram types, proposed Klinker et al. (2021, DOI 10.1177/20499361211011373), supported new antibiogram() function: Traditional Antibiogram (TA) e.g, susceptibility Pseudomonas aeruginosa piperacillin/tazobactam (TZP) Combination Antibiogram (CA) e.g, sdditional susceptibility Pseudomonas aeruginosa TZP + tobramycin versus TZP alone Syndromic Antibiogram (SA) e.g, susceptibility Pseudomonas aeruginosa TZP among respiratory specimens (obtained among ICU patients ) Weighted-Incidence Syndromic Combination Antibiogram (WISCA) e.g, susceptibility Pseudomonas aeruginosa TZP among respiratory specimens (obtained among ICU patients ) male patients age >=65 years heart failure section, show use antibiogram() function create antibiogram types. starters, included example_isolates data set looks like:","code":"example_isolates #> # A tibble: 2,000 × 46 #> date patient age gender ward mo PEN OXA FLC AMX #> #> 1 2002-01-02 A77334 65 F Clinical B_ESCHR_COLI R NA NA NA #> 2 2002-01-03 A77334 65 F Clinical B_ESCHR_COLI R NA NA NA #> 3 2002-01-07 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 4 2002-01-07 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 5 2002-01-13 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 6 2002-01-13 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 7 2002-01-14 462729 78 M Clinical B_STPHY_AURS R NA S R #> 8 2002-01-14 462729 78 M Clinical B_STPHY_AURS R NA S R #> 9 2002-01-16 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 10 2002-01-17 858515 79 F ICU B_STPHY_EPDR R NA S NA #> # ℹ 1,990 more rows #> # ℹ 36 more variables: AMC , AMP , TZP , CZO , FEP , #> # CXM , FOX , CTX , CAZ , CRO , GEN , #> # TOB , AMK , KAN , TMP , SXT , NIT , #> # FOS , LNZ , CIP , MFX , VAN , TEC , #> # TCY , TGC , DOX , ERY , CLI , AZM , #> # IPM , MEM , MTR , CHL , COL , MUP , …"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"traditional-antibiogram","dir":"Articles","previous_headings":"Analysing the data > Generate antibiograms","what":"Traditional Antibiogram","title":"How to conduct AMR data analysis","text":"create traditional antibiogram, simply state antibiotics used. antibiotics argument antibiogram() function supports (combination) previously mentioned antibiotic class selectors: Notice antibiogram() function automatically prints right format using Quarto R Markdown (page), even applies italics taxonomic names (using italicise_taxonomy() internally). also uses language OS either English, Chinese, Czech, Danish, Dutch, Finnish, French, German, Greek, Italian, Japanese, Norwegian, Polish, Portuguese, Romanian, Russian, Spanish, Swedish, Turkish, Ukrainian. next example, force language Spanish using language argument:","code":"antibiogram(example_isolates, antibiotics = c(aminoglycosides(), carbapenems())) #> ℹ For aminoglycosides() using columns 'GEN' (gentamicin), 'TOB' #> (tobramycin), 'AMK' (amikacin), and 'KAN' (kanamycin) #> ℹ For carbapenems() using columns 'IPM' (imipenem) and 'MEM' (meropenem) antibiogram(example_isolates, mo_transform = \"gramstain\", antibiotics = aminoglycosides(), ab_transform = \"name\", language = \"es\") #> ℹ For aminoglycosides() using columns 'GEN' (gentamicin), 'TOB' #> (tobramycin), 'AMK' (amikacin), and 'KAN' (kanamycin)"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"combined-antibiogram","dir":"Articles","previous_headings":"Analysing the data > Generate antibiograms","what":"Combined Antibiogram","title":"How to conduct AMR data analysis","text":"create combined antibiogram, use antibiotic codes names plus + character like :","code":"antibiogram(example_isolates, antibiotics = c(\"TZP\", \"TZP+TOB\", \"TZP+GEN\"))"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"syndromic-antibiogram","dir":"Articles","previous_headings":"Analysing the data > Generate antibiograms","what":"Syndromic Antibiogram","title":"How to conduct AMR data analysis","text":"create syndromic antibiogram, syndromic_group argument must used. can column data, e.g. ifelse() calculations based certain columns:","code":"antibiogram(example_isolates, antibiotics = c(aminoglycosides(), carbapenems()), syndromic_group = \"ward\") #> ℹ For aminoglycosides() using columns 'GEN' (gentamicin), 'TOB' #> (tobramycin), 'AMK' (amikacin), and 'KAN' (kanamycin) #> ℹ For carbapenems() using columns 'IPM' (imipenem) and 'MEM' (meropenem)"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"weighted-incidence-syndromic-combination-antibiogram-wisca","dir":"Articles","previous_headings":"Analysing the data > Generate antibiograms","what":"Weighted-Incidence Syndromic Combination Antibiogram (WISCA)","title":"How to conduct AMR data analysis","text":"create WISCA, must state combination therapy antibiotics argument (similar Combination Antibiogram), define syndromic group syndromic_group argument (similar Syndromic Antibiogram) cases predefined based clinical demographic characteristics (e.g., endocarditis 75+ females). next example simplification without clinical characteristics, just gives idea WISCA can created:","code":"wisca <- antibiogram(example_isolates, antibiotics = c(\"AMC\", \"AMC+CIP\", \"TZP\", \"TZP+TOB\"), mo_transform = \"gramstain\", minimum = 10, # this should be >= 30, but now just as example syndromic_group = ifelse(example_isolates$age >= 65 & example_isolates$gender == \"M\", \"WISCA Group 1\", \"WISCA Group 2\")) wisca"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"plotting-antibiograms","dir":"Articles","previous_headings":"Analysing the data > Generate antibiograms","what":"Plotting antibiograms","title":"How to conduct AMR data analysis","text":"Antibiograms can plotted using autoplot() ggplot2 packages, since AMR package provides extension function: calculate antimicrobial resistance sensible way, also correcting results, use resistance() susceptibility() functions.","code":"autoplot(wisca)"},{"path":"https://msberends.github.io/AMR/articles/AMR.html","id":"resistance-percentages","dir":"Articles","previous_headings":"Analysing the data","what":"Resistance percentages","title":"How to conduct AMR data analysis","text":"functions resistance() susceptibility() can used calculate antimicrobial resistance susceptibility. specific analyses, functions proportion_S(), proportion_SI(), proportion_I(), proportion_IR() proportion_R() can used determine proportion specific antimicrobial outcome. functions contain minimum argument, denoting minimum required number test results returning value. functions otherwise return NA. default minimum = 30, following CLSI M39-A4 guideline applying microbial epidemiology. per EUCAST guideline 2019, calculate resistance proportion R (proportion_R(), equal resistance()) susceptibility proportion S (proportion_SI(), equal susceptibility()). functions can used : can used conjunction group_by() summarise(), dplyr package: Author: Dr. Matthijs Berends, 26th Feb 2023","code":"our_data_1st %>% resistance(AMX) #> [1] 0.4288348 our_data_1st %>% group_by(hospital) %>% summarise(amoxicillin = resistance(AMX)) #> # A tibble: 3 × 2 #> hospital amoxicillin #> #> 1 A 0.342 #> 2 B 0.564 #> 3 C 0.372"},{"path":"https://msberends.github.io/AMR/articles/EUCAST.html","id":"introduction","dir":"Articles","previous_headings":"","what":"Introduction","title":"How to apply EUCAST rules","text":"EUCAST rules? European Committee Antimicrobial Susceptibility Testing (EUCAST) states website: EUCAST expert rules tabulated collection expert knowledge intrinsic resistances, exceptional resistance phenotypes interpretive rules may applied antimicrobial susceptibility testing order reduce errors make appropriate recommendations reporting particular resistances. Europe, lot medical microbiological laboratories already apply rules (Brown et al., 2015). package features latest insights intrinsic resistance unusual phenotypes (v3.1, 2016). Moreover, eucast_rules() function use purpose can also apply additional rules, like forcing ampicillin = R isolates amoxicillin/clavulanic acid = R.","code":""},{"path":"https://msberends.github.io/AMR/articles/EUCAST.html","id":"examples","dir":"Articles","previous_headings":"","what":"Examples","title":"How to apply EUCAST rules","text":"rules can used discard impossible bug-drug combinations data. example, Klebsiella produces beta-lactamase prevents ampicillin (amoxicillin) working . words, practically every strain Klebsiella resistant ampicillin. Sometimes, laboratory data can still contain strains ampicillin susceptible ampicillin. antibiogram available identification available, antibiogram re-interpreted based identification (namely, Klebsiella). EUCAST expert rules solve , can applied using eucast_rules(): convenient function mo_is_intrinsic_resistant() uses guideline, allows check one specific microorganisms antibiotics: EUCAST rules can used correction, can also used filling known resistance susceptibility based results antimicrobials drugs. process called interpretive reading, basically form imputation, part eucast_rules() function well:","code":"oops <- data.frame( mo = c( \"Klebsiella\", \"Escherichia\" ), ampicillin = \"S\" ) oops #> mo ampicillin #> 1 Klebsiella S #> 2 Escherichia S eucast_rules(oops, info = FALSE) #> mo ampicillin #> 1 Klebsiella R #> 2 Escherichia S mo_is_intrinsic_resistant( c(\"Klebsiella\", \"Escherichia\"), \"ampicillin\" ) #> [1] TRUE FALSE mo_is_intrinsic_resistant( \"Klebsiella\", c(\"ampicillin\", \"kanamycin\") ) #> [1] TRUE FALSE data <- data.frame( mo = c( \"Staphylococcus aureus\", \"Enterococcus faecalis\", \"Escherichia coli\", \"Klebsiella pneumoniae\", \"Pseudomonas aeruginosa\" ), VAN = \"-\", # Vancomycin AMX = \"-\", # Amoxicillin COL = \"-\", # Colistin CAZ = \"-\", # Ceftazidime CXM = \"-\", # Cefuroxime PEN = \"S\", # Benzylenicillin FOX = \"S\", # Cefoxitin stringsAsFactors = FALSE ) data eucast_rules(data)"},{"path":"https://msberends.github.io/AMR/articles/MDR.html","id":"type-of-input","dir":"Articles","previous_headings":"","what":"Type of input","title":"How to determine multi-drug resistance (MDR)","text":"mdro() function takes data set input, regular data.frame. tries automatically determine right columns info isolates, name species columns results antimicrobial agents. See help page info set right settings data command ?mdro. WHONET data (data), settings automatically set correctly.","code":""},{"path":"https://msberends.github.io/AMR/articles/MDR.html","id":"guidelines","dir":"Articles","previous_headings":"","what":"Guidelines","title":"How to determine multi-drug resistance (MDR)","text":"mdro() function support multiple guidelines. can select guideline guideline parameter. Currently supported guidelines (case-insensitive): guideline = \"CMI2012\" (default) Magiorakos AP, Srinivasan et al. “Multidrug-resistant, extensively drug-resistant pandrug-resistant bacteria: international expert proposal interim standard definitions acquired resistance.” Clinical Microbiology Infection (2012) (link) guideline = \"EUCAST3.2\" (simply guideline = \"EUCAST\") European international guideline - EUCAST Expert Rules Version 3.2 “Intrinsic Resistance Unusual Phenotypes” (link) guideline = \"EUCAST3.1\" European international guideline - EUCAST Expert Rules Version 3.1 “Intrinsic Resistance Exceptional Phenotypes Tables” (link) guideline = \"TB\" international guideline multi-drug resistant tuberculosis - World Health Organization “Companion handbook guidelines programmatic management drug-resistant tuberculosis” (link) guideline = \"MRGN\" German national guideline - Mueller et al. (2015) Antimicrobial Resistance Infection Control 4:7. DOI: 10.1186/s13756-015-0047-6 guideline = \"BRMO\" Dutch national guideline - Rijksinstituut voor Volksgezondheid en Milieu “WIP-richtlijn BRMO (Bijzonder Resistente Micro-Organismen) (ZKH)” (link) Please suggest (country-specific) guidelines letting us know: https://github.com/msberends/AMR/issues/new.","code":""},{"path":"https://msberends.github.io/AMR/articles/MDR.html","id":"custom-guidelines","dir":"Articles","previous_headings":"Guidelines","what":"Custom Guidelines","title":"How to determine multi-drug resistance (MDR)","text":"can also use custom guideline. Custom guidelines can set custom_mdro_guideline() function. great importance custom rules determine MDROs hospital, e.g., rules dependent ward, state contact isolation variables data. familiar case_when() dplyr package, recognise input method set rules. Rules must set using R considers ‘formula notation’: row/isolate matches first rule, value first ~ (case ‘Elderly Type ’) set MDRO value. Otherwise, second rule tried . maximum number rules unlimited. can print rules set console overview. Colours help reading console supports colours. outcome function can used guideline argument mdro() function: rules set (custom object case) exported shared file location using saveRDS() collaborate multiple users. custom rules set imported using readRDS().","code":"custom <- custom_mdro_guideline( CIP == \"R\" & age > 60 ~ \"Elderly Type A\", ERY == \"R\" & age > 60 ~ \"Elderly Type B\" ) custom #> A set of custom MDRO rules: #> 1. If CIP is \"R\" and age is higher than 60 then: Elderly Type A #> 2. If ERY is \"R\" and age is higher than 60 then: Elderly Type B #> 3. Otherwise: Negative #> #> Unmatched rows will return NA. #> Results will be of class 'factor', with ordered levels: Negative < Elderly Type A < Elderly Type B x <- mdro(example_isolates, guideline = custom) table(x) #> x #> Negative Elderly Type A Elderly Type B #> 1070 198 732"},{"path":"https://msberends.github.io/AMR/articles/MDR.html","id":"examples","dir":"Articles","previous_headings":"","what":"Examples","title":"How to determine multi-drug resistance (MDR)","text":"mdro() function always returns ordered factor predefined guidelines. example, output default guideline Magiorakos et al. returns factor levels ‘Negative’, ‘MDR’, ‘XDR’ ‘PDR’ order. next example uses example_isolates data set. data set included package contains full antibiograms 2,000 microbial isolates. reflects reality can used practise AMR data analysis. test MDR/XDR/PDR guideline data set, get: Frequency table Class: factor > ordered (numeric) Length: 2,000 Levels: 4: Negative < Multi-drug-resistant (MDR) < Extensively drug-resistant … Available: 1,729 (86.45%, NA: 271 = 13.55%) Unique: 2 another example, create data set determine multi-drug resistant TB: column names automatically verified valid drug names codes, worked exactly way: data set now looks like : can now add interpretation MDR-TB data set. can use: shortcut mdr_tb(): Create frequency table results: 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","code":"library(dplyr) # to support pipes: %>% library(cleaner) # to create frequency tables example_isolates %>% mdro() %>% freq() # show frequency table of the result #> Warning: in mdro(): NA introduced for isolates where the available percentage of #> antimicrobial classes was below 50% (set with pct_required_classes) # random_sir() is a helper function to generate # a random vector with values S, I and R my_TB_data <- data.frame( rifampicin = random_sir(5000), isoniazid = random_sir(5000), gatifloxacin = random_sir(5000), ethambutol = random_sir(5000), pyrazinamide = random_sir(5000), moxifloxacin = random_sir(5000), kanamycin = random_sir(5000) ) my_TB_data <- data.frame( RIF = random_sir(5000), INH = random_sir(5000), GAT = random_sir(5000), ETH = random_sir(5000), PZA = random_sir(5000), MFX = random_sir(5000), KAN = random_sir(5000) ) head(my_TB_data) #> rifampicin isoniazid gatifloxacin ethambutol pyrazinamide moxifloxacin #> 1 I I S I I I #> 2 I I I I I I #> 3 I S R R S S #> 4 I I R I R I #> 5 I S I I R I #> 6 S I R R R I #> kanamycin #> 1 I #> 2 I #> 3 S #> 4 R #> 5 S #> 6 S mdro(my_TB_data, guideline = \"TB\") my_TB_data$mdr <- mdr_tb(my_TB_data) #> ℹ No column found as input for col_mo, assuming all rows contain #> Mycobacterium tuberculosis. freq(my_TB_data$mdr)"},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/PCA.html","id":"transforming","dir":"Articles","previous_headings":"","what":"Transforming","title":"How to conduct principal component analysis (PCA) for AMR","text":"PCA, need transform AMR data first. example_isolates data set package looks like: Now transform data set resistance percentages per taxonomic order genus:","code":"library(AMR) library(dplyr) glimpse(example_isolates) #> Rows: 2,000 #> Columns: 46 #> $ date 2002-01-02, 2002-01-03, 2002-01-07, 2002-01-07, 2002-01-13, 2… #> $ patient \"A77334\", \"A77334\", \"067927\", \"067927\", \"067927\", \"067927\", \"4… #> $ age 65, 65, 45, 45, 45, 45, 78, 78, 45, 79, 67, 67, 71, 71, 75, 50… #> $ gender \"F\", \"F\", \"F\", \"F\", \"F\", \"F\", \"M\", \"M\", \"F\", \"F\", \"M\", \"M\", \"M… #> $ ward \"Clinical\", \"Clinical\", \"ICU\", \"ICU\", \"ICU\", \"ICU\", \"Clinical\"… #> $ mo \"B_ESCHR_COLI\", \"B_ESCHR_COLI\", \"B_STPHY_EPDR\", \"B_STPHY_EPDR\",… #> $ PEN R, R, R, R, R, R, R, R, R, R, R, R, R, R, R, R, R, R, R, R, S,… #> $ OXA NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ FLC NA, NA, R, R, R, R, S, S, R, S, S, S, NA, NA, NA, NA, NA, R, R… #> $ AMX NA, NA, NA, NA, NA, NA, R, R, NA, NA, NA, NA, NA, NA, R, NA, N… #> $ AMC I, I, NA, NA, NA, NA, S, S, NA, NA, S, S, I, I, R, I, I, NA, N… #> $ AMP NA, NA, NA, NA, NA, NA, R, R, NA, NA, NA, NA, NA, NA, R, NA, N… #> $ TZP NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ CZO NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, R, NA,… #> $ FEP NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ CXM I, I, R, R, R, R, S, S, R, S, S, S, S, S, NA, S, S, R, R, S, S… #> $ FOX NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, R, NA,… #> $ CTX NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, S, S, NA, S, S… #> $ CAZ NA, NA, R, R, R, R, R, R, R, R, R, R, NA, NA, NA, S, S, R, R, … #> $ CRO NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, S, S, NA, S, S… #> $ GEN NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ TOB NA, NA, NA, NA, NA, NA, S, S, NA, NA, NA, NA, S, S, NA, NA, NA… #> $ AMK NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ KAN NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ TMP R, R, S, S, R, R, R, R, S, S, NA, NA, S, S, S, S, S, R, R, R, … #> $ SXT R, R, S, S, NA, NA, NA, NA, S, S, NA, NA, S, S, S, S, S, NA, N… #> $ NIT NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, R,… #> $ FOS NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ LNZ R, R, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, R, R, R, R, R, N… #> $ CIP NA, NA, NA, NA, NA, NA, NA, NA, S, S, NA, NA, NA, NA, NA, S, S… #> $ MFX NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ VAN R, R, S, S, S, S, S, S, S, S, NA, NA, R, R, R, R, R, S, S, S, … #> $ TEC R, R, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, R, R, R, R, R, N… #> $ TCY R, R, S, S, S, S, S, S, S, I, S, S, NA, NA, I, R, R, S, I, R, … #> $ TGC NA, NA, S, S, S, S, S, S, S, NA, S, S, NA, NA, NA, R, R, S, NA… #> $ DOX NA, NA, S, S, S, S, S, S, S, NA, S, S, NA, NA, NA, R, R, S, NA… #> $ ERY R, R, R, R, R, R, S, S, R, S, S, S, R, R, R, R, R, R, R, R, S,… #> $ CLI R, R, NA, NA, NA, R, NA, NA, NA, NA, NA, NA, R, R, R, R, R, NA… #> $ AZM R, R, R, R, R, R, S, S, R, S, S, S, R, R, R, R, R, R, R, R, S,… #> $ IPM NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, S, S, NA, S, S… #> $ MEM NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ MTR NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ CHL NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ COL NA, NA, R, R, R, R, R, R, R, R, R, R, NA, NA, NA, R, R, R, R, … #> $ MUP NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA… #> $ RIF R, R, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, R, R, R, R, R, N… resistance_data <- example_isolates %>% group_by( order = mo_order(mo), # group on anything, like order genus = mo_genus(mo) ) %>% # and genus as we do here summarise_if(is.sir, resistance) %>% # then get resistance of all drugs select( order, genus, AMC, CXM, CTX, CAZ, GEN, TOB, TMP, SXT ) # and select only relevant columns head(resistance_data) #> # A tibble: 6 × 10 #> # Groups: order [5] #> order genus AMC CXM CTX CAZ GEN TOB TMP SXT #> #> 1 (unknown order) (unknown ge… NA NA NA NA NA NA NA NA #> 2 Actinomycetales Schaalia NA NA NA NA NA NA NA NA #> 3 Bacteroidales Bacteroides NA NA NA NA NA NA NA NA #> 4 Campylobacterales Campylobact… NA NA NA NA NA NA NA NA #> 5 Caryophanales Gemella NA NA NA NA NA NA NA NA #> 6 Caryophanales Listeria NA NA NA NA NA NA NA NA"},{"path":"https://msberends.github.io/AMR/articles/PCA.html","id":"perform-principal-component-analysis","dir":"Articles","previous_headings":"","what":"Perform principal component analysis","title":"How to conduct principal component analysis (PCA) for AMR","text":"new pca() function automatically filter rows contain numeric values selected variables, now need : result can reviewed good old summary() function: Good news. first two components explain total 93.3% variance (see PC1 PC2 values Proportion Variance. can create -called biplot base R biplot() function, see antimicrobial resistance per drug explain difference per microorganism.","code":"pca_result <- pca(resistance_data) #> ℹ Columns selected for PCA: \"AMC\", \"CAZ\", \"CTX\", \"CXM\", \"GEN\", \"SXT\", #> \"TMP\", and \"TOB\". Total observations available: 7. summary(pca_result) #> Groups (n=4, named as 'order'): #> [1] \"Caryophanales\" \"Enterobacterales\" \"Lactobacillales\" \"Pseudomonadales\" #> Importance of components: #> PC1 PC2 PC3 PC4 PC5 PC6 PC7 #> Standard deviation 2.1539 1.6807 0.6138 0.33879 0.20808 0.03140 9.577e-17 #> Proportion of Variance 0.5799 0.3531 0.0471 0.01435 0.00541 0.00012 0.000e+00 #> Cumulative Proportion 0.5799 0.9330 0.9801 0.99446 0.99988 1.00000 1.000e+00 #> Groups (n=4, named as 'order'): #> [1] \"Caryophanales\" \"Enterobacterales\" \"Lactobacillales\" \"Pseudomonadales\""},{"path":"https://msberends.github.io/AMR/articles/PCA.html","id":"plotting-the-results","dir":"Articles","previous_headings":"","what":"Plotting the results","title":"How to conduct principal component analysis (PCA) for AMR","text":"can’t see explanation points. Perhaps works better new ggplot_pca() function, automatically adds right labels even groups: can also print ellipse per group, edit appearance:","code":"biplot(pca_result) ggplot_pca(pca_result) ggplot_pca(pca_result, ellipse = TRUE) + ggplot2::labs(title = \"An AMR/PCA biplot!\")"},{"path":"https://msberends.github.io/AMR/articles/WHONET.html","id":"import-of-data","dir":"Articles","previous_headings":"","what":"Import of data","title":"How to work with WHONET data","text":"tutorial assumes already imported WHONET data e.g. readxl package. RStudio, can done using menu button ‘Import Dataset’ tab ‘Environment’. Choose option ‘Excel’ select exported file. Make sure date fields imported correctly. example syntax look like : package comes example data set WHONET. use analysis.","code":"library(readxl) data <- read_excel(path = \"path/to/your/file.xlsx\")"},{"path":"https://msberends.github.io/AMR/articles/WHONET.html","id":"preparation","dir":"Articles","previous_headings":"","what":"Preparation","title":"How to work with WHONET data","text":"First, load relevant packages yet . use tidyverse analyses. . don’t know yet, suggest read website: https://www.tidyverse.org/. transform variables simplify automate analysis: Microorganisms transformed microorganism codes (called mo) using Catalogue Life reference data set, contains ~70,000 microorganisms taxonomic kingdoms Bacteria, Fungi Protozoa. tranformation .mo(). function also recognises almost WHONET abbreviations microorganisms. Antimicrobial results interpretations clean valid. words, contain values \"S\", \"\" \"R\". exactly .sir() function . errors warnings, values transformed succesfully. also created package dedicated data cleaning checking, called cleaner package. freq() function can used create frequency tables. let’s check data, couple frequency tables: Frequency table Class: character Length: 500 Available: 500 (100%, NA: 0 = 0%) Unique: 38 Shortest: 11 Longest: 40 (omitted 28 entries, n = 57 [11.4%]) Frequency table Class: factor > ordered > sir (numeric) Length: 500 Levels: 3: S < < R Available: 481 (96.2%, NA: 19 = 3.8%) Unique: 3 Drug: Amoxicillin/clavulanic acid (AMC, J01CR02) Drug group: Beta-lactams/penicillins %SI: 78.59%","code":"library(dplyr) # part of tidyverse library(ggplot2) # part of tidyverse library(AMR) # this package library(cleaner) # to create frequency tables # transform variables data <- WHONET %>% # get microbial ID based on given organism mutate(mo = as.mo(Organism)) %>% # transform everything from \"AMP_ND10\" to \"CIP_EE\" to the new `sir` class mutate_at(vars(AMP_ND10:CIP_EE), as.sir) # our newly created `mo` variable, put in the mo_name() function data %>% freq(mo_name(mo), nmax = 10) # our transformed antibiotic columns # amoxicillin/clavulanic acid (J01CR02) as an example data %>% freq(AMC_ND2)"},{"path":"https://msberends.github.io/AMR/articles/WHONET.html","id":"a-first-glimpse-at-results","dir":"Articles","previous_headings":"","what":"A first glimpse at results","title":"How to work with WHONET data","text":"easy ggplot already give lot information, using included ggplot_sir() function:","code":"data %>% group_by(Country) %>% select(Country, AMP_ND2, AMC_ED20, CAZ_ED10, CIP_ED5) %>% ggplot_sir(translate_ab = \"ab\", facet = \"Country\", datalabels = FALSE)"},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"microorganisms-full-microbial-taxonomy","dir":"Articles","previous_headings":"","what":"microorganisms: Full Microbial Taxonomy","title":"Data sets for download / own use","text":"data set 52 171 rows 23 columns, containing following column names:mo, fullname, status, kingdom, phylum, class, order, family, genus, species, subspecies, rank, ref, oxygen_tolerance, source, lpsn, lpsn_parent, lpsn_renamed_to, gbif, gbif_parent, gbif_renamed_to, prevalence, snomed. data set R available microorganisms, load AMR package. last updated 14 July 2023 08:49:06 UTC. Find info structure data set . Direct download links: Download original R Data Structure (RDS) file (1.2 MB) Download tab-separated text file (11.7 MB) Download Microsoft Excel workbook (5.2 MB) Download Apache Feather file (5.5 MB) Download Apache Parquet file (2.6 MB) Download SAS transport (XPT) file (48.6 MB) Download IBM SPSS Statistics data file (17.8 MB) Download Stata DTA file (48.6 MB) NOTE: exported files SAS, SPSS Stata contain first 50 SNOMED codes per record, file size otherwise exceed 100 MB; file size limit GitHub. file structures compression techniques inefficient. Advice? Use R instead. ’s free much better many ways. tab-separated text file Microsoft Excel workbook contain SNOMED codes comma separated values.","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source","dir":"Articles","previous_headings":"microorganisms: Full Microbial Taxonomy","what":"Source","title":"Data sets for download / own use","text":"data set contains full microbial taxonomy five kingdoms List Prokaryotic names Standing Nomenclature (LPSN) Global Biodiversity Information Facility (GBIF): Parte, AC et al. (2020). List Prokaryotic names Standing Nomenclature (LPSN) moves DSMZ. International Journal Systematic Evolutionary Microbiology, 70, 5607-5612; . Accessed https://lpsn.dsmz.de December 11th, 2022. GBIF Secretariat (2022). GBIF Backbone Taxonomy. Checklist dataset . Accessed https://www.gbif.org December 11th, 2022. Public Health Information Network Vocabulary Access Distribution System (PHIN VADS). US Edition SNOMED CT 1 September 2020. Value Set Name ‘Microorganism’, OID 2.16.840.1.114222.4.11.1009 (v12). URL: https://phinvads.cdc.gov","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"example-content","dir":"Articles","previous_headings":"microorganisms: Full Microbial Taxonomy","what":"Example content","title":"Data sets for download / own use","text":"Included (sub)species per taxonomic kingdom: Example rows filtering genus Escherichia:","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"antibiotics-antibiotic-antifungal-drugs","dir":"Articles","previous_headings":"","what":"antibiotics: Antibiotic (+Antifungal) Drugs","title":"Data sets for download / own use","text":"data set 483 rows 14 columns, containing following column names:ab, cid, name, group, atc, atc_group1, atc_group2, abbreviations, synonyms, oral_ddd, oral_units, iv_ddd, iv_units, loinc. data set R available antibiotics, load AMR package. last updated 20 October 2023 12:51:48 UTC. Find info structure data set . Direct download links: Download original R Data Structure (RDS) file (44 kB) Download tab-separated text file (0.1 MB) Download Microsoft Excel workbook (72 kB) Download Apache Feather file (0.1 MB) Download Apache Parquet file (0.1 MB) Download SAS transport (XPT) file (1.5 MB) Download IBM SPSS Statistics data file (0.4 MB) Download Stata DTA file (0.5 MB) tab-separated text file Microsoft Excel workbook, SAS, SPSS Stata files contain ATC codes, common abbreviations, trade names LOINC codes comma separated values.","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-1","dir":"Articles","previous_headings":"antibiotics: Antibiotic (+Antifungal) Drugs","what":"Source","title":"Data sets for download / own use","text":"data set contains EARS-Net ATC codes gathered WHONET, compound IDs PubChem. also contains brand names (synonyms) found PubChem Defined Daily Doses (DDDs) oral parenteral administration. ATC/DDD index Collaborating Centre Drug Statistics Methodology (note: may used commercial purposes, freely available CC website personal use) PubChem US National Library Medicine WHONET software 2019 LOINC (Logical Observation Identifiers Names Codes)","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"antivirals-antiviral-drugs","dir":"Articles","previous_headings":"","what":"antivirals: Antiviral Drugs","title":"Data sets for download / own use","text":"data set 120 rows 11 columns, containing following column names:av, name, atc, cid, atc_group, synonyms, oral_ddd, oral_units, iv_ddd, iv_units, loinc. data set R available antivirals, load AMR package. last updated 20 October 2023 12:51:48 UTC. Find info structure data set . Direct download links: Download original R Data Structure (RDS) file (6 kB) Download tab-separated text file (17 kB) Download Microsoft Excel workbook (16 kB) Download Apache Feather file (16 kB) Download Apache Parquet file (13 kB) Download SAS transport (XPT) file (73 kB) Download IBM SPSS Statistics data file (32 kB) Download Stata DTA file (78 kB) tab-separated text file Microsoft Excel workbook, SAS, SPSS Stata files contain trade names LOINC codes comma separated values.","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-2","dir":"Articles","previous_headings":"antivirals: Antiviral Drugs","what":"Source","title":"Data sets for download / own use","text":"data set contains ATC codes gathered compound IDs PubChem. also contains brand names (synonyms) found PubChem Defined Daily Doses (DDDs) oral parenteral administration. ATC/DDD index Collaborating Centre Drug Statistics Methodology (note: may used commercial purposes, freely available CC website personal use) PubChem US National Library Medicine LOINC (Logical Observation Identifiers Names Codes)","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"clinical_breakpoints-interpretation-from-mic-values-disk-diameters-to-sir","dir":"Articles","previous_headings":"","what":"clinical_breakpoints: Interpretation from MIC values & disk diameters to SIR","title":"Data sets for download / own use","text":"data set 29 747 rows 12 columns, containing following column names:guideline, type, method, site, mo, rank_index, ab, ref_tbl, disk_dose, breakpoint_S, breakpoint_R, uti. data set R available clinical_breakpoints, load AMR package. last updated 20 October 2023 13:07:11 UTC. Find info structure data set . Direct download links: Download original R Data Structure (RDS) file (59 kB) Download tab-separated text file (2.3 MB) Download Microsoft Excel workbook (1.3 MB) Download Apache Feather file (1.2 MB) Download Apache Parquet file (96 kB) Download SAS transport (XPT) file (7.9 MB) Download IBM SPSS Statistics data file (4.6 MB) Download Stata DTA file (7.8 MB)","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-3","dir":"Articles","previous_headings":"clinical_breakpoints: Interpretation from MIC values & disk diameters to SIR","what":"Source","title":"Data sets for download / own use","text":"data set contains interpretation rules MIC values disk diffusion diameters. Included guidelines CLSI (2011-2023) EUCAST (2011-2023). Clinical breakpoints package validated imported WHONET, free desktop Windows application developed supported Collaborating Centre Surveillance Antimicrobial Resistance. can read website. developers WHONET AMR package contact sharing work. highly appreciate development WHONET software. CEO CLSI chairman EUCAST endorsed work public use AMR package (consequently use breakpoints) June 2023, future development distributing clinical breakpoints discussed meeting CLSI, EUCAST, , developers WHONET AMR package. NOTE: AMR package (WHONET software well) contains internal methods apply guidelines, rather complex. example, breakpoints must applied certain species groups (case package available microorganisms.groups data set). important considered using breakpoints use.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"intrinsic_resistant-intrinsic-bacterial-resistance","dir":"Articles","previous_headings":"","what":"intrinsic_resistant: Intrinsic Bacterial Resistance","title":"Data sets for download / own use","text":"data set 134 634 rows 2 columns, containing following column names:mo ab. data set R available intrinsic_resistant, load AMR package. last updated 16 December 2022 15:10:43 UTC. Find info structure data set . Direct download links: Download original R Data Structure (RDS) file (78 kB) Download tab-separated text file (5.1 MB) Download Microsoft Excel workbook (1.3 MB) Download Apache Feather file (1.2 MB) Download Apache Parquet file (0.2 MB) Download SAS transport (XPT) file (9.5 MB) Download IBM SPSS Statistics data file (7.4 MB) Download Stata DTA file (9.5 MB)","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-4","dir":"Articles","previous_headings":"intrinsic_resistant: Intrinsic Bacterial Resistance","what":"Source","title":"Data sets for download / own use","text":"data set contains defined intrinsic resistance EUCAST bug-drug combinations, based ‘EUCAST Expert Rules’ ‘EUCAST Intrinsic Resistance Unusual Phenotypes’ v3.3 (2021).","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"example-content-4","dir":"Articles","previous_headings":"intrinsic_resistant: Intrinsic Bacterial Resistance","what":"Example content","title":"Data sets for download / own use","text":"Example rows filtering Enterobacter cloacae:","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"dosage-dosage-guidelines-from-eucast","dir":"Articles","previous_headings":"","what":"dosage: Dosage Guidelines from EUCAST","title":"Data sets for download / own use","text":"data set 503 rows 9 columns, containing following column names:ab, name, type, dose, dose_times, administration, notes, original_txt, eucast_version. data set R available dosage, load AMR package. last updated 22 June 2023 13:10:59 UTC. Find info structure data set . Direct download links: Download original R Data Structure (RDS) file (3 kB) Download tab-separated text file (43 kB) Download Microsoft Excel workbook (25 kB) Download Apache Feather file (21 kB) Download Apache Parquet file (9 kB) Download SAS transport (XPT) file (0.1 MB) Download IBM SPSS Statistics data file (64 kB) Download Stata DTA file (0.1 MB)","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-5","dir":"Articles","previous_headings":"dosage: Dosage Guidelines from EUCAST","what":"Source","title":"Data sets for download / own use","text":"EUCAST breakpoints used package based dosages data set. Currently included dosages data set meant : (), ‘EUCAST Clinical Breakpoint Tables’ v11.0 (2021), ‘EUCAST Clinical Breakpoint Tables’ v12.0 (2022).","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"example_isolates-example-data-for-practice","dir":"Articles","previous_headings":"","what":"example_isolates: Example Data for Practice","title":"Data sets for download / own use","text":"data set 2 000 rows 46 columns, containing following column names:date, patient, age, gender, ward, mo, PEN, OXA, FLC, AMX, AMC, AMP, TZP, CZO, FEP, CXM, FOX, CTX, CAZ, CRO, GEN, TOB, AMK, KAN, TMP, SXT, NIT, FOS, LNZ, CIP, MFX, VAN, TEC, TCY, TGC, DOX, ERY, CLI, AZM, IPM, MEM, MTR, CHL, COL, MUP, RIF. data set R available example_isolates, load AMR package. last updated 21 January 2023 22:47:20 UTC. Find info structure data set .","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-6","dir":"Articles","previous_headings":"example_isolates: Example Data for Practice","what":"Source","title":"Data sets for download / own use","text":"data set contains randomised fictitious data, reflects reality can used practise AMR data analysis.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"example_isolates_unclean-example-data-for-practice","dir":"Articles","previous_headings":"","what":"example_isolates_unclean: Example Data for Practice","title":"Data sets for download / own use","text":"data set 3 000 rows 8 columns, containing following column names:patient_id, hospital, date, bacteria, AMX, AMC, CIP, GEN. data set R available example_isolates_unclean, load AMR package. last updated 27 August 2022 18:49:37 UTC. Find info structure data set .","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-7","dir":"Articles","previous_headings":"example_isolates_unclean: Example Data for Practice","what":"Source","title":"Data sets for download / own use","text":"data set contains randomised fictitious data, reflects reality can used practise AMR data analysis.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"microorganisms-groups-species-groups-and-microbiological-complexes","dir":"Articles","previous_headings":"","what":"microorganisms.groups: Species Groups and Microbiological Complexes","title":"Data sets for download / own use","text":"data set 521 rows 4 columns, containing following column names:mo_group, mo, mo_group_name, mo_name. data set R available microorganisms.groups, load AMR package. last updated 14 July 2023 08:49:06 UTC. Find info structure data set . Direct download links: Download original R Data Structure (RDS) file (6 kB) Download tab-separated text file (49 kB) Download Microsoft Excel workbook (20 kB) Download Apache Feather file (19 kB) Download Apache Parquet file (13 kB) Download SAS transport (XPT) file (0 kB) Download IBM SPSS Statistics data file (64 kB) Download Stata DTA file (81 kB)","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-8","dir":"Articles","previous_headings":"microorganisms.groups: Species Groups and Microbiological Complexes","what":"Source","title":"Data sets for download / own use","text":"data set contains species groups microbiological complexes, used clinical_breakpoints data set.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"microorganisms-codes-common-laboratory-codes","dir":"Articles","previous_headings":"","what":"microorganisms.codes: Common Laboratory Codes","title":"Data sets for download / own use","text":"data set 4 957 rows 2 columns, containing following column names:code mo. data set R available microorganisms.codes, load AMR package. last updated 8 July 2023 15:30:05 UTC. Find info structure data set . Direct download links: Download original R Data Structure (RDS) file (22 kB) Download tab-separated text file (0.1 MB) Download Microsoft Excel workbook (91 kB) Download Apache Feather file (85 kB) Download Apache Parquet file (57 kB) Download SAS transport (XPT) file (0 kB) Download IBM SPSS Statistics data file (0.1 MB) Download Stata DTA file (0.1 MB)","code":""},{"path":"https://msberends.github.io/AMR/articles/datasets.html","id":"source-9","dir":"Articles","previous_headings":"microorganisms.codes: Common Laboratory Codes","what":"Source","title":"Data sets for download / own use","text":"data set contains commonly used codes microorganisms, laboratory systems WHONET.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/articles/resistance_predict.html","id":"needed-r-packages","dir":"Articles","previous_headings":"","what":"Needed R packages","title":"How to predict antimicrobial resistance","text":"many uses R, need additional packages AMR data analysis. package works closely together tidyverse packages dplyr ggplot2. tidyverse tremendously improves way conduct data science - allows natural way writing syntaxes creating beautiful plots R. AMR package depends packages even extends use functions.","code":"library(dplyr) library(ggplot2) library(AMR) # (if not yet installed, install with:) # install.packages(c(\"tidyverse\", \"AMR\"))"},{"path":"https://msberends.github.io/AMR/articles/resistance_predict.html","id":"prediction-analysis","dir":"Articles","previous_headings":"","what":"Prediction analysis","title":"How to predict antimicrobial resistance","text":"package contains function resistance_predict(), takes input functions AMR data analysis. Based date column, calculates cases per year uses regression model predict antimicrobial resistance. basically easy : function look date column col_date set. running commands, summary regression model printed unless using resistance_predict(..., info = FALSE). text printed summary - actual result (output) function data.frame containing year: number observations, actual observed resistance, estimated resistance standard error estimation: function plot available base R, can extended packages depend output based type input. extended function cope resistance predictions: fastest way plot result. automatically adds right axes, error bars, titles, number available observations type model. also support ggplot2 package custom function ggplot_sir_predict() create appealing plots:","code":"# resistance prediction of piperacillin/tazobactam (TZP): resistance_predict(tbl = example_isolates, col_date = \"date\", col_ab = \"TZP\", model = \"binomial\") # or: example_isolates %>% resistance_predict( col_ab = \"TZP\", model = \"binomial\" ) # to bind it to object 'predict_TZP' for example: predict_TZP <- example_isolates %>% resistance_predict( col_ab = \"TZP\", model = \"binomial\" ) predict_TZP #> # A tibble: 32 × 7 #> year value se_min se_max observations observed estimated #> * #> 1 2002 0.2 NA NA 15 0.2 0.0562 #> 2 2003 0.0625 NA NA 32 0.0625 0.0616 #> 3 2004 0.0854 NA NA 82 0.0854 0.0676 #> 4 2005 0.05 NA NA 60 0.05 0.0741 #> 5 2006 0.0508 NA NA 59 0.0508 0.0812 #> 6 2007 0.121 NA NA 66 0.121 0.0889 #> 7 2008 0.0417 NA NA 72 0.0417 0.0972 #> 8 2009 0.0164 NA NA 61 0.0164 0.106 #> 9 2010 0.0566 NA NA 53 0.0566 0.116 #> 10 2011 0.183 NA NA 93 0.183 0.127 #> # ℹ 22 more rows plot(predict_TZP) ggplot_sir_predict(predict_TZP) # choose for error bars instead of a ribbon ggplot_sir_predict(predict_TZP, ribbon = FALSE)"},{"path":"https://msberends.github.io/AMR/articles/resistance_predict.html","id":"choosing-the-right-model","dir":"Articles","previous_headings":"Prediction analysis","what":"Choosing the right model","title":"How to predict antimicrobial resistance","text":"Resistance easily predicted; look vancomycin resistance Gram-positive bacteria, spread (.e. standard error) enormous: Vancomycin resistance 100% ten years, might remain low. can define model model parameter. model chosen generalised linear regression model using binomial distribution, assuming period zero resistance followed period increasing resistance leading slowly resistance. Valid values : vancomycin resistance Gram-positive bacteria, linear model might appropriate: model also available object, attribute:","code":"example_isolates %>% filter(mo_gramstain(mo, language = NULL) == \"Gram-positive\") %>% resistance_predict(col_ab = \"VAN\", year_min = 2010, info = FALSE, model = \"binomial\") %>% ggplot_sir_predict() example_isolates %>% filter(mo_gramstain(mo, language = NULL) == \"Gram-positive\") %>% resistance_predict(col_ab = \"VAN\", year_min = 2010, info = FALSE, model = \"linear\") %>% ggplot_sir_predict() model <- attributes(predict_TZP)$model summary(model)$family #> #> Family: binomial #> Link function: logit summary(model)$coefficients #> Estimate Std. Error z value Pr(>|z|) #> (Intercept) -200.67944891 46.17315349 -4.346237 1.384932e-05 #> year 0.09883005 0.02295317 4.305725 1.664395e-05"},{"path":"https://msberends.github.io/AMR/authors.html","id":null,"dir":"","previous_headings":"","what":"Authors","title":"Authors and Citation","text":"Matthijs S. Berends. Author, maintainer. Christian F. Luz. Author, contributor. Dennis Souverein. Author, contributor. Erwin E. . Hassing. Author, contributor. Casper J. Albers. Thesis advisor. Peter Dutey-Magni. Contributor. Judith M. Fonville. Contributor. Alex W. Friedrich. Thesis advisor. Corinna Glasner. Thesis advisor. Eric H. L. C. M. Hazenberg. Contributor. Gwen Knight. Contributor. Annick Lenglet. Contributor. Bart C. Meijer. Contributor. Dmytro Mykhailenko. Contributor. Anton Mymrikov. Contributor. Andrew P. Norgan. Contributor. Sofia Ny. Contributor. Jonas Salm. Contributor. Rogier P. Schade. Contributor. Bhanu N. M. Sinha. Thesis advisor. Anthony Underwood. Contributor. Anita Williams. Contributor.","code":""},{"path":"https://msberends.github.io/AMR/authors.html","id":"citation","dir":"","previous_headings":"","what":"Citation","title":"Authors and Citation","text":"Berends MS, Luz CF, Friedrich AW, Sinha BNM, Albers CJ, Glasner C (2022). “AMR: R Package Working Antimicrobial Resistance Data.” Journal Statistical Software, 104(3), 1–31. doi:10.18637/jss.v104.i03.","code":"@Article{, title = {{AMR}: An {R} Package for Working with Antimicrobial Resistance Data}, author = {Matthijs S. Berends and Christian F. Luz and Alexander W. Friedrich and Bhanu N. M. Sinha and Casper J. Albers and Corinna Glasner}, journal = {Journal of Statistical Software}, year = {2022}, volume = {104}, number = {3}, pages = {1--31}, doi = {10.18637/jss.v104.i03}, }"},{"path":"https://msberends.github.io/AMR/index.html","id":"the-amr-package-for-r-","dir":"","previous_headings":"","what":"Antimicrobial Resistance Data Analysis","title":"Antimicrobial Resistance Data Analysis","text":"Generates antibiograms - traditional, combined, syndromic, even WISCA Provides full microbiological taxonomy data antimicrobial drugs Applies recent CLSI EUCAST clinical breakpoints MICs disk zones Corrects duplicate isolates, calculates predicts AMR per antibiotic class Integrates WHONET, ATC, EARS-Net, PubChem, LOINC SNOMED CT Works Windows, macOS Linux versions R since R-3.0 completely dependency-free, highly suitable places limited resources https://msberends.github.io/AMR https://doi.org/10.18637/jss.v104.i03","code":""},{"path":"https://msberends.github.io/AMR/index.html","id":"introduction","dir":"","previous_headings":"","what":"Introduction","title":"Antimicrobial Resistance Data Analysis","text":"AMR package free open-source R package zero dependencies simplify analysis prediction Antimicrobial Resistance (AMR) work microbial antimicrobial data properties, using evidence-based methods. aim provide standard clean reproducible AMR data analysis, can therefore empower epidemiological analyses continuously enable surveillance treatment evaluation setting. Many different researchers around globe continually helping us make successful durable project! work published Journal Statistical Software (Volume 104(3); DOI 10.18637/jss.v104.i03) formed basis two PhD theses (DOI 10.33612/diss.177417131 DOI 10.33612/diss.192486375). installing package, R knows ~52,000 distinct microbial species (updated December 2022) ~600 antibiotic, antimycotic antiviral drugs name code (including ATC, EARS-Net, ASIARS-Net, PubChem, LOINC SNOMED CT), knows valid SIR MIC values. integral clinical breakpoint guidelines CLSI EUCAST included, even epidemiological cut-(ECOFF) values. supports can read data format, including WHONET data. package works Windows, macOS Linux versions R since R-3.0 (April 2013). designed work setting, including limited resources. created routine data analysis academic research Faculty Medical Sciences University Groningen, collaboration non-profit organisations Certe Medical Diagnostics Advice Foundation University Medical Center Groningen.","code":""},{"path":"https://msberends.github.io/AMR/index.html","id":"used-in-over-175-countries-translated-into-20-languages","dir":"","previous_headings":"Introduction","what":"Used in over 175 countries, translated into 20 languages","title":"Antimicrobial Resistance Data Analysis","text":"Since first public release early 2018, R package used almost countries world. Click map enlarge see country names. help contributors corners world, AMR package available English, Czech, Chinese, Danish, Dutch, Finnish, French, German, Greek, Italian, Japanese, Norwegian, Polish, Portuguese, Romanian, Russian, Spanish, Swedish, Turkish, Ukrainian. Antimicrobial drug (group) names colloquial microorganism names provided languages.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/index.html","id":"filtering-and-selecting-data","dir":"","previous_headings":"Practical examples","what":"Filtering and selecting data","title":"Antimicrobial Resistance Data Analysis","text":"One powerful functions package, aside calculating plotting AMR, selecting filtering based antibiotic columns. can done using -called antibiotic class selectors work base R, dplyr data.table: defined row filter Gram-negative bacteria intrinsic resistance cefotaxime (mo_is_gram_negative() mo_is_intrinsic_resistant()) column selection two antibiotic groups (aminoglycosides() carbapenems()), reference data microorganisms antibiotics AMR package make sure get meant: base R equivalent : base R code work version R since April 2013 (R-3.0). Moreover, code works identically data.table package, starting :","code":"# AMR works great with dplyr, but it's not required or neccesary library(AMR) library(dplyr) example_isolates %>% mutate(bacteria = mo_fullname()) %>% # filtering functions for microorganisms: filter(mo_is_gram_negative(), mo_is_intrinsic_resistant(ab = \"cefotax\")) %>% # antibiotic selectors: select(bacteria, aminoglycosides(), carbapenems()) library(AMR) example_isolates$bacteria <- mo_fullname(example_isolates$mo) example_isolates[which(mo_is_gram_negative() & mo_is_intrinsic_resistant(ab = \"cefotax\")), c(\"bacteria\", aminoglycosides(), carbapenems())] example_isolates <- data.table::as.data.table(example_isolates)"},{"path":"https://msberends.github.io/AMR/index.html","id":"generating-antibiograms","dir":"","previous_headings":"Practical examples","what":"Generating antibiograms","title":"Antimicrobial Resistance Data Analysis","text":"AMR package supports generating traditional, combined, syndromic, even weighted-incidence syndromic combination antibiograms (WISCA). used inside R Markdown Quarto, table printed right output format automatically (markdown, LaTeX, HTML, etc.). combination antibiograms, clear combined antibiotics yield higher empiric coverage: Like many functions package, antibiogram() comes support 20 languages often detected automatically based system language:","code":"antibiogram(example_isolates, antibiotics = c(aminoglycosides(), carbapenems())) antibiogram(example_isolates, antibiotics = c(\"TZP\", \"TZP+TOB\", \"TZP+GEN\"), mo_transform = \"gramstain\") antibiogram(example_isolates, antibiotics = c(\"cipro\", \"tobra\", \"genta\"), # any arbitrary name or code will work mo_transform = \"gramstain\", ab_transform = \"name\", language = \"uk\") # Ukrainian"},{"path":"https://msberends.github.io/AMR/index.html","id":"calculating-resistance-per-group","dir":"","previous_headings":"Practical examples","what":"Calculating resistance per group","title":"Antimicrobial Resistance Data Analysis","text":"manual approach, can use resistance susceptibility() function: use antibiotic class selectors select series antibiotic columns:","code":"example_isolates %>% # group by ward: group_by(ward) %>% # calculate AMR using resistance() for gentamicin and tobramycin # and get their 95% confidence intervals using sir_confidence_interval(): summarise(across(c(GEN, TOB), list(total_R = resistance, conf_int = function(x) sir_confidence_interval(x, collapse = \"-\")))) library(AMR) library(dplyr) out <- example_isolates %>% # group by ward: group_by(ward) %>% # calculate AMR using resistance(), over all aminoglycosides and polymyxins: summarise(across(c(aminoglycosides(), polymyxins()), resistance)) out # transform the antibiotic columns to names: out %>% set_ab_names() # transform the antibiotic column to ATC codes: out %>% set_ab_names(property = \"atc\")"},{"path":"https://msberends.github.io/AMR/index.html","id":"what-else-can-you-do-with-this-package","dir":"","previous_headings":"","what":"What else can you do with this package?","title":"Antimicrobial Resistance Data Analysis","text":"package intended comprehensive toolbox integrated AMR data analysis. package can used : Reference taxonomy microorganisms, since package contains microbial (sub)species List Prokaryotic names Standing Nomenclature (LPSN) Global Biodiversity Information Facility (GBIF) (manual) Interpreting raw MIC disk diffusion values, based CLSI EUCAST guideline (manual) Retrieving antimicrobial drug names, doses forms administration clinical health care records (manual) Determining first isolates used AMR data analysis (manual) Calculating antimicrobial resistance (tutorial) Determining multi-drug resistance (MDR) / multi-drug resistant organisms (MDRO) (tutorial) Calculating (empirical) susceptibility mono therapy combination therapies (tutorial) Predicting future antimicrobial resistance using regression models (tutorial) Getting properties microorganism (like Gram stain, species, genus family) (manual) Getting properties antibiotic (like name, code EARS-Net/ATC/LOINC/PubChem, defined daily dose trade name) (manual) Plotting antimicrobial resistance (tutorial) Applying EUCAST expert rules (manual) Getting SNOMED codes microorganism, getting properties microorganism based SNOMED code (manual) Getting LOINC codes antibiotic, getting properties antibiotic based LOINC code (manual) Machine reading EUCAST CLSI guidelines 2011-2021 translate MIC values disk diffusion diameters SIR (link) Principal component analysis AMR (tutorial)","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/index.html","id":"latest-official-version","dir":"","previous_headings":"Get this package","what":"Latest official version","title":"Antimicrobial Resistance Data Analysis","text":"package available official R network (CRAN). Install package R CRAN using command: downloaded installed automatically. RStudio, click menu Tools > Install Packages… type “AMR” press Install. Note: functions website may available latest release. use functions data sets mentioned website, install latest development version.","code":"install.packages(\"AMR\")"},{"path":"https://msberends.github.io/AMR/index.html","id":"latest-development-version","dir":"","previous_headings":"Get this package","what":"Latest development version","title":"Antimicrobial Resistance Data Analysis","text":"Please read Developer Guideline . latest unpublished development version can installed GitHub two ways: Manually, using: Automatically, using rOpenSci R-universe platform, adding R-universe address list repositories (‘repos’): , can install update AMR package like official release (e.g., using install.packages(\"AMR\") RStudio via Tools > Check Package Updates…).","code":"install.packages(\"remotes\") # if you haven't already remotes::install_github(\"msberends/AMR\") options(repos = c(getOption(\"repos\"), msberends = \"https://msberends.r-universe.dev\"))"},{"path":"https://msberends.github.io/AMR/index.html","id":"get-started","dir":"","previous_headings":"","what":"Get started","title":"Antimicrobial Resistance Data Analysis","text":"find conduct AMR data analysis, please continue reading get started click link ‘’ menu.","code":""},{"path":"https://msberends.github.io/AMR/index.html","id":"partners","dir":"","previous_headings":"","what":"Partners","title":"Antimicrobial Resistance Data Analysis","text":"development package part , related , made possible following non-profit organisations initiatives:","code":""},{"path":"https://msberends.github.io/AMR/index.html","id":"copyright","dir":"","previous_headings":"","what":"Copyright","title":"Antimicrobial Resistance Data Analysis","text":"R package free, open-source software licensed GNU General Public License v2.0 (GPL-2). nutshell, means package: May used commercial purposes May used private purposes May used patent purposes May modified, although: Modifications must released license distributing package Changes made code must documented May distributed, although: Source code must made available package distributed copy license copyright notice must included package. Comes LIMITATION liability Comes warranty","code":""},{"path":"https://msberends.github.io/AMR/reference/AMR-deprecated.html","id":null,"dir":"Reference","previous_headings":"","what":"Deprecated Functions — AMR-deprecated","title":"Deprecated Functions — AMR-deprecated","text":"functions -called 'Deprecated'. removed future release. Using functions give warning name function replaced (one).","code":""},{"path":"https://msberends.github.io/AMR/reference/AMR-deprecated.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Deprecated Functions — AMR-deprecated","text":"","code":"NA_rsi_ as.rsi(x, ...) facet_rsi(...) geom_rsi(...) ggplot_rsi(...) ggplot_rsi_predict(...) is.rsi(...) is.rsi.eligible(...) labels_rsi_count(...) n_rsi(...) random_rsi(...) rsi_df(...) rsi_predict(...) scale_rsi_colours(...) theme_rsi(...)"},{"path":"https://msberends.github.io/AMR/reference/AMR-deprecated.html","id":"format","dir":"Reference","previous_headings":"","what":"Format","title":"Deprecated Functions — AMR-deprecated","text":"object class rsi (inherits ordered, factor) length 1.","code":""},{"path":"https://msberends.github.io/AMR/reference/AMR-options.html","id":null,"dir":"Reference","previous_headings":"","what":"Options for the AMR package — AMR-options","title":"Options for the AMR package — AMR-options","text":"overview package-specific options() can set AMR package.","code":""},{"path":"https://msberends.github.io/AMR/reference/AMR-options.html","id":"options","dir":"Reference","previous_headings":"","what":"Options","title":"Options for the AMR package — AMR-options","text":"AMR_custom_ab Allows use custom antimicrobial drugs package. explained add_custom_antimicrobials(). AMR_custom_mo Allows use custom microorganisms package. explained add_custom_microorganisms(). AMR_eucastrules Used setting default types rules eucast_rules() function, must one : \"breakpoints\", \"expert\", \"\", \"custom\", \"\", defaults c(\"breakpoints\", \"expert\"). AMR_guideline Used setting default guideline interpreting MIC values disk diffusion diameters .sir(). Can guideline name (e.g., \"CLSI\") name year (e.g. \"CLSI 2019\"). default latest implemented EUCAST guideline, currently \"EUCAST 2023\". Supported guideline currently EUCAST (2011-2023) CLSI (2011-2023). AMR_ignore_pattern regular expression ignore (.e., make NA) match given .mo() mo_* functions. AMR_include_PKPD logical use .sir(), indicate PK/PD clinical breakpoints must applied last resort - default TRUE. AMR_ecoff logical use .sir(), indicate ECOFF (Epidemiological Cut-) values must used - default FALSE. AMR_include_screening logical use .sir(), indicate clinical breakpoints screening allowed - default FALSE. AMR_keep_synonyms logical use .mo() mo_* functions, indicate old, previously valid taxonomic names must preserved corrected currently accepted names. default FALSE. AMR_cleaning_regex regular expression (case-insensitive) use .mo() mo_* functions, clean user input. default outcome mo_cleaning_regex(), removes texts brackets texts \"species\" \"serovar\". AMR_locale language use AMR package, can one supported language names ISO-639-1 codes: English (en), Chinese (zh), Czech (cs), Danish (da), Dutch (nl), Finnish (fi), French (fr), German (de), Greek (el), Italian (), Japanese (ja), Norwegian (), Polish (pl), Portuguese (pt), Romanian (ro), Russian (ru), Spanish (es), Swedish (sv), Turkish (tr), Ukrainian (uk). default current system language (supported). AMR_mo_source file location manual code list used .mo() mo_* functions. explained set_mo_source().","code":""},{"path":"https://msberends.github.io/AMR/reference/AMR-options.html","id":"saving-settings-between-sessions","dir":"Reference","previous_headings":"","what":"Saving Settings Between Sessions","title":"Options for the AMR package — AMR-options","text":"Settings R saved globally thus lost R exited. can save options .Rprofile file, user-specific file. can edit using: file, can set options : add Portuguese language support antibiotics, allow PK/PD rules interpreting MIC values .sir().","code":"utils::file.edit(\"~/.Rprofile\") options(AMR_locale = \"pt\") options(AMR_include_PKPD = TRUE)"},{"path":"https://msberends.github.io/AMR/reference/AMR-options.html","id":"share-options-within-team","dir":"Reference","previous_headings":"","what":"Share Options Within Team","title":"Options for the AMR package — AMR-options","text":"global approach, e.g. within data team, save options file remote file location, shared network drive. work way: Save plain text file e.g. \"X:/team_folder/R_options.R\" fill preferred settings. user, open .Rprofile file using utils::file.edit(\"~/.Rprofile\") put : Reload R/RStudio check settings getOption(), e.g. getOption(\"AMR_locale\") set value. Now team settings configured one place, can maintained .","code":"source(\"X:/team_folder/R_options.R\")"},{"path":"https://msberends.github.io/AMR/reference/AMR.html","id":null,"dir":"Reference","previous_headings":"","what":"The AMR Package — AMR","title":"The AMR Package — AMR","text":"Welcome AMR package. AMR package free open-source R package zero dependencies simplify analysis prediction Antimicrobial Resistance (AMR) work microbial antimicrobial data properties, using evidence-based methods. aim provide standard clean reproducible AMR data analysis, can therefore empower epidemiological analyses continuously enable surveillance treatment evaluation setting. Many different researchers around globe continually helping us make successful durable project! work published Journal Statistical Software (Volume 104(3); doi:10.18637/jss.v104.i03 ) formed basis two PhD theses (doi:10.33612/diss.177417131 doi:10.33612/diss.192486375 ). installing package, R knows ~52 000 microorganisms (updated December 2022) ~600 antibiotic, antimycotic antiviral drugs name code (including ATC, EARS-Net, ASIARS-Net, PubChem, LOINC SNOMED CT), knows valid SIR MIC values. integral clinical breakpoint guidelines CLSI EUCAST included, even epidemiological cut-(ECOFF) values. supports can read data format, including WHONET data. package works Windows, macOS Linux versions R since R-3.0 (April 2013). designed work setting, including limited resources. created routine data analysis academic research Faculty Medical Sciences public University Groningen, collaboration non-profit organisations Certe Medical Diagnostics Advice Foundation University Medical Center Groningen. AMR package available English, Chinese, Czech, Danish, Dutch, Finnish, French, German, Greek, Italian, Japanese, Norwegian, Polish, Portuguese, Romanian, Russian, Spanish, Swedish, Turkish, Ukrainian. Antimicrobial drug (group) names colloquial microorganism names provided languages.","code":""},{"path":"https://msberends.github.io/AMR/reference/AMR.html","id":"source","dir":"Reference","previous_headings":"","what":"Source","title":"The AMR Package — AMR","text":"cite AMR publications use: Berends MS, Luz CF, Friedrich AW, Sinha BNM, Albers CJ, Glasner C (2022). \"AMR: R Package Working Antimicrobial Resistance Data.\" Journal Statistical Software, 104(3), 1-31. doi:10.18637/jss.v104.i03 BibTeX entry LaTeX users :","code":"@Article{, title = {{AMR}: An {R} Package for Working with Antimicrobial Resistance Data}, author = {Matthijs S. Berends and Christian F. Luz and Alexander W. Friedrich and Bhanu N. M. Sinha and Casper J. Albers and Corinna Glasner}, journal = {Journal of Statistical Software}, year = {2022}, volume = {104}, number = {3}, pages = {1--31}, doi = {10.18637/jss.v104.i03}, }"},{"path":"https://msberends.github.io/AMR/reference/AMR.html","id":"reference-data-publicly-available","dir":"Reference","previous_headings":"","what":"Reference Data Publicly Available","title":"The AMR Package — AMR","text":"data sets AMR package (microorganisms, antibiotics, SIR interpretation, EUCAST rules, etc.) publicly freely available download following formats: R, MS Excel, Apache Feather, Apache Parquet, SPSS, SAS, Stata. also provide tab-separated plain text files machine-readable suitable input software program, laboratory information systems. Please visit website download links. actual files course available GitHub repository.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/reference/AMR.html","id":"author","dir":"Reference","previous_headings":"","what":"Author","title":"The AMR Package — AMR","text":"Maintainer: Matthijs S. Berends m.s.berends@umcg.nl (ORCID) Authors: Christian F. Luz (ORCID) [contributor] Dennis Souverein (ORCID) [contributor] Erwin E. . Hassing [contributor] contributors: Casper J. Albers (ORCID) [thesis advisor] Peter Dutey-Magni (ORCID) [contributor] Judith M. Fonville [contributor] Alex W. Friedrich (ORCID) [thesis advisor] Corinna Glasner (ORCID) [thesis advisor] Eric H. L. C. M. Hazenberg [contributor] Gwen Knight (ORCID) [contributor] Annick Lenglet (ORCID) [contributor] Bart C. Meijer [contributor] Dmytro Mykhailenko [contributor] Anton Mymrikov [contributor] Andrew P. Norgan (ORCID) [contributor] Sofia Ny (ORCID) [contributor] Jonas Salm [contributor] Rogier P. Schade [contributor] Bhanu N. M. Sinha (ORCID) [thesis advisor] Anthony Underwood (ORCID) [contributor] Anita Williams (ORCID) [contributor]","code":""},{"path":"https://msberends.github.io/AMR/reference/WHOCC.html","id":null,"dir":"Reference","previous_headings":"","what":"WHOCC: WHO Collaborating Centre for Drug Statistics Methodology — WHOCC","title":"WHOCC: WHO Collaborating Centre for Drug Statistics Methodology — WHOCC","text":"antimicrobial drugs official names, ATC codes, ATC groups defined daily dose (DDD) included package, using Collaborating Centre Drug Statistics Methodology.","code":""},{"path":"https://msberends.github.io/AMR/reference/WHOCC.html","id":"whocc","dir":"Reference","previous_headings":"","what":"WHOCC","title":"WHOCC: WHO Collaborating Centre for Drug Statistics Methodology — WHOCC","text":"package contains ~550 antibiotic, antimycotic antiviral drugs Anatomical Therapeutic Chemical (ATC) codes, ATC groups Defined Daily Dose (DDD) World Health Organization Collaborating Centre Drug Statistics Methodology (WHOCC, https://www.whocc.) Pharmaceuticals Community Register European Commission (https://ec.europa.eu/health/documents/community-register/html/reg_hum_atc.htm). become gold standard international drug utilisation monitoring research. WHOCC located Oslo Norwegian Institute Public Health funded Norwegian government. European Commission executive European Union promotes general interest. NOTE: WHOCC copyright allow use commercial purposes, unlike info package. See https://www.whocc./copyright_disclaimer/.","code":""},{"path":"https://msberends.github.io/AMR/reference/WHOCC.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"WHOCC: WHO Collaborating Centre for Drug Statistics Methodology — WHOCC","text":"","code":"as.ab(\"meropenem\") #> Class 'ab' #> [1] MEM ab_name(\"J01DH02\") #> [1] \"Meropenem\" ab_tradenames(\"flucloxacillin\") #> [1] \"culpen\" \"floxacillin\" \"floxacillin sodium\" #> [4] \"floxapen\" \"floxapen sodium salt\" \"fluclox\" #> [7] \"flucloxacilina\" \"flucloxacillin\" \"flucloxacilline\" #> [10] \"flucloxacillinum\" \"fluorochloroxacillin\" \"staphylex\""},{"path":"https://msberends.github.io/AMR/reference/WHONET.html","id":null,"dir":"Reference","previous_headings":"","what":"Data Set with 500 Isolates - WHONET Example — WHONET","title":"Data Set with 500 Isolates - WHONET Example — WHONET","text":"example data set exact structure export file WHONET. files can used package, example data set shows. antibiotic results example_isolates data set. patient names created using online surname generators place practice purposes.","code":""},{"path":"https://msberends.github.io/AMR/reference/WHONET.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Data Set with 500 Isolates - WHONET Example — WHONET","text":"","code":"WHONET"},{"path":"https://msberends.github.io/AMR/reference/WHONET.html","id":"format","dir":"Reference","previous_headings":"","what":"Format","title":"Data Set with 500 Isolates - WHONET Example — WHONET","text":"tibble 500 observations 53 variables: Identification number ID sample Specimen number ID specimen Organism Name microorganism. analysis, transform valid microbial class, using .mo(). Country Country origin Laboratory Name laboratory Last name Fictitious last name patient First name Fictitious initial patient Sex Fictitious gender patient Age Fictitious age patient Age category Age group, can also looked using age_groups() Date admissionDate hospital admission Specimen dateDate specimen received laboratory Specimen type Specimen type group Specimen type (Numeric) Translation \"Specimen type\" Reason Reason request Differential Diagnosis Isolate number ID isolate Organism type Type microorganism, can also looked using mo_type() Serotype Serotype microorganism Beta-lactamase Microorganism produces beta-lactamase? ESBL Microorganism produces extended spectrum beta-lactamase? Carbapenemase Microorganism produces carbapenemase? MRSA screening test Microorganism possible MRSA? Inducible clindamycin resistance Clindamycin can induced? Comment comments Date data entryDate data entered WHONET AMP_ND10:CIP_EE 28 different antibiotics. can lookup abbreviations antibiotics data set, use e.g. ab_name(\"AMP\") get official name immediately. analysis, transform valid antibiotic class, using .sir().","code":""},{"path":"https://msberends.github.io/AMR/reference/WHONET.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Data Set with 500 Isolates - WHONET Example — WHONET","text":"Like data sets package, data set publicly available download following formats: R, MS Excel, Apache Feather, Apache Parquet, SPSS, SAS, Stata. Please visit website download links. actual files course available GitHub repository.","code":""},{"path":"https://msberends.github.io/AMR/reference/WHONET.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Data Set with 500 Isolates - WHONET Example — WHONET","text":"","code":"WHONET #> # A tibble: 500 × 53 #> `Identification number` `Specimen number` Organism Country Laboratory #> #> 1 fe41d7bafa 1748 SPN Belgium National … #> 2 91f175ec37 1767 eco The Netherlands National … #> 3 cc4015056e 1343 eco The Netherlands National … #> 4 e864b692f5 1894 MAP Denmark National … #> 5 3d051fe345 1739 PVU Belgium National … #> 6 c80762a08d 1846 103 The Netherlands National … #> 7 8022d3727c 1628 103 Denmark National … #> 8 f3dc5f553d 1493 eco The Netherlands National … #> 9 15add38f6c 1847 eco France National … #> 10 fd41248def 1458 eco Germany National … #> # ℹ 490 more rows #> # ℹ 48 more variables: `Last name` , `First name` , Sex , #> # Age , `Age category` , `Date of admission` , #> # `Specimen date` , `Specimen type` , #> # `Specimen type (Numeric)` , Reason , `Isolate number` , #> # `Organism type` , Serotype , `Beta-lactamase` , ESBL , #> # Carbapenemase , `MRSA screening test` , …"},{"path":"https://msberends.github.io/AMR/reference/ab_from_text.html","id":null,"dir":"Reference","previous_headings":"","what":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","title":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","text":"Use function e.g. clinical texts health care records. returns list antimicrobial drugs, doses forms administration found texts.","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_from_text.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","text":"","code":"ab_from_text( text, type = c(\"drug\", \"dose\", \"administration\"), collapse = NULL, translate_ab = FALSE, thorough_search = NULL, info = interactive(), ... )"},{"path":"https://msberends.github.io/AMR/reference/ab_from_text.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","text":"text text analyse type type property search , either \"drug\", \"dose\" \"administration\", see Examples collapse character pass paste(, collapse = ...) return one character per element text, see Examples translate_ab type = \"drug\": column name antibiotics data set translate antibiotic abbreviations , using ab_property(). default FALSE. Using TRUE equal using \"name\". thorough_search logical indicate whether input must extensively searched misspelling faulty input values. Setting TRUE take considerably time using FALSE. default, turn TRUE input elements contain maximum three words. info logical indicate whether progress bar printed - default TRUE interactive mode ... arguments passed .ab()","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_from_text.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","text":"list, character collapse NULL","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_from_text.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","text":"function also internally used .ab(), although searches first drug name throw note drug names returned. Note: .ab() function may use long regular expression match brand names antimicrobial drugs. may fail systems.","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_from_text.html","id":"argument-type","dir":"Reference","previous_headings":"","what":"Argument type","title":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","text":"default, function search antimicrobial drug names. text elements searched official names, ATC codes brand names. uses .ab() internally, correct misspelling. type = \"dose\" (similar, like \"dosing\", \"doses\"), text elements searched numeric values higher 100 resemble years. output numeric. supports unit (g, mg, IE, etc.) multiple values one clinical text, see Examples. type = \"administration\" (abbreviations, like \"admin\", \"adm\"), text elements searched form drug administration. supports following forms (including common abbreviations): buccal, implant, inhalation, instillation, intravenous, nasal, oral, parenteral, rectal, sublingual, transdermal vaginal. Abbreviations oral ('po', 'per os') become \"oral\", values intravenous ('iv', 'intraven') become \"iv\". supports multiple values one clinical text, see Examples.","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_from_text.html","id":"argument-collapse","dir":"Reference","previous_headings":"","what":"Argument collapse","title":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","text":"Without using collapse, function return list. can convenient use e.g. inside mutate()):df %>% mutate(abx = ab_from_text(clinical_text)) returned AB codes can transformed official names, groups, etc. ab_* functions ab_name() ab_group(), using translate_ab argument. using collapse, function return character:df %>% mutate(abx = ab_from_text(clinical_text, collapse = \"|\"))","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_from_text.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Retrieve Antimicrobial Drug Names and Doses from Clinical Text — ab_from_text","text":"","code":"# mind the bad spelling of amoxicillin in this line, # straight from a true health care record: ab_from_text(\"28/03/2020 regular amoxicilliin 500mg po tid\") #> [[1]] #> Class 'ab' #> [1] AMX #> ab_from_text(\"500 mg amoxi po and 400mg cipro iv\") #> [[1]] #> Class 'ab' #> [1] AMX CIP #> ab_from_text(\"500 mg amoxi po and 400mg cipro iv\", type = \"dose\") #> [[1]] #> [1] 500 400 #> ab_from_text(\"500 mg amoxi po and 400mg cipro iv\", type = \"admin\") #> [[1]] #> [1] \"oral\" \"iv\" #> ab_from_text(\"500 mg amoxi po and 400mg cipro iv\", collapse = \", \") #> [1] \"AMX, CIP\" # \\donttest{ # if you want to know which antibiotic groups were administered, do e.g.: abx <- ab_from_text(\"500 mg amoxi po and 400mg cipro iv\") ab_group(abx[[1]]) #> [1] \"Beta-lactams/penicillins\" \"Quinolones\" if (require(\"dplyr\")) { tibble(clinical_text = c( \"given 400mg cipro and 500 mg amox\", \"started on doxy iv today\" )) %>% mutate( abx_codes = ab_from_text(clinical_text), abx_doses = ab_from_text(clinical_text, type = \"doses\"), abx_admin = ab_from_text(clinical_text, type = \"admin\"), abx_coll = ab_from_text(clinical_text, collapse = \"|\"), abx_coll_names = ab_from_text(clinical_text, collapse = \"|\", translate_ab = \"name\" ), abx_coll_doses = ab_from_text(clinical_text, type = \"doses\", collapse = \"|\" ), abx_coll_admin = ab_from_text(clinical_text, type = \"admin\", collapse = \"|\" ) ) } #> Loading required package: dplyr #> #> Attaching package: ‘dplyr’ #> The following objects are masked from ‘package:stats’: #> #> filter, lag #> The following objects are masked from ‘package:base’: #> #> intersect, setdiff, setequal, union #> # A tibble: 2 × 8 #> clinical_text abx_codes abx_doses abx_admin abx_coll abx_coll_names #> #> 1 given 400mg cipro and 5… CIP|AMX Ciprofloxacin… #> 2 started on doxy iv today DOX Doxycycline #> # ℹ 2 more variables: abx_coll_doses , abx_coll_admin # }"},{"path":"https://msberends.github.io/AMR/reference/ab_property.html","id":null,"dir":"Reference","previous_headings":"","what":"Get Properties of an Antibiotic — ab_property","title":"Get Properties of an Antibiotic — ab_property","text":"Use functions return specific property antibiotic antibiotics data set. input values evaluated internally .ab().","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_property.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get Properties of an Antibiotic — ab_property","text":"","code":"ab_name(x, language = get_AMR_locale(), tolower = FALSE, ...) ab_cid(x, ...) ab_synonyms(x, ...) ab_tradenames(x, ...) ab_group(x, language = get_AMR_locale(), ...) ab_atc(x, only_first = FALSE, ...) ab_atc_group1(x, language = get_AMR_locale(), ...) ab_atc_group2(x, language = get_AMR_locale(), ...) ab_loinc(x, ...) ab_ddd(x, administration = \"oral\", ...) ab_ddd_units(x, administration = \"oral\", ...) ab_info(x, language = get_AMR_locale(), ...) ab_url(x, open = FALSE, ...) ab_property(x, property = \"name\", language = get_AMR_locale(), ...) set_ab_names( data, ..., property = \"name\", language = get_AMR_locale(), snake_case = NULL )"},{"path":"https://msberends.github.io/AMR/reference/ab_property.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Get Properties of an Antibiotic — ab_property","text":"x (vector ) text can coerced valid antibiotic drug code .ab() language language returned text - default current system language (see get_AMR_locale()) can also set package option AMR_locale. Use language = NULL language = \"\" prevent translation. tolower logical indicate whether first character every output transformed lower case character. lead e.g. \"polymyxin B\" \"polymyxin b\". ... case set_ab_names() data data.frame: columns select (supports tidy selection column1:column4), otherwise arguments passed .ab() only_first logical indicate whether first ATC code must returned, giving preference J0-codes (.e., antimicrobial drug group) administration way administration, either \"oral\" \"iv\" open browse URL using utils::browseURL() property one column names one antibiotics data set: vector_or(colnames(antibiotics), sort = FALSE). data data.frame columns need renamed, character vector column names snake_case logical indicate whether names -called snake case: lower case spaces/slashes replaced underscore (_)","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_property.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Get Properties of an Antibiotic — ab_property","text":"integer case ab_cid() named list case ab_info() multiple ab_atc()/ab_synonyms()/ab_tradenames() double case ab_ddd() data.frame case set_ab_names() character cases","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_property.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Get Properties of an Antibiotic — ab_property","text":"output translated possible. function ab_url() return direct URL official website. warning returned required ATC code available. function set_ab_names() special column renaming function data.frames. renames columns names resemble antimicrobial drugs. always makes sure new column names unique. property = \"atc\" set, preference given ATC codes J-group.","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_property.html","id":"source","dir":"Reference","previous_headings":"","what":"Source","title":"Get Properties of an Antibiotic — ab_property","text":"World Health Organization () Collaborating Centre Drug Statistics Methodology: https://www.whocc./atc_ddd_index/ European Commission Public Health PHARMACEUTICALS - COMMUNITY REGISTER: https://ec.europa.eu/health/documents/community-register/html/reg_hum_atc.htm","code":""},{"path":"https://msberends.github.io/AMR/reference/ab_property.html","id":"reference-data-publicly-available","dir":"Reference","previous_headings":"","what":"Reference Data Publicly Available","title":"Get Properties of an Antibiotic — ab_property","text":"data sets AMR package (microorganisms, antibiotics, SIR interpretation, EUCAST rules, etc.) publicly freely available download following formats: R, MS Excel, Apache Feather, Apache Parquet, SPSS, SAS, Stata. also provide tab-separated plain text files machine-readable suitable input software program, laboratory information systems. Please visit website download links. actual files course available GitHub repository.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/reference/ab_property.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Get Properties of an Antibiotic — ab_property","text":"","code":"# all properties: ab_name(\"AMX\") #> [1] \"Amoxicillin\" ab_atc(\"AMX\") #> [1] \"J01CA04\" ab_cid(\"AMX\") #> [1] 33613 ab_synonyms(\"AMX\") #> [1] \"actimoxi\" \"amoclen\" \"amolin\" #> [4] \"amopen\" \"amopenixin\" \"amoxibiotic\" #> [7] \"amoxicaps\" \"amoxicilina\" \"amoxicillin\" #> [10] \"amoxicillin hydrate\" \"amoxicilline\" \"amoxicillinum\" #> [13] \"amoxiden\" \"amoxil\" \"amoxivet\" #> [16] \"amoxy\" \"amoxycillin\" \"amoxyke\" #> [19] \"anemolin\" \"aspenil\" \"atoksilin\" #> [22] \"biomox\" \"bristamox\" \"cemoxin\" #> [25] \"clamoxyl\" \"damoxy\" \"delacillin\" #> [28] \"demoksil\" \"dispermox\" \"efpenix\" #> [31] \"flemoxin\" \"hiconcil\" \"histocillin\" #> [34] \"hydroxyampicillin\" \"ibiamox\" \"imacillin\" #> [37] \"lamoxy\" \"largopen\" \"metafarma capsules\" #> [40] \"metifarma capsules\" \"moksilin\" \"moxacin\" #> [43] \"moxatag\" \"ospamox\" \"pamoxicillin\" #> [46] \"piramox\" \"promoxil\" \"remoxil\" #> [49] \"robamox\" \"sawamox pm\" \"tolodina\" #> [52] \"topramoxin\" \"unicillin\" \"utimox\" #> [55] \"vetramox\" ab_tradenames(\"AMX\") #> [1] \"actimoxi\" \"amoclen\" \"amolin\" #> [4] \"amopen\" \"amopenixin\" \"amoxibiotic\" #> [7] \"amoxicaps\" \"amoxicilina\" \"amoxicillin\" #> [10] \"amoxicillin hydrate\" \"amoxicilline\" \"amoxicillinum\" #> [13] \"amoxiden\" \"amoxil\" \"amoxivet\" #> [16] \"amoxy\" \"amoxycillin\" \"amoxyke\" #> [19] \"anemolin\" \"aspenil\" \"atoksilin\" #> [22] \"biomox\" \"bristamox\" \"cemoxin\" #> [25] \"clamoxyl\" \"damoxy\" \"delacillin\" #> [28] \"demoksil\" \"dispermox\" \"efpenix\" #> [31] \"flemoxin\" \"hiconcil\" \"histocillin\" #> [34] \"hydroxyampicillin\" \"ibiamox\" \"imacillin\" #> [37] \"lamoxy\" \"largopen\" \"metafarma capsules\" #> [40] \"metifarma capsules\" \"moksilin\" \"moxacin\" #> [43] \"moxatag\" \"ospamox\" \"pamoxicillin\" #> [46] \"piramox\" \"promoxil\" \"remoxil\" #> [49] \"robamox\" \"sawamox pm\" \"tolodina\" #> [52] \"topramoxin\" \"unicillin\" \"utimox\" #> [55] \"vetramox\" ab_group(\"AMX\") #> [1] \"Beta-lactams/penicillins\" ab_atc_group1(\"AMX\") #> [1] \"Beta-lactam antibacterials, penicillins\" ab_atc_group2(\"AMX\") #> [1] \"Penicillins with extended spectrum\" ab_url(\"AMX\") #> Amoxicillin #> \"https://www.whocc.no/atc_ddd_index/?code=J01CA04&showdescription=no\" # smart lowercase transformation ab_name(x = c(\"AMC\", \"PLB\")) #> [1] \"Amoxicillin/clavulanic acid\" \"Polymyxin B\" ab_name(x = c(\"AMC\", \"PLB\"), tolower = TRUE) #> [1] \"amoxicillin/clavulanic acid\" \"polymyxin B\" # defined daily doses (DDD) ab_ddd(\"AMX\", \"oral\") #> [1] 1.5 ab_ddd_units(\"AMX\", \"oral\") #> [1] \"g\" ab_ddd(\"AMX\", \"iv\") #> [1] 3 ab_ddd_units(\"AMX\", \"iv\") #> [1] \"g\" ab_info(\"AMX\") # all properties as a list #> $ab #> [1] \"AMX\" #> #> $cid #> [1] 33613 #> #> $name #> [1] \"Amoxicillin\" #> #> $group #> [1] \"Beta-lactams/penicillins\" #> #> $atc #> [1] \"J01CA04\" #> #> $atc_group1 #> [1] \"Beta-lactam antibacterials, penicillins\" #> #> $atc_group2 #> [1] \"Penicillins with extended spectrum\" #> #> $tradenames #> [1] \"actimoxi\" \"amoclen\" \"amolin\" #> [4] \"amopen\" \"amopenixin\" \"amoxibiotic\" #> [7] \"amoxicaps\" \"amoxicilina\" \"amoxicillin\" #> [10] \"amoxicillin hydrate\" \"amoxicilline\" \"amoxicillinum\" #> [13] \"amoxiden\" \"amoxil\" \"amoxivet\" #> [16] \"amoxy\" \"amoxycillin\" \"amoxyke\" #> [19] \"anemolin\" \"aspenil\" \"atoksilin\" #> [22] \"biomox\" \"bristamox\" \"cemoxin\" #> [25] \"clamoxyl\" \"damoxy\" \"delacillin\" #> [28] \"demoksil\" \"dispermox\" \"efpenix\" #> [31] \"flemoxin\" \"hiconcil\" \"histocillin\" #> [34] \"hydroxyampicillin\" \"ibiamox\" \"imacillin\" #> [37] \"lamoxy\" \"largopen\" \"metafarma capsules\" #> [40] \"metifarma capsules\" \"moksilin\" \"moxacin\" #> [43] \"moxatag\" \"ospamox\" \"pamoxicillin\" #> [46] \"piramox\" \"promoxil\" \"remoxil\" #> [49] \"robamox\" \"sawamox pm\" \"tolodina\" #> [52] \"topramoxin\" \"unicillin\" \"utimox\" #> [55] \"vetramox\" #> #> $loinc #> [1] \"101498-4\" \"15-8\" \"16-6\" \"16365-9\" \"17-4\" \"18-2\" #> [7] \"18861-5\" \"18862-3\" \"19-0\" \"20-8\" \"21-6\" \"22-4\" #> [13] \"25274-2\" \"25310-4\" \"3344-9\" \"55614-2\" \"55615-9\" \"55616-7\" #> [19] \"6976-5\" \"6977-3\" \"80133-2\" #> #> $ddd #> $ddd$oral #> $ddd$oral$amount #> [1] 1.5 #> #> $ddd$oral$units #> [1] \"g\" #> #> #> $ddd$iv #> $ddd$iv$amount #> [1] 3 #> #> $ddd$iv$units #> [1] \"g\" #> #> #> # all ab_* functions use as.ab() internally, so you can go from 'any' to 'any': ab_atc(\"AMP\") #> [1] \"J01CA01\" \"S01AA19\" ab_group(\"J01CA01\") #> [1] \"Beta-lactams/penicillins\" ab_loinc(\"ampicillin\") #> [1] \"101477-8\" \"101478-6\" \"18864-9\" \"18865-6\" \"20374-5\" \"21066-6\" #> [7] \"23618-2\" \"27-3\" \"28-1\" \"29-9\" \"30-7\" \"31-5\" #> [13] \"32-3\" \"33-1\" \"3355-5\" \"33562-0\" \"33919-2\" \"34-9\" #> [19] \"43883-8\" \"43884-6\" \"6979-9\" \"6980-7\" \"87604-5\" ab_name(\"21066-6\") #> [1] \"Ampicillin\" ab_name(6249) #> [1] \"Ampicillin\" ab_name(\"J01CA01\") #> [1] \"Ampicillin\" # spelling from different languages and dyslexia are no problem ab_atc(\"ceftriaxon\") #> [1] \"J01DD04\" ab_atc(\"cephtriaxone\") #> [1] \"J01DD04\" ab_atc(\"cephthriaxone\") #> [1] \"J01DD04\" ab_atc(\"seephthriaaksone\") #> [1] \"J01DD04\" # use set_ab_names() for renaming columns colnames(example_isolates) #> [1] \"date\" \"patient\" \"age\" \"gender\" \"ward\" \"mo\" \"PEN\" #> [8] \"OXA\" \"FLC\" \"AMX\" \"AMC\" \"AMP\" \"TZP\" \"CZO\" #> [15] \"FEP\" \"CXM\" \"FOX\" \"CTX\" \"CAZ\" \"CRO\" \"GEN\" #> [22] \"TOB\" \"AMK\" \"KAN\" \"TMP\" \"SXT\" \"NIT\" \"FOS\" #> [29] \"LNZ\" \"CIP\" \"MFX\" \"VAN\" \"TEC\" \"TCY\" \"TGC\" #> [36] \"DOX\" \"ERY\" \"CLI\" \"AZM\" \"IPM\" \"MEM\" \"MTR\" #> [43] \"CHL\" \"COL\" \"MUP\" \"RIF\" colnames(set_ab_names(example_isolates)) #> [1] \"date\" \"patient\" #> [3] \"age\" \"gender\" #> [5] \"ward\" \"mo\" #> [7] \"benzylpenicillin\" \"oxacillin\" #> [9] \"flucloxacillin\" \"amoxicillin\" #> [11] \"amoxicillin_clavulanic_acid\" \"ampicillin\" #> [13] \"piperacillin_tazobactam\" \"cefazolin\" #> [15] \"cefepime\" \"cefuroxime\" #> [17] \"cefoxitin\" \"cefotaxime\" #> [19] \"ceftazidime\" \"ceftriaxone\" #> [21] \"gentamicin\" \"tobramycin\" #> [23] \"amikacin\" \"kanamycin\" #> [25] \"trimethoprim\" \"trimethoprim_sulfamethoxazole\" #> [27] \"nitrofurantoin\" \"fosfomycin\" #> [29] \"linezolid\" \"ciprofloxacin\" #> [31] \"moxifloxacin\" \"vancomycin\" #> [33] \"teicoplanin\" \"tetracycline\" #> [35] \"tigecycline\" \"doxycycline\" #> [37] \"erythromycin\" \"clindamycin\" #> [39] \"azithromycin\" \"imipenem\" #> [41] \"meropenem\" \"metronidazole\" #> [43] \"chloramphenicol\" \"colistin\" #> [45] \"mupirocin\" \"rifampicin\" colnames(set_ab_names(example_isolates, NIT:VAN)) #> [1] \"date\" \"patient\" \"age\" \"gender\" #> [5] \"ward\" \"mo\" \"PEN\" \"OXA\" #> [9] \"FLC\" \"AMX\" \"AMC\" \"AMP\" #> [13] \"TZP\" \"CZO\" \"FEP\" \"CXM\" #> [17] \"FOX\" \"CTX\" \"CAZ\" \"CRO\" #> [21] \"GEN\" \"TOB\" \"AMK\" \"KAN\" #> [25] \"TMP\" \"SXT\" \"nitrofurantoin\" \"fosfomycin\" #> [29] \"linezolid\" \"ciprofloxacin\" \"moxifloxacin\" \"vancomycin\" #> [33] \"TEC\" \"TCY\" \"TGC\" \"DOX\" #> [37] \"ERY\" \"CLI\" \"AZM\" \"IPM\" #> [41] \"MEM\" \"MTR\" \"CHL\" \"COL\" #> [45] \"MUP\" \"RIF\" # \\donttest{ if (require(\"dplyr\")) { example_isolates %>% set_ab_names() # this does the same: example_isolates %>% rename_with(set_ab_names) # set_ab_names() works with any AB property: example_isolates %>% set_ab_names(property = \"atc\") example_isolates %>% set_ab_names(where(is.sir)) %>% colnames() example_isolates %>% set_ab_names(NIT:VAN) %>% colnames() } #> [1] \"date\" \"patient\" \"age\" \"gender\" #> [5] \"ward\" \"mo\" \"PEN\" \"OXA\" #> [9] \"FLC\" \"AMX\" \"AMC\" \"AMP\" #> [13] \"TZP\" \"CZO\" \"FEP\" \"CXM\" #> [17] \"FOX\" \"CTX\" \"CAZ\" \"CRO\" #> [21] \"GEN\" \"TOB\" \"AMK\" \"KAN\" #> [25] \"TMP\" \"SXT\" \"nitrofurantoin\" \"fosfomycin\" #> [29] \"linezolid\" \"ciprofloxacin\" \"moxifloxacin\" \"vancomycin\" #> [33] \"TEC\" \"TCY\" \"TGC\" \"DOX\" #> [37] \"ERY\" \"CLI\" \"AZM\" \"IPM\" #> [41] \"MEM\" \"MTR\" \"CHL\" \"COL\" #> [45] \"MUP\" \"RIF\" # }"},{"path":"https://msberends.github.io/AMR/reference/add_custom_antimicrobials.html","id":null,"dir":"Reference","previous_headings":"","what":"Add Custom Antimicrobials — add_custom_antimicrobials","title":"Add Custom Antimicrobials — add_custom_antimicrobials","text":"add_custom_antimicrobials() can add custom antimicrobial drug names codes.","code":""},{"path":"https://msberends.github.io/AMR/reference/add_custom_antimicrobials.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Add Custom Antimicrobials — add_custom_antimicrobials","text":"","code":"add_custom_antimicrobials(x) clear_custom_antimicrobials()"},{"path":"https://msberends.github.io/AMR/reference/add_custom_antimicrobials.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Add Custom Antimicrobials — add_custom_antimicrobials","text":"x data.frame resembling antibiotics data set, least containing columns \"ab\" \"name\"","code":""},{"path":"https://msberends.github.io/AMR/reference/add_custom_antimicrobials.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Add Custom Antimicrobials — add_custom_antimicrobials","text":"Important: Due R works, add_custom_antimicrobials() function run every R session - added antimicrobials stored sessions thus lost R exited. two ways circumvent automate process adding antimicrobials: Method 1: Using package option AMR_custom_ab, preferred method. use method: Create data set structure antibiotics data set (containing least columns \"ab\" \"name\") save saveRDS() location choice, e.g. \"~/my_custom_ab.rds\", remote location. Set file location package option AMR_custom_ab: options(AMR_custom_ab = \"~/my_custom_ab.rds\"). can even remote file location, https URL. Since options saved R sessions, best save option .Rprofile file loaded start-R. , open .Rprofile file using e.g. utils::file.edit(\"~/.Rprofile\"), add text save file: Upon package load, file loaded run add_custom_antimicrobials() function. Method 2: Loading antimicrobial additions directly .Rprofile file. Note definitions stored user-specific R file, suboptimal workflow. use method: Edit .Rprofile file using e.g. utils::file.edit(\"~/.Rprofile\"). Add text like save file: Use clear_custom_antimicrobials() clear previously added antimicrobials.","code":"# Add custom antimicrobial codes: options(AMR_custom_ab = \"~/my_custom_ab.rds\") # Add custom antibiotic drug codes: AMR::add_custom_antimicrobials( data.frame(ab = \"TESTAB\", name = \"Test Antibiotic\", group = \"Test Group\") )"},{"path":[]},{"path":"https://msberends.github.io/AMR/reference/add_custom_antimicrobials.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Add Custom Antimicrobials — add_custom_antimicrobials","text":"","code":"# \\donttest{ # returns NA and throws a warning (which is suppressed here): suppressWarnings( as.ab(\"testab\") ) #> Class 'ab' #> [1] # now add a custom entry - it will be considered by as.ab() and # all ab_*() functions add_custom_antimicrobials( data.frame( ab = \"TESTAB\", name = \"Test Antibiotic\", # you can add any property present in the # 'antibiotics' data set, such as 'group': group = \"Test Group\" ) ) #> ℹ Added one record to the internal antibiotics data set. # \"testab\" is now a new antibiotic: as.ab(\"testab\") #> Class 'ab' #> [1] TESTAB ab_name(\"testab\") #> [1] \"Test Antibiotic\" ab_group(\"testab\") #> [1] \"Test Group\" ab_info(\"testab\") #> $ab #> [1] \"TESTAB\" #> #> $cid #> [1] NA #> #> $name #> [1] \"Test Antibiotic\" #> #> $group #> [1] \"Test Group\" #> #> $atc #> [1] NA #> #> $atc_group1 #> [1] NA #> #> $atc_group2 #> [1] NA #> #> $tradenames #> [1] NA #> #> $loinc #> [1] NA #> #> $ddd #> $ddd$oral #> $ddd$oral$amount #> [1] NA #> #> $ddd$oral$units #> [1] NA #> #> #> $ddd$iv #> $ddd$iv$amount #> [1] NA #> #> $ddd$iv$units #> [1] NA #> #> #> # Add Co-fluampicil, which is one of the many J01CR50 codes, see # https://www.whocc.no/ddd/list_of_ddds_combined_products/ add_custom_antimicrobials( data.frame( ab = \"COFLU\", name = \"Co-fluampicil\", atc = \"J01CR50\", group = \"Beta-lactams/penicillins\" ) ) #> ℹ Added one record to the internal antibiotics data set. ab_atc(\"Co-fluampicil\") #> [1] \"J01CR50\" ab_name(\"J01CR50\") #> [1] \"Co-fluampicil\" # even antibiotic selectors work x <- data.frame( random_column = \"some value\", coflu = as.sir(\"S\"), ampicillin = as.sir(\"R\") ) x #> random_column coflu ampicillin #> 1 some value S R x[, betalactams()] #> ℹ For betalactams() using columns 'coflu' (co-fluampicil) and #> 'ampicillin' #> coflu ampicillin #> 1 S R # }"},{"path":"https://msberends.github.io/AMR/reference/add_custom_microorganisms.html","id":null,"dir":"Reference","previous_headings":"","what":"Add Custom Microorganisms — add_custom_microorganisms","title":"Add Custom Microorganisms — add_custom_microorganisms","text":"add_custom_microorganisms() can add custom microorganisms, non-taxonomic outcome laboratory analysis.","code":""},{"path":"https://msberends.github.io/AMR/reference/add_custom_microorganisms.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Add Custom Microorganisms — add_custom_microorganisms","text":"","code":"add_custom_microorganisms(x) clear_custom_microorganisms()"},{"path":"https://msberends.github.io/AMR/reference/add_custom_microorganisms.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Add Custom Microorganisms — add_custom_microorganisms","text":"x data.frame resembling microorganisms data set, least containing column \"genus\" (case-insensitive)","code":""},{"path":"https://msberends.github.io/AMR/reference/add_custom_microorganisms.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Add Custom Microorganisms — add_custom_microorganisms","text":"function fill missing taxonomy , specific taxonomic columns missing, see Examples. Important: Due R works, add_custom_microorganisms() function run every R session - added microorganisms stored sessions thus lost R exited. two ways circumvent automate process adding microorganisms: Method 1: Using package option AMR_custom_mo, preferred method. use method: Create data set structure microorganisms data set (containing least column \"genus\") save saveRDS() location choice, e.g. \"~/my_custom_mo.rds\", remote location. Set file location package option AMR_custom_mo: options(AMR_custom_mo = \"~/my_custom_mo.rds\"). can even remote file location, https URL. Since options saved R sessions, best save option .Rprofile file loaded start-R. , open .Rprofile file using e.g. utils::file.edit(\"~/.Rprofile\"), add text save file: Upon package load, file loaded run add_custom_microorganisms() function. Method 2: Loading microorganism directly .Rprofile file. Note definitions stored user-specific R file, suboptimal workflow. use method: Edit .Rprofile file using e.g. utils::file.edit(\"~/.Rprofile\"). Add text like save file: Use clear_custom_microorganisms() clear previously added microorganisms.","code":"# Add custom microorganism codes: options(AMR_custom_mo = \"~/my_custom_mo.rds\") # Add custom antibiotic drug codes: AMR::add_custom_microorganisms( data.frame(genus = \"Enterobacter\", species = \"asburiae/cloacae\") )"},{"path":[]},{"path":"https://msberends.github.io/AMR/reference/add_custom_microorganisms.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Add Custom Microorganisms — add_custom_microorganisms","text":"","code":"# \\donttest{ # a combination of species is not formal taxonomy, so # this will result in \"Enterobacter cloacae cloacae\", # since it resembles the input best: mo_name(\"Enterobacter asburiae/cloacae\") #> [1] \"Enterobacter cloacae cloacae\" # now add a custom entry - it will be considered by as.mo() and # all mo_*() functions add_custom_microorganisms( data.frame( genus = \"Enterobacter\", species = \"asburiae/cloacae\" ) ) #> ℹ Added Enterobacter asburiae/cloacae to the internal microorganisms data #> set. # E. asburiae/cloacae is now a new microorganism: mo_name(\"Enterobacter asburiae/cloacae\") #> [1] \"Enterobacter asburiae/cloacae\" # its code: as.mo(\"Enterobacter asburiae/cloacae\") #> Class 'mo' #> [1] CUSTOM1_ENTRB_ASB/ # all internal algorithms will work as well: mo_name(\"Ent asburia cloacae\") #> [1] \"Enterobacter asburiae/cloacae\" # and even the taxonomy was added based on the genus! mo_family(\"E. asburiae/cloacae\") #> [1] \"Enterobacteriaceae\" mo_gramstain(\"Enterobacter asburiae/cloacae\") #> [1] \"Gram-negative\" mo_info(\"Enterobacter asburiae/cloacae\") #> $mo #> [1] \"CUSTOM1_ENTRB_ASB/\" #> #> $kingdom #> [1] \"Bacteria\" #> #> $phylum #> [1] \"Pseudomonadota\" #> #> $class #> [1] \"Gammaproteobacteria\" #> #> $order #> [1] \"Enterobacterales\" #> #> $family #> [1] \"Enterobacteriaceae\" #> #> $genus #> [1] \"Enterobacter\" #> #> $species #> [1] \"asburiae/cloacae\" #> #> $subspecies #> [1] \"\" #> #> $status #> [1] \"accepted\" #> #> $synonyms #> NULL #> #> $gramstain #> [1] \"Gram-negative\" #> #> $oxygen_tolerance #> [1] NA #> #> $url #> [1] \"\" #> #> $ref #> [1] \"Self-added, 2023\" #> #> $snomed #> [1] NA #> #> $lpsn #> [1] NA #> #> $gbif #> [1] NA #> # the function tries to be forgiving: add_custom_microorganisms( data.frame( GENUS = \"BACTEROIDES / PARABACTEROIDES SLASHLINE\", SPECIES = \"SPECIES\" ) ) #> ℹ Added Bacteroides/Parabacteroides to the internal microorganisms data #> set. mo_name(\"BACTEROIDES / PARABACTEROIDES\") #> [1] \"Bacteroides/Parabacteroides\" mo_rank(\"BACTEROIDES / PARABACTEROIDES\") #> [1] \"genus\" # taxonomy still works, even though a slashline genus was given as input: mo_family(\"Bacteroides/Parabacteroides\") #> [1] \"Bacteroidaceae\" # for groups and complexes, set them as species or subspecies: add_custom_microorganisms( data.frame( genus = \"Citrobacter\", species = c(\"freundii\", \"braakii complex\"), subspecies = c(\"complex\", \"\") ) ) #> ℹ Added Citrobacter braakii complex and Citrobacter freundii complex to the #> internal microorganisms data set. mo_name(c(\"C. freundii complex\", \"C. braakii complex\")) #> [1] \"Citrobacter freundii complex\" \"Citrobacter braakii complex\" mo_species(c(\"C. freundii complex\", \"C. braakii complex\")) #> [1] \"freundii complex\" \"braakii complex\" mo_gramstain(c(\"C. freundii complex\", \"C. braakii complex\")) #> [1] \"Gram-negative\" \"Gram-negative\" # }"},{"path":"https://msberends.github.io/AMR/reference/age.html","id":null,"dir":"Reference","previous_headings":"","what":"Age in Years of Individuals — age","title":"Age in Years of Individuals — age","text":"Calculates age years based reference date, system date default.","code":""},{"path":"https://msberends.github.io/AMR/reference/age.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Age in Years of Individuals — age","text":"","code":"age(x, reference = Sys.Date(), exact = FALSE, na.rm = FALSE, ...)"},{"path":"https://msberends.github.io/AMR/reference/age.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Age in Years of Individuals — age","text":"x date(s), character (vectors) coerced .POSIXlt() reference reference date(s) (default today), character (vectors) coerced .POSIXlt() exact logical indicate whether age calculation exact, .e. decimals. divides number days year--date (YTD) x number days year reference (either 365 366). na.rm logical indicate whether missing values removed ... arguments passed .POSIXlt(), origin","code":""},{"path":"https://msberends.github.io/AMR/reference/age.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Age in Years of Individuals — age","text":"integer (decimals) exact = FALSE, double (decimals) otherwise","code":""},{"path":"https://msberends.github.io/AMR/reference/age.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Age in Years of Individuals — age","text":"Ages 0 returned NA warning. Ages 120 give warning. function vectorises x reference, meaning either can length 1 argument larger length.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/reference/age.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Age in Years of Individuals — age","text":"","code":"# 10 random pre-Y2K birth dates df <- data.frame(birth_date = as.Date(\"2000-01-01\") - runif(10) * 25000) # add ages df$age <- age(df$birth_date) # add exact ages df$age_exact <- age(df$birth_date, exact = TRUE) # add age at millenium switch df$age_at_y2k <- age(df$birth_date, \"2000-01-01\") df #> birth_date age age_exact age_at_y2k #> 1 1950-06-04 73 73.38356 49 #> 2 1953-03-27 70 70.57260 46 #> 3 1968-11-29 54 54.89589 31 #> 4 1997-09-02 26 26.13699 2 #> 5 1981-09-03 42 42.13425 18 #> 6 1974-02-03 49 49.71507 25 #> 7 1947-03-27 76 76.57260 52 #> 8 1972-02-11 51 51.69315 27 #> 9 1959-12-06 63 63.87671 40 #> 10 1935-10-25 87 87.99178 64"},{"path":"https://msberends.github.io/AMR/reference/age_groups.html","id":null,"dir":"Reference","previous_headings":"","what":"Split Ages into Age Groups — age_groups","title":"Split Ages into Age Groups — age_groups","text":"Split ages age groups defined split argument. allows easier demographic (antimicrobial resistance) analysis.","code":""},{"path":"https://msberends.github.io/AMR/reference/age_groups.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Split Ages into Age Groups — age_groups","text":"","code":"age_groups(x, split_at = c(12, 25, 55, 75), na.rm = FALSE)"},{"path":"https://msberends.github.io/AMR/reference/age_groups.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Split Ages into Age Groups — age_groups","text":"x age, e.g. calculated age() split_at values split x - default age groups 0-11, 12-24, 25-54, 55-74 75+. See Details. na.rm logical indicate whether missing values removed","code":""},{"path":"https://msberends.github.io/AMR/reference/age_groups.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Split Ages into Age Groups — age_groups","text":"Ordered factor","code":""},{"path":"https://msberends.github.io/AMR/reference/age_groups.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Split Ages into Age Groups — age_groups","text":"split ages, input split_at argument can : numeric vector. value e.g. c(10, 20) split x 0-9, 10-19 20+. value 50 split x 0-49 50+. default split young children (0-11), youth (12-24), young adults (25-54), middle-aged adults (55-74) elderly (75+). character: \"children\" \"kids\", equivalent : c(0, 1, 2, 4, 6, 13, 18). split 0, 1, 2-3, 4-5, 6-12, 13-17 18+. \"elderly\" \"seniors\", equivalent : c(65, 75, 85). split 0-64, 65-74, 75-84, 85+. \"fives\", equivalent : 1:20 * 5. split 0-4, 5-9, ..., 95-99, 100+. \"tens\", equivalent : 1:10 * 10. split 0-9, 10-19, ..., 90-99, 100+.","code":""},{"path":[]},{"path":"https://msberends.github.io/AMR/reference/age_groups.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Split Ages into Age Groups — age_groups","text":"","code":"ages <- c(3, 8, 16, 54, 31, 76, 101, 43, 21) # split into 0-49 and 50+ age_groups(ages, 50) #> [1] 0-49 0-49 0-49 50+ 0-49 50+ 50+ 0-49 0-49 #> Levels: 0-49 < 50+ # split into 0-19, 20-49 and 50+ age_groups(ages, c(20, 50)) #> [1] 0-19 0-19 0-19 50+ 20-49 50+ 50+ 20-49 20-49 #> Levels: 0-19 < 20-49 < 50+ # split into groups of ten years age_groups(ages, 1:10 * 10) #> [1] 0-9 0-9 10-19 50-59 30-39 70-79 100+ 40-49 20-29 #> 11 Levels: 0-9 < 10-19 < 20-29 < 30-39 < 40-49 < 50-59 < 60-69 < ... < 100+ age_groups(ages, split_at = \"tens\") #> [1] 0-9 0-9 10-19 50-59 30-39 70-79 100+ 40-49 20-29 #> 11 Levels: 0-9 < 10-19 < 20-29 < 30-39 < 40-49 < 50-59 < 60-69 < ... < 100+ # split into groups of five years age_groups(ages, 1:20 * 5) #> [1] 0-4 5-9 15-19 50-54 30-34 75-79 100+ 40-44 20-24 #> 21 Levels: 0-4 < 5-9 < 10-14 < 15-19 < 20-24 < 25-29 < 30-34 < ... < 100+ age_groups(ages, split_at = \"fives\") #> [1] 0-4 5-9 15-19 50-54 30-34 75-79 100+ 40-44 20-24 #> 21 Levels: 0-4 < 5-9 < 10-14 < 15-19 < 20-24 < 25-29 < 30-34 < ... < 100+ # split specifically for children age_groups(ages, c(1, 2, 4, 6, 13, 18)) #> [1] 2-3 6-12 13-17 18+ 18+ 18+ 18+ 18+ 18+ #> Levels: 0 < 1 < 2-3 < 4-5 < 6-12 < 13-17 < 18+ age_groups(ages, \"children\") #> [1] 2-3 6-12 13-17 18+ 18+ 18+ 18+ 18+ 18+ #> Levels: 0 < 1 < 2-3 < 4-5 < 6-12 < 13-17 < 18+ # \\donttest{ # resistance of ciprofloxacin per age group if (require(\"dplyr\") && require(\"ggplot2\")) { example_isolates %>% filter_first_isolate() %>% filter(mo == as.mo(\"Escherichia coli\")) %>% group_by(age_group = age_groups(age)) %>% select(age_group, CIP) %>% ggplot_sir( x = \"age_group\", minimum = 0, x.title = \"Age Group\", title = \"Ciprofloxacin resistance per age group\" ) } #> Loading required package: ggplot2 # }"},{"path":"https://msberends.github.io/AMR/reference/antibiogram.html","id":null,"dir":"Reference","previous_headings":"","what":"Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram","title":"Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram","text":"Generate antibiogram, communicate results plots tables. functions follow logic Klinker et al. Barbieri et al. (see Source), allow reporting e.g. R Markdown Quarto well.","code":""},{"path":"https://msberends.github.io/AMR/reference/antibiogram.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram","text":"","code":"antibiogram( x, antibiotics = where(is.sir), mo_transform = \"shortname\", ab_transform = NULL, syndromic_group = NULL, add_total_n = TRUE, only_all_tested = FALSE, digits = 0, col_mo = NULL, language = get_AMR_locale(), minimum = 30, combine_SI = TRUE, sep = \" + \", info = interactive() ) # S3 method for antibiogram plot(x, ...) # S3 method for antibiogram autoplot(object, ...) # S3 method for antibiogram knit_print( x, italicise = TRUE, na = getOption(\"knitr.kable.NA\", default = \"\"), ... )"},{"path":"https://msberends.github.io/AMR/reference/antibiogram.html","id":"source","dir":"Reference","previous_headings":"","what":"Source","title":"Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram","text":"Klinker KP et al. (2021). Antimicrobial stewardship antibiograms: importance moving beyond traditional antibiograms. Therapeutic Advances Infectious Disease, May 5;8:20499361211011373; doi:10.1177/20499361211011373 Barbieri E et al. (2021). Development Weighted-Incidence Syndromic Combination Antibiogram (WISCA) guide choice empiric antibiotic treatment urinary tract infection paediatric patients: Bayesian approach Antimicrobial Resistance & Infection Control May 1;10(1):74; doi:10.1186/s13756-021-00939-2 M39 Analysis Presentation Cumulative Antimicrobial Susceptibility Test Data, 5th Edition, 2022, Clinical Laboratory Standards Institute (CLSI). https://clsi.org/standards/products/microbiology/documents/m39/.","code":""},{"path":"https://msberends.github.io/AMR/reference/antibiogram.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram","text":"x data.frame containing least column microorganisms columns antibiotic results (class 'sir', see .sir()) antibiotics vector antibiotic name code (evaluated .ab(), column name x, (combinations ) antibiotic selectors aminoglycosides() carbapenems(). combination antibiograms, can also set values separated \"+\", \"TZP+TOB\" \"cipro + genta\", given columns resembling antibiotics exist x. See Examples. mo_transform character transform microorganism input - must \"name\", \"shortname\", \"gramstain\", one column names microorganisms data set: \"mo\", \"fullname\", \"status\", \"kingdom\", \"phylum\", \"class\", \"order\", \"family\", \"genus\", \"species\", \"subspecies\", \"rank\", \"ref\", \"oxygen_tolerance\", \"source\", \"lpsn\", \"lpsn_parent\", \"lpsn_renamed_to\", \"gbif\", \"gbif_parent\", \"gbif_renamed_to\", \"prevalence\", \"snomed\". Can also NULL transform input. ab_transform character transform antibiotic input - must one column names antibiotics data set: \"ab\", \"cid\", \"name\", \"group\", \"atc\", \"atc_group1\", \"atc_group2\", \"abbreviations\", \"synonyms\", \"oral_ddd\", \"oral_units\", \"iv_ddd\", \"iv_units\", \"loinc\". Can also NULL transform input. syndromic_group column name x, values calculated split rows x, e.g. using ifelse() case_when(). See Examples. add_total_n logical indicate whether total available numbers per pathogen added table (default TRUE). add lowest highest number available isolate per antibiotic (e.g, E. coli 200 isolates available ciprofloxacin 150 amoxicillin, returned number \"150-200\"). only_all_tested (combination antibiograms): logical indicate isolates must tested antibiotics, see Details digits number digits use rounding col_mo column name names codes microorganisms (see .mo()) - default first column class mo. Values coerced using .mo(). language language translate text, defaults system language (see get_AMR_locale()) minimum minimum allowed number available (tested) isolates. isolate count lower minimum return NA warning. default number 30 isolates advised Clinical Laboratory Standards Institute (CLSI) best practice, see Source. combine_SI logical indicate whether susceptibility determined results either S , instead S (default TRUE) sep separating character antibiotic columns combination antibiograms info logical indicate info printed - default TRUE interactive mode ... used R Markdown Quarto: arguments passed knitr::kable() (otherwise, use) object antibiogram() object italicise logical indicate whether microorganism names knitr table made italic, using italicise_taxonomy(). na character use showing NA values","code":""},{"path":"https://msberends.github.io/AMR/reference/antibiogram.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram","text":"function returns table values 0 100 susceptibility, resistance. Remember filter data let contain first isolates! needed exclude duplicates reduce selection bias. Use first_isolate() determine data set one four available algorithms. types antibiograms listed can plotted (using ggplot2::autoplot() base R plot()/barplot()). antibiogram object can also used directly R Markdown / Quarto (.e., knitr) reports. case, knitr::kable() applied automatically microorganism names even printed italics default (see argument italicise). can also use functions specific 'table reporting' packages transform output antibiogram() needs, e.g. flextable::as_flextable() gt::gt().","code":""},{"path":"https://msberends.github.io/AMR/reference/antibiogram.html","id":"antibiogram-types","dir":"Reference","previous_headings":"","what":"Antibiogram Types","title":"Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram","text":"four antibiogram types, proposed Klinker et al. (2021, doi:10.1177/20499361211011373 ), supported antibiogram(): Traditional Antibiogram Case example: Susceptibility Pseudomonas aeruginosa piperacillin/tazobactam (TZP) Code example: Combination Antibiogram Case example: Additional susceptibility Pseudomonas aeruginosa TZP + tobramycin versus TZP alone Code example: Syndromic Antibiogram Case example: Susceptibility Pseudomonas aeruginosa TZP among respiratory specimens (obtained among ICU patients ) Code example: Weighted-Incidence Syndromic Combination Antibiogram (WISCA) Case example: Susceptibility Pseudomonas aeruginosa TZP among respiratory specimens (obtained among ICU patients ) male patients age >=65 years heart failure Code example: Note combination antibiograms, important realise susceptibility can calculated two ways, can set only_all_tested argument (default FALSE). See example two antibiotics, Drug Drug B, antibiogram() works calculate %SI:","code":"antibiogram(your_data, antibiotics = \"TZP\") antibiogram(your_data, antibiotics = c(\"TZP\", \"TZP+TOB\", \"TZP+GEN\")) antibiogram(your_data, antibiotics = penicillins(), syndromic_group = \"ward\") library(dplyr) your_data %>% filter(ward == \"ICU\" & specimen_type == \"Respiratory\") %>% antibiogram(antibiotics = c(\"TZP\", \"TZP+TOB\", \"TZP+GEN\"), syndromic_group = ifelse(.$age >= 65 & .$gender == \"Male\" & .$condition == \"Heart Disease\", \"Study Group\", \"Control Group\")) -------------------------------------------------------------------- only_all_tested = FALSE only_all_tested = TRUE ----------------------- ----------------------- Drug A Drug B include as include as include as include as numerator denominator numerator denominator -------- -------- ---------- ----------- ---------- ----------- S or I S or I X X X X R S or I X X X X S or I X X - - S or I R X X X X R R - X - X R - - - - S or I X X - - R - - - - - - - - --------------------------------------------------------------------"},{"path":"https://msberends.github.io/AMR/reference/antibiogram.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Generate Antibiogram: Traditional, Combined, Syndromic, or Weighted-Incidence Syndromic Combination (WISCA) — antibiogram","text":"","code":"# example_isolates is a data set available in the AMR package. # run ?example_isolates for more info. example_isolates #> # A tibble: 2,000 × 46 #> date patient age gender ward mo PEN OXA FLC AMX #> #> 1 2002-01-02 A77334 65 F Clinical B_ESCHR_COLI R NA NA NA #> 2 2002-01-03 A77334 65 F Clinical B_ESCHR_COLI R NA NA NA #> 3 2002-01-07 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 4 2002-01-07 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 5 2002-01-13 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 6 2002-01-13 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 7 2002-01-14 462729 78 M Clinical B_STPHY_AURS R NA S R #> 8 2002-01-14 462729 78 M Clinical B_STPHY_AURS R NA S R #> 9 2002-01-16 067927 45 F ICU B_STPHY_EPDR R NA R NA #> 10 2002-01-17 858515 79 F ICU B_STPHY_EPDR R NA S NA #> # ℹ 1,990 more rows #> # ℹ 36 more variables: AMC , AMP , TZP , CZO