This section shows how to run the norovirus model with two levels of vaccination-derived protection against symptomatic disease.

using Epidemics
using Plots

# all arguments have appropriate defaults
data = epidemic_norovirus()

# plot exposed in unvaccinated group
plot(data, vars=(0, 5:8))
Example block output
# plot exposed in single-vaccinated group
plot(data, vars=(0, 25:28))
Example block output
# plot exposed in double-vaccinated group
plot(data, vars=(0, 45:48))
Example block output

Benchmarking

using Epidemics
using BenchmarkTools

@benchmark epidemic_norovirus()
BenchmarkTools.Trial: 11 samples with 1 evaluation.
 Range (min … max):  394.349 ms … 613.455 ms  ┊ GC (min … max): 14.66% … 18.40%
 Time  (median):     432.277 ms               ┊ GC (median):    14.98%
 Time  (mean ± σ):   479.798 ms ±  91.825 ms  ┊ GC (mean ± σ):  15.11% ±  2.16%

  █  ▁ ▁▁   ▁ ▁                                      ▁▁   ▁   ▁  
  █▁▁█▁██▁▁▁█▁█▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁██▁▁▁█▁▁▁█ ▁
  394 ms           Histogram: frequency by time          613 ms <

 Memory estimate: 647.82 MiB, allocs estimate: 6271431.