Abstract
The current COVID-19 pandemic has caused several damages to the world, especially in public health sector. This study considered a simple deterministic SIR (Susceptible-Infectious-Recovered) model to characterize and predict future course of the pandemic in the West African countries. We estimated specific characteristics of the disease’s dynamics such as its initial conditions, reproduction numbers, true peak, reported peak with their corresponding times, final epidemic size and time-varying attack ratio. Our findings revealed a relatively low proportion of susceptible individuals in the region and in the different countries (1.2% across West Africa). The detection rate of the disease was also relatively low (0.9% for West Africa as a whole) and < 2% for most countries, except for Cape-Verde (9.5%), Mauritania (5.9%) and Ghana (4.4%). The reproduction number varied between 1.15 (Burkina-Faso) and 4.45 (Niger) and the peak time of the pandemic was between June and July for most countries. Most generally, the reported peak time came a week (7–8 days) after the true peak time. The model predicted 222,100 actual active cases in the region at peak time while the final epidemic size accounted for 0.6% of the West African population (2,526,700 individuals). Results obtained showed that the COVID-19 pandemic has not severely affected West Africa as noticed in other regions of the world, but current control measures and standard operating procedures should be maintained over time to ensure trends observed and even accelerate the declining trend of the pandemic.
Introduction
The COVID-19 pandemic is a severe acute respiratory syndrome caused by the new coronavirus, SARS-CoV-2, which has emerged from Wuhan, Hubei Province (China) in the last days of 2019 [25,39]. It is currently the most important threat to global public health. By August 15th, 2020, about 21,026,758 total confirmed cases and 755,786 deaths were recorded worldwide [44]. The disease was rapidly spread around the world (about 212 countries) [27] including the 54 African countries. By mid August 2020, The World Health Organisation (WHO) reported 936,062 and 152,483 confirmed cases and 18,286 and 2,351 deaths across Africa and the West-African region, respectively [44]. Healthcare services in the region have particularly faced critical time sensitive decisions regarding patients and their treatment [28]. It is clear that the COVID-19 pandemic has severely affected people’s life, health and economy. Actually, it led to an important increase in the demand of hospital beds and artificial respirators (mechanic and non-invasive). According to the WHO global health observatory data, most countries in West Africa have less than 5 hospital beds and 2 medical doctors per 10,000 of the population, while 50% of the countries have health expenditures lower than US$50 per capita [46]. In contrast, European countries such as Italy and Spain, 34 and 35 hospital beds, respectively, per 10,000 of the population are noted, 41 medical doctors per 10,000 of the population, and US$2,840 and US$2,506 per capita expenditure [34]. Moreover, medical staff in the World were directly exposed to infections [26]. Since vaccines are still under development, and antiviral drugs are not available for effective curative treatment of COVID-19 infections, the actual cure practice is hospitalization and intensive care unit management [32]. Prevention measures used are essentially non-pharmaceutical interventions such as regular hand washing with soap, mask wearing and social distancing. To be efficient, these non-pharmaceutical measures required a good understanding of the dynamic of the spread of the disease.
Mathematical and statistical models can be extremely helpful tools to make decisions in public health. They are also important to ensure optimal use of resources to reduce the morbidity and mortality associated with epidemics, through estimation and prediction [27–30,33]. The prediction of essential epidemiological parameters such as the peak time, the duration and the final size of the outbreak is crucial and important for policy makers and public health authorities to make appropriate decisions for the control of the pandemic [32]. Therefore, modelling and forecasting the numbers of confirmed and recovered COVID-19 cases play an important role in designing better strategies for the control of the COVID-19’s spread in the world [27, 29]. Since the appearance of the first COVID-19 case in the world, several studies have been conducted to model the dynamics of the disease. The main methods used were: deterministic modelling techniques (SIR, eSIR, SEIR, SEIRD, etc. compartmental models) [24,29,30,33,41–43], autoregressive time series models based on two-piece scale mixture normal distributions [27], stochastic modelling methods [49, 50], machine learning techniques [47, 48], growth models [43,52,53] and bayesian method [48].
Among these modeling techniques, deterministic models are the most considered because of their simplicity. However, they fail to provide accurate results due to non-identifiability problem when the number of compartments and the number of parameters are high [39]. Actually, complex deterministic models have been showed to be less reliable than simpler model such as SIR model framework [39], which performs better in describing trends in epidemiological data. This under-performance may be worse when meta-population confirmed-cases data are considered. However, only few studies related to COVID-19 in Africa used mathematical models and prevalence data to study the dynamics, analyze the causes and key factors of the outbreak [24,25,51]. Recently, [24] assessed the current pattern of COVID-19 spread in West Africa using a deterministic compartmental SEIR-type model.
In this study, we used a simple deterministic SIR model to characterize and predict future trend of the spread of the pandemic in West Africa. Specifically, we aimed to estimate specific characteristics of COVID-19 dynamics (initial conditions of the pandemic, reproduction numbers, true peak, reported peak and their times and dates, final epidemic size and time-varying attack ratio). The originality of this work is that it focuses on the sixteen West African countries and the whole region as well. It is the first study dealing with the dynamic of the pandemic in each of the West African countries.
Methods
Model description
Problems of non-identifiability in parameters estimation in deterministic compartment models (especially complex models) are common in epidemiological modelling studies, which often imply biased estimations of parameters. [39] recommended the use of simpler models which overperformed complex models in estimating reliable parameters. Hence, in this study, the simple SIR model [31] was considered with the particularity of two removal rates and illustrated in the system below [17]:
with initial conditions,
In (1), S = S(t), I = I(t) and R = R(t) represent the number of susceptible, infected and removed individuals at time t, respectively while N is defined as the total population size for the disease transmission. The parameters β,ν1 and ν2 are the transmission rate, the removal rate of reported infected individuals (detected) and the removal rate of infected individuals due to all other unreported causes (mortality, recovery or other reasons), respectively. We considered the removal rate ν2 as constant with value ν2 = 1/10 [54].
From the second differential equation of (1), one can notice that ν1I0 represents the daily confirmed cases (Ir0) at time 0 of the outbreak. Thus, the relationship between the initial number of infected individuals and the detection rate, ν1 is as follow and used in the estimation process:
Data consideration and parameter estimation procedure
For each country, the data considered for the modeling spans the period from the date of detection of the first case of COVID-19 in the country and August, 12, 2020. Data considered were the daily numbers of reported cases that have been assimilated to ν1I.
These data were downloaded from the Global Rise of Education website [55]. Table 1 presents the demographic patterns [11], initial dates of the pandemic [55] and testing efforts (identification of new cases) of the countries [11]. We fitted the model (1) to the observed daily cases to study the dynamics of COVID-19 pandemic in the sixteen West African countries.
To improve the prediction power of (1), we used a cross-validation procedure of parameter estimation, where 90% of the observations were considered to estimate values of the six unknown model parameters (S0,I0,R0,β,ν1,ν2) and the remaining observations were used to validate the model. The Root Mean Square Error (RMSE) statistics was used as the measure of estimation precision:
where 
and θ are the predicted and observed number of daily cases respectively; k is the number of observations considered. We considered as RMSE1, the Root Mean Square Error computed on the 90% of the observations and RMSE2, computed on the remaining observations (10%)
The solutions of (1) were obtained using the built-in function ODE45 of Matlab [56]. Then, the non linear least square estimate technique was performed to estimate the six parameters in (1) given starting values, using the built-in function fminsearchbnd of Matlab [56].
Afterwards, we simulated 2,000 different starting values of the six parameters using a resampling method (function resample of Matlab) for S0, I0 and R0 and the uniform distribution (function rand of Matlab) for β and ν1. Then, we estimated for each of the starting points, values of the six parameters using the non-linear least square technique described above. The final values considered for the these parameters are those related to both smallest values of RMSE2 and RMSE1 to guarantee a good fit of the model and a good predictive power. At the end of this process, we obtained reliable estimates of the model parameters with 95 % confidence intervals. Curves were plotted to show evolution trends of predicted daily new COVID-19, daily reported cases and the attack ratio for the 16 countries. With the estimated values of the six parameters from (1), the COVID-19 dynamic in each country was characterized computing the following parameters with their 95% confidence intervals. Table 1 presents demographic patterns and testing efforts in the West African region.
– Reproduction number,
[1, 16]: It is the average number of new cases of infection, caused by an average infected individual (during his infectious period), in a fully constituted population:
– Running reproductive number,
[1]: measures the number of secondary infections caused by a single infected individual in the population at time t.
– True peak size, npp and True peak time, Tpp. The true peak size indicates the largest daily number of new infectious cases in the population:
while the true peak time, Tpp represents the time at which the largest daily new infected cases is obtained. These two parameters were determined numerically.– Reported peak size, nrp and Reported peak time, Trp. The reported peak size indicates the largest number of daily reported cases:
while the reported peak time is the associated time to nrp. They were determined numerically.– Maximum number of active cases, Imax: since I0,R0 << S0, we assumed the number of initial susceptible individuals to be approximately equal to N (S0 ≈ N). Thus, Imax can be approximated as follow [16]:
– Final epidemic size, Itotal [16]: it is the total number of cases over the course of the epidemic wave.
S∞ can be approximated considering the entire population as initially susceptible (S0 ≈ N); hence, following [16]:
For each country, the equation (7) was solved numerically to determine S∞ through an iterative process.
– Attack ratio, Ar [17]: is the fraction of susceptible population that becomes infected. It is calculated along the epidemic wave as follow:
For the West Africa as a whole, S0, I0 and R0 were first computed by summing the corresponding estimated values of the 16 countries. Afterwards, the model (1) was fitted to daily reported cases of the region using initial conditions computed. This allowed the estimation of β and ν1 and computation of the characteristics of COVID-19 dynamics across the region.
Results
Current patterns of COVID-19 transmission in West Africa
Table 1 reveals a great heterogeneity in the region in terms of population density and testing efforts. Countries like Cape-Verde, Mauritania and Ghana have put relatively more effort into identifying infected individuals while Niger, Nigeria and Guinea are countries with the lowest number of tests per 1 million people. Combining both the testing effort and the mean number of reported cases per test indicate a relatively less testing effort to identify many infected individuals (Guinea and Gambia). On the other hand, countries like Cape-Verde, Benin and Togo put much effort to find few COVID-19 cases (Table 1).
Results obtained from the estimation of initial conditions of COVID-19 pandemic in West Africa revealed the relatively low proportion of the susceptible individuals in most countries (about 1% of their total populations). However, countries such as Guinea-Bissau, Gambia and Cape-Verde showed relatively large proportion of susceptible individuals to COVID-19 with 17.5%, 6.0% and 2.4% of their total populations, respectively. The proportion of the susceptible individuals across West Africa was also relatively low (1.2%) (Tables 1 and 2). Moreover, before the detection of the first cases, infected individuals were present in the population of all the countries with some already recovered individuals. The detection rate of infected individuals was relatively low (less than 1%) for Benin, Burkina, Mali, Niger, Nigeria, Sierra-Leone and West Africa as a whole. However, some countries like Cape-Verde (9.5%) Mauritania (5.9%) and Ghana(4.4%) recorded the highest detection rates, respectively. (Table 3).
Estimates with 95% confidence intervals of the initial parameters of the SIR-model; ν2 = 0.10 is constant for all countries
In most countries, the model estimated an average of 1 new case of infection caused by an infected individual during his infectious period (
) except for Sierra-Leone, Nigeria, and Côte d’Ivoire with 
≈ 2 and Niger, which recorded the highest reproduction number (
) (Table 3).
Long term dynamics of COVID-19 in West Africa
We analyzed the long term dynamic of COVID-19 in West Africa by first focusing on the true peak of the pandemic. In general, the estimated reported peak time came a week (7–8 days) after the true peak time in all the countries while their estimated reported peak sizes accounted in average for 21% of the estimated true peak size (Table 3). Most countries have already experienced the peak of the epidemic wave. The true peak time was estimated in June for Sierra-Leone (14th), Mauritania (19th), Benin (27th), Côte d’Ivoire (24th) and Ghana (13th), while it was estimated in July for Liberia (4th) and Nigeria (5th), Senegal (18th), Guinea (23th) and Cape-Verde (22th) and October 4th for Togo. Niger recorded the earliest true peak time (April 8th) while the latest true peak time was on December 10th, 2020 for Gambia (Figure 1 and Table 3). The true peak time across the region was July 1st with 25,267 new cases (Table 3 and Figure 4a and 4b). Margin error of the estimated true peak time varies from 1 day to 16 days with average value of 5 days in the region. Half of the countries (8 out of 16) recorded true peak size less than 1000 new cases cases at peak time. The highest numbers of new cases at peak time were estimated at 19,021 and 17,703 for Niger and Nigeria respectively (Table 3 and Figure 1). The estimates of the reported peak size was 1,891 daily cases across the region (Table 3 and Figures 4a-b).
Evolution trend of the COVID-19 daily cases per country.
The final epidemic size would account for 0.6% of the population of West Africa. This estimate is generally low (< 1% of the population size) for more than half of the countries though Gambia, Guinea-Bissau and Cape-Verde would record the highest final epidemic sizes (> 9% of their populations). The estimates of the maximum number of actual daily active cases at peak time are for most countries greater than 1,000 cases, though, it is 107,200 and 195,000 for Niger and Nigeria respectively (Table 3).
The running reproduction number helped assess the evolution trend of the disease. It decreased over time in all countries from the beginning of the outbreak (1.2–4.5) to a stability point, which varied according to countries (0.50–0.82, Figure 2). As expected, the fraction of susceptible individuals being infected (attack ratio) increased over time from 0% to 40%–70%, depending on countries. These evolving trends in reproduction number and attack ratio are similar to those noted for West Africa as a whole (Figure 4b-c).
Running reproduction number per country in West Africa.
Evolution trend of the attack ratio per country in West Africa.
Trend in COVID-19 dynamics across West Africa. (a) Prediction of true peak and reported peak date and size of COVID-19. (b). Evolution trend of the reproduction number. (c) Evolution trend of the attack ratio.
Discussion
In epidemiology, understanding the dynamics of an epidemic outbreak and predicting its future course is a major research question, which is often studied using modelling techniques [2, 18–23]. Estimation and prediction rely on mathematical and statistical models, which inform public health decisions and ensure optimal use of resources to reduce the morbidity and mortality associated with epidemics [14,27–30,33]. For instance, estimation of epidemiological parameters and prediction on the Influenza outbreak dynamics in Canada was done using the Richard’s model [14], while a three-parameter logistic growth model was used to study and forecast the final epidemic size in real-time of the Zika virus outbreaks in Brazil from 2015 to 2016 [2].
In this study, we used deterministic SIR model to understand COVID-19 dynamics in West African countries and estimated the overall number of susceptible individuals that accounted for 1.2% of West African population and 1% for most countries, except Guinea-Bissau, and Gambia where the susceptible individuals account for more than 9%. In general, small countries with relatively high population density are those with high proportion of susceptible individuals, indicating how high population density with small area can affect epidemics dynamics [57, 58]. Our findings however, revealed a great disparity between countries in terms of testing rate of COVID-19. Countries like Guinea and Gambia and in less extent, Côte d’Ivoire and Nigeria, showed a relatively less testing effort to identify many infected individuals. This suggests that there may not be enough tests being carried out to properly monitor the outbreak [55]. In contrast, countries like Cape-Verde, Benin and Togo, which have recorded less than or equal to 50 positive cases per 1,000 tests, seem to be effectively controlling the pandemic according to the WHO criteria [45]. Compared to relatively wealthier countries like Australia, South Korea and Uruguay, it takes hundreds of tests to find one case [55].
The detection rate considered in model (1) is a better indicator of testing effort since it represents the proportion of active cases in the population that are identified daily. Our results revealed a relatively low detection rate of COVID-19 in West Africa with less than 2% in most countries except Cape-Verde, Mauritania and Ghana (> 5%). These three countries are also the ones with the highest testing rates (see table 1) confirming the link between detection rate and testing effort [5]. Thus, fairly low detection rates in most West African countries demonstrate low testing effort and may be explained by a number of factors, including the availability of testing kits and qualified healthcare workers and low ability to control the disease due to their low GDP. For instance the average detection rate of COVID-19 in the world in April 2020 was estimated at 6% [59]. It is also useful to note that the estimated average detection rate hides a great variability in the testing effort over time. Indeed, it is generally accepted that the testing rate is relatively low at the very beginning of an epidemic outbreak but can increase rapidly over time when a better response mechanisms are put in place [60].
The dynamics of the COVID-19 pandemic in West Africa shows a reproduction number greater than 1 in all countries (from 1.2 in Burkina-Faso to 4.4 in Niger) while it was 1.4 for the whole region. This value is relatively low compared to that estimated by [24] (1.6) for the same region in June in a recent study using a modified SEIR model and thus reveals either a declining trend in the pandemic over time or the result of using a different modeling approach to estimate the reproduction number. The COVID-19 pandemic appears more serious than the Ebola outbreak in Africa given the basic reproduction number. Indeed, [8] estimated the basic reproduction number at 1.1, 1.2 and 1.2 for Guinea, Liberia and Sierra-Leone, respectively against 1.4, 1,3 and 1.6 for the same countries as far as COVID-19 is concerned. These comparisons indicate that COVID-19 is on average 1.29 times more reproducible than Ebola in these countries.
The trend in the running reproduction number reveals a rapid declining in all countries except Burkina-Faso and reveals some efficiency of the control measures put in place and being implemented in the different countries. Most countries have already experienced the peak in the new COVID-19 cases in June and July. In general, the reported peak was very low compared with the true peak.
However, current control measures should be maintained overtime to avoid possible second wave as observed in other parts of the world.
Conclusion
Our study shows that the novel COVID-19 pandemic, although highly contagious has not seriously impacted West Africa in terms of prevalence, compared to other parts of the world, in particular the Europe and the USA. Actually, the total number of susceptible individuals and final epidemic size account for 1.2% and 0.6% of the total population size of West Africa, respectively. But the relatively low reported cases are related to very low testing effort in the West African countries. The study also indicates a relatively low detection rate and for most countries in the region, the dates of the true peak of infection have passed (June-July 2020). Nevertheless, the pandemic is still ongoing in the region and it is important that the non-pharmaceutical measures currently in place continue over time to help reduce its spread dynamics, pending adequate effective treatment.
Data Availability
I declare that all data referred in the manuscript are available.
Authors’ contributions
RGK conceived the ideas with SHH, designed the methodology, wrote the codes and supervised the work; SHH gathered the data, run the codes for some countries and wrote most sections of the manuscript; HBT run the codes for some countries and write some sections of the manuscript. All authors reviewed the drafts and gave final approval for publication.
Conflicts of Interest
The authors declares no conflict of interest.
Acknowledgments
HBT acknowledges the support of IMU-CDC through GRAID program. RGK acknowledges the support from the African German Network of Excellence in Sciences (AGNES) and the Alexander von Humboldt Foundation (AvH).
