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")
Run the code above in your browser using DataLab