# NOT RUN {
## load data on smallpox in Kosovo in 1972
data("Smallpox1972")
## estimate the instantaneous reproduction number (method "NonParametricSI")
EstimateR(Smallpox1972$Incidence, T.Start=27:51, T.End=33:57, method="NonParametricSI",
SI.Distr=Smallpox1972$SI.Distr, plot=TRUE, leg.pos=xy.coords(1,15))
# the second plot produced shows, at each each day,
# the estimate of the instantaneous reproduction number over the 7-day window finishing on that day.
# }
Run the code above in your browser using DataLab