Loading [MathJax]/jax/output/HTML-CSS/jax.js

Question 6.1: Simulating an SIR process

In the simulated outbreaks, the overall incidence is much too low, and the outbreak dies out after only a few weeks. To attempt to simulate data for which the observed data is a more plausible realization, we might try increasing the force of infection.

Taking it farther….

While this increases the overall incidence, the epidemic is now peaking too quickly. To counteract this, we might try reducing the recovery rate

Additionally, we might have a look at the effects of changing the initial susceptible fraction, η. It’s possible to get something not too awful to contemplate by just manipulating η, in fact:

Question 6.2: Simulating an SEIR process

The existing code may be modified as follows:

Using the rough estimate that the latent period in measles is 810 days, we take μEI0.8wk1 and μIR1.3wk1 (so as to have roughly the same generation time as before).

Again one can increase the force of infection:

Now increase β to, by trial and error, 37.

Notice that it seems that the simulated time series peak too late. We can decrease η back to 0.03, and increase β to 73.

Despite some uncertainties, this seems like a reasonable model for the time series. So we can choose β=73,μEI=0.6,μIR=1,ρ=0.7,η=0.03 for the SEIR model.


References

1. King, A.A., Ionides, E.L., and Lin, Q. (2021). Simulation-based inference for epidemiological dynamics. Summer Institute in Statistics and Modeling in Infectious Diseases. Available at: https://kingaa.github.io/sbied/.