# NOT RUN {
## load data on measles in Hallegoch in 1861
data("Measles1861")
## estimate the reproduction number (method "non_parametric_si")
res <- estimate_R(Measles1861$incidence, method="non_parametric_si",
config = make_config(list(
t_start = seq(17, 42),
t_end = seq(23, 48),
si_distr = Measles1861$si_distr)))
plot(res)
## the second plot produced shows, at each each day,
## the estimate of the reproduction number
## over the 7-day window finishing on that day.
# }
Run the code above in your browser using DataLab