This section shows how to run the Daedalus model.

using Epidemics
using Plots

# all arguments have appropriate defaults
data = epidemic_daedalus(time_end=600.0)

# plot exposed group
plot(data, vars=(0, 50:99))
Example block output
# plot exposed among vaccinated - should start at vax_time = 200.0
plot(data, vars=(0, 394:442))
Example block output

Benchmarking

# benchmark for a typical daedalus run of 600 days
using Epidemics
using BenchmarkTools

@benchmark epidemic_daedalus(time_end = 600.0)
BenchmarkTools.Trial: 25 samples with 1 evaluation.
 Range (minmax):  140.335 ms834.293 ms   GC (min … max): 16.82% … 83.49%
 Time  (median):     154.677 ms                GC (median):    17.91%
 Time  (mean ± σ):   201.168 ms ± 136.237 ms   GC (mean ± σ):  28.36% ± 15.21%

  █     ▃▅▃▁▆▃▃▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▃ ▁
  140 ms           Histogram: frequency by time          834 ms <

 Memory estimate: 263.37 MiB, allocs estimate: 1464186.