From e5d32cafe04bf136e5570dbb0da929fec179010f Mon Sep 17 00:00:00 2001
From: "Matthijs S. Berends" <m.s.berends@umcg.nl>
Date: Sun, 12 Aug 2018 17:44:06 +0200
Subject: [PATCH] new unit tests for ggplot, small fixes

---
 DESCRIPTION                      |   8 ++--
 R/classes.R                      |  11 +++--
 R/ggplot_rsi.R                   |  11 +++--
 R/n_rsi.R                        |   2 +-
 R/portion.R                      |  52 +++++++++++------------
 README.md                        |  68 +++++++++++++++++++------------
 man/figures/rsi_example.png      | Bin 12592 -> 0 bytes
 man/figures/rsi_example1.png     | Bin 0 -> 9651 bytes
 man/figures/rsi_example2.png     | Bin 0 -> 7874 bytes
 man/figures/rsi_example3.png     | Bin 0 -> 7160 bytes
 man/ggplot_rsi.Rd                |   6 +--
 tests/testthat/test-ggplot_rsi.R |  35 ++++++++++++++++
 tests/testthat/test-portion.R    |  17 +++++++-
 13 files changed, 144 insertions(+), 66 deletions(-)
 delete mode 100755 man/figures/rsi_example.png
 create mode 100644 man/figures/rsi_example1.png
 create mode 100644 man/figures/rsi_example2.png
 create mode 100644 man/figures/rsi_example3.png
 create mode 100644 tests/testthat/test-ggplot_rsi.R

diff --git a/DESCRIPTION b/DESCRIPTION
index d6d8fd22..9dbc4f8b 100755
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: AMR
 Version: 0.2.0.9023
