Skip to main content

Estimation of SIR model’s parameters of COVID-19 in Algeria

Abstract

Background

The coronavirus disease 2019 (COVID-19) is reported in Algeria on February 25th, 2020. Since then, the number is still increasing leading to a total number of 36,699 cases and 1333 deaths on August 12th, 2020. Thus, comprehension of the epidemic curve is very important to predict its evolution and subsequently adapt the best prevention strategies. In this way, the current study was conducted to estimate the parameters of the classical SIR model and to predict the peak of the COVID-19 epidemic in Algeria using data from February 25th, 2020 to August 12th, 2020.

Results

Results showed that the peak of the epidemic will be reached on September 8th, 2021 and the total infected persons will exceed 800,000 cases at the end of the epidemic. Also, more than 15 million persons will be susceptible. The reproduction number (R0) is estimated at 1.23254.

Conclusion

These results may be helpful for the Algerian authorities to adapt their strategies and may be taken into consideration in the future phase of discontainment.

Background

The coronavirus disease (COVID-19) is a respiratory viral disease first reported in Wuhan City, in China on December 2019 and turned epidemic worldwide (Lounis 2020). This disease has affected more than 21 million persons through the world leading to more than 6.000 deaths (JHUM 2020). COVID-19 has been declared in Algeria on February 25th, 2020 in the South department of Ouargla, and it continues to spread in the country afflicting a total number of 36,699 cases on August 12th, 2020. The total positive cases, deaths and recovered persons in Algeria are shown in Fig. 1 (Algerian Heath Ministry 2020).

Fig. 1
figure 1

Total positive cases, deaths and recovered cases in Algeria

Regarding daily reported cases as shown in Fig. 2, we remark that the number of cases seemed to decrease after a peak in late April to middle May. This may be due to the measures implemented since March 12th, 2020. These interventions consisted of border banning, school closure, social distancing and cancellation of all public gatherings (Lounis 2020).

Fig. 2
figure 2

COVID-19 cases of Algeria as per time series from 25th February 2020

However, since middle June a considerable increase in daily reported cases was observed. Increased numbers have been reported due to different factors including the increase in daily performed tests but also the remarkable relaxation and lighting in prevention measures associated with the alleviation of the containment measures. This number is still increasing despite the strengthening of these measures in some departments. In this way, the prediction of COVID-19 evolution is of great importance for providing information for the public authorities that could help them to adopt or readapt their prevention measures. Multiple forecasting methods were used to analyze and predict the future trends of COVID-19 such as the logistic growth model (Zou et al. 2020), stochastic susceptible–infectious–removed (SIR) model (Bagal et al. 2020), susceptible–exposed–infectious–removed (SEIR) model (Al-Sayed et al. 2020) and natural growth model (Huang and Qiao 2020).

In this work, the classical SIR model has been utilized to predict the number of case counts, the rate of transmission, the reproduction number (R0), the size and the extinction date of the COVID-19 epidemic in Algeria.

Methods

Estimation of SIR model’s parameters

In the current study, the standard susceptible–infectious–recovered (SIR) model was adopted. In this model, we divide the total population (N) into three categories: Susceptible (S), Infected (I) and Recovered (R). The susceptible is the part of the total population which is vulnerable and is at risk of being infected (if previously unexposed to the pandemic disease); the infected is the fraction that has been infected (currently colonized by the pandemic disease), and the recovered is the fraction of the total population having recovered/removed (either by death or recovery). The basic structure of SIR mathematical model is described in Fig. 3 (Bagal et al. 2020).

Fig. 3
figure 3

Scheme of the fraction of the SIR model

The specific description of the SIR model is as follows:

$$\frac{{{\text{d}}S}}{{{\text{d}}t}} = - \frac{\beta }{N}{\text{SI }}$$
(1)
$$\frac{{{\text{d}}I}}{{{\text{d}}t}} = \frac{\beta }{N}SI - \gamma I$$
(2)
$$\frac{{{\text{d}}R}}{{{\text{d}}t}} = \gamma I$$
(3)

where β represents the effective transmission rate and γ represents the removal or the recovery rate. γ is defined as the inverse of the duration of recovery d (γ = 1/d). Knowing that the total population size (N = S + I + R) is time independent. The population of Algeria is estimated at 43,851,044 inhabitants in the middle of the year according to United Nations reports (Worldmeter 2020) which is the (N) value for the SIR modeling.

Initially, in the absence of infection we have I + R = 0 and S≈N. We obtain the following equation from Eq. (2):

$$\frac{{{\text{d}}I}}{{{\text{d}}t}}\sim I\left( {\beta - \gamma } \right)$$
(4)

Then, integration of Eq. (4) gives the Equation:

$$I = I_{0} e^{{\left( {\beta - \gamma } \right)t}}$$
(5)

