set up a simulation as demonstrated by John C Mallery

demo_simulation(
  gt.type = "gamma",
  gt.val = c(3, 1.5),
  begin = 1,
  end = 30,
  epid.nb = 1,
  epid.length = 30,
  family = "poisson",
  R0 = 1.67,
  peak.value = 500
)

Arguments

gt.type: character(1) type setting for R0:generation.time
gt.val: numeric() val setting for R0:generation.time
begin: numeric(1) passed to R0::estimate.R
end: numeric(1) passed to R0::estimate.R
epid.nb: numeric(1) passed to R0:sim.epid
epid.length: numeric(1) passed to R0:sim.epid
family: character(1) passed to R0:sim.epid
R0: numeric(1) passed to R0:sim.epid
peak.value: numeric(1) passed to R0:sim.epid

Value

result of estimate.R

Examples

if (FALSE) {
d1 = demo_simulation()
# plot(d1) ## note: R0:::plot.R0.sR is very peculiar and runs dev.new()
plotfit(d1)  # same here but seems to have simpler behavior
# s.a <- sensitivity.analysis(res=d1$estimates$EG, begin=1:15, end=16:30, sa.type="time")
}