COVID-19 outbreak in Algeria: A mathematical Model to predict cumulative cases
==============================================================================

* Mohamed Hamidouche, Sr.

## Abstract

**Introduction** Since December 29, 2019 a pandemic of new novel coronavirus-infected pneumonia named COVID-19 has started from Wuhan, China, has led to 254 996 confirmed cases until midday March 20, 2020. Sporadic cases have been imported worldwide, in Algeria, the first case reported on February 25, 2020 was imported from Italy, and then the epidemic has spread to other parts of the country very quickly with 139 confirmed cases until March 21, 2020.

**Methods** It is crucial to estimate the cases number growth in the early stages of the outbreak, to this end, we have implemented the Alg-COVID-19 Model which allows to predict the incidence and the reproduction number R0 in the coming months in order to help decision makers.

The Alg-COVIS-19 Model initial equation 1, estimates the cumulative cases at t prediction time using two parameters: the reproduction number R0 and the serial interval SI.

**Results** We found R0=2.55 based on actual incidence at the first 25 days, using the serial interval SI= 4,4 and the prediction time t=26. The herd immunity HI estimated is HI=61%. Also, The Covid-19 incidence predicted with the Alg-COVID-19 Model fits closely the actual incidence during the first 26 days of the epidemic in Algeria Fig. 1.A. which allows us to use it.

According to Alg-COVID-19 Model, the number of cases will exceed 5000 on the 42th day (April 7th) and it will double to 10000 on 46th day of the epidemic (April 11th), thus, exponential phase will begin (Table 1; Fig.1.B) and increases continuously until reaching à herd immunity of 61% unless serious preventive measures are considered.

**Discussion** This model is valid only when the majority of the population is vulnerable to COVID-19 infection, however, it can be updated to fit the new parameters values.

Keywords
*   coronavirus
*   COVID-19
*   A mathematical model
*   Algeria outbreak

## Introduction

On December 29, 2019, Wuhan, the capital city of Hubei Province in Central China, has reported four cases of pneumonia with unknown etiology (unknown cause), the next day, the WHO China Country Office was informed (WHO.a 2020) about this pneumonia cases that were found to have a link with Huanan seafood and animal market in Wuhan, the Centers for Disease Control and Prevention (CDC) and Chinese health authorities determined and announced later that a novel coronavirus denoted as Wuhan (CoV) had caused the pneumonia outbreak (CDC 2020). Since then, the outbreak has rapidly spread over a short span of time and has received considerable global attention.

On January 7, 2020 the etiological agent of the outbreak was identified as a novel coronavirus (2019-nCoV) and its gene sequence was quickly submitted (GenBank 2019), the coronavirus was renamed COVID-19 by WHO on February 12, 2020. It has since been identified as a zoonotic coronavirus, similar to SARS and MERS coronaviruses. (Ying, et al. 2020)

On January 30, WHO announced the listing of this novel coronavirus-infected pneumonia (NCP) as a “public health emergency of international concern”, A total of 254 996 confirmed cases of infection with COVID-19, including 10444 deaths have been reported worldwide as on midday March 20th, 2020 (WHO.a 2020).

Sporadic cases have been imported to Europe, Africa and North and South America via returning travellers from China. In Algeria, the first case of COVID-19 was reported on February 25, 2020, when an Italian national tested positive in Ouargla region in the south of the country, a few days later, on March 1, 2020, two cases were reported in Blida region in the North of Algeria, following their contacts with two Algerian nationals who came from France for holidays, they were detected positive after their return to France, since then, a COVID-19 outbreak has started in this region (Blida) that form a cluster of more than 5,4 million inhabitants with the surrounding cities (Algiers, Boumerdes, Tipaza) (Algerian Ministry of Health 2020, ONSA 2020), now, the epidemic is spreading to other parts of the country, until March 22, 2020, the Algerian authorities have declared 200 confirmed cases with a fatality rate of 8,5% (17 deaths) (Algerian Ministry of Health 2020).

