Plot sequestration rate in addition to probability

vignettes/articles/model-description.Rmd

@@ -744,13 +744,24 @@ We denote the number of iRBCs sequestered in a time-step as \(I_{c \to q}\):
 | \(\zeta_{\mathrm{sl}}\) | Slope parameter  | \(`r p$zeta_sl`\)        |
 | \(\zeta_{50}\)          | Half-maximal age | \(`r p$zeta_a50`\) hours |
 
-```{r zeta, class.source = 'fold-hide', echo = ! knitr::is_latex_output(), fig.cap = 'The age-specific probability of iRBC sequestration.'}
+```{r zeta, class.source = 'fold-hide', echo = ! knitr::is_latex_output(), fig.cap = 'The age-specific rate \\(\\zeta(a)\\) of iRBC sequestration.'}
 p_pv <- baseline_parameters("Pv")
 df_zeta <- data.frame(
   age = rep(1:p$T_irbc, 2),
   prob = 1 - exp(- c(p$zeta, p_pv$zeta)),
+  rate = c(p$zeta, p_pv$zeta),
   species = rep(c("Pf", "Pv"), each = p$T_irbc)
 )
+ggplot(df_zeta, aes(age, rate, colour = species)) +
+  geom_line() +
+  geom_vline(xintercept = p$zeta_a50, linetype = "dashed") +
+  expand_limits(y = 0) +
+  xlab("iRBC age (hours)") +
+  ylab("Sequestration rate") +
+  scale_colour_hue(name = "Species")
+```
+
+```{r zeta-prob, class.source = 'fold-hide', echo = ! knitr::is_latex_output(), fig.cap = 'The age-specific probability of iRBC sequestration.'}
 ggplot(df_zeta, aes(age, prob, colour = species)) +
   geom_line() +
   geom_vline(xintercept = p$zeta_a50, linetype = "dashed") +