library(daedalus)
This vignette is intended as a visual check on daedalus features in development.
output <- daedalus2(100, n_strata = 3) plot(colSums(output$E))
plot(colSums(output$I))
This should be a curve and not monotonically increasing.
plot(colSums(output$cases_inc))