In this paper, we propose a susceptible-exposed-infection-asymptomatic-hospitalized-recovered (SEIAHR) model with parameters on retrospective social distancing and masking. We estimated the model parameters from information published on the World Health Organization (WHO) website. We found that the actual reproduction number Rt varies over the period from 03 March to 07 June 2020 and moreover, effective control over contacts and the frequency of population movement would reduce the evolution of the epidemic (control c ≥ 50%). And the contact check has an influence on the base reproduction number R 0.
COVID-19 is an infectious disease caused by the latest coronavirus (SARS-CoV -2) which outbroke on November 17, 2019 in Hubei province (in central China), more precisely in the city of Wuhan. On March 11, 2020, the COVID-19 epidemic is declared a pandemic by the World Health Organization (WHO), which calls for essential protective measures to prevent the saturation of intensive care services and also to strengthen preventive hygiene measures (elimination of physical contact, kisses and handshakes, etc.).
In the Central African Republic, surveillance for the COVID-19 epidemic began from March 3, 2020. This surveillance concerns travelers coming from countries (France, Italy and China) affected by this epidemic. After the confirmation of the first case imported from Italy on March 14, 2020, and some sporadic cases, the government through the Ministry of Health decided that all people coming by the different routes of entry and exit of the country are systematically screened.
School closures followed, removal of gatherings of more than 15 people, and barrier measures applied. It should also be noted that the population was very resentful about the barrier measures due to the habit of over a decade of greeting with handshake, kisses and gatherings during major celebrations (marriage, mourning, baptism, etc.). The government has therefore taken great measures to educate the population about the COVID-19 epidemic.
Pioneers Kermack and McKendrick made mathematical models gain much more attention today. Usually, models included three conditions in the total population: S for susceptible, I for infected, and R for cured, for example, L. Allen in [
Some pathologies take into account exposure to the disease, which means that the population is in this case divided into four (4) compartments (S for susceptible, E for exposed, I for infectious and R for recovered), more details in [
Adam Kucharski in [
Abhijit Majumder1et al. [
Jonathan M. Readet al. [
On the other hand, certain diseases such as hepatitis B, hepatitis C, AIDS and today COVID-19, the exposed host E plays a very important role when dynamic behaviors are expected to be discussed.
Recent literature [
In the world, in order to understand the dynamics and this pandemic, several mathematical models in epidemiology have been developed, for more details on these different models in [
The spread and control of infectious diseases have been described by mathematical models, in order to study the dynamic properties of the disease.
We used an extension of the model (SEIR) by including asymptomatic infected and symptomatic infected compartments after being infected. The use of this extension comes from the reality of the dynamics of transmission of COVID-19 between 03 March 2020 and 07 June 2020 in the Central African Republic. We have found that when a person is exposed, they can become infected and transmit the disease whether asymptomatic or symptomatic. However, we do not take into account the delay in onset of symptoms since it varies from person to person.
In our work, we propose an SEIR model extended to nonlinear incidence with evaluation of social distancing on the dynamics of COVID-19 in order to estimate the effective number of Rt reproduction in the Central African Republic.
Compared to the epidemiological reality of the Central African Republic concerning this pandemic disease of COVID-19, our model is as follows:
{ ∂ S ( t ) ∂ t = Λ − β S ( t ) q E ( t ) + I ( t ) + A ( t ) N ( t ) ∂ E ( t ) ∂ t = β S ( t ) q E ( t ) + I ( t ) + A ( t ) N ( t ) + ε E ( t ) ∂ I ( t ) ∂ t = p ε E ( t ) − ( γ 1 + τ ) I ( t ) ∂ A ( t ) ∂ t = ( 1 − p ) ε E ( t ) − γ 2 A ( t ) ∂ H ( t ) ∂ t = p ε E ( t ) ∂ R ( t ) ∂ t = ( γ 1 + τ ) I ( t ) + γ 2 A ( t ) (1)
Biological interpretation of parameters:
· N ( t ) = S ( t ) + E ( t ) + I ( t ) + H ( t ) + R ( t ) the total sum of the population at time t;
· β Average transmission rate;
· q = ( 1 − c ) transmissibility reduction factor of the exposed class where c corresponds to the control intensity;
· Λ recruitment rate of the population;
· ε exposure rate;
· γ 1 infection rate of symptomatic;
· γ 2 infection rate of asymptomatic;
· τ COVID-19 death rate.
The sensitive population is increased by the net influx of recruitment Λ of individuals and is decreased by natural death (which we neglect over the study period and therefore is not part of our model).
The susceptible population also decreases following an infection, acquired by contact between a susceptible person and an infected person, which can be symptomatic, asymptomatic. The transmission coefficients for these classes of infected individuals are respectively β and q. In this paper, we use a basic transmission coefficient, β (which models both the infectivity of COVID-19 and contact rates), with modifying factors for asymptomatic individuals q.
These are individuals exposed to the COVID-19 virus, but there is no evidence that they are infected and can be infected asymptomatic or symptomatic.
These are individuals exposed to the virus and who subsequently became infected and confirmed in the laboratory. They are therefore either asymptomatic and symptomatic. The symptomatic population is generated after the development of clinical symptoms of COVID-19 by members of the exposed class.
0 < p ≤ 1 proportion of exposed individuals progresses to infectious class I ( t ) at a rate e and the rest ( 1 − p ) goes to the asymptomatic class.
These are laboratory confirmed and infected individuals, but have not developed clinical signs related to COVID-19. This is the proportion ( 1 − p ) of infected presentations that progressed to the asymptomatic class.
This is the proportion of exposed symptomatic infected patients requiring hospitalization on the notified doctor’s opinion.
It is the symptomatic and asymptomatic infected who are recovered.
This subsection is provided to prove the positivity and boundedness of solution of the system (1) with initial ( S ( 0 ) , E ( 0 ) , I ( 0 ) , A ( 0 ) , H ( 0 ) , R ( 0 ) ) T ∈ R 6 .
This research process for points of stability is much more detailed in [
For the above system (1), we find a region of attraction which is given by Lemma 1.
Lemma 1 [
Suppose Ω ⊂ R × C n is open, f i ∈ ( Ω , R ) , i = 1 , 2 , 3 , ⋯ , n .
If f i | x i ( t ) = 0 , X t ∈ C n ≥ 0 , X t = ( x 1 t , x 2 t , ⋯ , x n t ) T , i = 1 , 2 , 3 , ⋯ , n , then C n { ϕ = ( ϕ 1 , ⋯ , ϕ n ) : ϕ ∈ C ( [ − τ ,0 ] , R + n ) } is the invariant domain of the following equations
d x i ( t ) d t = f i ( t , X t ) , t ≥ σ , i = 1 , 2 , ⋯ , n ; (2)
where R + n = { ( x 1 , ⋯ , x n ) : i = 1 , ⋯ , n } .
Proposition 1 The system (1) is invariant in R + 6 .
Proof. We have by the system (1):
d X d t = F ( X ( t ) ) , X ( 0 ) = X 0 (3)
F ( X ( t ) ) = ( F 1 ( X ) , F 2 ( X ) , ⋯ , F 6 ( X ) ) T
We note that
d S d t | S = 0 = Λ ≥ 0 , d E d t | E = 0 = β S I + A N ( t ) ≥ 0 , d I d t | I = 0 = ( 1 − p ) ε ≥ 0
d A d t | A = 0 = p ε E ≥ 0 , d H d t | H = 0 = ( 1 − p ) ε E ≥ 0 , d R d t | R = 0 = ( γ 1 + τ ) I ( t ) + γ 2 A ( t ) ≥ 0
Then it follows from the Lemma 1 that R + 6 is an invariant set.
Lemma 2 The system (1) is bounded in the region
Ω = { ( S , E , I , A , H , R ) ∈ R + 6 | S + E + I + A + H + R ≤ Λ } (4)
The proof of Lemma 2 is obvious since we observed from the system that
d N d t ≤ Λ (5)
which ⇒ l i m t → + ∞ sup N ( t ) ≤ Λ
The DFE of our model (1) is given by:
Ω 0 = ( S 0 , E 0 , I 0 , A 0 , H 0 , R 0 ) = ( Λ , 0 , 0 , 0 , 0 , 0 )
The local stability of Ω 0 can be established on the system (1) by using the next generation operator method.
The basic reproduction number R0, the expected number of secondary infections produced by an index case in a completely susceptible population, is a key factor used in estimating the transmissibility of infectious disease.
For our model (1), the basic reproduction number is:
R 0 = β ( q ε + p γ 1 + 1 − p γ 2 ) (6)
For more details on the calculation method of R0, [
The threshold quantity R0 is the basic reproduction number of the disease. This represents the average number of secondary cases generated by an infected person in a fully susceptible population. The epidemiological significance is that when R0 is less than unity, a low number of infected individuals in the population will not cause major epidemics and so much disease to go away over time.
It should also be noted that R0 is by no means a characteristic of a given disease, but the one of a given disease in a given host population. As a result, the same disease could have two different R0 in two different host populations and two different diseases could also have two different R0 in the same host population.
The effective reproduction number Rt is defined as the expected number of secondary infections per infectious person in a population composed of both susceptible and non-susceptible hosts [
The formula for calculating Rt is as follows:
R t = R 0 × s ^ ( t ) (7)
where s ^ ( t ) is the fraction of the host population sensitive to time t.
If:
· R t > 1 : the number of new cases will be increased;
· R t = 1 : the disease will become endemic;
· R t < 1 : there will be a decrease in the number of new cases.
Lemma 3 The DFE, Ω 0 of the model (1) is locally-asymptotically stable if Ω whenever R t > 1 , and unstable if R t < 1 .
The data were observed over the period from March 03, 2020 to June 07, 2020, 14,259 samples were taken and tested with real-time PCR in the laboratory.
We have:
· 1570 positives, therefore a national prevalence over the period of 11%;
· 5 epidemic peaks observed (31 to 24 April, 49 to 09 May, 80 to 16 May, 119 to 30 May, 163 to 06 June);
· The dynamics of the epidemic are nolinear.
The first case of COVID-19 confirmed on March 14, 2020 in the country was a case imported from Italy. As a result, other cases started to be reported. The Ministry of Public Health has therefore decided to systematically test all people arriving by land in Central African territory since the airport was already closed during the period.
It is therefore from these data that we estimated certain biological parameters in our mathematical model which is a non-linear model of the dynamics and circulation control of COVID-19 in the Central African Republic.
The general population is not quarantined during this period compared to other countries given the country’s very weak economic context. Only barrier measures and social distancing have been adopted by the Ministry of Public Health and the population has been widely made aware of these various measures.
The data that we used to build our mathematical model are the data for the period from March 03, 2020 to June 07, 2020 retrieved from the World Health Organization website (http://www.who.int/fr/emergencies/diseases/novel-coronavirus-2019).
In
In order to fit the model to time series data from confirmed COVID-19 cases, it is necessary to define certain parameter values. In our model, birth and natural death rates were assumed to ignore over our study period. It is also assumed that the case distribution profile in the different provinces is similar across the country.
We can clearly see in
| Data | Values |
|---|---|
| Total number samples | 14,259 |
| Number infected | 1570 |
| Number revered | 29 |
| Number of death | 5 |
We used the maximum likelihood method to estimate the parameters of our model [
The force of infection β = − q ∗ log ( 1 − p ) with p the probability of infection during infectious contact (intrinsic property of the disease).
And q = 1 − c where c control of contact between individuals. We therefore evaluate the different values of R0 for different values of c (weak, medium and strong).
We used R software with R0 package to estimate reproduction numbers for epidemic outbreaks.
Assuming that there is no difference between cure rates in symptomatic and asymptomatic we have:
If we assume that there is a difference between the cure rates in symptomatic and asymptomatic patients, we have:
| Parameters | Biological Meaning | Estimated Values (95% CI) |
|---|---|---|
| Λ | Recruitment rate of population | 0.87 (0.86 - 0.89) |
| τ | COVID-19 death rate | 0.0032 (0.0004 - 0.006 ) |
| γ 1 | symptomatic infection rate | 0.17 (0.15 - 0.19 ) |
| γ 2 | asymptomatic infection rate | 0.83 (0.81 - 0.85 ) |
| c: control | q = 1 − c | β | R0 |
|---|---|---|---|
| 0.25 | 0.75 | 0.22 | 1.45 |
| 0.5 | 0.5 | 0.15 | 0.89 |
| 0.85 | 0.15 | 0.04 | 0.23 |
| c: control | q = 1 − c | β | R0 |
|---|---|---|---|
| 0.25 | 0.75 | 0.22 | 0.72 |
| 0.5 | 0.5 | 0.15 | 0.4 |
| 0.85 | 0.15 | 0.04 | 0.09 |
COVID-19 is a pandemic that strikes all countries in the world. Each country controls this disease according to the means it has but also internationally recognized barrier gestures. During the study period, the Central African Government had also taken decisions going in the direction of social distancing and masking. We must recognize that if we really have to evaluate these different measures in the population, especially in Bangui, we cannot even reach 25%.
This means that very few of the people of Bangui respect these social distancing and masking. This leads, according to
On the other hand, if the respect of social distancing and masking by the population reaches at least 50%, the value of β = 0.15 and R 0 = 0.85 < 1 , the disease will die out over time. Indeed, our study goes well with the reality over the study period. At the start of the pandemic in the country, the social distancing and masking decided by the Government were not respected at all. It took a lot of awareness, and the squaring of major roads by the police by controlling buses and motorcycle taxis made it possible to control the pandemic well despite the limited resources of the Government.
In
It must also be recognized that the low density unevenly spread and the low frequency of displacement of the Central African population played in its favor. Since over this period, there is a low frequency of displacement, unlike in countries with high density and frequency of displacement of the population.
In this paper, our objective is to estimate the number of reproduction in real time Rt in retrospective and to study the measures of social distancing and masking controls. We have found that indeed by contact control through social distancing and masking remains the best possibility to control the COVID-19 pandemic. In addition, the high population density and population dynamics are an important factor in the spread of the disease.
We have shown that a strong contact control and the frequency of displacement of the population would reduce the evolution of the epidemic.
Therefore, for countries with low income and low density and frequency of population displacement, social distancing and masking would be the best possibility to control COVID-19.
Sincere thanks to the AM members for their professional performance, and a special thank you to the editorial members for a rare high quality attitude.
The authors declare that they have no conflict of interest.
Yambiyo, B.M., Manirakiza, A. and N’Guérékata, G.M. (2021) Estimation of Reproduction Numbers for the Period of the COVID-19 Epidemic Outbreak in the Central African Republic. Applied Mathematics, 12, 477-488. https://doi.org/10.4236/am.2021.126033