# NOT RUN {
## load data on pandemic flu in Baltimore in 1918
data("Flu1918")
## estimate the reproduction number (method "non_parametric_si")
res <- estimate_R(Flu1918$incidence,
method = "non_parametric_si",
config = make_config(list(si_distr = Flu1918$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