| Title: | Impact Study of Vaccination Campaigns |
|---|---|
| Description: | Tools to estimate the impact of vaccination campaigns at population level (number of events averted, number of avertable events, number needed to vaccinate). Inspired by the methodology proposed by Foppa et al. (2015) <doi:10.1016/j.vaccine.2015.02.042> and Machado et al. (2019) <doi:10.2807/1560-7917.ES.2019.24.45.1900268> for influenza vaccination impact. |
| Authors: | Yohann Mansiaux [aut, cre], Alexandre Blake [aut], James Humphreys [aut], Baltazar Nunes [aut] |
| Maintainer: | Yohann Mansiaux <[email protected]> |
| License: | MIT + file LICENSE |
| Version: | 0.1.0 |
| Built: | 2026-06-07 08:16:23 UTC |
| Source: | https://github.com/epiconcept-paris/vaccinationimpact |
Compute events averted by increasing the final vaccine coverage
compute_events_avertable_by_increasing_coverage( number_of_events, cumulative_coverage, vaccine_coverage_increase, vaccine_effectiveness )compute_events_avertable_by_increasing_coverage( number_of_events, cumulative_coverage, vaccine_coverage_increase, vaccine_effectiveness )
number_of_events |
number of events |
cumulative_coverage |
cumulative vaccination coverage |
vaccine_coverage_increase |
percentage increase in final vaccine coverage (between 0 and 1) |
vaccine_effectiveness |
vaccine effectiveness |
a list with the new vaccine coverage ("new_vaccine_coverage") and the estimated number of events averted ("nabe")
data(coverage_and_incidence_mock_data) data(ve_mock_data) coverage <- coverage_and_incidence_mock_data$coverage_data incidence <- coverage_and_incidence_mock_data$incidence_data vaccine_effectiveness <- ve_mock_data$ve nabe <- compute_events_avertable_by_increasing_coverage( number_of_events = incidence$events, cumulative_coverage = coverage$cumulative_coverage, vaccine_coverage_increase = 0.1, # 10% increase in final coverage vaccine_effectiveness = vaccine_effectiveness ) plot(nabe$new_vaccine_coverage, type = "l", xlab = "Time", ylab = "Vaccine coverage with 10% increase") plot(nabe$nabe, type = "l", xlab = "Time", ylab = "Events averted")data(coverage_and_incidence_mock_data) data(ve_mock_data) coverage <- coverage_and_incidence_mock_data$coverage_data incidence <- coverage_and_incidence_mock_data$incidence_data vaccine_effectiveness <- ve_mock_data$ve nabe <- compute_events_avertable_by_increasing_coverage( number_of_events = incidence$events, cumulative_coverage = coverage$cumulative_coverage, vaccine_coverage_increase = 0.1, # 10% increase in final coverage vaccine_effectiveness = vaccine_effectiveness ) plot(nabe$new_vaccine_coverage, type = "l", xlab = "Time", ylab = "Vaccine coverage with 10% increase") plot(nabe$nabe, type = "l", xlab = "Time", ylab = "Events averted")
Compute events averted by vaccination
compute_events_averted_by_vaccination( number_of_events, cumulative_coverage, vaccine_effectiveness )compute_events_averted_by_vaccination( number_of_events, cumulative_coverage, vaccine_effectiveness )
number_of_events |
number of events |
cumulative_coverage |
cumulative vaccination coverage |
vaccine_effectiveness |
vaccine effectiveness |
The number of events averted by vaccination is calculated as described by Machado et al. (2019) https://doi.org/10.2807/1560-7917.ES.2019.24.45.1900268.
estimated number of events averted
data(coverage_and_incidence_mock_data) data(ve_mock_data) coverage <- coverage_and_incidence_mock_data$coverage_data incidence <- coverage_and_incidence_mock_data$incidence_data vaccine_effectiveness <- ve_mock_data$ve nae <- compute_events_averted_by_vaccination( number_of_events = incidence$events, cumulative_coverage = coverage$cumulative_coverage, vaccine_effectiveness = vaccine_effectiveness ) plot(nae, type = "l", xlab = "Time", ylab = "Events averted")data(coverage_and_incidence_mock_data) data(ve_mock_data) coverage <- coverage_and_incidence_mock_data$coverage_data incidence <- coverage_and_incidence_mock_data$incidence_data vaccine_effectiveness <- ve_mock_data$ve nae <- compute_events_averted_by_vaccination( number_of_events = incidence$events, cumulative_coverage = coverage$cumulative_coverage, vaccine_effectiveness = vaccine_effectiveness ) plot(nae, type = "l", xlab = "Time", ylab = "Events averted")
Compute the number of individuals needed to vaccinate to prevent one event according to Machado et al. method
compute_number_needed_to_vaccinate_machado( number_of_events, number_of_events_averted, population_size, vaccine_effectiveness )compute_number_needed_to_vaccinate_machado( number_of_events, number_of_events_averted, population_size, vaccine_effectiveness )
number_of_events |
number of events |
number_of_events_averted |
number of events averted |
population_size |
population size |
vaccine_effectiveness |
vaccine effectiveness |
The number of individuals needed to vaccinate to prevent one event is calculated as described by Machado et al. (2019) https://doi.org/10.2807/1560-7917.ES.2019.24.45.1900268.
The number of individuals needed to vaccinate to avert one event
data(coverage_and_incidence_mock_data) data(ve_mock_data) coverage <- coverage_and_incidence_mock_data$coverage_data incidence <- coverage_and_incidence_mock_data$incidence_data vaccine_effectiveness <- ve_mock_data$ve nae <- compute_events_averted_by_vaccination( number_of_events = incidence$events, cumulative_coverage = coverage$cumulative_coverage, vaccine_effectiveness = vaccine_effectiveness ) nnv_machado <- compute_number_needed_to_vaccinate_machado( number_of_events = incidence$events, number_of_events_averted = nae, population_size = 1234, vaccine_effectiveness = vaccine_effectiveness ) nnv_machadodata(coverage_and_incidence_mock_data) data(ve_mock_data) coverage <- coverage_and_incidence_mock_data$coverage_data incidence <- coverage_and_incidence_mock_data$incidence_data vaccine_effectiveness <- ve_mock_data$ve nae <- compute_events_averted_by_vaccination( number_of_events = incidence$events, cumulative_coverage = coverage$cumulative_coverage, vaccine_effectiveness = vaccine_effectiveness ) nnv_machado <- compute_number_needed_to_vaccinate_machado( number_of_events = incidence$events, number_of_events_averted = nae, population_size = 1234, vaccine_effectiveness = vaccine_effectiveness ) nnv_machado
Compute the number of individuals needed to vaccinate to prevent one event according to Tuite and Fisman method
compute_number_needed_to_vaccinate_tuite_fisman( number_of_vaccinated, number_of_events_averted )compute_number_needed_to_vaccinate_tuite_fisman( number_of_vaccinated, number_of_events_averted )
number_of_vaccinated |
number of vaccinated individuals |
number_of_events_averted |
number of events averted |
The number of individuals needed to vaccinate to prevent one event is calculated as described by Tuite and Fisman (2013) https://doi.org/10.1016/j.vaccine.2012.11.097.
The number of individuals needed to vaccinate to avert one event
data(coverage_and_incidence_mock_data) data(ve_mock_data) coverage <- coverage_and_incidence_mock_data$coverage_data incidence <- coverage_and_incidence_mock_data$incidence_data vaccine_effectiveness <- ve_mock_data$ve nae <- compute_events_averted_by_vaccination( number_of_events = incidence$events, cumulative_coverage = coverage$cumulative_coverage, vaccine_effectiveness = vaccine_effectiveness ) nnv_tuite_fisman <- compute_number_needed_to_vaccinate_tuite_fisman( number_of_vaccinated = coverage$number_of_vaccinated, number_of_events_averted = nae ) nnv_tuite_fismandata(coverage_and_incidence_mock_data) data(ve_mock_data) coverage <- coverage_and_incidence_mock_data$coverage_data incidence <- coverage_and_incidence_mock_data$incidence_data vaccine_effectiveness <- ve_mock_data$ve nae <- compute_events_averted_by_vaccination( number_of_events = incidence$events, cumulative_coverage = coverage$cumulative_coverage, vaccine_effectiveness = vaccine_effectiveness ) nnv_tuite_fisman <- compute_number_needed_to_vaccinate_tuite_fisman( number_of_vaccinated = coverage$number_of_vaccinated, number_of_events_averted = nae ) nnv_tuite_fisman
Coverage and incidence mock data. Coverage values are computed considering a sample size of 1234 individuals.
coverage_and_incidence_mock_datacoverage_and_incidence_mock_data
A list with two data frames:
data.frame with weekly incidence data
data.frame with weekly coverage data
Simulated coverage and incidence data
Vaccine effectiveness data.
ve_mock_datave_mock_data
A data frame with 52 rows and 2 variables:
Date
numeric: weekly vaccine effectiveness
Simulated vaccine effectiveness data