What to do
The tasks below are described in a way that assumes that 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.
Task 1:
- Set the model parameters such that it corresponds to the following setting:
- 1000 initially susceptible individuals, S0, 1 initially infected host, I0, no pathogen in the environment (E0 = 0) , simulation duration, tmax, approximately 5 years.
- We’ll first look at direct transmission. Set transmission rate of bD = 0.001, environmental transmission bE = 0.
- Set environmental pathogen shedding, p, and clearance, c, to 0.
- Assume that the duration of the infectious period, g, is about 2 months long.
- Turn off births, m, and deaths, n, for now.
- Run the simulation. You should get the usual direct transmission dynamics and a single outbreak.
- Record the number and fraction of susceptible/infected/recovered remaining at the end of the outbreak. You should be left with around 202 susceptibles at the end of the simulation. Remember that in a model like this, where the only flow from the susceptible class is outflow through becoming infected, everyone who is not in the susceptible compartment at the end has been infected. Thus the total number of infected during the outbreak is the different between susceptibles at the beginning and at the end.
Task 2:
- Now try various values for the rate of shedding into the environment, p, and environmental clearance c. Leave everything else unchanged.
- As those parameters move away from 0, what do you expect to see?
- Think about your expectations, run the simulation and compare your expectations with the results.
- Anything surprising happening? Do you understand why you see what you see?
Task 3:
- Now set p = 10 and c = 10 . Turn off direct transmission. Run to make sure you don’t get an outbreak.
- Turn on environmental transmission, set it to bE = 0.001. Leave everyhing else unchanged.
- What do you expect to see? Run simulations and compare with the direct transmission results.
Task 4:
- Now also turn on direct transmission by setting bD = bE = 0.001. Leave everyhing else unchanged.
- What do you expect to see? Run simulations and see how your expectations hold up.
Task 5:
- Now set the initial number of infected to 0 and initial pathogen in the environment to 1. Leave everyhing else unchanged.
- What do you expect to see? Run simulations and see how your expectations hold up.
Task 6:
- Try different combinations of direct and environmental transmission only and having only infected hosts or only pathogen in the environment at the start.
- What do you expect to see? Run simulations and see how your expectations hold up.
Task 7:
- Keep exploring by trying different parameters and transmission settings and see how they influence results.
- You can also turn on births/deaths and see what you get.
- As you continue your exploration, think about potential real infectious diseases that might be approximated by either one of the transmission types, and what settings for other model parameters would be for those ID.