What to do
The tasks below are described in a way that assumes everything is in units of MONTHS (rate parameters, therefore, have units of inverse months). If any quantity is not given in those units, you need to convert it first (e.g. if it says a year, you need to convert it to 12 months).
Some of the simulations might take a few seconds to run. Be patient.
Task 1:
- Set the model parameters such that it corresponds to the following setting:
- 1000 susceptible and 1 initially symptomatic host.
- Simulation duration approximately 5 years.
- Assume that only symptomatic individuals transmit, at rate 0.002.
- Assume that the duration of the symptomatic period is 1 month long, the duration of the presymptomatic period is half a month long.
- Assume immunity does not wane and that there are no births and deaths.
- Set the number of simulations to 1.
- With parameters set to correspond to the scenario just described, run the simulation.
- If you rerun the simulation without changing anything, what do you expect to see? Test your expectations by running the simulation several times.
- To speed things up, set the number of simulations to 20, run all 20.
Task 2:
- In the previous task, you should have seen that occasionally, no outbreak occurs (because the infected person recovers/dies before they can infect others). Explore how the initial number of infected influences the chance of (not) getting outbreaks. Set the initial number of infected to 5, everything as above.
- Run 20 outbreaks at a time, repeat a few times. Observe the number of times you do (not) get an outbreak.
- Repeat with the initial number of infected at 2,3, and 10.
Task 3:
- Set everything as in task 1, with 1 initial infected. Now explore how the rate of transmission (and with it the reproductive number) affects outbreaks.
- Try values for the transmission rate of 0.001, 0.0015 and 0.004. Compare the results.
Task 4:
- Set susceptibles to 400, infected to 100, the rate of immunity loss to 0.2, the rate of symptomatic transmission to 0.005, and no presymptomatic transmission.
- Set the duration of the symptomatic period to 1 month and the duration of the presymptomatic period to half a month.
- Simulation time about 5 years.
- No births and deaths.
- Set the number of simulations to 1. Run the model several times.
- You should see the stochastic version of a steady state, namely S-E-I-R values that fluctuate around some level.
Task 5:
- Record the fraction of susceptibles at the end of the simulation for the setting you just had. They will change.
- Now change the number of simulations to 5, 10, 20 (be patient). See how the average number of susceptibles at the end of the simulation settles down to some number. It will still fluctuate a bit, but as you take the average of more simulations, results become more stable (deterministic).
Task 6:
- Set simulation number to 5, rest as before.
- Run the model a few times for around 50 years (be patient). Observe if/when you get extinctions.
- Now change the rate of transmission to 0.003, again run multiple simulations, then repeat with transmission at 0.0025.
- Contemplate why you see what you see.
Task 7:
- Keep exploring. Turn on births and deaths. Play around with anything else you can think of.