Determination of transmission rate (β) and the m (m = βγ) value

When the infection occurs, the susceptible population is almost equal to the total population (S ≈ N). Subsequently, the number of infectious individuals I(t) first increases exponentially.

$$\frac{{{\text{d}}I}}{{{\text{d}}t}}\sim mI$$
(6)

m is the difference between transmission and recovery rates (m = βγ.

$$I\left( t \right) = I_{0} e^{mt}$$
(7)
$${\text{In }}I = mt + {\text{In}} I_{0}$$
(8)

We can calculate the value of m from the log-plot data and for the best-line fit we can use, for example, the least squares.

In the current work, the used data file was considered from the JHU-CSSE (Johns Hopkins University Center for Systems Science and Engineering) 2019 Novel Coronavirus Visual Dashboard database (JHU-CSSE 2020). The simulation of the SIR model’s Python code (SMPC 2020) based on the dataset of Algeria was done on the Google Colab platform, and the online date calculator was used to estimate the futuristic dates (DC 2020).

Determination of the recovery rate (γ)

If we suppose that I is a constant (I(t) = I0), we obtain the equation:

$$\frac{{{\text{d}}R}}{{{\text{d}}t}} = \gamma I_{0}$$
(9)

By integrating Eq. (9), then it is found the following Equation:

$$R\left( t \right) = \gamma tI_{0}$$
(10)

If time of recoveryt = T days, and \(R\left(T\right)={I}_{0}\), or \(\gamma T=1\).

Therefore, we obtain the next Equation:

$$\gamma \approx \frac{1}{T}$$
(11)

when time change dt = α, we obtain the following Equation from Eq. (3):

$$\frac{{R\left( {t + \alpha } \right)}}{\alpha } = \gamma I$$
(12)

or

$$\gamma \approx \frac{{R \left( {t + 1} \right) - R\left( t \right)}}{I\left( t \right)}$$
(13)

Estimation of Imax(percentage infected persons at the peak of the epidemic) and Sinf (percentage of susceptible people remaining after the end of the epidemic):

when we divide Eq. (2) by Eq. (1), we get the equation that follows:

$${\raise0.7ex\hbox{${\left( {\frac{{{\text{d}}I}}{{{\text{d}}t}}} \right)}$} \!\mathord{\left/ {\vphantom {{\left( {\frac{{{\text{d}}I}}{{{\text{d}}t}}} \right)} {\left( {\frac{{{\text{d}}S}}{{{\text{d}}t}}} \right)}}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{${\left( {\frac{{{\text{d}}S}}{{{\text{d}}t}}} \right)}$}} = - 1\frac{\gamma }{\beta } N \frac{1}{S}$$
(14)

or

$$\frac{{{\text{d}}I}}{{{\text{d}}S}} = - 1 + \frac{\gamma }{\beta } N\frac{1}{S}$$
(15)

If we integrate both sides, we obtain the Equation:

$$I = - S + \frac{\gamma }{\beta } N\log S + C$$
(16)

C is a constant.

Knowing that at the beginning of the infection, S≈N and the number of infected (I) is extremely low, this mean that At: \(t=0, I\sim 0\) and \(S\sim N\)

Therefore, if we substitute these values in Eq. (16), we get:

$$0 = - N + \frac{\gamma }{\beta }N\;In\;N + C$$
(17)

or

$$C = {\text{N}}\left( {1 + \frac{{\upgamma }}{{\upbeta }}{\text{InN}}} \right)$$
(18)

From Eq. (18), we take the value of C and we change it into Eq. (16). We will obtain Eq. (19):

$$I = N - S + \frac{\gamma }{\beta } N {\text{In}}\frac{S}{N}$$
(19)

Equation (19) is valid for all times. In general, the number of infected persons (I) increases exponentially from the beginning of the infection to reach a peak, before gradually shrinking to zero. We will estimate the percentage infected persons at the peak of COVID-19 (Imax) and remained susceptible people after the end of the epidemic.

Here, we use Eq. (19) to resolve Eq. (2) related to infection rate.

For simplification, we suppose \(\mathrm{S}=\mathrm{Ns},\mathrm{ I}=\mathrm{Ni},\) \(\mathrm{S}=\mathrm{Ns},\mathrm{I}=\mathrm{Ni},\) and \(\mathrm{R}=Nr\) \(\mathrm{S}=\mathrm{Ns}\) \(\mathrm{S}=\mathrm{Ns}\).

where s is the fraction of total susceptible population.

i is the fraction of total infected population.

r is the fraction of total recovered/removed population.

Therefore, following equations are obtained:

$$\frac{{{\text{d}}t}}{{{\text{d}}t}} = \beta is - \gamma i = i\left( {\beta s - \gamma } \right)$$
(20)

and

$$i = 1 - s + \frac{\gamma }{\beta }{\text{In}} S$$
(21)

At the peak of the infection when d \(\mathrm{di}/\mathrm{dt}=0\) \(\mathrm{di}/\mathrm{dt}=0\) \(\mathrm{di}/\mathrm{dt}=0\) \(\mathrm{di}/\mathrm{dt}=0\), s is calculated by the equation that follows:

$$s = \frac{\gamma }{\beta }$$
(22)

if we substitute the value of s from Eq. (22) into Eq. (21), we get Eq. (23):

$$i_{{\max}} = 1 + \frac{\gamma }{\beta }\left( {{\text{In}}\frac{\gamma }{\beta } - 1} \right)$$
(23)

We, now, need to find out \({S}_{inf}=\underset{t\to \infty }{\mathrm{lim}}S(t)\) the percentage of remaining susceptible individuals after the end of the infection. It is noted that at the end of infection t tends to infinity (∞), and thus, i = 0. We can rewrite Eq. (21) as follows:

$$1 - s + \frac{\gamma }{\beta }{\text{In}} s = 0$$
(24)

We can solve Eq. (24) numerically to obtain the value of s.

Determination of the basic reproduction number (R0)

R0 or the basic reproduction number is defined as the average number of people infected by a single individual. Mathematically, it represents the ratio of transmission and recovery rates.

$$R_{0} = \frac{\beta }{\gamma }$$
(25)

To calculate R0, we used the following Equations:

$$i_{{\max}} = 1 - \frac{1}{{R_{0} }}\left( {1 - {\text{In}}\frac{1}{{R_{0} }}} \right)$$
(26)

and

$$R_{0} = 1 - \frac{{{\text{In}}S_{{\inf}} }}{{S_{{\inf}} - 1}}$$
(27)

Results

In this work, we applied a SIR model to estimate different parameters related to COVID-19 epidemic curve for Algeria.

The values of the rate of transmission (β) and the constant m obtained are shown in Fig. 4.

Fig. 4
figure 4

a Total confirmed cases, deaths and recovered cases in Algeria with “Gamma,” “Beta” and “R0” parameters of SIR model. b Plot for “m” parameter of SIR model in case of Algeria

By approximating directly from the dataset of Algeria, we obtained the value of gamma plotted in Fig. 5.

Fig. 5
figure 5

Plot for “γ” parameter of SIR model in case of Algeria

The value of Imax was estimated at 1.91%. Taking into account the total number of the Algerian population, the number of sick peoples will be 8,37,555 persons (1.91% of 43,851,044).

The values of Sinfobtained for Algeria are 35.23% (Fig. 6).

Fig. 6
figure 6

Estimated values of tmax, Imax and Sinf in Algeria

Thus, the number of susceptible people remaining after the end of the COVID-19 epidemic in Algeria will be 15,448,723 persons (35.23% of 43,851,044).

According to Fig. 6, we can also see that the peak of the epidemic tmax will be reached after approximately 561 days, i.e., 8th September 2021 (predicted date) from its emergence (February 25th, 2020) in Algeria.

At last, from the dataset of Algeria, we obtained the values of R0 = 1.23 which is shown in Fig. 4.

Discussion

We estimated the different SIR parameters for COVID-19 in Algeria using data reported by the Algeria health ministry and the Johns Hopkins Institute dataset. This model which allows a prediction of a long period showed that the peak of the epidemic will be reached on September 2021. The basic reproduction number is estimated at 1.23 which is lower than that reported in India using the same model by Bagal et al. 2020.

At last, these parameters are estimated using only RT-PCR confirmed cases knowing that the number of tests do not exceed 2500 tests / day and the hospitalized patients diagnosed with CT scan is higher which could be considered as a limitation of the epidemic predictions.

Conclusion

In this study, we estimated the parameters of SIR models for the COVID-19 epidemic in Algeria basing on the daily reported/confirmed cases by the Algerian Ministry of Heath from February 25th, 2020 to August 12th, 2020.The SIR model showed that the peak of the epidemic will be reached after an estimated period of 561 days from the beginning of the epidemic. The model estimated that 1.91% of the total population will be affected and the percentage of susceptible people remaining after the end of the infection will be about 35.23%

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

SIR model:

Suspected–infected–recovered model

COVID-19:

coronavirus disease 2019

References

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

M L concepted the work, collected data, wrote, read and approved the final manuscript.

DK B did the simulations, read and approved the final manuscript.

Corresponding author

Correspondence to Mohamed Lounis.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lounis, M., Bagal, D.K. Estimation of SIR model’s parameters of COVID-19 in Algeria. Bull Natl Res Cent 44, 180 (2020). https://doi.org/10.1186/s42269-020-00434-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42269-020-00434-5

Keywords