From united nations- world health organization website https roberts
Modelling Ebola
with the SIR Model1 Modelling Ebola with the SIR Model
When I was younger, my mother would come home detailing the mathematical concepts behind her current projects at the CDC (Center for Disease Control). But as a child I did not realize the concepts, let alone math behind her explanations. It was not until my highschool biology classes where I started to understand
epidemiology- the study of diseases- and the importance of the mathematical models my mom analyzes and creates. My biological interest in the field rose along with my understanding of calculus, until I finally found a fascination with such graphs and charts displaying the progression of diseases globally.Specifically, I began to study the SIR Model which I will focus on in this investigation to show the behavior of disease and how to avoid infection.
Rt = the population of those who have recovered
from the illness meaning they cannot spread
or contract it once again. This includes those
who are deceasedBecause the SIR Model only accounts for a proportion of the entire population as it only models a closed system (a constant population), only a fraction of those susceptible to infection can be modelled. This will be denoted by P and in relation to the set of variables above are as follows:
rt = P Rt Part of the population
who is recovered
or removed from the disease
These assumptions together lessen the external probabilities of the unpredictable factors pertaining to the spread of disease.
3. The Math Behind the Model
product of | St and | S |
---|
The susceptible population with regards to time is:
st+1 = st − β t tTherefore the derivation of the differential equation for the susceptible class is:
st+1 − st = st − β t t − st
dt= d i |
|
---|
3. Deriving the Differential Equation for the Recovered/ Removed Class The recovered/ removed population is the last category of the SIR model and accounts for those who have fully recovered from the disease in addition to those who have died from the disease. But both groups have no possibility of spreading or contracting the disease. Therefore this class has two components: those who are already in the population and those who are newly recovered as they have just passed the infection stage. The
recovered/ removed population with regards to time is:
rt+1 = rt + Φ t
And the respective differential equations are:
dt= d s |
|
||
---|---|---|---|
dt= d i | |||
dt d r | = |
population as the net population seen here is zero meaning no person in the community enters or exits.
4. Modelling SIR with an Ebola in Guinea
= N = 12,000,000 people S0
There no are people initially infected because the disease is introduced: I0 =0 people
Therefore the Ebola virus in Guinea is modeled by the equations:
dt= d s |
|
---|
derived differential equations. From these equations, a graph modelling the Ebola virus through SIR was developed as seen in Figure 1 from data in the appendix.
see how the susceptible population is decreasing while the infected population increases
showing more people are becoming exposed to the virus and the spread of the virus is
In all, I gained a greater appreciation for the concepts highlighted in my exploration and applying math to real life situations such as the spread of disease. I also learned how to calculate and quantify gathered data and interpret critical points on graphs displaying
10 such data. Modeling SIR for diseases such as the present Coronavirus and Influenza will definitely help in efforts to eradicate.
12
7. Appendix 2: Raw Data for SIR Model of Ebola
Time in Days | Susceptible (st) | Infected (it) | Recovered (rt) |
---|---|---|---|
1 | 1.000 | 0.008 | 0.008 |
2 | 1.000 | 0.009 | 0.009 |
3 | 1.000 | 0.008 | 0.008 |
4 | 1.000 | 0.009 | 0.009 |
5 | 1.000 | 0.009 | 0.009 |
6 | 1.000 | 0.009 | 0.009 |
7 | 1.000 | 0.009 | 0.009 |
8 | 1.000 | 0.009 | 0.009 |
9 | 1.000 | 0.009 | 0.009 |
10 | 1.000 | 0.009 | 0.009 |
11 | 1.000 | 0.009 | 0.009 |
12 | 1.000 | 0.009 | 0.009 |
13 | 1.000 | 0.009 | 0.009 |
14 | 1.000 | 0.009 | 0.009 |
15 | 1.000 | 0.009 | 0.009 |
16 | 1.000 | 0.009 | 0.009 |
17 | 1.000 | 0.009 | 0.009 |
18 | 1.000 | 0.009 | 0.009 |
19 | 1.000 | 0.009 | 0.009 |
20 | 1.000 | 0.009 | 0.009 |
21 | 1.000 | 0.009 | 0.009 |
22 | 1.000 | 0.009 | 0.009 |
23 | 0.989 | 0.009 | 0.009 |
24 | 0.986 | 0.010 | 0.010 |
25 | 0.980 | 0.013 | 0.010 |
55 | 0.009 | 0.090 | 0.877 |
---|---|---|---|
56 | 0.009 | 0.083 | 0.886 |
57 | 0.009 | 0.079 | 0.894 |
58 | 0.009 | 0.078 | 0.902 |
59 | 0.009 | 0.073 | 0.910 |
60 | 0.009 | 0.065 | 0.919 |
61 | 0.009 | 0.059 | 0.933 |
62 | 0.009 | 0.055 | 0.942 |
63 | 0.009 | 0.051 | 0.952 |
64 | 0.009 | 0.042 | 0.961 |
65 | 0.009 | 0.038 | 0.962 |
66 | 0.009 | 0.035 | 0.962 |
67 | 0.009 | 0.034 | 0.963 |
68 | 0.009 | 0.032 | 0.964 |
69 | 0.009 | 0.029 | 0.969 |
70 | 0.009 | 0.026 | 0.977 |
71 | 0.009 | 0.024 | 0.977 |
72 | 0.009 | 0.021 | 0.977 |
73 | 0.009 | 0.018 | 0.978 |
74 | 0.009 | 0.016 | 0.979 |
75 | 0.009 | 0.014 | 0.979 |
76 | 0.009 | 0.013 | 0.977 |
77 | 0.009 | 0.013 | 0.980 |
78 | 0.009 | 0.012 | 0.980 |
79 | 0.009 | 0.012 | 0.980 |
80 | 0.009 | 0.012 | 0.982 |
81 | 0.009 | 0.011 | 0.982 |
82 | 0.009 | 0.011 | 0.983 |
83 | 0.009 | 0.010 | 0.984 |