# Simulation of the Final Size of the Evolution Curve of Coronavirus Epidemic in Morocco using the SIR Model

Authors: Ousama Ifguis; Mohamed El Ghozlani; Fouzia Ammou; Abdelaziz Moutcine; Zeroual Abdellah
Journal: Journal of Environmental and Public Health
 (2020)

---

## Abstract

Since the epidemic of COVID-19 was declared in Wuhan, Hubei Province of China, and other parts of the world, several studies have been carried out over several regions to observe the development of the epidemic, to predict its duration, and to estimate its final size, using complex models such as the SEIR model or the simpler ones such as the SIR model. These studies showed that the SIR model is much more efficient than the SEIR model; therefore, we are applying this model in the Kingdom of Morocco since the appearance of the first case on 2 March 2020, with the objective of predicting the final size of the epidemic.


---

## Body

## 1. Introduction

During Christmas 2019, COVID-19 caused an epidemic in the city of Wuhan, Hubei Province of China [1]. It spread to the other parts of China and subsequently to many other countries around the world. Morocco is one of the countries affected by COVID-19. On 2 March 2020, the country identified the first case [2, 3], and as of 30 March 2020, confirmed cases reached up to 516 with a death toll of 27 [4]. Also, as the number of infected cases is increasing, it is necessary for modellers to estimate the severity of the epidemic in terms of the total number of people infected, the total number of confirmed cases, the total number of deaths, and basic reproduction and to predict the duration of the epidemic, the arrival of its peak, and its final size. This information can help public health agencies make informed decisions.In this work, we used the SIR model to predict the development of the epidemic in the Kingdom of Morocco from the identification of the first case on 2 March 2020 in the city of Casablanca [2], given the reliability of the data and the definition of confirmed cases during this period and the simplicity of our forecasts and analyses. We were able to determine the detailed results of the SIR model calibration and the predictions of our model, including the distribution of the peak period, the prediction interval of future confirmed cases, and the total number of infected persons.


## 2. The SIR Model

In the SIR model, compartmentS refers to the sensitive population in Morocco, I refers to the infectious population, and R refers to confirmed cases [5]. The latency of COVID-19 infection is biologically realistic due to an incubation period of up to 14 days; newly infected persons may not be contagious during this time period as the virus will be incubating in the organism.Here, we highlight the difference between the latency period, the period between when an individual is infected and when he or she is infectious, and the incubation period, the period between when an individual is infected and when clinical symptoms, including fever and cough symptomatic of COVID-19, appear.The transfer diagram for the model is shown in Figure1. The biological significance of all model parameters is given in Table 1. A key assumption in both models is that deaths occurring in compartments S, E, and I are negligible during the model prediction period. The differential equation system for the SIR model is given as follows [5]:(1)dSdt=−βSI,(2)dIdt=βSI−γI,(3)dIdt=γI.Figure 1
Transfer diagrams for the SIR model for COVID-19 in Morocco.Table 1
Daily increased number of COVID-19 cases in Morocco [2].
DatePositive casesDeath casesReestablished casesNew cases2 March 202010013 March 202010004 March 202010005 March 202020016 March 202020007 March 202020008 March 202020009 March 2020200010 March 2020310111 March 2020610312 March 2020611013 March 2020811214 March 202018111015 March 202028111016 March 20203711917 March 20204421718 March 202054211019 March 20206322920 March 202086322321 March 202096331022 March 2020115431923 March 2020143452824 March 2020170562725 March 2020225675526 March 20202751085027 March 202034523117028 March 202040225125729 March 202047926137730 March 202055633157731 March 2020617362461We can estimate the nature of the disease in terms of the power of infection:(4)R0=βγ.It is called the basic reproduction number.R0 is the average number of people infected from another person. If it is high, the probability of the pandemic is also higher.


## 3. Simulation with COVID-19 Data in Morocco

For this study, we used the data accumulated since 2 March 2020 concerning Morocco [2]. Four parameters to be estimated in the SIR model from the data are as follows: basic reproduction number (R0), contact rate (β), removal rate (γ), and final number of cases (I).We do not consider the effect of the natural death or birth rate so that the total population remains constant.(i)
N=constant(ii)
N0 = total population in Morocco = 40,000,000(iii)
I0 = 1(iv)
R0 = 0


## 4. Results

To solve the ordinary differential equation of the SIR model (equations (1)–(3)), we must first simulate the daily data of the SIR model, compare them with the real data, and then execute the optimization algorithm (reduce the difference between the real data and the corresponding simulated data), so that it searches for the values of α and β, which minimize the difference between the real data and the simulated ones.In the optimal model, we supposed that the predicted infected cases should definitely be close to the actual number of infected cases; we found that the critical number of susceptible cases was about 1,812 (Table2).Table 2
The estimation of parameters of the SIR model.
CountryInitial number of cases (I0)Contact rate (β)Removal rate (γ)Basic reproduction number (R0)Critical number of susceptible casesFinal number of casesFinal number of susceptible casesMorocco10.499 (1/day)0.249 (1/day)2.0031,8121,446366The transfer diagram for the SIR model used in the simulation of the size of COVID-19 infection in the Moroccan state is provided in detail in Figure1.We note from Figure2 that the number of confirmed cases is not very high; the highest daily confirmed number was 70 on 28 March 2020. This can be explained by the strategy of Morocco, which declared a state of health emergency and confinement on 20 March 2020, Friday, at 6 p.m., to limit the displacement of the population as much as possible, the only way to keep the coronavirus under control.Figure 2
(a) Confirmed cases, (b) death cases, and (c) recovered cases from 2 March 2020 to 31 March 2020.
(a)(b)(c)Figure3, based on optimal SIR models, shows that the start of acceleration of the epidemic is around 21 March 2020, the regular growth will begin on 8 April 2020, and the end of the epidemic in Morocco would be around 26 April 2020, with a total of 1,446 infected cases and 366 final number of susceptible cases (Table 2).Figure 3
Presentation of the evolution of total infected cases (R0 = 2.003, β = 0.499, γ = 0.249, N = 1,812, Cend = 1,446, Send = 366, RMSE = 7).The peak of the reported cases will be around 28 March 2020 with 80 confirmed cases (Figure4).Figure 4
Daily estimated number of reported cases.The optimal estimates ofR0 ranges from 1.22 to 2.55 (Table 3); our study confirmed that the transmissibility of COVID-19 is R0 = 2.003 (Table 3), which was comparable to that of severe acute respiratory syndrome coronavirus (SARSCoV) (ranging from 2.2 to 3.7) [6] and much higher than that of Middle East respiratory syndrome coronavirus (MERS-CoV) (ranging from 0.47 to 0.91) [7].Table 3
Basic reproduction numberR0 (t) of the time-dependent SIR model of COVID-19 in Morocco based on the given data from 21 March 2020 to 31 March 2020.
DaysThe reproduction numberR021 March 20201,75122 March 202014923 March 20201,18224 March 20201,22925 March 20201,33926 March 20201,58627 March 20201,33728 March 20202,55929 March 20201,372


## 5. Conclusion

Our simulation study on the optimization of the final size of COVID-19 epidemic evolution in the Kingdom of Morocco, with the SIR model, has allowed us to accurately predict the peak of the infected and death cases (Table2), although the number of people tested is very low, about 3,079, until 31 March 2020.The Moroccan government should probably increase the number of cases tested on a daily basis in order to accurately identify the true size of the pandemic in Morocco.



---

*Source: 9769267-2020-06-02.xml*