-Date: 2018-08-11
+Date: 2018-08-12
 Title: Antimicrobial Resistance Analysis
 Authors@R: c(
     person(
@@ -57,14 +57,14 @@ Imports:
     Rcpp (>= 0.12.14),
     readr,
     rvest (>= 0.3.2),
-    tibble,
-    ggplot2
+    tibble
 Suggests:
     testthat (>= 1.0.2),
     covr (>= 3.0.1),
     rmarkdown,
     rstudioapi,
-    tidyr
+    tidyr,
+    ggplot2
 VignetteBuilder: knitr
 URL: https://github.com/msberends/AMR
 BugReports: https://github.com/msberends/AMR/issues
diff --git a/R/classes.R b/R/classes.R
index 4e2f7dba..0962e240 100755
--- a/R/classes.R
+++ b/R/classes.R
@@ -102,9 +102,14 @@ print.rsi <- function(x, ...) {
 #' @noRd
 summary.rsi <- function(object, ...) {
   x <- object
-  lst <- c('rsi', sum(is.na(x)), sum(x == "S"), sum(x %in% c("I", "R")), sum(x == "R"), sum(x == "I"))
-  names(lst) <- c("Mode", "<NA>", "Sum S", "Sum IR", "-Sum R", "-Sum I")
-  lst
+  c(
+    "Mode" = 'rsi',
+    "<NA>" = sum(is.na(x)),
+    "Sum S" = sum(x == "S", na.rm = TRUE),
+    "Sum IR" = sum(x %in% c("I", "R"), na.rm = TRUE),
+    "-Sum R" = sum(x == "R", na.rm = TRUE),
+    "-Sum I" = sum(x == "I", na.rm = TRUE)
+  )
 }
 
 #' @exportMethod plot.rsi
diff --git a/R/ggplot_rsi.R b/R/ggplot_rsi.R
index d18d4bfd..00cf1e4d 100644
--- a/R/ggplot_rsi.R
+++ b/R/ggplot_rsi.R
@@ -18,7 +18,7 @@
 
 #' AMR bar plots with \code{ggplot}
 #'
-#' Use these functions to create bar plots for antimicrobial resistance analysis. All functions rely on internal \code{\link{ggplot}} functions.
+#' Use these functions to create bar plots for antimicrobial resistance analysis. All functions rely on internal \code{\link[ggplot2]{ggplot}} functions.
 #' @param data a \code{data.frame} with column(s) of class \code{"rsi"} (see \code{\link{as.rsi}})
 #' @param position position adjustment of bars, either \code{"stack"} (default) or \code{"dodge"}
 #' @param x parameter to show on x axis, either \code{"Antibiotic"} (default) or \code{"Interpretation"}
@@ -28,13 +28,13 @@
 #' \strong{The functions}\cr
 #' \code{geom_rsi} will take any variable from the data that has an \code{rsi} class (created with \code{\link{as.rsi}}) using \code{\link{portion_df}} and will plot bars with the percentage R, I and S. The default behaviour is to have the bars stacked and to have the different antibiotics on the x axis.
 #'
-#' \code{facet_rsi} creates 2d plots (at default based on S/I/R) using \code{\link{facet_wrap}}.
+#' \code{facet_rsi} creates 2d plots (at default based on S/I/R) using \code{\link[ggplot2]{facet_wrap}}.
 #'
 #' \code{scale_y_percent} transforms the y axis to a 0 to 100% range.
 #'
 #' \code{scale_rsi_colours} sets colours to the bars: green for S, yellow for I and red for R.
 #'
-#' \code{theme_rsi} is a \code{\link{theme}} with minimal distraction.
+#' \code{theme_rsi} is a \code{\link[ggplot2]{theme}} with minimal distraction.
 #'
 #' \code{ggplot_rsi} is a wrapper around all above functions that uses data as first input. This makes it possible to use this function after a pipe (\code{\%>\%}). See Examples.
 #' @rdname ggplot_rsi
@@ -67,6 +67,11 @@
 ggplot_rsi <- function(data,
                        x = "Antibiotic",
                        facet = NULL) {
+
+  if (!"ggplot2" %in% rownames(installed.packages())) {
+    stop('this function requires the ggplot2 package.', call. = FALSE)
+  }
+
   p <- ggplot2::ggplot(data = data) +
     geom_rsi(x = x) +
     scale_y_percent() +
diff --git a/R/n_rsi.R b/R/n_rsi.R
index 72273e34..cac652cd 100644
--- a/R/n_rsi.R
+++ b/R/n_rsi.R
@@ -35,7 +35,7 @@
 #'             combination_n = n_rsi(cipr, gent))
 n_rsi <- function(ab1, ab2 = NULL) {
   if (NCOL(ab1) > 1) {
-    stop('`ab` must be a vector of antimicrobial interpretations', call. = FALSE)
+    stop('`ab1` must be a vector of antimicrobial interpretations', call. = FALSE)
   }
   if (!is.rsi(ab1)) {
     ab1 <- as.rsi(ab1)
diff --git a/R/portion.R b/R/portion.R
index a7e874f4..71ed9d9c 100755
--- a/R/portion.R
+++ b/R/portion.R
@@ -176,6 +176,32 @@ portion_S <- function(ab1,
            as_percent = as_percent)
 }
 
+#' @rdname portion
+#' @importFrom dplyr bind_cols summarise_if mutate
+#' @export
+portion_df <- function(data, translate = getOption("get_antibiotic_names", TRUE)) {
+  resS <- bind_cols(data.frame(Interpretation = "S", stringsAsFactors = FALSE),
+                    summarise_if(.tbl = data,
+                                 .predicate = is.rsi,
+                                 .funs = portion_S))
+  resI <- bind_cols(data.frame(Interpretation = "I", stringsAsFactors = FALSE),
+                    summarise_if(.tbl = data,
+                                 .predicate = is.rsi,
+                                 .funs = portion_I))
+  resR <- bind_cols(data.frame(Interpretation = "R", stringsAsFactors = FALSE),
+                    summarise_if(.tbl = data,
+                                 .predicate = is.rsi,
+                                 .funs = portion_R))
+
+  res <- bind_rows(resS, resI, resR) %>%
+    mutate(Interpretation = factor(Interpretation, levels = c("R", "I", "S"), ordered = TRUE)) %>%
+    tidyr::gather(Antibiotic, Percentage, -Interpretation)
+  if (translate == TRUE) {
+    res <- res %>% mutate(Antibiotic = abname(Antibiotic, from = "guess", to = "official"))
+  }
+  res
+}
+
 rsi_calc <- function(type,
                      ab1,
                      ab2,
@@ -244,29 +270,3 @@ rsi_calc <- function(type,
     found / total
   }
 }
-
-#' @rdname portion
-#' @importFrom dplyr bind_cols summarise_if mutate
-#' @export
-portion_df <- function(data, translate = getOption("get_antibiotic_names", TRUE)) {
-  resS <- bind_cols(data.frame(Interpretation = "S", stringsAsFactors = FALSE),
-                    summarise_if(.tbl = data,
-                                 .predicate = is.rsi,
-                                 .funs = portion_S))
-  resI <- bind_cols(data.frame(Interpretation = "I", stringsAsFactors = FALSE),
-                    summarise_if(.tbl = data,
-                                 .predicate = is.rsi,
-                                 .funs = portion_I))
-  resR <- bind_cols(data.frame(Interpretation = "R", stringsAsFactors = FALSE),
-                    summarise_if(.tbl = data,
-                                 .predicate = is.rsi,
-                                 .funs = portion_R))
-
-  res <- bind_rows(resS, resI, resR) %>%
-    mutate(Interpretation = factor(Interpretation, levels = c("R", "I", "S"), ordered = TRUE)) %>%
-    tidyr::gather(Antibiotic, Percentage, -Interpretation)
-  if (translate == TRUE) {
-    res <- res %>% mutate(Antibiotic = abname(Antibiotic, from = "guess", to = "official"))
-  }
-  res
-}
diff --git a/README.md b/README.md
index a8a55fa6..3bbdcace 100755
--- a/README.md
+++ b/README.md
@@ -136,43 +136,61 @@ guess_bactid("VRSA") # Vancomycin Resistant S. aureus
 This package contains two new S3 classes: `mic` for MIC values (e.g. from Vitek or Phoenix) and `rsi` for antimicrobial drug interpretations (i.e. S, I and R). Both are actually ordered factors under the hood (an MIC of `2` being higher than `<=1` but lower than `>=32`, and for class `rsi` factors are ordered as `S < I < R`). 
 Both classes have extensions for existing generic functions like `print`, `summary` and `plot`.
 
-```r
-# Transform values to new classes
-mic_data <- as.mic(c(">=32", "1.0", "8", "<=0.128", "8", "16", "16"))
-rsi_data <- as.rsi(c(rep("S", 474), rep("I", 36), rep("R", 370)))
-```
 These functions also try to coerce valid values.
 
-Quick overviews with `summary`:
+#### RSI
+The `septic_patients` data set comes with antimicrobial results of more than 40 different drugs. For example, columns `amox` and `cipr` contain results of amoxicillin and ciprofloxacin, respectively.
 ```r
-summary(rsi_data)
-#  Mode  :rsi  
-#  <NA>  :0    
-#  Sum S :474  
-#  Sum IR:406  
-#  -Sum R:370  
-#  -Sum I:36
+summary(septic_patients[, c("amox", "cipr")])
+#      amox          cipr     
+#  Mode  :rsi    Mode  :rsi   
+#  <NA>  :1002   <NA>  :596   
+#  Sum S :336    Sum S :1108  
+#  Sum IR:662    Sum IR:296   
+#  -Sum R:659    -Sum R:227   
+#  -Sum I:3      -Sum I:69  
+```
+
+You can use the `plot` function from base R:
+```r
+plot(septic_patients$cipr)
+```
+
+![example_1_rsi](man/figures/rsi_example1.png)
+
+Or use the `ggplot2` and `dplyr` packages to create more appealing plots:
+```r
+septic_patients %>%
+  select(amox, cipr) %>%
+  ggplot_rsi()
+```
+
+![example_2_rsi](man/figures/rsi_example2.png)
+
+```r
+septic_patients %>%
+    select(amox, cipr) %>%
+    ggplot_rsi(x = "Interpretation", facet = "Antibiotic")
+```
+
+![example_3_rsi](man/figures/rsi_example3.png)
+
+
+#### MIC
+
+```r
+# Transform values to new class
+mic_data <- as.mic(c(">=32", "1.0", "8", "<=0.128", "8", "16", "16"))
 
 summary(mic_data)
 #  Mode:mic      
 #  <NA>:0        
 #  Min.:<=0.128  
 #  Max.:>=32 
-```
 
-A plot of `rsi_data`:
-```r
-plot(rsi_data)
-```
-
-![example1](man/figures/rsi_example.png)
-
-A plot of `mic_data` (defaults to bar plot):
-```r
 plot(mic_data)
 ```
-
-![example2](man/figures/mic_example.png)
+![example_mic](man/figures/mic_example.png)
 
 Other epidemiological functions:
 
diff --git a/man/figures/rsi_example.png b/man/figures/rsi_example.png
deleted file mode 100755
index b1500c2e49e268c468e2e9795e8b3faf8f0df404..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 12592
zcmd6NcT`hpyYI$^ih_=Sh!h>ANEMOZ#DYqdF1;y84WWk?97P;jMrlHTr~x8Giu4kK
zLlrPY2oNBUD4|y&kN~0F&CGoFyZ5a7$64#Fb<g=Di@ZB~Z}#)P&wig@d47LjeD~yW
z!Q%h`oV<VUwiy5%Kmq_q_t7KZk+?N;ZtxFBpxND<KuM3#EI2sip=+cI0G|+C+b)N}
z@iG5<Hh}=Z-Ms(d;FHvk0f4J*_iyW3ggGvbHQJ6Wh;sf!?YTYWZqz)Wd0?Thv+8_u
zdC>e&F;{ZE7S=qY`Sh35VQRcRvGAo+mv4lbUvkw!O?R2UynpXzsyxfiSUeK|Vg$EK
zH_ifpV@}yI;D+yWBLMKa<Z}RcaNgw%5HoTVeB^$;p+mr1<z{gJs5-680endMZ#~ro
z4LwmJPhVKLY;5dKvdMOE>~OZ1SvBeEF;M4ht619(4@vGRsqip=);+A&AS2sgJfK`u
zqfv$9DfRKUS%}CA8QvpT&ki|J?C^u)vs#4crraL(eP)9dE3=d}#9|9S#7ub$dx`rr
z+s@wlaBg~xBEX$sQ|EF5(g3AuCXDEC($Fgx<Vdi(C>eL5&*v9fzO#H+B<qrB+OzX2
z8)=Q0;brTU>#RB(&5gT_aWs{%i!Nsd1)SgsOU)Q+xyxv-==;|BvLj7Bm%0xc2ttE(
ze($`DTq-U1Abw&i7qNaQopinE?e3=&G!mJ_TApl&i^ktAo=GZT2f{5`R2;(13A5`h
zJSXAPEQszHcG|EC+CtC9nP%zj%v8mT%BF=CmtH&HI{DG|XL#YYFrwGiinHftg+;S%
z)z3=(0BFOrG_U5^s2&4T{rT*CSle%zl?$8ZM;nqD!{M+iKYenPEMFN=9JSiViWIS|
ziyGC#^_Fdfe5sWjeKcoNg{;JI1nmlSfkI1jl4&PbxTz=1uk@)AduV->obx;-6@jKQ
z=7t{L={eFlD6M2@&Aj(c9gVg1@tPKO71miYJm2Q6&{#^qKtfw)lrg1MOW88ntXPIx
zle<&+0^7#xz#mlNNeRykZ|9w5N<X9Il4{IC%b5?~mRm3H#NQ4Q-*Kg;m_lCP`D1&!
znEK&OooCK&p82Qk4`iwI(kv^dQ2fq7A#QASzSyw1ZP9OJD#**o9a-Y2><hWxd0=A=
zhRx8TR7S$2mJJ4!t)`19aa2;>Xb{y3ZZ+w-xt8TkozEmJyJSsG%afNaWE7EI-6bOI
z?Wz!s_Oa5^ITnOAU(wa>o9c+&igHj774xoX2q-BkSnR-~`jiKa^U0$tf)7KBrzyjm
z>c6#|QFcjc<G~K*k>*Xd&;<CO<hZP&GMnSib9l)GdigxBV*M+wHa(2C;7((Z?S~MU
zlMt>i3Edg2ON%tabp<Q9v%HT6vI$19qzqM$^=YRjmsU1pwlTl}BdXLAW$Y!cd?6-(
zx(<)V8&Kv_t1UOF7PuwHj#}cnWpE<<W0RgIynDK2YczEl!efvy-5BIM*@TwYVrt1c
z5J|p?<-_?=k@Z%rx+8rB8DtYC<%m#L*)i|5-M8{4HXmjlqx`zzJy;!q*ONes&Bg`e
zi;o#`E%m+_o-r%tQ$?qy@v7qto|#SMuBSh*V0eoh1LxzMW)DW-7reqAj)8ApwQndv
zJC!Go_fR^bw2!pPh@$fT*7eTAg+0j%C2%3i<+*koJ^G}aJWRU%)I=%C(9(JAcoh?;
zjkRbpF}KvA1m$KFYb|#gAfim>5f!}Q(mlH6n}S1;+#c1w2f2-qMi0G5pw6i~Rm9zw
z@`6-77VaFO0>&<ERKmBm`7{hOdyoXr25dajr%-DMx7<1{$s+pNk9F)2EGy67A-`Cm
zy^&R4;HStGz>u|Z8x4IW*fUO)ww9X=h7Tt;oz~$QX*&m738&<Gwq^_;ldUR^oU^3P
zwXyQ)*lSL^zv6k>qZ<LINwPfjC1QP!uWU%RUu0KueV1xULP!o_3fA_76S)9QjBs$`
zbFANe;jIj7WWi@~67p15RBCbsLy2GhQ%EGCZ>lK6{-9#%r;}mlo5g*4pU_&PsCL1;
z=zg^i8L`B1k4?^DIW1#-4u3hWioJEW!M%z5%-^9DErdID<=q?WkS)(Z#BeWt{F{RO
zw~gxw1}$Mid^{-huRi@QBF8<1Yo;2z?9MtVqiP;S3hQ_XUwbjOpe6ZRcm+GiI4%_K
zc*@GmBh@{~V=TcgXCZuUnr!)4ys(U|pv~9{U+r7JyVBoLGM*Jx<&Et#U`2`M(?>hI
zOP`X5kKlw|&J;a*=QfB9qDJ6OC+TUFl$8He%e)O|T1OelPv$(LG~Cw>;dec&Efm+g
zW7Dg_{INV$M%v!yeW9^~D>Q5i5^YG7G#81Vy3-=eEn0z?P<bHLtAC3xsT8igF4gpg
z&)o0oZeWep2ak(u<Urcb6;v&9@)=c9oTW0n-rrM#oyFxW1SS&e@Li6T(TY>C8CdcH
z*1>Kbw>A`7y{T%AYghUyNQqIo>~k~CjBlolEh|^?*=G0D{aeYJX)$RZiv>rmMpU^z
z-8tsP99CWD?GsdxKWcm2Z}Y*ud>1QyHl8;02likFF%E89e1v9_Ly?oVR@_U3n4rEe
zq*m73S3Ih>`H3Qd6+yz;={)rot79WsxLjzWU1_OoZq@#pn-keiKk_^90|j(@;nH=>
zlBd>n$=?@`A|DhA(;_OqzbT?-F{B23ymE<HOKXa!SblXBUDw`PjMP{A(ql5@5#3a?
zM@J#I2t_eTYVxB`Px@>^@rytFM_vTkqMf`X(qms?6jVv785M316?2`n@9mJrGZ+fH
z`W+ji)gcW&5!OiNyHR|=;0xaxsWq6VY0#F<)%!p4G#~3eW|k9W336%bv>yGzzV(%2
z3!|o6IjhDUwbScW<N5O*BeaONXlgrbHp{^6Q4e*`YMD)4-m@w#TV_k{`4zY;k;n#g
zzg&;N^aa>S3&;7=AmKleox3j&r0F+Wkc|zrZg)32lMN?D#uIL#WlWIa{?>+EjV9p1
z3=FtMkGgBi3cJDHWMEL%teG(gxYs=-nb!u0x~J1>@^_rv+?J<q4qTg$cp%T+q*&FW
zWLdSKFF!k$3}ww@x2l+p+SpX$a93rJ{XB89$hRB6HpHr8SK*0De+)C()qb|2r&~BA
zho4uOLDo!zA!|N*5=viCqcZ7T*GtBu`|W1le4AT)-&F7;nmvUNU>93S*m+brJZ2ZX
z4~@jjqu-3Cvhsvg3X5&4FL5PIvCB&u-I-Z1C3<Ai?5aV>(5}cZCv4?%iikF!aB)hq
zMDWm1`J#jGXxu^++f<#Rk(Ql~{bn68N&ITSYapf{k?CkY8JHcqG~+?GNf!ywNXuFZ
zNn|@ZbvxEazt%GlSSKyt-h7kbUENN(g>T0%qoWk#qnRSwjrQv83sHm+m(E=BjmfHA
zTk>gKQc9|7LdR#FC(#BtUS&w`*<C8qBFVO&6tb1_@XHU;wkH3?&v9=bo!NM3M(M{W
zm(LqqL5%!$@S@KF(c8ToMKBSQ7M{W0&55_Wv2L1WX4CH%Nj92EiAO$)foD9B=5pDo
zg)hE<EqCo-cK0n|9|xM2|Bpb6$=7|7>0A3qaL_12K{EXi@ZBGTwc{(M0#2wq6&i_F
z5MGvblIzmzKObjj)BGH#V$D7r22?9qB9qB7J_WxJ+(#XJcwWII8yjBs?pXGpRxuya
zEYIc!BCmG*aX0y+U=%N+<t534LU*H+JX?P%*c}97ZmE@*4JQ(dU&b}j$2Y_2O+2TI
ziNk^bAk_6LmY^3|_oS_ilGtse6>}tm6q0k+$*b@tC*Z7dm#q85>r{rsG}KPhu6J%=
z&?-<5X#{-#mV6N*<T<3N2|{kw>%z|gU;}fmUl@enLJmWQbdTI@!;TTdw)Ql0opDXB
z91k)jqWB>20CkBAN#6*pnOxGn&lcBfv^uf8`^Xyr{3Y-5E7RN`u;z+3+l~P-i(sU9
z!rP@ZONPN20H6rl3|;@A0en_a_`p`oFA5UN20mXDz9q6hJca~5zXQLk)xSdkfU3y9
zK|cWutzY3o_-OZBr>`^zz&-n=*<Hg=96eLjM>k%R7;tm$!61E3f@Ed#UIZFJ-}-oN
zj1+C?$9+b^ymQtxR~`7?qWVHJ0C!9Dq#S)>TUmYQzR<6b<NN>d7*g~NNGgOz6ednL
z_Wrk1jQ`97+!a*rfJTM<A|Wwf&@I4gZ~*XEy?a)<gY04R#;1L(O;SAN{(3(`m3VV{
zc)j)^7nnAx<z;HgI-xC3*7cz90m&&#u*82L$Nv~Z{XK18nx`a1;!a~w%7+Ly|DWkb
z>6Y$*i~eehs!#$~ClyS9q>X^yc$b=U8OU(pq#Sx!SZ{vZpc<r@aRJxvYk4L2T#z&A
zp;w1hzaK~d6TDbgaXnn8%bXL46xdJj>pmj6`xBp42FvR^U&Dpr22}Ihjw0`>^tfNE
zl{xu6SA}XXJQ3xF?2+w5R{I_RoEy+wNY4i>Mx0|}{5LX!A!DzMju?@T{WH=c^jZBu
zKxhce6xRnVGkFG+Y|PxVTEo|16ayAF{2qh?G0!y0%j*0dut*PcAMcM_{u#&pA7h&n
zbB;wKqT9x3Uv4xTJ9_C!RXzAN)Px2JKNqxS_3A9@Rm=oOixj$jlMLr6M1T9J60Rr$
zHyq#9?>_{nUZ(qKn;)+-76$z9?Z;U12xPK!aJ&1TNNOXV6Ljr!vs>?H^^~|}Pub_T
z4hM%H0sw3X7xG!R%%j~zBKm$}c(JSdZFYeuW~|^dxbfx-H&3aSz^QszmXp0wyp74o
ziWV2DWvNtv7ZePf-Q&vC;yYYgXEAAcnck7=vLYP7=l8<R?uyzqzl!WyS*P(t)ATB}
zJq-RieyJon$aQFkzuQO=%wK5kV?qxC>vIJF9tK_O<VDonYwRUn>P{B}N$8UFNWh@4
z5ZQw;!<~SUR_{fP-fjr9itCq#PlH+Y;F;T(OAYIK#45VD#+#rUjM>F!?ovitYM$UL
z@x;?fe}4<p*bMIk^JGxOg`p+TfogE63`f_}QUfR|o2QIg8gw?p-hz8+JGV>4^*8MF
z*>dIzLYb8)-yPA(Bl>sZgRyA%ZRKp6N}1Y!EXnOzUL!GZOZ8Jl^8_t=H%OS*UULGf
zW#BXgJz&Y0mRLEQ9T=$5AuUQLa8SLA>-}0Y60X0D3=F(7_x)Pdm`kM`enzZmAnWyf
z{uxbQX$|(~aMzbZmenPS_~eRVVuIb<o_u3xc53=^!}h#E4a7k;dArB1g6@<vxH#Vx
zvdVf=6UYrL6S!RhZG*=)ZO3XF>>YOaEy+n(HpZ;Ki05bI)~1?XO6OEaj$zVX`z`IF
z8-=avMk)-JMr^w~>-ts&YTX1)85N2)bM^{x!;kge8CB}$kmib`tsX8%A>N-Idy-Bf
zf42@YVVjO5CAE`jOxxXo$A4Zsd9eBd`&zZXF{4lv7IlL;SsdVmtxh4!J(}L81YF>Z
zVL{@`>gwj|?HzJ-qR&`r(G5tAP@Pca@I;sF_AU3aXB%nyzemt2Djc7{5suyhV@~t&
zr#W9C)}8Qeo^+|*2&u!w#J$$BGQP7;eFsRt%rE>Z-E1Lc1b<0oIi%2oM9EPXi)a_*
ztk1aSudEYmjcx^6N)Xm5_op?vCcI750GzqTbHec3CBgvn&7I(Na78>#Ru;~{Mf8*V
z_3wMOt(i_Zol?eoker*WudI9z<^SNW8{lWV8g&+4%}1Qf&(aAIwUUN68h@>rs`UIk
z>HO*U8xEEWuo$29_bR^4IO>KvUY3*fnn1}svWh9=2JX7Xm1TBf53-?IhkF9PtE(Y!
zUP^p@C*(*M=I8Z&9IdGZixbxpcw+-7O_AX??CP$*zWiCVfR_7@=^q}mqrS@%Y1TG8
zb`hzX=_lz5kiZo&Iop=n?i3G&JWFbO2+fR;ssxcf$XMJYZrEGUr+pKMWc#zyWE0p(
zfb1m4wf8KVWu%mQvEs%6AsK&0r^sxLn=-8S@!;o(*HNcL)1oEH^T?TxE{1ybBP-u4
z%_hZ9*X`z|2-{A3z9WY|$CNa$Y2LxATtVT3lI;>Cf{JLJ2o0k&9aJYRWq^S?jE@|{
zY_<h#Dtoa1aC9qG@e4nZz&2Na5VRY=K6(9k2t)2HCFBRHSy!PHHT*E?Zov09U{N0~
ze#97iQ*}7o`qotbot#N2vL`XjGMU2b>y$B4TisOPex_+ktV7X6*ijD(C4ZdJb$nvq
zhLEw5O})%-`Gu=|a*44>V7{|2vMD(rpXOR>s>GaY#<gQqWkvAmxaTV2D>S+F9q%mu
z%Nh`j>I;V#y9Z)et{r7{LrmY`5s@Q~B+c82ZX-IL!PkwR)i>Q>1o9{N(ODxZ!6V6`
z=4pOx_g{8;uGSOA1wFOf5jcquX(?(+p_v`w&q>&5$QJC>VyGWI(xHqHK@z=DI1&-}
zL9KWW;?Sx;I2xz{`1^z$M{aV$u^BuPsU!G&7rWKI52o@7+_SSrkK^LV<mA%A)`csF
z`RIOANx<22{X6-}&zmpv<lSL+E17-IPJdo-r9p4ixwd4Q?ON!bONbY~Csyp;q>9t~
z{W1Rr)3-uxH!wpu>WZO<IrUZKP?9?{?31DyRFjX6U#bxXRJFl8Z{z90wfZaHCan4^
zL*7P_5Iu`9<~vPUBVYTSvyE|#He*?sJS_i?9p6Q7v2Uqyw?=9W>o_+;U}d(f!_G93
z83RX)D*KWir_3d3mZv&(e!P71WL3QUW3HWv_&#z`(Jpy}8^Bp>)C@;%#Tba$=9?9?
zXC%Z<uS<`do1)iaUcc?<(&jsq7T&eDn-%rMI_{9ZyRUS{T(^1JvU$VM@HH+cNwX`D
zJTv+3R3rr49NReTT*tXpc__llRjc?nM{E7Y+J%NxPKK;lTSHj!BJlbA$VDHylXJP&
zcG|6fvL!93Ir(p>&KvoNUm<qq3iE1@&JC}z=f`WaHm_kwUSqS@2w~f%j^bPi9+Uo=
z_)}l?-zEh2R9san+C6Z!kxg_hoWdQ_u#^sc-{kUk7nile40Y_=5ix><lzD#N9y3I~
zu&87oiV0n~%;ix?+16uF_>3v;1yDJ7w8*5(rM=;?JBhoW8kl*fDq9NaMtQIMBfpVb
z#!AXV_L>zAl$C<@M<bSS`S4f9%-Q@GBp!S};$mMhEB@4${Ap!c#<l^GjzHKIY1L*^
z8MNYu9Fm(4TvF~dPer5@tzrh)qVkXe`OsOvdm|SQmd!kSPLLH1Dr^-i9b?z%tHk~T
z))z+4N9A}<b5frL{h4VvFAk_auM$JT&J#S*UyE%UJw4vR?{EE#hR`k_I<9Q`3(wCL
zT!&`)!sNiro}^2P+hzX)kusu}ea!9*O$clTb8b9LzUZ?w<X$0AsiE<9&NF=e-uInf
zu5+k9&(T-LEUN4gR<4upUt~8%g<+z83w~}CgP7Kudi5AoW;GFU@Uk=C#5;aY9RNfr
z3KOq#;<VMHW|ML_fQ{*eQ_9GdIzA}d4r>Si%nCcm+P0i-P46tc)s90O;S~gO4{Zx1
zfk1EAi^z!s#eZ;oZzwPRA1warEbY8701OUa@+p{Gb3MhXm>0mPWk^6na`qub(O9$l
z+=Wz{Xw)xgN9>z(z4+M<&zdX#AkgvzAxQE(s-|~t=PcH(zJ?ybuN^n$<Y#>hu*(zt
z(`NDe=$!03>*LUNGPK&FK~=?JFJPKyJSqb&OrNjlQ^;*=zBm<YTY)Mz=3>{EE%b$~
z+73;^ghs#$3o^3xXZ0%DIs6D9rmFmpApEsTgX7U4AhGK`Z?u{$Qg{S9<lEk9-rSYG
zIhf@0bYU7jb4Mj>a!qN2Z;Z6tbL&iUY13lWw)Z{Vieas!=#a=+0f`;bjIzPIvM9xz
zrmO4dl<2`Li{9G2^cGT`?XYnrU7WzK7W$%=_cx@7YZ`!1zK1fAq^)e(!+-^-wtMw`
zH5)?Jo-fc0-m{1ps;x&U+cfNG6)Szs@N8?#sJw9rZ$Ba0b>~gjZ!<)h$S9hU8uz$$
zs6b9BMD5{^H6^>oWyZx{C^R%#P>)648Xx$(mP%`j5rRG%C>7h>TMSLFPN~j0ym}Hj
zm9%(D`8*VS&ZzFKYBQM}d|@6f=Wo#VR$nCBwsB_iPOw-Kwifk~diAlb_r@vjF4)wn
z89|hfFghr(-E(kM)7w~M&yia{Cx$WIhH$#0;aBd4qroH;>GK=aLJtz0@%E&VZ0eN1
z;MYLXjcB&{>{P@jJI;5wF!o(dWs~=rjhi;BiS8i^$_%rpP`(2(Rjc1qqa1fb`B3pT
zl!>YrP(fIUSLE*rgIAY!^JcCm>)4EYIYf5;I)Km5Y75bDn&v3+kI-H`lwDL{UO2b7
znH7SVtVT<dmVCt+Mce&B&ICM9bg<i-*l|&5jw(J>Mv)mdi!sS~5&8;sjoYuLEE5^6
zXyap|AkX9k(g`DWR>9>C=qD64`2p|N0eIuf{$xhpxv`lxxThOC-K+gVfh^7HJU$Iu
zaNbF-{GLi;NAKF^Z|^UMC~$o(5Wc3{TPMM1ZyBm-9KR!z;U$jl@F^X%kvi-b?Fc`C
ztJ4Cha2W!32*S@D`YJT$em-KS&1*W~OM{6Kxhp?wQtD^7S#UsA5qA?8LlE`8x}?cO
z_$gRMGh(OHa@R9)Mesqr^~`v;hh*CGnySDd`*S7lHbPtzjerjiinM80l=*Zk(!z5f
z-Y6@n0U`U`4qtv_^Oy5r!NW;?7Na+<JkZp{dBXoFDq8aFpeH=+%z}ANs^3XDNX^e_
zc0Zmq&C@6_Nh>(MS(7sRc8T%h91!zpds^<rFnnz4W4{7hh7i#Ik#=hnttk4e{Ad<M
zqfaK?D;?WTaWapINbDM+lO;bz<|`xXx{8*@!&a3GEv>Na9dMe-^Md`CZmZpzrD{UA
zd3=_6Q~oM&;iSuT2vIes@2-2gq+6<1E}x@u$6FyiOp_a8v612>7j8PR8aU#^&k>w&
zF9msKvvdx1Mc+xfgW02@U%a#@dGZK;bk>s;5*7#+O#<C#OOFE}rdn5Fdaj{`?}pec
zJ*6glx=dBpLKs8jgl>GM%-aG%M%%lWebBtuqeaFyuv$GBsTZzaHOwv7ylAz;)cvPa
zceJcW{@GRY8sqdjKjAPi<YhaM>osU@!dM*WBQ@-SAy-4Tdh4`^WVNQn+fAbrzV{so
z(J>=mn1st)7}C^1lACus9Wn5F)RSZI1v6t2DKII_NLZW*J!30tN~kd(SW0=OnweL=
z?2zaB(_N5NBs1nXPZQdebeL!WC5|vrGvYeq<?daPOKd0A_PTvG6lcT_+NeTkT+$33
zv8Y8wSk4XXF&M6xDV_y^)i2xwGxS<*#n#h6<e?<&-g}n+ZdXlSr2G+RX3=1%uwT5&
zr(2y<O_oB9(qj8sNmRDiTIAc|R24|jAt~_2qUKv{Zc8eCdY7|M&*V;@C+`SFmB=kA
z+&@>j!+cokZCmBkC=V8u9S-1a7f7mvRcu_13t>BJA!50k4*bSh&`^b<A%yGPU6I5p
z#lwK3fz&rwIJhjVM?&zWf@d^8T=z|FZfofmp>X|!Uk(622bo8ah3E&SMh>GBI%A4E
zz6LyZUoGiZbuVHBQO0DNC_mtAQpil*S>#tX=X{rfT7g~fG&teW=fGv$Ku|f6jfCXb
zC^trYuD1)}0Dkg^97R^JGj^?^n71L3CIbUwPLTeK9xcU+-R1tHhBdvo4cBa0fLY$h
z^ZAjKc&^7^ghsuzB>EP}uNF=cBLE}p8nz&4I%eVh&``;~aSiT|D%K7l&lP!FA8x?2
z37_D%1c2o+F#gllI_e4w(Fv0UnVtgguOIT5?&vrknwFpcmpZl>!uf7=tvby!j8$i~
zrpMy*RpS8-6Ayx4HvZ3ICG+p2oBbR>|95cZ_WA3F<2YoHS3_IDze?uV3+-@~5h!VS
z0@3$*X?8-+5bc=vG*f+NSd$e1lG3~+m%woEFEEX^q?VJ7X-p&jUorHEo+BBY!1pNI
z^UCJmd_pdTqhW<VS~fL~W#56U?Ru(YNq9&LX8>S&4T4lk>x;p`I+c>~CP!&MIOsaI
zGL7odrnhT02Wtq20Ozyb&F<?rmX=qd>6;gpNv-E3XiM>R0D$AXJE;34soiqaM6Wl&
zg((WoM_b2HWSyS<m<c6gP^%qOb{jTt>>K~Azde8Nr3-7(OZ^PTx@V^V|0nxu%0fum
za@@bPx`CPhaqR*Bw{{np=Mo5t6%x>Za=-dS-A|W6VJtAnB``87J+>3neiL1<Wm=f*
zyky7(<vJh|sz6ok`QNgt%D=Q@$5^HjD8L9%cRB|F;QIqmO-&w$3P86!Z+(-wuAQTI
z;vZ20ypFIj8XvaN(1;7zW?g*0uJ*sx;+OJgiv>cE(7Hm-Z~syp{?=+GQu@(TGbhGZ
zCp7&MJfv+KK>|J?c-_2|h1LuPKF2{%z$WgsO#S#vnLUQQE7b^X@Sg+KMD-VCupqjR
zIO&t}^s`#)Dq6N21NY|`>E9~dmD|b<pqwhRF?p%BqqAM)7O1K;gW?>>#rN7G4Ot<h
zhI+jbs)&>d?`O3K0o)}}HzuE)qrSwl@2{CO^zy250MsNfEF@<zgvR^SJWU<=+~EJn
zB~bYTyR=FeKBBTSGnzZSulo5v1l8>i&|pzC^;UN2Z0v6TiwaB4!@zgv#EMPBj;ken
zz|T+TJiyd~gZ}b=wbSHtoMNr!?<pW0Jxi<+bm|!He0(AvW3swAT-{e$-Uziv<+M(A
z>CZ2Kat-JN7_oFvJ>&Ird8r6RBQ6cKm)*hh91Ayw4r!-@z5tCM5&Al|&|T$4TSc9W
zQ1rasnw~;t;$M=>t)agp7br<C7ovdQttPLQUOY_L7i?1QhnyJh*GcjVyrSwpJwIXB
z3a@x9zCRxs2Rt<V`od|<qsxF;v-r=akN$+B@-lxL?1WjADD0jIZLzYC^n7F@xH9>k
zZ{=srQG1(HJDaB(mm9=8bhqcY(Tw$zDefy-#Z;Dz9RJ8ah*0KU0lhY9uq|WZsK>}k
zV}|hj;<67{Y^5WAHOlJX(6h2*kplF26qGV|r_}uG(B;DfnK#vDosZ!eeER#i52|GH
z1bH&r<7xF7#>Y$$a;0sp(QMY&=dC`~(A}^7&NkAhJN022)K?WUN59_kI$aHcON||3
zN`NL%%ZB=@?ebp9xVeHinrY<di&IVn-Ss=(0-<m-oZ!BMXL_I(Qs~~5RdYJ0H2<e^
zaKx^nHlkpFjCSOo>IwBRue-WP&AnFIMQolV;SnNk0+^M;PXggYDdnM!AVCFO!gdg_
zpbM@@j19YMFGSvVwCJ?;-ukcg1eq8m^!dgOdEZD-cM-VVl~qh6jBnawE1Qz!&-Qk`
zr_rBmMe5R}mQW)WawaSaOAI3C{|f<U8?~DSS~J6Y5c>8+Xk^p~RZ?!L+|5zW-8b5?
z2@ae;iu~N$ZLr|(@>4187IkO9tkAvgdf9sk_G^*}mR)_dbmfY1LR@KyWLwnn*IZij
zAiI#Y%5>V9k*OXwjNggaz&<%MgeW^ERAaWjDSek<nQ!(lq2e#x7CH97fZp_~i-#8X
z>peM(IE>fY_zy$Jc-Z--;G+i|pMqCbt)LnN-1mZ!4COdj+?tC@)-atsA5`%WR9FNd
za3hW-1u5@0PUzOvv=N^v{ed4ajjvU<rlt&@G8LCIs$QgAQU&t`>B_P_*pMI@vF$~-
zT^h3r>`f>_p-Ao(M?XlrXe{$utR+I=1QE+rXYF7isa%OW$9ydVffyr@NEW9@ysdGw
zeWNDZ{WStvqkT76_xpG`M4ambjMN?bBOXLsKXTpOAjTd2qd3EhJa6?ip6(n&e%86b
z^}s^{tMxJT`d^&hon25%>bYN3y?k3|`G$pTsk_sXKjX-VO3Dn%7<tAJJBnP{aD&!j
zaE~@R=*`zc&!W>VzFW<mFHU$U`1*`9HsyX%oC{<YQEtw;Aii#|gP3V@t=dJZGk4+;
z^8%|c%Qg+0#mf-e+LH(!*MLM%*(A~?BQQ}pm3koVqEECl_#?_9O{chBy9Af1uSCj8
z4iV542@7nV;uHs!|2Xl#&e1g48ha|WThP1<T5dB*+<)&w)WiAIp9<gx`&X;mcZYB8
zXSButYJZza+BFPE9|YnGRMY`r@A`l2iTlR{daV*)oGa({K;{_m`3LCP2|1x2yOf^;
zgT%2cC;#gvkQ;o!l9BL_cjwbSZPnUBTbxY$m+{EPe#~;p0)U|<(2@gbW~D1ge<uq>
zZanLL5)#*2zwx>lWb5yQW*4P6p|AjvnzjnQV#ma#Qm>O)g;Ea3B>px^M$hkL{v`qI
z+mU16?t;~V=RXGgZ$B{U9U`|=R+=M*^}X?~no9~#_rnjk&=OOs23fLHR|%Z1|75e^
zkzh`+yHNkX+*p{x^9b}-C7|V39oz$d<eam2Cr{bfBNC?l{Xu+fx64*T_1J&J#guh+
z`~eR`-CdJV$oz_lo>uOO`fzZefaJ7sKFegs_(TJtc8!qZ@r&~2sUuAghR$o1MOVFd
zpa=h#HomKWK+}L7!3@1V&a7!dYa+@zIj<%+@wBz@CDqmw@tN|N2S7;wOB2kn4iCJd
zsk%%DoSO?=sXFSz<RylxDe>UXu6X%1466L^-w^Mj&ZL;i+g|b4R`MdN9sQmUYWwds
z5@#YUZ?!UxA|r)B77zF8&mHk9=+V6rG<V@pSxT(@0;w=lj9?5e?-LgZ5%XEq<?~gl
z_&&UR0!^{#I#y;`rxIQyD>Zd23-`RNbLG5KhYfQv!Xc3Gmn9dsfekw299mfYLC2Q|
zl5OO8b7*2?zg=q3Ma`I=lD-JAJ!!YQkS$9pgjZq8U;uCX&+jaRL!TsHA7XF%R)*7V
z$)5VWyI=oXuh!Gjck#{1&yBQS^q-G<y`w7>g0}G+jov-g*n4@Yl-K)!rilg~V>(vR
zR7&H8teyZN&v5~;&lGG+_@C5bDs@uMWPD=ca!*M|TecSHQ3pKGOxCq`r)25ClgUHv
z-1}P9dT?%q{Q1Goo`KFzn}x`tvptI%>)BC}>)uDgc%B<MTk_XNQ>)*TJlo3H((ECc
zH9v0;B%$Lz=#;0_dY&~wO)6Ed*mBp&cTY35!h5s_p#)$VECnA0_UaNqLx#M3wl|)#
zQ5<wE228=HVxYmV0_H+e2;IQG@>f6WzNwG#fm@IwKNEAu$xQ28D38-hcFUSbcu0;f
zv^o+WP*@ZW-NAL?H-?*)0}UvWiKWEcbOHu*xNp@18kJM@fkc$6P>aQlgyW4Ln?G$z
zV1KGLT*rUS68Ty%iCT8AUTh0YfnZlhjb%}B13M1tXE&?MVb6R!Y@)cyc<PHN`-o^M
za-ntylGbo?np-vd*bSQ|#owP8g`;OlrlnOsBa3uqAHlT;<yCfI7}TC*#ad)P^3l0!
z*PZL7DY6c^DiL@bp(qL|N9yw6We?36<wj>ubkx*C6p*H$AQdpuqBxe_x13u2%BBV~
z8JsJV&2WWV^$^C1-4q*Gu@bw%UoJPfqOv%6+I?1NX)Am(YIJY=qo3P+1bA!aD93)g
z%fqW1QOGD|r`C-2H5h&FZjdVoRl3CPnHN*xP$j!1gQ5V@_I==R7CVoD(cA_NA?gqZ
zrQHkB!Es!aH&u)V6!8o9q`_BI5`Amhx4P#;KmS!b34@zZcHp~bJ<TEKtyi2N=IjC8
zZM^TDWj5~VkHDV0!%NfE<C>4y5h!8O>aJXFLSbxRFaCF#?8g{aN5M>)k{b3rp|O1v
zMvY+gu{griT+X~Xk4e_;{_$Sese;HTx&YlO<aH24jaP=%Vq@32EE!V!ZXuIcJ$W`-
zd-qIX?B_0Z0#W+uBUolTCMbw!B~&TU!|Y2VDz4OGfU@9j7m1H0(}MiKK3!F6e8O<e
zv&`)`{LjJ%>WVabemJBy*~!RMbbt@hya`*1U?tA~66OFx)JG!<l83veZZ_0<O8j6A
znSS6yA=_Fq9v-{QM<LuIgmd`{5HD?^sqK?aWYC}4WT%;E#MVvcHtjf-Rx>UcHaCDU
zG1jKzZM0W+<8@G$l>IbfI`hr3FuYPyEC)NZGy{!#1LE|WrK?K9GCgFYuhNOMo6_Ov
zRkjjp?wtA)A{oor{)h|x904|;n?G~v6_%^UBE36Uj^6pcqq=-zwmbvtBLnf&b?NyD
z-GEQ25Lo64O{On(q+IR`dAXImT1eqFb~0ysQwx}|cM(;xO999+J|e7<tf~_M-k1HW
z&w)DMi*AFx?~Bmi)FTf6*=Tg|XBoSgH)g+$!%T4Zi4ACj=OxP!`utU&%{q_gh%%I8
zGh-LI<<Ob>Q_=RtL1mrAmDV!VnT7L0Uu!h9cxt8-LN8=-D^{wOcCA8<P1>^PUl(`7
zP1v;^n>KKK2z)t)b!X_%s#j0H!V|YIPwek7ZBS7|%K>-a*<jvJHLE6+7Vk4A$J7y?
z{T>QLXF)xkX_N;ITU`-i1ozP2so2_JNDVI=hC?cDq!h)F{i<h*wP{nt4b80EbCQh$
zFPvKD>KwFq3UhoVwkD0$ycdGAV;4(nDiInvDB+S!VVcJGD7DPIBiDg#yS^K9pop<_
vN{++vX6L_~`9T)|@IwFp+lrNb-_TJCEG~(n&w_j!0Pf#0zFl(D`M3WCEBEsr

diff --git a/man/figures/rsi_example1.png b/man/figures/rsi_example1.png
new file mode 100644
index 0000000000000000000000000000000000000000..73fdc0c71b079d39b35d3855d5a7f9b08e899534
GIT binary patch
literal 9651
zcmc(FWmKF^w`BvtEkJ<a!7WJe;4Y23TOhd8NN@rKcPGIuI3Y-A+zIXkNJEeW3*C65
zVH)zi-@P+y)|xv%?#vHXQ+8_Csj7W;RXx!fY6`em<X8Xz09Q#-Rto??h5`Ucy%;En
zoJhHL7DUFahKi0HB8h~Agp7=gf`ans(IZq;R5Ua+baZqK3=B+6Oe`!cY;0^C92{I+
zTs%BHe0==Jj~^2d5D*d)5)lzSdGdsqn3#lwgp`!@>C>lVWMt&z<P;PXl$4ZIR8-W|
z)HF0S&z?P_rKP2#qob#%XJBApWMpJwVq#`yW?^ArWo2b!V`FD$fByVA2L}fyCnpyd
z7dJQeix)3=czAeudHMMG`1$z-1Ox;H1%-r!UcP)OEG#S{A|fg(Dkdf-E-o%1At5O#
zDJ3N(EiEl0BO@y-D<>x>FE6j4prELzsHCI>1Ok<nl~q(!R8>{g)YR0~)ipFUG&MD~
zw6wIfwRLoKbai$0^z`)g^$iRR3=Iv9jEszpjZI8UOifMA%*@Qq%`Ge}EG;cxy?SM3
zWo2z`{rdH58yg#2TU$FjJ9~S32L}g7M@J_oCue787Z(>-S64STH+Oe;4-XGdPfrjC
z<mKh%?d|R3<Kye=`{vCXKR-WzfB%4hfWW}OprD}O;NXyukkHW3u&}W3@bHL;h<ESa
zMMg$OMMXtNN5{m(#Ky+P#l^+P$0sBtBqk;%B_$;%C#R&Oq^72(rKP2(r)OkjWM*b&
zWo2b&XXoVPynp}x!-o&Kxw(0HdHMPI1qB6#g@r{$Ma9L%B_$=LrKM$MW##4NA3uJq
zsHmu{tgNc4s;;iCsi~>0t*xu8tFN#B^yyPWL&N9KpBoz+o0^)Mo10r&TE2Yw(%RYz
z27}w$+S=ROJ32Z#J3G6&y1Ki&dwP0$dwcu(`uh9(2L=WP2M33ShK7fSzkdBXGBPqc
zIyyEsHa<Q+F)=YYIXN{oH9b8&Gcz+gJ3BWwH$OkWu&}VWxcKecx22_}<>lp-m6h+`
zze6C9)z#ItwY48Veyp#rZ)|LAZf<UEZT<ZD^VhFm+uPeaJ3G6(yL)?k`}_O9fB!x>
zIDkT-Fc|Fc@bKv9==k{f<mBY^^z`iP?EL&34u@Y{TwGpWUR_;XUtizc+}z&Y{`vFg
z?(WW@#xMi10WjSZ4Ltz>yuSNCq#@Uj)&RgWfRe0~j_-$qTuk35y5M)}287S3q*I?&
zHH10!k)O*+fgWQ{c*)lZ^X>`7=qr#A6DE@pXj;#?xhFZN$)z#=e*f_qIr3ut%aCt-
zj~{;pAk#LWHk`e2InaevKuc!kceP)xU5o1I>)&qIbyb=Pmsho^6JP-ZlN3n-l0T6k
zND7iZ-2_uqZjV6dl>gmHi(KqF1@K#}4FgDB$(gIs{(|P`6B+gSyvMJ#975eIFikNo
zVGY{GW~GN-Yev?lWH0nQmuH|uEngdo1)RYi?$Yr#A6(_E!DcwL%5@jMM$J%<G)#qv
zgCR5Q=FlV%jJ1TT4*E{)lkCy~@ZdI{xdi3TEp$?AhgvO2T@iN4Zf-9-tW9g3JUM;b
z>KUk>*;z2hMq@^(eqGrhXc3bp%)+nK;VJjJ|H9?t?+*(r>B|wU;W$jinAyKR&s~+#
zhiZFuW%oxzvS1kYn(LIPc_iV!oKxeYnb&!`fuY>x^67&sIed&;ofhvgirN~JpKbYt
z2cD$CTr(&03<`dGGx@Kj%n7R;W9)7b*?G+-_f&v~i1ceOz64zYb=l-<&idZ>H@0(W
zf8bD~(14Y#wD?9WRr}`*oSE*W_f4$je_L;LJYAci`MKESGO7-rgYg{+pY4yaXyOwd
z2lrNrL>XVx=gAFkb(So=W)p0^((%`yB<{UgOl<a?x4-RHV@y?7?htM@k!d0vfZ>1y
z$#XcL<S;qwzxm1+p5Dql2w9j8!$rMnKKWYxJD%L~<Fkw}iG(=<&hIA6=K4yVWwW~e
z$bRo~7TIk6p#y<2_21fEVP2+@d*##n>=rLzY`=;Um580euT=I7PVV!IK7rI;@2_eF
z<be82;LhPC8h$TxL@VaG^OA0PF)uY(3tOIiNrpyi0u9)iUo`XSF#i-3&1C76F*}s#
z1XJ{mtP=IPP)yw;YCIYLi>LsfHj#b#sX!?;rx*RkjaBb7c0MvO!jCjCHmdD@i6GTL
zFIOdpYEhmR@?GGr1ja7h-Va@Cb7SBU&(?dIKeEGWtYcz6*#+k!;<%*)G<?j~W)&@!
zGbJBi*emQBu|}x4TY1FlPg(M%dp2f$c6dtr^|Yn82sgVT{5zu{|2@dp{Gyq@nMo5X
zGZaL}5r^SP=OeGgsww@8V9iG%WN&Y=)5?Rg=Q~&f4J6_T@ospp2TI<l?xJhfu~e=M
zf8fr&!`je+*iR1;T`iC#_{S#(dex~rOo%JB)T}7(0;Cf(!5!-P5E{aUnaw5rFj4jh
z%#N|J-ss~+)19c_%+?{T&{Nbxd}SJD)tGrK@5elOGZPF~q%ViMSj1En6~LC$MV^fo
zu?$snb5k=Kxs4--Ca<2+DzKW$3_o?4GsR~g%G+sk7BzeKLkktC-8HtinkHTiQcyTb
z|8*3QuUernVvr15BPDUGhSH(#!=C){yV?^M)Dt3h9;AG1L!>0bc<3rC6(rWiQloOX
zC^l*PsS6qz(Y3oexb@O`KbVO@2G}7xK?|yGPm$hUeN%7DH|@;Ri11qtk7%Z15^#di
znwx*GHZt8x{S+eb3g?_lBj@V%!g<-8C9$3pV{~}Ql@m}*YcPg0!rQMeMB6L)<SaM}
zKf;BuDlggl*m(ZCp$}F$j{|w>nW7b=Uy2{wibMHU*$g-~mLjPi0G@RznM5jVsnyf?
zp_jU<+VoU2spF{scBqjHr}nAQ-u&09(``b`l4R}<MaAI%a7h1S6R@QFUZx?0!fUl$
z1OFwPDIu@Zi0@-u;umOPjOnqQ40~S=50#x~LT8LIU#gv$&KyDhEMfY2$?;XNca#&Z
zvIdG<D=u7zQ|cIFm0L2FniKJ81QKgcrqy4BoPv#ve2zGdR>DQ)mK_#F1S<WK$XmtV
zV!(>JN>8JDh&xIPh3(VVTl)bOHB;GdQpmj)`$j%7U%61=XVHGElPt|OxAAgK!;mQy
z9BMg!EM=ch?<5*!BqxyP`re|x#N$YmR<J8sygYU;m;}s6FHSC(Mp`*xocQN+&IR-k
zR~qbU3`v*+2UOy<UQc}zEUr>dSChGYv$0aqn&$n?dd=%VWj?>F8kH87K4Y3tWLuTW
za@?TT6Bh1tj)74Zb9YKmB1jwj)nqn1ig7shEaTl?+Rd_O^viJ2QgbB%ZN7Ff)X(7M
zx$K(rCKB%>3hcR+%h!5v`PeBeDwW}(Eg{;LODCX6l!M$9ME=;Hg{!6`!bb(WKiluD
zWV21uP~3cHcbIf`m=0(=>WRXoTrOp_S{F5QM>mL_JU%ecLY#Sz3RH4<k!#dFFI>`e
z3cRDD?F)SLMt50fD5u}FnsdR3yfRgfEKm6fermo%6;tYq@bnJmcMt&!3;_~5yY2Ue
z?2-xGI`vv^IeaO#oevaJL$~9?jCz6&p@h8JNs+ef+-k~jm$@sYT#NVNuX;0Uu`BJp
z9w&Yzh)ZhDvU@H<K6sW@8t@ME$0z#0?D2LYM333B(7P;%a<`MeDI^1dF=fIN-}!nJ
zGn--3i(!6}43#wpQlYlUx)^cNsbgB*GFO$VPga@Ek_YKB>nn7u$ceiWWcOE(Ypb!k
zhv~GPih#M~FM~SUk9|7N{Nx(37fg)grgAaag2+3Wx5h_Y$Z(YXUzG<g@eW6pk<Ml7
ze8<@Q84((jm$Ls7Kjou#`U;N_!ECM{^444EnC7&0gW?7DpB11B*@uSnts;fZtPUo$
zt2gK%+k}{Xv%tfze@m=uy(`lA>L`f{{nOW%kKQ8GhYO|?8W<%Z1R6qLAo2ra2Tut3
z|6{UV4ldKBy?Yx^bqT=0Z{Epk1D2Y;dq@cQ*j)X|@3;ii6aAM!B1}z21Lyoij*%Wy
zAMX9n<BA1umoV_xEY)AfG9y<<LPK6x?{*6QniwF9<jG(fqT<F!I0TU={)3=Sl-Dxq
z4FH7x+)CAFQTX+*QSdLMDXRK7eu#ew8!vC_7)pVN1(_qbu1-K;KVUr~FtrFJzz++8
zG-}Y!lbTGaj8WkQn*i(=B5I8;LOL=VPTGS;xhHVk&;%0cGQP|w1OUR(W5G03vprO^
z0Y5&BA}tj2q8JS3K>at&E*)yGZwF9T@_7G&WS7FN)D}jP=r=410Hb7lx_B`=)NAc(
z3e@N4`B$AF#~f}YSDD>YrhlN?whEdr&0Z0MznE42YcTMXT}~~}0HoV-cDcXR=J%nD
z62y;;aS_;mEUJXiQ`Y)RP3FreyEW@ib?Fs2^rS&YEcnyRDApL(Ep@k2mZ&z-V&5d5
z`os?x&aJ@FYX*#wR(Tv5-UQeRoX1hrXOcH_A0240p-uFXK=2*bvX#CO0`ns1=u=pz
zRy>saW=uzBcA+iJYH=dK4rl5m{u>vdzE~+g3;5D^Dm>zI^pGtRy~E71ME};)X1CuN
z8yY~g1ca$zCH*dfr6B!}t(*ccrS1D_(ri*i?WLj}Q@L}6-~CFW1ViLv!M6(7tKS9f
zYTvw&Z8mwYM9qHV^)&C*cj^43Kat1^Jvz4X*5Az|l=^bh@@Pp*Q6-YudrzzI8d+YP
z&B+ArBFo&pt`<^1Uxm4}_C8;1rB|>XZ6edp;5JON{@&YvRG5K7iQ)+MCpRT5BWerE
zc%$FYUQ$q&<(c<>Mx97$zwa`aYj0Xi51VMlcJ+LET7Ue6AkqC7;Y%(R-W>9l`Ud?d
zMwIRD94NgRX*uZzX2p3&gPqC#{Po=;75x}-E?Xh5t%hsDGQ_3cI4g+{s9t46bZ}66
z{9OU=q{X6wKhNhvTWtAYBMX{6J>v-fajl+v9d=b|eCSd2#S$i%nZM=Ld9G@lauIz~
zJ!okueXeclnD{nlN*^xj5jygA-}v!Q^bg%x&_oEcn&;P-z)59?%IDN95^8x0*Al*L
z{F*_5>fdFEllg0I59cAq$~SeN8r$5lvvH-_$%HVPL0#1EImJ!r)-da4WiP+Zxn^h>
z@z^Es1(!L}dO2ru&P0KXm2W}#tJ#I{81rxHYWyyiUhpI!n%Y|qhy=53>5yiZQEU6C
zo#ihisr`yGGi4fxy!YKMaasap-lZMI<>I<}m+W5yBE@TCzMRl6me`<UUfOh=508w&
z>h_Ehzr+^`daPP?G8$P*HTxpB5)D54XJWW2D}b-4k~bG{OEQFcrDM(9!{8JxlIZL2
z=x@3+pn9xM0c>+mB=|WSiwS>Ve5eu072u!43HF~)5Dww3SAd&yvN7m3KrbO4GRjFS
z9BX;N5Tg3+34fP#@}Dq_y4F5C7a9_yYs;4R91zPPO{CuhtrIHbmYJ{5sfoL>2xe~$
z$`wVZ!t-o@i#-n^pewr`QDmn|n>_fNYV!#wNyqEu^H0RGKy?I{TWk7RIZ)_#IU^xY
z`8?1}L(+`<c*cQp6;ZbjHaPz#&UDJFuE=1agWISb2=@pvAGiF0+s%{VC85;6hU@H>
zjgC0p@sN1st7Ahp1gpzHfBEbAo>~rMp-8Y424;HyP-Ou$`Nw)9+`4cj@4_wQQwxH*
zvSUnr;MZmgt*9r&M#1(ns(^4zV{x0SUrOokXjYQ(Lk=xI!B(#o4*?(>z^K7oploUv
zDtM>_aguy)bq2m3aE~{e#HIw?^W7a`sIBan>+B>N#sh!utJ(ju`c%^Nz&P3j5Z3Xv
zqm;ji7Z4O=L=&NaAGm;OMB{%BsNo!7*z6#^hid~p4Pd>81j6Mgar`^WUO^ntiz3YX
ze~MBdNQdiL+gDzDFT34)S~i6@?oAR*@>*-6rucHn?@jmv45Su2&SKMi+Y?cr?(?yT
zuzPHa0obTa;ma%i6cu%KL|Vf4Bi#YRZzSh#VGtQ0Is)hp*mW-;j?#ptjjaT9Pyg%E
zTjwi88`rKF98~wNI0E~t0kvjxaND|k>AL-P@E!=<3TpDfoyE%D<2%v`0FfuZXxdq9
zr}HW^`x2=bU34SemI>#@18ajIyg90?P5!`9by7hZ3hw15wd~7=+6QkH@8><oy}p=%
z?QAyzk>%iux;8dUe<vLbJ(>x$X$J2SA8QxvoXro~<Z`vwyJ~?#)Yb1Jjww9BHELE+
zc*^x_ddNE2*H6;uT?-u|C=Wpo1lexOFc~)8Gja>qsSD#jG#i;Xx#AEhUiEbis^aR^
zWB{C7Tp%C5wF&17SiJ7J-L7W3)N4C2D|T+yLl3FBF!7nZUApsH$w#pgKbe3%USXdi
zrJ2q~dcJx0Lc1s0V%o71_jpX%*A*7G<MlVlNWmScLmAzBidITBme<8Be+}gQEvUWD
zczUEgzc+~%>y4e8WpH$jLSg!^BL75Y&NCsnXX0&75sXa69#+~rR&|kuDfHMBF1&9)
zVtZQL(^J&tC6U+E{Igb>eIM<@NE9(E2*-U%3Od!#AuGr%I=r35e3fWB{OfNlHietF
zd^>18`#pOw;*cD?;pE<(em}9b*NbQ7Wt&B>7Q)l;5CQguNmP5D!tOv-VQCmb6yQj}
zgy_}ai`cC9(_Xw2G_XlawlSn)A5E;7QVMy)MeVbSf>XTH8O<WA`cX2)ab;cv=C$nX
z=mq!|6Ku;QJB4sh7*ZXQ!c$P!UNXjCGy9#QM$nGPhYqI?T-;&1Q?pWtJq?CJsRT;y
z`vpOQ1=rp^EWKKPy?n##BRygr`<tPteq&vYh#oJ*{Z8|anFpP5Avy>h$7(cCx67_n
z3;Q0-eZuR$+lPPdx5amTVz)=oDlDPi<v?BS29}N&P4QLc1))LrOU4nd%ZM1=qSI1m
zqcBpygE}(xAq8>8k#qsiJ<J~#@?gPrq9bc=kdW`Gv|mb;m8Vm8(#P)59Wd;pK-!Eq
zf6Jrzyn8k(sL)0>sSAG?z%_$;@#)nL0HuX}x@5_j!%zFF*0JGVE1a4=FxbLtY9G3k
zKbD>J9JE42!}*uSKoF4*dbPD(gFL;~PIb@Ktiq`qPYOqrZX3k<oj1;V-va(qjm2<D
zD_8op{5%~k1oE_#;L?Rin(z<xA=(j|UmivyS}s{PkN6^l>jUeieCn|xA9Vxcl<UpC
zFZ$B66Gswil@$-yE*@Oo{56<UL|CL^1^%iyq-*Wn19mu*>KKbB!`ru{8yV6*KXjfa
z%GNja(CZh7xG^mR=iKHOwji=q*tiL$j^d`vk!$Yf%!w>o5*@LGdK4dp&~p<a{lGAm
zr~5Xm&FE4bSc+sJ*c>P87T5*M6nj8zKoMP_@-R)#3~zZt5$(Bu@IOd$>G_VsEnJbA
zdw2~{{rqaF%G)EYod<?9gm3RVcec0v%$HWbC#qp^OBS0-Mw`sPQVONmseJ<N&JhTE
z)&pOKEjL1;1u<5!(1E>pcte|vPd9fA6y__M$I-9d(`FMKEdoPQDA-XOA=fPaP=L^(
zX!(`M|2`l5UkK;dQ9m$c{t?G*3#`263;h4mr>pe3h#O{W>W82)8yML1$#5H+>-YO7
zR6l05@0oKWJv(m8Kw7cJ-v`YA|30*ILZKaLx!b|Y?^6-DZN4}3m4N{mD9r(5AQc25
z%})LvMwj%huk<&g=D13&I*>sA{u$?f(icDyylJxW0p*jm_!|;@|3dJ`4Q6EvtvDwx
ztveU^M>Wq10|k?+(p-9mRJubAX2t%X0J~_I#sYjlKdezJRrtDTP||AgFMP7s|H6K|
z!(knFSpV65Eb}i_4PYl<8hc~Sin#`mT8JSSBFvl@>EG@%22IXE+es-l#$5my%!vK!
zE7@_fTIaSc*3G>$2KaZq9XzMPg6&CS2cKsB(DC#~C}0%Nl&!+aE-KxwTXL@@^gyeX
zvUCl<k}gMCywcT93A=UI(hN0AN#}}jyG>Q~01?r%_sEb?GK7XMxrRfT7gVdgHO#)m
ziGp8n{NdPhJ{hiE_t6ILv2t|o{Mzi9so*%ve$eG&-pcMz;eGf!cq&<>GS-m%Ad@T{
za(XURa%k!4tpH!-*3<!~(cFQAQP_B3*e`A4JzhG}_ml)>fz9Do_+z2a?r=|qgkdnV
zH6~Hwao9By3@erE(Ik<!y4YMMC-<rOV9bd%>=v5AK6Z@glq}gjn7Ev1DAJx;B?Q_{
zsRR#bs8x3wLVZWGktr$^^B+<=*1wwM{b_bfAn0@+fLWR`NTyD|J2=nozhoZptU!cY
zVdnOceXI4I#6H(lMjqw7kg+d|u&%W2M8PU$(;CF0<r_o%%&cr{Tjt$(Vp_(#Y=Syv
zd@|g=4tC-5fao6u=ZwfZX1Ym@V#?l_vv6UyQLOTLKgo{Oti{V^L1(26=Ht)gkpoQE
z^WMnRs9hO;F=KX@b$x|7>#M|6#-O!`x+q;S2Itr+bCNQDtDkcESAvss3+n$_Nf!&y
z+Q41rjL5Y_z4LFfmY(6vqXTluEE1{@e(3gY;g<T`_Sdv7@9~;qS2(tTlI2TQjP}ag
z)NSP*E6TqhwDx2*lzOn^@P08^I}0ubI?pO^_)g+)O;I%cQ3e}v)&%*XmeNPsiZ9Ct
zQCy!Y&#C)ZAsX8&9=&|Yj;<C|M(8uL3YzK|DX?p+v~_~>AJ4*Y>K2zel*3y}=}m7F
zhNt(@W}fAe)s)f>r0I=}NTvJk+wb80E~r1Mz-qp!a4)P$e>9K-Wi{-jDXVZcxot*8
z2)%Y5X9rY@KL@(T`xpuQkm#vh^8bh(c=WI$)eq1V(?ahhQ#Eu0Pga~ZM?-k8-BAK6
z7y3oaBaqM2PvM5b8KxbtpYI20Q4*xPJm^{5^%^4k3Y1xKoNGr-_+{f1M-WJ$Gz_s5
zeMrLS4nW*Nu6KWo=5qPkYUEsjYhJ4N+jwH5^CLoeXYE;Dl$G#9K1*ct(i-kGGmBIH
zAl09OSQfy{@H_(ZsnvB|<rZ*S=sZbOA!CF$==8*W8Hw3J{qanp^xfr^e^n;pa0mK_
zzu*gNf`276i3soK;pBpc<9sZbv|rNahe>s(ZmtjKj>C=_?Q^_(gVs$pJw~_JOzqZ-
zx@&WpW6-C92BM-1=mt7N_H$Y~dtYL7F?Xo-uu>Y=8dQC{n!H|ppT@$e{}?~}vFzh_
zcj79-DY5GIDc$*w{wAd?H7e~;B$9b|@}?*0BDJq{k-?FvHMKAuNHCp#_^v(`8+Uh+
zuv&3Ok;jm3HSk-JzdU>Ta{5J;$N5W$Ma6mF&D)#}_8!reGbEBGck;6*>B0yjh(mO2
zgD!Y;-n86zhyI;A%V3jP$6YYw#oRF<_DU$;FJ^>N{}Ws2gvv$M^q&|h_}fh50bKvC
z)|>v3|6BnH(l=yxNenKxoZbjME>1@x2^_Ld)7JU*cZ}j{(RGv7qnMXV{~QiookQfD
zKgf&8y-uhyNDqNp$`fAv`}Fs#p>b;Fe=#3t<E{xTeuw{cjksQ*yGUHEIIoE3ex#B<
zX!WAquS);gd2ZnqG61APU`B!%KI<}eS&|ZNg&g<kDiNcDBfIy=yk9tj(|2}fM-o~I
zcf<M^ro<klu)D_5j($iu?Xe^@26Yu8RBv#mnD>?^2D6CV^6s<e@u$=GC2KX1H<L}%
z2?2QP_{n<0Pc=iKzef3Y#&kaPIU`o$<JGwPdOOs}J8dGkh@V*PZ$CXs@bOarH;-;4
zUXw~n4o06X;iT%1{)V>MK-aO{SYEu>&s7sdEdB|K-So$G<2&KBf5Qqt#B@I`n`@t(
z^v?L&t(X<sj{X#JK_oVpylTTlMw&m<`TF+4unh~PI|y-;)w%qw`Xqw%Ts3A;JuDRE
zlkaOLdc*c{4@5Y9SNast8~bw~x~jbNX4OpxS@PSbODiUT!T+;cz#`yIr;X&TUq@z|
zEI4lJQ0Q96H)@z;a$5nMR3LQybo6;;F4%G$`%|49IAO}I9lF0zV=UGu#yS)6Xkcl`
zplp-QJ%@V#cc#zKKBw4oCxUh1x0}DMo@@K33XRjil}oTOqbm!qLxl-y>4hf52YFhW
zE)V5P95Hecxy#i5<VxT%dO`EfJWEC~2!{5dN4}kgD}I)U#xb80KLIC-B40VWMuCV6
z*L@+3gm+ZC(r-5qovd>ZmXqeGu*88FX%Qt6(QpYkPArI4%KiLB#xU*nK&1I_g`z$o
zp{|q<(1$<TjjHFqgr#$r_{dUOlv7szap{fg(kgay6+6wzSJ$;$H;)cMdQ$~E9bAZd
zeFv$NzP{Rwf<O%HEbp`V_o^u9r4tNwWgf|b{lEG3wnbjy*lahVUwZ+WIM7-&|M6M0
zD#|wMcTu>ck8*&_E7-AV!!Ga(h?F-IJ703vF9aFX(w}y*JI{6ToM0ZjAwlj~fyxG+
z6sIi%XBsbS-FYIr<oOr05w{;uOGy3M{XD$PZXdE5GY$@ul=NYzWLUAU6B&5f01xsR
zCGIT_DS3Lik>>|#5Z^Z}-^5#b3A>r7?cBIo-4btPsIhWNoG+S0zitC_mU|2@cDZPj
zFOuvb9=UD6t@^pmpW8i1@GsXBz^3{djOAc8EFEq8vyaAe$4FR^n&g@Ktcj1Ds}zxc
z)q}3zJL?b0Ua*GWn^H%)3r?Li>D93B7fV$Ny4T4mW{JUychH<x{-!INp6#+a*9+x2
zKh`aI=~&eLJ4(-SL`eNbm75Xp!RVHQ?e$TTcb&7ZfJ5w!iZL^qLq((C9EZatAOf0k
z*i5Os`+{y&Ik`mB&uGTvwk2kU?5=AZ@uiH${kH(QTjps{`s24{G@5>1J4QezEL@@t
z3}-ELOmQwko58b8>0`51)!o<pp;0AmFHR1I!b>O5inT~uk1lABdHaO#bb2M6WBtSd
z5RPQVH5K<S+{)+LVS5{0u5vUTF%8syp)9WUN@%4LsV0`3^qK~&oVPq?FK*g?vWhA1
zx{Q5z=6uNXGImBtGswnZJ0Eya?yNb9UeLyv=Xh8%w5dew&zmA44)#y;$%%p5D&~Dx
zG;!BU;rg8S>US;L@y5|Xd4OMS8gy{%2wjZM!LPzyjB(I5SDkG=zgbR<CIqU-(t0_$
z5f=HJB@()N)qCB)e^5T^YD^^lc%fZL(xBBK>y1m+6OT{JhQP4>Rr8&n{`uJ@W#$ak
zenCluRwmo;fqP>t#^F%;0{_de3>&&n_MH7ArHVcpqN{k72iM&t^1~NI85Abs{Iz{I
zzbF|jZN%eF>aB1_n8d$f@{~*5Q^p$eBIR|B3iQ0SVlY+B7qZ$E+O4U4T3kbfe5T_d
zwntLfMzCWRG&l%JqRy}LR?D<Akv8>&D9^amKmtF%GWxngwFgS&l%*znvgq6~<E{Gs
zo78n)yohCOTUK`MZ=<$by2`S<N_D42P~~z}EfGr7zpX+kDTCO7B;5ocyeTnjw@0}s
xAbc`Pscj8kx*1|2>i^R7#eexhX5$WtF6|CQX}9ql@$m^jNls0+R@&n2{{XD%#T5Vm

literal 0
HcmV?d00001

diff --git a/man/figures/rsi_example2.png b/man/figures/rsi_example2.png
new file mode 100644
index 0000000000000000000000000000000000000000..d9240011e779f57055c1e4314a0c02e0a8fdc6cf
GIT binary patch
literal 7874
zcmdT}2{hDwzgNoALfJ}UEJadDvV}wmC9-Erq7lh9nX#+Th%7zI&X6MOWQoQyBeE4E
zp~8&GzR!$xFvi^fjJD@F&w1~8@B7|+?&}=q%>0(`@4I}y>pZ+{ptEz^zHKZlEIajd
z&l$6@Y$md>Y^vj61y<nN)rS}hm-S7xfv=l3ZQ8tfGb=0W)~#FF*x1<F**Q2kI5|1D
zZQHhe`}Q3>cI@1_lZ%Uso11&ru3fu#@8;p*;pOG!<KyG!=NAwV*t2KP-o1MT1qJu*
z+b1L>w15A8VPRnr5fKmwBq}N@CMG5>E-oP<ap1s#g9i^vN=iyeNgX<LNLpI@@ZrNU
zGBQVw965UQ=&@tRWMyTKA3rW9Cnqm2fAZu>1qB5~MMWhgB^4EwQ>RX;s;a7~si~{0
zpFVwBLqp@tnKPQ2np#>~XV0EJckZ0Fwzi(0-h~SnE?&H-udjdU(j@}}14BbYBO@bY
zV`CE&6H`-DGcz-Db8`y|iz`>ISXx?ISy@?ITie*!*xK4&y?WKo&hFZ^Yxef`*RNl9
zaBy&Rbi8rn1{e%>a&o$P^QN=2vx|$%ty{NTU0rYAzU}7b=I-w9;o;%w>FMR=1%W{B
z+_~fJ?G1%OVKA7FkB_gfub-b^KtMoXU|>*CP;hYY-Me@1-Me@H{{4`UkkHW3u&}TP
z4<5kb@bK{PhYue{L_|bJMm~D<C@Lx{IyyQgCMGsE_A&6+el0F8?#YuUPoF-GkB?7C
zNJvafOiD^hPEJlqNl8sjO-oBlPftf65E&U6nVFf-o;}OT%6k6%`HL4Xva_>ua&lh2
ze3_e@o0pfDpPye)P*7M{SX5MmL?U0kdR1IpTvAd}T3T9GR`&Y!>+<sQii(OiZ{DC#
zsLIO9s;a8$>gt-B8Z;V>!C-1@Yu~<oTUS?CUtj<3-MfZ{hQ`LmrlzLm=H`}`miO=9
zfB5jBwY9aatqqICwzs!`{P?k>qvO-3Po15eU0q$>-QAx*fBy32%h#`8dwP0$dwcu(
z`fxZL9*^(u?<Wul-@bht7#J8F92^=N8Xg`V85tQJ9UU7R8y_E^n3$NHoSd4Pnx3AX
znVBIHi6jz<OeW9H&d$xv&Cky-EG#T8E-o!CQ7Dw<<>i%?mDSZ%DwVnxWx2MtMx)Wz
z*VpNE44qEj%z5(@fOt+1T}ucH3s*hke^ZnDOGg$KAwRuynx;OfGkE`i?!~~SS%k?1
z-?=U`&(~Vpi-cZd%`2{!U-d}IO)%QEn@9HLornlsekm?(*R9?zX6DshZGJiLwM~x=
z9`xB{TKhskd$-tT9xc(>+CMr;!BiSn&&$ATu4N^iK!Kb_xX<}(fUJ+DE~~ast3@Vu
zhZAa94xG=Rak4aP5@SW1g?YB1&fXV=>Fyfl;=GR&(~MO;z~Keo@=t8W9)q~Q<glAY
zsDoRlX@yAhN5a%{xo*0|!o((z$-@+1i?rS70JjT)v?r)S<a>>)cJ{_F($m#7mc8Cj
zsl64?a7{>gFrFXSvqK(HinJ&X)Tf!>74J?^yLI79J3b4i^MKGoN$*>A(e_=5=(X&L
z`TDBsJ^dgn+pChy=X$Ug9z=;cEalP7Kz<Yks;#}!mlSie-(8z(zw7>{f|Ci#E>;6L
zS$cC~Ok)LlTcM5Pdmt7>R0_?&H96SGHUqdZeRr42k)F>M1M<nLCfNbBrsf*ny|s`q
zh<md7(lhpCD?7n?0wQU=Lc#*N!VP^+JI?WzC7pNGL2`q}ey_36O=@tQiakra7|Z~;
z+MLGAM*e@X90nn%IG8pB&p^~t^q4A9vrwj_fO7NGXm5sk^n1BOcJ26vszb%RR04ld
zo^;KwQ7|t9v5D=ijp)~ehCvW(DeX%Asym$IdumzmtFGD2<+=5;iSo&(QgVy25li)=
zTUf6+AGy(CTe57?O=~Pj4=qh{y~fhckEIn=Eq&Ja*~+|{E|cny@7kSkOn1OlQV7t|
zvV-LRm8#BhlXDznbEQuYV;_o~sX8lmLZJaLR|o4&*r9XT*~_u2QhqEmTw%>!d?2r}
zbKDvD%w0*zgJ+4U<D1G<`NU030sxDu--W`p99<%~QIKcsABmU4^)Jl==3Zdw72fvv
z&%~V`Q<FcyV}e0|=oyT&x{ocxOPg!1`$9&-lv7&Rn6x5VwNkZ-Ljp5aK7!tgo$Ko6
zz++m>*Q@p|B^)0wvK(I5w6yFwmD=VQ!ekszca?<9)c1j^(@&*vN|8pGoaj<CeA{>A
zD=0{KD8}t%io(YqwmA}q^sStoo%+k7lEMDd&b(M~i6%E0(Q?_8ZyUKu|9>j!|Ee1Q
z!-z-XK8XLtxIley3p@w_o{19a6Wcw_4Z+4=+_JD>7a&FfP-Jct=VziERrXweO8rA;
zC(|+j1Y<=2!#VHUlWmaNR(30n9{}h_24de8J+OT-g<?s;SrGQ4YkTwsm?$*{4C(5H
zczrzJUTM-1=E>(T5Aj_YPOrRXz0BhzZ@Ga)BcnsLAbHRK4#NLoy8rijw#flA6*)JY
z)4btukm&;Sve`G#neV}1kWbHzxV65G`f_|m&w3k;Rqk}*-k*M(ql<&cu5t6NNH;tO
zf(X+Um9=?BL`>81?%y0FIRxcNs_@k1JIl8|tQS`w_8A$LXi{|Scdj`|GjgZBu_KYp
z!2Mit?)!sVi!Dx;+rA8`L*)oOQXPx;`xykZ!5vpmQB1}P3YC3etNx3FXPpZ71=mJO
zDtq?!nZ0e1ZKJ*KmZTXh#UC!sA?a;Tusbd~k5^x#KMdI;E2t*zpSiSc41GuI60NZG
zS8Pbq_MgC!6u{P3Zu601+_nXKj!#G?ir~v@6IB9DCSWBMrb)eh1o_dZ3<ZiQBCJsj
z6ln7ytXC$L`tEh(w3SX;o1@+jKuecLLO-F9CnmhIF2@w?Af(H^+k(`o=-yYV{`HR_
zDgTShM!7d&sW5$OVs}T=3`aR2X$=l%rxvJ-=>DSQy&P`TZgX;N^mbQ|zL4GPXFqiL
zwfj|ZBrd-MMTts}VTZpU6uoz>mY&~)s(2)F(zj!-9H+#t8kM_KvBEaIm6pqDkRn(>
z3Uq+S?1d6?x-m)mX>DbYkl$S?4T=heIWI4*b{|E-4Z-PC_|%oVsn%$Y(@NPVsE*l`
z?TtsK;|`}%V`ZMfKzqLqNF3XVAbvrF09o#blbv!&ph~2)mWHHT(yozOj-t4)o}}mV
zTAj|AiPur=rX<e|tcktVnFtoIo}J{4KPDIuC7ntYBb_4^t<Q6TX1`7+_rh}+d2u65
z9FB<gf?$Tj>mT{9J(m>@K5KH3Xt?v#+#e4<-?w&qczW{Ty&8DY%D~-)zJaKnLsW_(
z{k&r)4x=vPN6C#-4S0vPbik$+^Zv{}Mn-&D5o6rpUPUJ;A5#v8ya-9ZhH_&E5?3ZO
zJ4yHI4YlKstadyO#88iXxS{%qc$tHlqb~3m4%Z_~ZP<bz-u;;YWu+*2W!&hS$KK+B
zcofRaOb2!0W<)!!Tb|^wzOSV7KY(D*4YxEPlJ@&|-g;WGUy0uFcI*`)3P_~Pm{ljj
zRho*!B&=6$q9EL^PIb@9@+f!yc>7J&Kzy7(sk`F?Wi=W9)(Hp8`hf!vSp88-Oj3|=
zP4$>HxZz&ZU-qOZEAW0(OW#g0LyYg*%Ji?BWSzk~4QeP13GF&=w|||srxho(&xcem
z3L6$t;>nR2ZJ|13$FPz6H6M6Dtv0v>J?!M4A0A#`y}RIEMGvQNod_-M+(^{xx~WBp
zQc@HQh8*P(Z^30*@6Rq+*W=98D@X7F1B=$UIBKSQ;=DkLZ_e!BQ0*seg-{nq>rYq1
zV#QL^1eh^ngG&_@efc>he&^~pwseeh&%Iv^=`D@2`Y)auN7_Lv;W!*&+TY*h&c`X~
zitSJXW!0lh7!AZ6@SnfeI9)>*GN61UDv49N_RVY8^euAy#$QiLJ?>K4kjYHOPzpkr
z2z`s?bjk^1T6FEPqdOOb;Qe6&xE@bM=jDmghJ-Qdbu;(Vc(DfCnfO?iXz35QwrdM}
z_0%2i+bDNe+CcO-<G->2S;p7QCgdpN&#uQMecJW=<^Ad~q<ZUwNUs7~b-pvtZO?S-
zpTs-_gT!{K509=_g~+#qamT)@3{00X9WPYlkLk7vT#PqE1MH?`x-L}a*B;Z9Ei^L4
zS2pe9FO|HTi179D(#Sta-&B0^ucvF~x~S&fOUS%Vpv@%K_vkBbWL8H3fyc_Wa)wfL
z%NNn06U6!KZ18W<*nT$yAAIY9iP`oR+$oi;KmD%QmM@y|;x<HE+03Ly7k*X-sIVvb
z&41#jYG0nO&@run(j;+(F=HqB+w|IqAKr1@qJ9Z$y=LNc<yM7q%0#I`<8&U7?!Pl-
z$|u}1$|wyym!&Ro_rrgWM}M@X0fopBv`^QJ|E~zoudP#p8gi9%V#32S;+!HZwQcwj
z)7S=&ImX;3AFVh@N9So0oALKuiHP(QU5Qzk6)}>q<4-p?c+A&KJChUXG(3EDDdcry
zH9amE{KiG(XK^CU(}Ux-+h1cdvOHvz_S2VB#UHoPrnYA#%wrloy48L}o^X=`B`8CU
zAh$N`7(1@{uRBVFr4i9_)i7FDP52M34Dvm6B{n`Uzil9AQ_A}%LKW)=n~Z^Kdn0r)
zp9E#g*~BpF`R}P%JcU}qn}dm_7yQBdeaTo9#wDPU|6cP)%wlMUV*E1%ri$;z*32W!
zZe;#W#($4S{-plzS73uqW|qGxs(I94t>EzXADVg2H{(&VY-z^zgJoUm-hohr*~mt}
zk)ykeOk@Yx06G)^t_F|Ec5>W71V-Bd_(Q7r0Oq;~v&i|eX?TmcI7FW#+gtl0)6Oy(
z0sPuCB9tECYXw1c1!F(CbC1+OxtI)aIOQ}*^Q_B^I=2PZ9&s>dQ}vNef58+8@B(FR
z<8Kj&RL`yi_-(IBcsY(f?d`BH#M^C8mNPTJGau!9Gsqom)<=%)Oe1C<rHHN?Ngk#c
z27xya29V)P3Trsulfq4{OSIjij1Ms*`WLCye`Dq~fH$E0jf^T$LafjTBxohNrG?f4
zjEk<>J885%SsjtFCjyw=%oI0b<{23vB;-cLCY-<RSqXO~S_*48_9HJPs(LDR#={CX
z^bD{CRoK-A^=p8FitP_{LSVgc3#IcYAoGb)Y?^CAQ4}tOovV6D7fQogV1<KO1Owsc
zqLI;++fQm%>~w;P7ax?S65r0r?hYm3Jr(1&`mr!1;ScKu7=R^ENU%|Gj{mE6;XmCS
z26D)<%ypn^0z7$;kCV<nFUn{Ie_DXqQZhAL7u=x6>>mH#jNFWdQxu5%|5C&G@2}P~
zajNP!^2V5XjM3?GjudUoJ%FzG_}wu;z8&9GhOwNiK3phB9ZUK@Tt4CzOfP{n7J8{g
z=u;VU88x6oKKL|ONoIeHcxTYwP+(qQIgBoGCRryK={#lx!HvN}&(*q&D)1t%W;(~L
z%gC`?rr@T=4z1&NpH3Qg+r*j?cN!XGM2nYxAISivXi=3hxXEJ(%9q3>;8Wy5%$=*C
zs4Mo&9y+lxu1Rf(OUj#-K;}b#as$@@5B|R$;V_3uz@UJ%-lr0$$Y>^+!zqG#zzkYP
z>gRX5Po<nfR&Hq$G{Euu@N4=IV=UuG><=Iu0rE~1(yD{0YZ;-7W)Mw96ef6twT;oS
zvXKG$?KUW}fQLZu$Q<4TnyjQ}XG;QUj+K?6McMc;EF<{<wr7kB{%HvObAb65PnWjT
zh^-=NfNr(}x)MOOS4}=59|RifyG-OUrjCH1(Ay=?n3y()J}nyCT##O+!-yQWuef~I
zpX;GZ9|Q^@O+bqc$omgI{0C?R%Lrh#Ea|1cJ^(d>vwIXj(<@#XPb*zK{TWr2*0v?~
zb><eazG5t6pDuD&d`4F9y<!0brkcW8%f4~s8|xGfj|Fcg-}?xX-?Yr~uTGH}^XE@(
z5cmB{BwL(S{QUVo`XV3bs%k)M$fdMkZA?|h6taxH?Q*JHI(3e-esbIz^xjvZbM@>3
zJ51jCysZ1$bY9abdPI=<7w-;rrw8KR^Eh{Z=gEsuP){lPwGd>in%t0_atc$N6Ikq8
z2_r^jP1>R50&JJY6X75c)j??!CCHtX8E<gUaMw~~Ih@;frp#!l-<j^1J^fThHaUm(
zu&T<@denP?UtQj>Zj{_TIB$(!bG~8m0MqXl`QB3kX>-ksGA=03m+L>a-gu|fX_-ha
z_FwEQL0!skSCl99SCELR`jelx(r2v}USfR?yq>HxrQ*gC$*fg1r70(nS;+B@QTCh%
zNU&K&z)BLb3@>1zmSP18>Q?bzcqACnRRcj730U!*hEh##T#_^cpZ~nC)ZvMoH^038
zLM%e)bSZKwmYOayi&@yCL%^pX+{4k&T;_XJ2J(B$Wi@C|%lpI^TS}_x?b>>>oLVnG
zcT4jh3McQjJ^p$I<uWQ(BrLcdY11HU+(DJ7bEi#FD%$Z2;Z44ykby|<^Vug*7Tkux
zMum1kmi+=S-0J8O*`jZzgBN_rxB^RSdpHnz)T!GEgchVs6JQkmxJtK%g|TOp9Q9;B
zCr5~3D{bBgsSx0|B9SB|p9^=_4>me#7c~08wW5J+*fvNGX8B1haFn;>gQqM&t5Ptl
zDU|NF)S8yPw03ELs-TqaKYO8*7oudw3pZ1%mZkBlpLEgkUX65FnUl8zR~-$~m@JW!
zB6J7Eh{&Dwtn?<v##g&IoF(ExUKI~;RT^tKONC;8^Lt0M;uiQcHZXVy?a`H&_t`ci
zkE=pA$nE29FW8+av9}H=l8&Zf@27a`^P|FUkJoS+62ruC_jYQPJcQ65Krfnfv<Pte
z0>|Ok#1d5tzhKEx8FSM?tRQ*KLo`P)Q7L*~l@#7*XxSmU=2`KgKxe>IMvw;0lsa#4
zXqY(q0KJngw5Av!q3n(AjR=9|VBk6^dfUbMXdwhPaPV8hMK;3|M=KHqo7(NxLQKju
zNTFIp1LI(Y#^(mh3FN|WR~^d36e>p1z^n$Y2zSq3J#un%1)UpjZLFi=)Y_YAyLLw9
z36~eF<ir~%eC1q-OP!cS(t;A{GwxPYkBI9wLP2ytZnx`g)D7#)Su{PVy<V^ra=Iz_
zIoW43koJ>k6%Iu27q%2@H$BYu3$ay^qsZuOcG+mCRS<C?F3E|i-G0ru#_l-S9#JkO
zA6P%8NCMGxW(?JnfU_1M+8z1Q>~3wN>&@ABhy<&X^fyzVFkXwC7ZY}%8_RIlvv=5#
zu^-(PWr9aJ!*I51Ycdn`wT4wP!Ni(52Kj`-C=b6xC}jt~un(ZGTZFoU-Q+Hrur00R
nc)v&Z82_C?|L)hVbx`PB7LQX&tupY_9E+Z|!MP%>Yxn;PO6o$t

literal 0
HcmV?d00001

diff --git a/man/figures/rsi_example3.png b/man/figures/rsi_example3.png
new file mode 100644
index 0000000000000000000000000000000000000000..36857cc5738418c4d01b58f4a3c5b4a4b146023b
GIT binary patch
literal 7160
zcmd5>cT`i~mZld4g3^&9DuN=?M34}gA}UQ(ique}bR{7oHB_aDh)U=n2uKx>me5fI
z=^cbiQH0PER3Jb=<_6#I&8(SOzW3I8v({v-{BiExXSeU$d!Lg$zHXqyen9X54Gj(Z
z)hm~7(9qBkX=rGhS?Iwh5!$tq-~;k?y<6Jgl9ra1j*gC=o}Q7Bk%@_knVFe|g@u)s
z^}vAx2M-=(V`F1yXXoJHfIuLeoScUa9XfpYFc%jWH#avA4-YRdFCQNtKR>^KfPkQ&
z;E^LogoK2Ig@r{#L_|eJj~+cLCMI_5*s<fskDoYkLR?(@<jIo~5)zV<lBZ6cI(_=I
zl$6w&GiRiwrO%!{D<dN#D=RA}Cnqm2ub`lC?%X*=MaA>y&nqh{U$}5VMMVV)g{rEm
z!eB5pH8pj0bqx&-O-)TLEv<_eFJ8KINn2Z6M@Q$%l`B`TUcGkhny#*{o}Qk*zP^Eh
zfuW({_3PJ-jErvFxN-C5&0Du_85<j$n3$NFnwpuJ-M)R>+}zy4!ot$h(#p!p+S=O2
z#s&_D+uGXN+1c6K+dDWoI668yIXNK^h&y-g+`W6(+1c5}#l_Xt)y>V#-QC^8!^6|l
z^WMFC_wV2L^729=kth_(+uIwBMn8D)z{khu;lqc1et!P`{s932fq{WRK|vS{CO9}a
zBqSs>G&C$MEId3sA|m3^qeqV)KaPxyjEahij*gCriFxwmNo;KF)2C14z)uhS*|TTy
z@$t`}KY#J!#mkp36A}`zSZrcqVp394a&mG?N=j;KYFb)adU|?hX6CC`ud=eTva_>a
zzkZ#Qlarg9o0pgO=FOY@{QQD~g2KYWqN1YW;^LB$lDBW)mX?;5m6es3m%n@WuA-u%
zva+(Os;auWx~8TEhr{9V_}bdqy1Kgh`uc{3hQ`LmrlzLm=H`}`miO=9fB5jBwYBx*
z$B%7oZSC#tpFVx+=;-L|?Ck35>hA9D>FMe1?d|L9>+kO$7#J8F92^=N8W|ZG9UUDT
z8~gnE^Z5Aq#KZ&u0F#rGQ&UsZ)6+9EGhe=Z`TF(i?Ck8fZ{OzT=H}<;2?WBz!ouR>
z;?mO6^78V^$_kN4B#}s~tE=C?fB*60$IqWX*Vfk7*Vi{THppc1=H}+s*4Fm+Hibgj
ziMhSAv$MOq`|H=Qy}gFLy*<VLsyvYTte#hHd(qHvv{1jaA3gFMXlVFzuU^u;^&n*>
z)!*-y@xzZLs1MFGao$o3-i+dl=Ywz9<P(Pz9LFmznVdQn$$vWPT3E)9ex8CN&QUAw
zDiq$%#4oer(Sd>>Ccd^;x*FOhK8a<VC*E^j(P}%0sYS~B{vy;VB+eQ<gL)A=Yv$7^
zDO2R~?umAo&wUbUi|85OJ}bvh^BKPBPeap2kE4$XB5EAZfcCNwIL^`i!zIfjVy@T8
zlVizl<Am+>;$=bfLSF7T#bA2OC0^p9T4fO7*RP_Mr4;hncW&=3Z3wkkWggS|sT<x)
z{q<bc*=K9|X6+9?BJRDD4-}ueR+s)+z<1X1XAbh(6ch!KIEt%^cD01H*Y;S;=X>_8
zC+w!AVtu~Zbg4ECPmJAL!(}-gbjf{Edr$Yv2fdAw`aqx~K(yMg-s^<7hxlsqi>iRT
zlSQ~3h2l43frtQAG5osHT)>{9Q%@|?wW+u6fuIZSUDNUhJf$q5DuKgZVZO4mJMqIY
zNUe%b(a>vwC}_aXk+@4hf}B^Tx8<M+zUr{qVQbHHAF@n@%GXHo_@}@WFKW5&(b&)9
z!WOc~K;`jS{c=9HxMgqNcBh((XyaIiLSlYI;!tx(mlOpjf>XM?B|~O^F{l6F!zn4u
z#lmHaIV(3)%RYz0Wl^ESopS2V9iVHN-nwt{<rho7p{;yd%}po?b4pI8z@#m#zSD1U
zYW!s$G;2ah$En_7a=A0+JyIQc_gv@$%dxt$<n_kztT`R7AIF6(jv}W#lYMT=y(`5F
zE6i0^bf>nfw7Mm2=q`O96pxdxu_bhyZM9Y8Zp2s%;js&nwMv3o2+LBl!)5q2vg?a2
zh?7=8sQtlEqC%h(F!>IL7e*-*3pobh*PF8n=aNF^zS_s1ZO9DOM8^wPeY3P071W2%
zv^A)6vu7$4o1I}D<B8cM{^W9!nA1_N)F5#52Nze|r3V)10AE^Qp=^vqe;x8qQS8rr
zypr!5wVul~o-l=;`mwPIJh?ld=2*DcATDg(r~!|R*CV51pD21v`RCVcy{Y()-;seB
z*v6%&zeM?A=rC^#<N>mu9jXMbZcjmcApNL%IajWq-mX#$$28uO*s8z})#1A&at6rn
zN#3*bq&%5LYoiK!A#bix+lE*Z+WgSd%i)8#DaXo{L96f3D841MUi;p;%2!U1y!%{c
zFQ#+~XcU9R=OikK3LX4XQjew*{VL<WC8+E+`E_Z+o#V~Picu{y+NAfS&cKX7UfV6^
zhD+pa&@p{yprs>sclDd+&Xdx;?a15eRxL_4{+Px=CK(0Vf2~OWoWcL)?Tw|$)z0iw
z4X!P+2i!Taxa?v7naI}dZys=dMkN+Ks~bxJrV94@19H*9M;LO>zCoq$_4h^cwL;mm
z@-|tX*RU>UP$TCFo<$it11@Q4o_k5{1}&0;pEOgG8y9!f%o*=+=9E|{*BiSbTt`&d
zvZk3ZZ*-l3^uS+8z+TkuFEf*yYgUjtzNX)9Ogo;;IMgeAhXKhfNDE*tYoPz_%LpPK
zVIVmgWy`&DRMv5(*Z))Hlm0jSmzyATH+eaLTo2dNX*1`dUu6l09>PLZ7}}RaH~=y@
zG*-}pmyN9HHt?*j4b<4b9qZIT$N>bdEq*sm^Ax<?XUU#4!7w~a?mMkfe$)qO$)s#P
zG(q!)8?A;kJ`cM{uJNd<|E~P^W#k2}Ib$#n>d)yYXroR4V{)tVTvwR>?;-~^XQmT<
zc&V|zZd@CH32~4Dcv~blGr6$=I0fD<qnH8N$l=SqUrxT|ZHOHfZbnw|x$zuzf9SbT
zzl9245bj<}nT$L4Hb|eunS-yp<C8O-!yZeVS#X*c3OAY_tBXaD9WXOvR_h3|)xoYK
zAm|a-;mt@g<&xbh!TOf#B*-rHP`Uv_a@(%roOFyW_l*pyZTc9@KNv~fJL!^_M7Y}+
z2t+od1TzxNgo~%gLFs{ZS|AleoDCv&YoPCIpb;6+3mJHY-c+GpK!9Mt@jt}HZjzts
zw?>{%%PKD^tBkoGX&6L2DiE-cc4NZ2pm4?K;|+g0Lt|@DP4H{qMZ@VgpAI4=B3yx(
z{(B+eMGoPt_DjWTi1*`jA1f!YqVj4&iVw_nmo@zikwX$yY@cwM=PKP!PnmmaZWIm|
zYr&|>`%w+vU1_HfsnKet9x7FL5KRK2vi6r|a>eWiOGX0h&qt{dCb{f#nXHX!SAJ!d
zBG5k@E?h`HKi|1@_BjiZx&3J|<zVy>cG&#KWNqWlG(oP@%yq>#wp5Y+cpo$1t5q1V
zz!&YN^aZE0&8M(YUr}6s{g^kP&L2M>IKnj}B6W<UpI)7p=)%v4tnxgB-M%v7FA<gG
z@>r#Ma3bJjf`CPl|9Z5r(i3e_BuUOxyHXWRVzk1U-d-}`AF>KNT)OaUSYXO`SQ)!s
z0yBB8YBI9(^;oejcQ4Sa4PQKNh1Tsg0GQ23zP1Jtyd`t)Il5x}duC^TD!DK9?>_2l
zG^@LpU2y%_dB6*C8)73eZ2utz7OhAvr#?%!i&)|i=2g2sTt7gx2pYQJF7vrFm$EKC
zx6x#kic98GAUUDl^Nwx|3hD?*t<+=NKR24)+q>cn^lYWBdya^b&E)Q`LnO`=&LmBW
z>7i@5PRlzsBsa>|wcLBC6HJskJqB$tB|Q{(3CJa0I5HCWv$UJh%Gc!9d9v?f<F;ki
znQei&?y-dDuHA_pVAC^>TVE|Qf|*3ycvmmJ-U2(C%4av0fjZdr&;brT|9$`B%*o6>
zqK%_ROjG=4;%5wiDM7W(rkB-CslJ8A1>ue{ll}v}gvy%o*YLKd1jw<Q3-ucljIS@Z
zuNuQlRFVC)4vB4w(OOud|Hng|mJB&+&B&frWi7ja&^sk%-{+d>O9FG}BVTkZE!s%W
z3SOoKq^UxoC{l9eD+<vk2lICz^IlOfQDG@h^vTzs`&V$7nh1eKkoO9Ibh+WhK>O{R
zF8ybifObt&@|2^@rrW}+2D80~-S5vuX`=7<4()~$3?poWS^~yP=4B3h((MZtS`Dx8
z(T`pITCf)RTPokitZMJVvlV59u2Yj>T_I66fyChce#ZIQ-S1(0SSLBK%y77J>8!oA
zS`QDnOS<dTa(=L8^Z@a^-X+JJvdjg(5aMh`0cAL|PZU3u67YM6jS`B&IDX9woomB}
ziw`L#&buSX?siM{34arSd~STH*2J-9*Z`_ywgNnAzC=m4kYR~UC262<lG&8R-EIFB
zi-z`Pq$@5^)$^`HSs$rrjT&BFHq<`XruIzvK<VDmAGcuk6&xPP-{^4h7o~yQ_Dl(%
zc`x^Lm2~S;RY#cPl%xgY#}?kBJPw9M6UL#|Q1x9qFxJxx?)<F+=?HI?dw|pz52CEd
zwfzm(a0<iQ@h<YvRBJwkL>izy^Cht`K*K559=+qR4CL<zd|eN{QPFAeF^963!h=t`
zy1&yIhXA#muZb|YN|qkG2eujRek+FaK+kyc-QQ`M>RW<5s$tyE$m@}_<BvUzHR+J3
z*uCkx`De}GrO6RH*#rrMFJ3f7OXP^MM?hukLPUQx4o18s`u8{;h)sw6THx9r#Z=;|
zULUMsE$k*`p(12~740MMRhEMX%M+cQ_B_6NsiVM=_C1c>78l4Z`Duc|yZ0|SsZf_r
z-Ni9IUKClGJ?n}FmY9~3ucu6ZQ`8pc<yWM>7B)I?6_&hE07NV5H?48FQz?o>p|M@J
z<YumKzMwi869qa|Oynp}9L-@qDGdu-5Z+m3{+rk=9Y~vIT>0hJfb?H^f&>Fj98R>}
z>70F__ZDU!LMpy=VETL~(HpXrWPDZ5Up#Rub4NR^Zb<GE&&JXfC~lD4#|z}9-^^!R
zI5}ctN1ToiEQa2pI{Wq!cKhj@cMX8@ASw-^_Hjwa<Me*h8}`%2E_k1uQ`}DGzj<hJ
z`OIV_IX{?z3XKrS_#5*+7>VGp(KI91PZcaueaGRTwWdRi`xRmM!`j#Y#FDT`T(vf)
z+HD)OqI>h7l69spM8Q3YINi+#*2#hc&p<`de;-?%Lv=ONJv~QkQf;Xmc!ldZQkFBD
z-zUHe*3Hp;O7F2&H3b6)6VSbjVgMv&H@+mcq;nbh1-I8{(rO;F>o22%lxSJwP+Dcm
zQ9G~@FkrZ3@ykcgDVJ*mcsu!NHiPm;r&V5$<VlN$FJ?U@(mM;d;MB_pW_n>i;vrqc
z3v5}+^*!OT0GTn;-kq*sgK`G|o2hT#SARTKi#*=2TM&PjKTRWn=n<|ZmQ`JXRw$w{
zD&Xg3qhC2kfb=1;a7O?P-lKofr(oM#E{rIxTdd}LqbHX)0l6@#AIrn+RAml;Zd(6t
zdWK{d@j}T4i{1*n0%`X4bK8f^GH2;-iiRP3<V!dTn88BvE6GnXG{H{JFm*Q$>WBe*
zGsMq}5~Q|eWHss|7_aB=c;y4zGR9*ayj3<+1)j2nA^B**srknv0ku2P0GjfDJS+U+
zXt9ox7>1F=c=?f#`{ddPGk7no4frWyhtVV-MfdnVVq^kBCjKbUnn{@y;L0VPknh&G
zM&1M~>lCIv<RY0__qeE`*nfN2)10fGII}jg#T7&}`b&1{GmKREbiig_eWQk2!bF~C
z975)^ud2zj#?vsUlHux!w`_>R4m1fwX%DqIaJS3GC}$2uq=a)A=G!;fBKeHL<l9({
zPp|}Lh=o%L23`u*{6bT)jd!9IIEU*_*tAgR`*GRF*~74Bdw0h_ZCz>=Xd|S7b**z5
zet4!D%gQsc`+z6z9@7onM^W%#!Kx?nAo0|{n2CRws{aoHAYAhh!xU2t*u;RZ5Iz5F
zhLf>_xb&2_xhv}zz-|Y?_%Lr`qZ>~=GbhXNDVQh;3_z>#^T%Y@;|^gjRTl%$mpm6Z
zvPRv(6VG^^mKk|j8v<UVo+9SwSkoAWeJ38Dd4Qavv~>uBz1!gV7-x=eqAAN#THT$Z
z<VNSe!O5jHsL#v+n5y4!pR{|ZL)Tu<0#1VC;<7&3!*%j|x0Hoi*Eop88Wuf&2Zobz
zgSfz_2zlQB>d316&8PTz*n7GVq9!qo%z2Ov$6Vt(dGU-aL1dQb11W9H@cpL+0=yxZ
zt4uh`7>wbh0D*`gb7CDb@ae!s_%jUg|MF$H4w-kU=ZM`y=0MP9M`LzVCIzSSF{x@F
z2!=U%SszbA{XAEJFURi_+>$l^cQW}gH3thXKEf?DqE_5dJ+}PXK9AXL0jci?Y=Lq<
zfOPWhE=ysysschG$O0W<%dP+7d?>1<Un&&UWBsm4AI{EFSa!<ZZ6M47tw*|YbSps<
zewIF-v5E*bsaL)%osRIb>0#jks3j7DsbSOOeURARLll#3*FdY|BN|W81G~#2bn$FR
za=9;yo_sIxZ+I!8`d0msR$rq$BNd4mum#Ohtx2fL64HdIsDv^z0r6DJX{#W*_o)bR
ziSQ#0IjCMV(}xh-sc0EK0Grqasz<C4Hk<_WKc@`r!?FhwfS=-2WA1=lbh9Jbj}_b>
zv7f%87;5_Zey8t)K3qDD>7*Xoklv+BW3+R67)V$UZpO_o$VN^MK%X(pt*kwPGnk|m
z7YVc!;ag-Q2T&Pf`tUCBYCyN<?<hg<4)}x{(HvpX)2v3#3tXmkJu1sMGn|UxXMBrk
z<i14c(4)y?|2N#A+lKs?h0#0kZTBX{4=<sW(`i?uVxx)@)mZe-JVPN0D>lVu6t#bB
zF3we&?&-UBt!cvT;)<bARPJUHU1D$B=)xYSg@QjP)+~D|U|YYbhanhms_cOe44lp9
z+t<<4(Scvdm`<|ms3`Jrm{3@^bIF0waYd~m$QxCnJk~d5hjgFUaBG8l16L+hHU0E4
z(dx}+yTIS%{D3itYTl_QA9MwmY-G598}sYkRf9N3eqHm5m*5?K{M+9L#yo%gJq3^!
zvvA;XQ%$$EadbN(*7uP?8AEh0UFw!7_J2_UBVh^lbSWC>Du-g|Dm#J5tBcM;PGs0h
zCE7NR;tF~MzE^IQN-NDpL_>7g(}oZGa!zqst+z{ul(|dBxFvt8;<bDVE2t>oe9C}C
z>)s-}j4UVKA8sdjt><_}dGn@xr&YMj9Uv%wdH3Opl2bV$ga{M!cAn<>FisQ|^z~NS
zOO7wU>(yOL)Y9PS1*DAt8+c*syWwZ#ImS6=omkseJC`Zn=pq6B$NRVFb=?NIEB??3
ztiJ(UN-jdBh6M>_;4NO^NDf9`(nib)N#i}U^7aK4C26H!TT7|t=JMNK&J*i0bjhb(
z0M;xV=qf$vYTnJy)kONDrF-@oNGSI}wR~k`J>tje){!dILcsew@WM&ueq)w4S{(j$
zE^#CQJI;(oExKEG)+P@fvVf1@bP;kad`Gc6=|XS+VfT2i6<F-jthkaq2NIi2-O`J~
ztyJ8+tL0wTM-@Kg%59^u3DO(-C9jw=0s;#~*xV?)fH)gM1p6om#joitdVb%pk>xUe
z@)K_;7A?%`qlcBMy<#hzhu#+U;XMQef%E7}ZZ?6#yP7AE2?vmTS6S}mvu__iOej_+
o#9raBTqYk1%Kr1uX*=pRvX!vS2W~mw{}?n^wGA#6Yr%v615`U3+5i9m

literal 0
HcmV?d00001

diff --git a/man/ggplot_rsi.Rd b/man/ggplot_rsi.Rd
index e2e7cffb..95c6bd6a 100644
--- a/man/ggplot_rsi.Rd
+++ b/man/ggplot_rsi.Rd
@@ -31,7 +31,7 @@ theme_rsi()
 \item{position}{position adjustment of bars, either \code{"stack"} (default) or \code{"dodge"}}
 }
 \description{
-Use these functions to create bar plots for antimicrobial resistance analysis. All functions rely on internal \code{\link{ggplot}} functions.
+Use these functions to create bar plots for antimicrobial resistance analysis. All functions rely on internal \code{\link[ggplot2]{ggplot}} functions.
 }
 \details{
 At default, the names of antibiotics will be shown on the plots using \code{\link{abname}}. This can be set with the option \code{get_antibiotic_names} (a logical value), so change it e.g. to \code{FALSE} with \code{options(get_antibiotic_names = FALSE)}.
@@ -39,13 +39,13 @@ At default, the names of antibiotics will be shown on the plots using \code{\lin
 \strong{The functions}\cr
 \code{geom_rsi} will take any variable from the data that has an \code{rsi} class (created with \code{\link{as.rsi}}) using \code{\link{portion_df}} and will plot bars with the percentage R, I and S. The default behaviour is to have the bars stacked and to have the different antibiotics on the x axis.
 
-\code{facet_rsi} creates 2d plots (at default based on S/I/R) using \code{\link{facet_wrap}}.
+\code{facet_rsi} creates 2d plots (at default based on S/I/R) using \code{\link[ggplot2]{facet_wrap}}.
 
 \code{scale_y_percent} transforms the y axis to a 0 to 100% range.
 
 \code{scale_rsi_colours} sets colours to the bars: green for S, yellow for I and red for R.
 
-\code{theme_rsi} is a \code{\link{theme}} with minimal distraction.
+\code{theme_rsi} is a \code{\link[ggplot2]{theme}} with minimal distraction.
 
 \code{ggplot_rsi} is a wrapper around all above functions that uses data as first input. This makes it possible to use this function after a pipe (\code{\%>\%}). See Examples.
 }
diff --git a/tests/testthat/test-ggplot_rsi.R b/tests/testthat/test-ggplot_rsi.R
new file mode 100644
index 00000000..e46e24d6
--- /dev/null
+++ b/tests/testthat/test-ggplot_rsi.R
@@ -0,0 +1,35 @@
+context("ggplot_rsi.R")
+
+test_that("ggplot_rsi works", {
+
+  skip_if_not("ggplot2" %in% rownames(installed.packages()))
+
+  library(dplyr)
+  library(ggplot2)
+
+  # data should be equal
+  expect_equal(
+    (septic_patients %>% select(amcl, cipr) %>% ggplot_rsi())$data %>%
+      summarise_all(portion_IR) %>% as.double(),
+    septic_patients %>% select(amcl, cipr) %>%
+      summarise_all(portion_IR) %>% as.double()
+  )
+
+  expect_equal(
+    (septic_patients %>% select(amcl, cipr) %>% ggplot_rsi(x = "Interpretation", facet = "Antibiotic"))$data %>%
+      summarise_all(portion_IR) %>% as.double(),
+    septic_patients %>% select(amcl, cipr) %>%
+      summarise_all(portion_IR) %>% as.double()
+  )
+
+  expect_equal(
+    (septic_patients %>% select(amcl, cipr) %>% ggplot_rsi(x = "Antibiotic", facet = "Interpretation"))$data %>%
+      summarise_all(portion_IR) %>% as.double(),
+    septic_patients %>% select(amcl, cipr) %>%
+      summarise_all(portion_IR) %>% as.double()
+  )
+
+  expect_error(geom_rsi(x = "test"))
+  expect_error(facet_rsi(facet = "test"))
+
+})
diff --git a/tests/testthat/test-portion.R b/tests/testthat/test-portion.R
index fec33cc3..f55e4d43 100755
--- a/tests/testthat/test-portion.R
+++ b/tests/testthat/test-portion.R
@@ -1,6 +1,6 @@
 context("portion.R")
 
-test_that("resistance works", {
+test_that("portions works", {
   # amox resistance in `septic_patients`
   expect_equal(portion_R(septic_patients$amox), 0.6603, tolerance = 0.0001)
   expect_equal(portion_I(septic_patients$amox), 0.0030, tolerance = 0.0001)
@@ -46,6 +46,9 @@ test_that("resistance works", {
   expect_warning(portion_S(as.character(septic_patients$amcl)))
   expect_warning(portion_S(as.character(septic_patients$amcl,
                                              septic_patients$gent)))
+  expect_equal(n_rsi(as.character(septic_patients$amcl,
+                                  septic_patients$gent)),
+               1570)
 
 
   # check for errors
@@ -59,6 +62,9 @@ test_that("resistance works", {
   expect_error(portion_S("test", as_percent = "test"))
   expect_error(portion_S(septic_patients %>% select(amox, amcl)))
   expect_error(portion_S("R", septic_patients %>% select(amox, amcl)))
+  expect_error(n_rsi(septic_patients %>% select(amox, amcl)))
+  expect_error(n_rsi(septic_patients$amox, septic_patients %>% select(amox, amcl)))
+
 
   # check too low amount of isolates
   expect_identical(portion_R(septic_patients$amox, minimum = nrow(septic_patients) + 1),
@@ -102,6 +108,15 @@ test_that("old rsi works", {
                            combination_n = n_rsi(cipr, gent)) %>%
                  pull(combination_n),
                c(202, 482, 201, 499))
+
+  # portion_df
+  expect_equal(
+    septic_patients %>% select(amox) %>% portion_df(TRUE) %>% pull(Percentage),
+    c(septic_patients$amox %>% portion_S(),
+      septic_patients$amox %>% portion_I(),
+      septic_patients$amox %>% portion_R())
+  )
+
 })
 
 test_that("prediction of rsi works", {