loading package
library (R0)
## Loading required package: MASS
load data
data(Germany.1918)
plot graph
plot(Germany.1918)
Plot log
plot(log(Germany.1918), main = "Exponential growth rate", xlab = " Time", ylab = "log(Germany.1918)")
generation time distribution
GT = generation.time ("gamma", c(3,1.5))
estimate Ro using exponential growth method using start date of first day and end dat 0 to 30
est.R0.EG(Germany.1918,GT,begin=1,end=30)
## Waiting for profiling to be done...
## Reproduction number estimate using Exponential Growth method.
## R : 1.43944[ 1.416545 , 1.462916 ]
Description: Reproductive number by estimation of exponential growth of 1918 influenza pandemic in Germany is 1.44