Estimation of Reproduction Numbers for the Period of the COVID-19 Epidemic Outbreak in the Central African Republic

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 in-formation published on the World Health Organization (WHO) website. We found that the actual reproduction number R t 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 50% c ≥ ). And the contact check has an influence on the base reproduction number R 0 .


Introduction
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 confir-mation 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 [1] to study the positivity of solutions in the case of an SIS model. Likewise, Y. Cai et al. in [2] introduced the stochastic form of the SIS model incorporating media coverage to study the effects of environmental fluxes. Also, B. M.
Yambiyo et al. studied the SIR model using the fractional derivative of Caputo.
They studied the stability of solutions using the homotopy pertubation technique [3].
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 [4] [5] [6]. Adam Kucharski in [7], they developed the same type of model (SEIR). They took into account the delays in the onset and reporting of symptoms, thereby including compartments to reflect transitions and reporting states and disease states in Wuhan.
Abhijit Majumder1 et al. [8], used another type of epidemic model, notably the stochastic model, to determine the end point of disease extinction. They also adjusted the epidemic data for Italy and estimated the epidemic burden for that country.
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 [10] [11] has shown that a model of the Susceptible-Exposed-Infected-Recovered (SEIR) type is more suitable for studying the dynamics of COVID-19.
In the world, in order to understand the dynamics and this pandemic, several mathematical models in epidemiology have been developed, for more details on 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 R t reproduction in the Central African Republic.

Mathematical Model
Compared to the epidemiological reality of the Central African Republic concerning this pandemic disease of COVID-19, our model is as follows: 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.

Population of Susceptible S(t)
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.

Exposure Population E(t)
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.

Infected I(t)
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 1 p < ≤ proportion of exposed individuals progresses to infectious class ( ) I t at a rate e and the rest ( ) 1 p − goes to the asymptomatic class.

Asymptomatic A(t)
These are laboratory confirmed and infected individuals, but have not developed clinical signs related to COVID-19. This is the proportion ( ) presentations that progressed to the asymptomatic class.

Hospitalized Serious Symptomatic and Notified H(t)
This is the proportion of exposed symptomatic infected patients requiring hospitalization on the notified doctor's opinion.

Recovered R(t)
It is the symptomatic and asymptomatic infected who are recovered.

Positivity and Boundedness of the Solution for the Model(1)
This subsection is provided to prove the positivity and boundedness of solution This research process for points of stability is much more detailed in [16] [17].
For the above system (1), we find a region of attraction which is given by Lemma 1.
is the invariant domain of the following equations Proof. We have by the system (1): Then it follows from the Lemma 1 that 6 + R is an invariant set.
The proof of Lemma 2 is obvious since we observed from the system that

Local Stability of Disease-Free Equilibrium (DFE)
The DFE of our model (1) is given by: The local stability of 0 Ω can be established on the system (1) by using the next generation operator method.
The basic reproduction number R 0 , 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: For more details on the calculation method of R 0 , [19] [20] [21]. The threshold quantity R 0 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 R 0 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 R 0 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 R 0 in two different host populations and two different diseases could also have two different R 0 in the same host population.
The effective reproduction number R t is defined as the expected number of secondary infections per infectious person in a population composed of both susceptible and non-susceptible hosts [22].
The formula for calculating R t is as follows: where ( ) s t is the fraction of the host population sensitive to time t.

Description of the Epidemiological Situation in the Central African Republic
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 Figure 1, the cumulative curve weighs on the data for epidemiological weeks, this fact that we do not visualize very well the number of people who have died from COVID-19 (May 19, 2020, date of the first case of mortality linked to COVID-19, then May 29, 2020 and June 04, 2020) over the entire period of our study.

Estimation of Parameters
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. Table 1 describes the epidemiological dataset on COVID-19 in the Central African Republic. From these data, we estimated the parameters of our model (1).
We can clearly see in Table 1 that the prevalence of the COVID-19 pandemic is 11% in just 4 months. Figure 2 shows the number of new cases confirmed by the Rt-PCR technique in the laboratory. The very first COVID-19 case in the country was detected on March 16, 2020.   We used the maximum likelihood method to estimate the parameters of our model [23] for Table 2.

Analytical Results
We used R software with R 0 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:

Discussion
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 Table 3, for a weak control or non-compliance with the social distancing and masking, the value of 0.22 β = and 0 1.45 1 R = > . Then the epidemic will continue in the population.
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 0 0.85 1 R = < , 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 Table 4, we have made the scenario on the cure rate of symptomatic and asymptomatic, assuming that the cure rate in asymptomatic is higher than symptomatic. And there, we see that if the cure rate of asymptomatic patients is very high, the epidemic stabilizes. Figure 3 shows the incidence and the number of reproduction R t in real time over the study period.
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.

Conclusions
In this paper, our objective is to estimate the number of reproduction in real time R t 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.