Meanwhile, there is considerable uncertainty as to the extent of the epidemic and its parameters, the COVOD-19 reproduction number (R0) has been estimated in various studies, they found, 2,35 (95% CI 1.15–4.77) (Kucharski, et al. 2020), two studies used stochastic methods to estimate R0 have reported a range of 2.2 to 2.68 with an average of 2.44. (Joseph, et al. 2020, Riou et Althaus 2020), Six studies that used mathematical methods to estimate R0 produced a range from 1.5 to 6.49 with an average of 4.2. (Shen, et al. 2020, Read, et al. 2020, Chong, et al. 2020, Natsuko, et al. 2020, Tang, et al. 2020) The three studies using statistical methods such as exponential growth estimated an R0 ranging from 2.2 to 3.58, with an average of 2.67. (Liu, et al. 2020, Zhao, et al. 2020). Also, the COVID-19 fatality rate vary by region from 0.39% in Norway to 8,3% in Italy, in China the fatality rate is 4%, (Wilson, et al. 2020, Johns Hopkins 2020), however, the highest mortality rate remain in Algeria at 8,5% until March 22, 2020 (Algerian Ministry of Health 2020)

In the early stages of a new infectious disease outbreak, as COVID-19 in Algeria, it is crucial to estimate the transmission dynamics and inform predictions about potential future epidemic growth (Viboud, et al. 2018), it can provide insights into the epidemiological situation which help decision makers to adapt the health system capacities, thus, a prediction Model can help to do that and identify whether outbreak control measures are having a measurable effect or not (Funk, et al. 2017, Riley, et al. 2003) and guide the design of alternative interventions (Kucharski, et al. 2015), in addition, a prediction Model can be updated to help estimate risk to other countries (Cooper, et al. 2006).

To this end, we have implemented the Alg-COVID-19 Model which allows to predict the cumulative cases in the coming weeks, and to calculate the actual basic reproduction number (R0), consequently, this Model can show hospitals what to expect in terms of Covid-19 patients, the percentage who need to be in an intensive care unit (ICU) or on a ventilator and the future number of deaths based on a given data.

## Methods

### Epidemiological data

We retrieved information on cases number with confirmed COVID-19 infection based on official reports from governmental institutes in Algeria (Algerian Ministry of Health 2020).

### The mathematical Model

In order to predict the cumulative cases of COVID-19 Algerian epidemic in the coming weeks we used the mathematical model (Alg-COVID-19) defined as: ![Formula][1]</img> 

**I(t):** The incidence (cases number) at t time

**R0:** The reproduction number

**SI:** The serial interval

**t:** The prediction time

### Estimation of Basic Reproduction Number (R0) and the herd immunity

The basic reproduction number is defined as the (average) number of new infections generated by one infected individual during the entire infectious period in a fully susceptible population. It can be also understood as the average number of infections caused by a typical individual during the early stage of an outbreak when nearly all individuals in the population are susceptible to infection. The basic reproduction number reflects the ability of an infection spreading under no control, it has three components (Roy, et al. 1992): ![Formula][2]</img> 

CB: effective contact rate (C: contact frequency, B: contact efficacy)

N: susceptible population size

D: infectious phase duration

The approach implemented to estimate the basic reproduction numbers (R0) in this model is to calculate the average R0 using the actual chronological cumulative cases during the first 26 days of the COVID-19 Algeria epidemic between February 25 and March 22, 2020, so that, we used the equation (3) derived from the equation (1). The R0 will be used also to calculate the herd immunity (HI) needed to stop the epidemic spontaneously based on the equation 4. ![Formula][3]</img>  ![Formula][4]</img> 

### Estimation of the serial interval (SI)

The serial interval (SI) is the time between symptom onset of a primary and secondary case, a previous studies reported that (SI = 4.4±3.17) days (Chong, et al. 2020), so we used this value in our model to predict the incidence.

### Statistical analysis and software

This study was conducted using Excel 2013 and STATA/IC 15 software

## Results

### The initial parameters used in Alg-COVID-19 Model

Without considering the prevention measures and other factors, this paper focused on the predicted COVID-19 incidence at any time of the epidemic using Alg-COVID-19 Model (equation 1), based on three parameter combinations that create plausible epidemic curves, namely, the average basic reproduction number calculated with the actual cases number at the first 26 days R0=2.55 (95% CI 2.15–2.93) (equation 3), the serial interval (SI=4.4) and the prediction time (t) we are looking for. The herd immunity (HI) estimated using the (equation 4) is HI=61%.

### Model fitting

The Covid-19 incidence (cases number) predicted with the Alg-COVID-19 Model fits closely the actual incidence during the first 26 days of the epidemic in Algeria (Fig. 1.A), this allows us to use this Model to predict the Covid-19 incidence in the next months of the epidemic (Fig. 2.B).

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/03/27/2020.03.20.20039891/T1)

Table 1. 
Alg-COVID-19 Model results for next weeks

![Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/27/2020.03.20.20039891/F1.medium.gif)

