From 052ae7a0c06fc98a9d62751751b350b2b5b0227d Mon Sep 17 00:00:00 2001 From: Andree Valle Campos Date: Tue, 21 Jul 2026 23:53:16 +0100 Subject: [PATCH 1/2] Add callout defining vaccination rate before vaccination() model call Learners see `nu` used in epidemics::vaccination() without a definition; this adds the formula and a worked example (age 65+, Feb 2021) just before the object is passed into model_default(). Co-Authored-By: Claude Sonnet 5 --- episodes/modelling-interventions.Rmd | 17 +++++++++++++++++ 1 file changed, 17 insertions(+) diff --git a/episodes/modelling-interventions.Rmd b/episodes/modelling-interventions.Rmd index 2327cc40..0d01cf10 100644 --- a/episodes/modelling-interventions.Rmd +++ b/episodes/modelling-interventions.Rmd @@ -431,6 +431,23 @@ vaccinate <- epidemics::vaccination( ) ``` +::::::::::::::::::::::::::::::::::::: callout +### Vaccination rate + +The vaccination rate $\nu_{g_i}$ is the number of doses given daily to age group $g_i$ among the total population in that age group: + +$$ +\nu_{g_i} = \frac{\text{daily doses in group } g_i}{\text{population in group } g_i} +$$ + +For example, in February 2021, 425,000 daily doses were administered to the population aged 65+, with a population size of 12,730,000: + +$$ +\nu_{65+} = \frac{425{,}000}{12{,}730{,}000} = 0.0334 +$$ + +:::::::::::::::::::::::::::::::::::::::::::::::: + We pass our vaccination object into the model using the argument `vaccination = vaccinate`: ```{r output_vaccinate} From bc3503e67342aeec25ef1466401521f1b1e9f52d Mon Sep 17 00:00:00 2001 From: Andree Valle Campos Date: Wed, 22 Jul 2026 12:32:22 +0100 Subject: [PATCH 2/2] update values to reference Co-authored-by: Andree Valle Campos --- episodes/modelling-interventions.Rmd | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/episodes/modelling-interventions.Rmd b/episodes/modelling-interventions.Rmd index 0d01cf10..c5ab858d 100644 --- a/episodes/modelling-interventions.Rmd +++ b/episodes/modelling-interventions.Rmd @@ -440,12 +440,14 @@ $$ \nu_{g_i} = \frac{\text{daily doses in group } g_i}{\text{population in group } g_i} $$ -For example, in February 2021, 425,000 daily doses were administered to the population aged 65+, with a population size of 12,730,000: +For example, during the NHS England autumn 2023 COVID-19 booster campaign, an average of 98,000 daily doses were administered to the population aged 65+, with a population size of 12,730,000: $$ -\nu_{65+} = \frac{425{,}000}{12{,}730{,}000} = 0.0334 +\nu_{65+} = \frac{98{,}000}{12{,}730{,}000} = 0.0077 $$ +61.5% of people aged 65+ had received an autumn COVID-19 booster by 30 November 2023, since the campaign began on 11 September 2023 (80 days). 61.5% × 12,730,000 ≈ 7,829,000 doses ÷ 80 days ≈ 98,000 doses/day. ([NHS England, Vaccinations: COVID-19](https://www.england.nhs.uk/statistics/statistical-work-areas/covid-19-vaccinations/)) + :::::::::::::::::::::::::::::::::::::::::::::::: We pass our vaccination object into the model using the argument `vaccination = vaccinate`: