Name: Diogo Cosin Ayres de Oliveira
Course: Data Engineering Masters
First of all, let's see the SIR model time reponse with the following default parameters:
Parameter | Value |
---|---|
c | 1.0 |
a | 0.025 |
0.03 | |
d | 0.1 |
0.3 | |
0.15 | |
0.25 |
Plotting the first 60 days of the SRI model time response, we get the following graph.
The graph show that after approximately 20 days, almost all 100 individuals will be dead. This may be explained by the fact that the death rate of healthy individuals d
is 0.1 per day, that is, 10% of all healthy individuals will die. At the beggining, this is equivalent to 8 individuals.
This extreme scenario would be a permanent state, that is, no changes on individuals dying, getting infected or recovered and immigrating since all parameters are set to zero. To check, let's see the time response plot.
As expected, all individuals levels remain the same as the initial conditions.
By default the c
is set to 1, meaning that every susceptible will necessarily be in contact with an infected individual. Let's change this parameter to 0.5.
Now the populations takes slightly longer to decrease, but it is still almost completely eliminated as it happens using c =
1.0 and the default parameters.
Finally let's check the extreme case where no susceptible individual will be in contact with infected individuals, in other words, c =
0.
In this case, as susceptible people are not getting infected, they will die according to the d
parameters and not to
Let's check how the probability of transmission of the disease upon contact, a
, affects the SIR model time reponse. For this, let's set the extreme value 1 to it, meaning that every individual who contacted to an infected individual will certainly get infected as well. The other parameters are set to the default values.
As expected, right at the first day, all susceptible individuals get infected, elevating the infected population to more than 80 individuals. After that, with 25% of the infected population dying each day, it takes less than 25 days to almost completely eliminate the whole population.
Let's check what it the model's behaviour when all infected individuals get recovered, that is, d
parameter.
Indeed, there are more recovered individuals, but still the population is almost completely eliminated after 60 days.
Now let's check what happens when the individual cannot be recovered,
In this scenario, there are no recovered individuals besides the five intials ones. Also, the population is eliminated faster within 20 days.
Let's imagine that humanity has evolved to a level where only 0.1% of healthy people die. That way we would obtain the following time response in the SIR model.
That's an interisting result. After around 30 days, when there is no longer infected individuals as they all died, the population begins to grow again. This may explained by the fact that there a new people coming according to the d
was set to a low value.
Let's consider that we are modelling a disease which it is not possible to get infected twice, that is, let's set
Now that recovered individual can't be infected again, their population increases and dies according to the d
parameter.
In mass immigration cases,
In that case, it is possible to notice that immigrants are feeding the disease to a point that the number of individuals stabilize after 20 days with the majority infected. In other words, the individuals are dying, but, as there is always new individuals coming in, the populations is not eliminated.
Now let's play with an unrealistic situation: infected individuals acquire superpowers and, because of that, they do not die. In the SIR model, this means setting
Infected individuals do not die, but they are still recovered according to the c
parameter.
Let's not allow individuals to be recovered setting
Accordingly to what was previously discussed, now that there is no recover chance to infected individuals and they don't die, the infected population grows to infinite due to the immigated new susceptible individuals while the healthy population is completely eliminated in less than 10 days.