[Fig. 1.](http://medrxiv.org/content/early/2020/03/27/2020.03.20.20039891/F1)

Fig. 1. 
A: Alg-COVID-19 Model match to actual data (Cases number) of the first 26 days of epidemic. B: Alg-COVID-19 Model before and before (R0=2.55) and after mitigation (R0<2.55).

### Alg-COVIS-19 Model results

The estimations presented by this Model are based on the first 26 days data and cover the next months of COVID-19 epidemic in Algeria, according to that, the number of cases will exceed 1000 case on the 35th day of the epidemic (March 31, 2020), 5000 on the 42th day (April 7th) and it will double to 10000 on 46th day of the epidemic (April 11th), thus, exponential phase will increases continuously until reaching à herd immunity of 61% or a serious preventive measures are considered (Table 1; Fig.1.B).

## Discussion

Regarding the used parameters (R0 and SI) by Alg-CIVID-19 Model for the predictions, the estimated R0=2.55 (95% CI 2.15–2.93) until March 19, 2020, is very close to R0 that was found on the Chinese epidemic data R0=2,35 (95% CI 1·15–4·77) until February 4th, 2020, (Kucharski, et al. 2020), also, it is in the range of basic reproduction numbers of other studies mentioned below (Shen, et al. 2020, Read, et al. 2020, Chong, et al. 2020, Natsuko, et al. 2020, Tang, et al. 2020), on the other hand, the parameter that we cannot update and can affect the predictions is the serial interval (SI), we found two studies each one has different value that range from 4.4±3,17 (Chong, et al. 2020) to 7.5 days (95% CI, 5.3 to 19) (Li, et al. 2020), we used (SI=4.4) because of its standard deviation is not very large.

This model is valid only when the majority of the population is vulnerable to COVID-19 infection, since it does not take into account the herd immunity acquired ongoing epidemic. In case of R0 decreases over the next weeks of the epidemic, Alg-COVID-19 Model can be revalidated by recalculating a new effective reproductive number R1 (equation 3). However, this will happen only if the herd immunity rises spontaneously during the epidemic or the state act on one of the components of R0 by preventive measures (equation 2).

In the case of the Algerian COVID-19 epidemic, this Model can be valid for the first few months if the herd immunity remain low and susceptible population size is not changed, thus, it is applicable whenever a new regional outbreak begins with a naive population. On the other hand, if the state implement a good preventive policy according to the WHO guidelines (WHO.b 2020) and set up measures to reduce the effective contact rate, this will obviously reduce R0, so that, (R1) the effective reproduction number with control and preventive measures can be calculated to estimate the implemented preventive actions effectiveness, and Alg-COVID-19 Model can be updated with R1 values adapted to the new situation and its forecasts would correspond as closely as possible to the real epidemic growth (Fig.1.B). Also, it is obvious that the measures must be strong enough to reduce the basic reproduction number below 1 (R0 <1) to control the epidemic, otherwise, the epidemic will spread over time, which can also be a good strategy for maintaining the healthcare demand at a manageable pace.

According to Alg-COVID-19 Model results on Algeria data, the exponential phase 3 will probably announced after the 35th day of the epidemic (1000 case) (March 31, 2020), consequently, the response to the epidemic will change from trying to curb the spread of the COVID-19 on the territory to a strategy that lead to mitigate the effects of the epidemic wave, so that, after R0 and SI update (Fig. 1.B), the Model can be used to predict the hospitalization and ICU or on a under ventilator patients numbers, and also, predict the number of deaths that remain the highest in the wold (Algeria: 8.5% until March 22nd) since the start of the pandemic, without any concrete explanation.

## Data Availability

All data refferd to in the manuscript is available

## Conflict of Interest Statement

The authors have no conflicts of interest to declare

*   Received March 20, 2020.
*   Revision received March 23, 2020.
*   Accepted March 27, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  Algerian Ministry of Health. 2020. [http://covid19.sante.gov.dz/carte/](http://covid19.sante.gov.dz/carte/).
    
    

2.  CDC. 2020. [https://www.cdc.gov/coronavirus/novel-coronavirus-2019.html](https://www.cdc.gov/coronavirus/novel-coronavirus-2019.html).
    
    

3.  Chong, You, et al. “Estimation of the Time-Varying Reproduction Number of COVID-19 Outbreak in China.” medRxiv, 2020.
    
    

4.  Cooper, BS.,  RJ. Pitman,  WJ. Edmunds, and  NJ. Gay. “Delaying the international spread of pandemic influenza.” PLoS Med, 2006: 3: e212.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.0030212&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16640458&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F27%2F2020.03.20.20039891.atom) 

5.  Funk, S.,  I. Ciglenecki,  A. Tiffany, and et al. “The impact of control strategies and behavioural changes on the elimination of Ebola from Lofa County, Liberia.” Philos Trans R Soc Lond B Biol Sci, 2017: 372(1721).
    
    

6.  GenBank. 2019. [http://virological.org/t/novel-2019-coronavirus-genome/319](http://virological.org/t/novel-2019-coronavirus-genome/319).
    
    

7.  Johns Hopkins. 20 03 2020. [https://coronavirus.jhu.edu/](https://coronavirus.jhu.edu/).
    
    

8.  Joseph, T. Wu,  Leung Kathy,  M. Leung Gabriel, and et al. “Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study.” THE LANCET, 2020: Vol 395.
    
    

9.  Kucharski, Adam J., et al. “Early dynamics of transmission and control of COVID-19: a mathematical modelling study.” Lancet Infect Dis Published Online March 11, 2020, 2020: (20)30144–4.
    
    

10. Kucharski, AJ.,  A. Camacho,  F. Checchi, and et al. “Evaluation of the benefits and risks of introducing Ebola community care centers, Sierra Leone‥” Emerg Infect Dis., 2015: 21: 393–99.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2103.141892&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25694150&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F27%2F2020.03.20.20039891.atom) 

11. Li, Q., et al. “Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus–Infected Pneumonia.” NEJM, 2020.
    
    

12. Liu, T., et al. “Transmission dynamics of 2019 novel coronavirus (2019-nCoV).” bioRxiv, 2020.
    
    

13. Natsuko, Imai, et al. Report 3: Transmissibility of 2019-nCoV. COVID-19 Response Team, Imperial College London, 2020.
    
    

14. ONSA. Office national des statistiques d’Algérie 2020. [http://www.ons.dz/](http://www.ons.dz/).
    
    

15. Read, J.M.,  JRE. Bridgen,  DAT. Cummings,  A. Ho, and  CP. Jewell. “Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions.” medRxiv, 2020.
    
    

16. Riley, S.,  C. Fraser,  CA. Donnelly, and et al. “Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions.” Science, 2003: 300: 1961–66.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMDAvNTYyNy8xOTYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMjcvMjAyMC4wMy4yMC4yMDAzOTg5MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

17. Riou, J., and  C. L. Althaus. “Pattern of early human-to-human transmission of Wuhan 2019-nCoV.” bioRxiv, 2020.
    
    

18. Roy, M., Anderson, May Robert M., and et al. Infectious Diseases of Humans: Dynamics and Control. oxford science publications, 1992.
    
    

19. Shen, M.,  Z. Peng,  Y. Xiao, and  L. Zhang. “Modelling the epidemic trend of the 2019 novel coronavirus outbreak in China.” bioRxiv, 2020.
    
    

20. Tang, Biao, et al. “Estimation of the Transmission Risk of 2019-nCov and Its Implication for Public Health Interventions.” 2020.
    
    

21. Viboud, C.,  K. Sun,  R. Gaffey, and et al. “The RAPIDD Ebola forecasting challenge: synthesis and lessons learnt.” Epidemics, 2018: 22: 13–21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.EPIDEM.2017.08.002.&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28958414&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F27%2F2020.03.20.20039891.atom) 

22. WHO.a. 2020. [https://www.who.int/csr/don/05-january-2020-pneumonia-of-unkown-cause-china/en/](https://www.who.int/csr/don/05-january-2020-pneumonia-of-unkown-cause-china/en/).
    
    

23. WHO.b. 2020. [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance).
    
    

24. Wilson, N.,  A. Kvalsvig,  L. Telfar, and  MG. Baker. “Case-Fatality Risk Estimates for COVID-19 Calculated by Using a Lag Time for Fatality.” Emerg Infect Dis., 2020: 10 Volume 26, Number 6.
    
    

25. Ying, Liu,  A. Gayle Albert,  Wilder-Smith Annelies, and Rocklöv Joacim. “The reproductive number of COVID-19 is higher compared to.” Journal of Travel Medicine, 2020.
    
    

26. Zhao, S., et al. “Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A datadriven analysis in the early phase of the outbreak.” bioRxiv, 2020.

 [1]: /embed/graphic-1.gif
 [2]: /embed/graphic-2.gif
 [3]: /embed/graphic-3.gif
 [4]: /embed/graphic-4.gif