key: cord-261599-ddgoxape
authors: Nabi, Khondoker Nazmoon; Abboubakar, Hamadjam; Kumar, Pushpendra
title: Forecasting of COVID-19 pandemic: From integer derivatives to fractional derivatives
date: 2020-09-21
journal: Chaos Solitons Fractals
DOI: 10.1016/j.chaos.2020.110283
sha: 
doc_id: 261599
cord_uid: ddgoxape

In this work, a new compartmental mathematical model of COVID-19 pandemic has been proposed incorporating imperfect quarantine and disrespectful behavior of the citizens towards lockdown policies, which are evident in most of the developing countries. An integer derivative model has been proposed initially and then the formula for calculating basic reproductive number [Formula: see text] of the model has been presented. Cameroon has been considered as a representative for the developing countries and the epidemic threshold [Formula: see text] has been estimated to be ∼ 3.41 [Formula: see text] as of July 9, 2020. Using real data compiled by the Cameroonian government, model calibration has been performed through an optimization algorithm based on renowned trust-region-reflective (TRR) algorithm. Based on our projection results, the probable peak date is estimated to be on August 1, 2020 with approximately 1073 [Formula: see text] daily confirmed cases. The tally of cumulative infected cases could reach ∼ 20, 100 [Formula: see text] cases by the end of August 2020. Later, global sensitivity analysis has been applied to quantify the most dominating model mechanisms that significantly affect the progression dynamics of COVID-19. Importantly, Caputo derivative concept has been performed to formulate a fractional model to gain a deeper insight into the probable peak dates and sizes in Cameroon. By showing the existence and uniqueness of solutions, a numerical scheme has been constructed using the Adams-Bashforth-Moulton method. Numerical simulations enlightened the fact that if the fractional order α is close to unity, then the solutions will converge to the integer model solutions, and the decrease of the fractional-order parameter (0 < α < 1) leads to the delaying of the epidemic peaks.

In the late December 2019, the World Health Organization (WHO) was enlightened about sudden new cases of a unique pneumonia and the specific cause of this infection was an enigma for the health authorities of the city of Wuhan in China. A few days later, it has been unearthed that this pneumonia is caused by a virulent virus and it has been named officially as the novel coronavirus (2019-nCoV). The modes of transmission are similar to the viruses responsible for the previous epidemics of SARS and MERS [31] . However, this virus is more contagious and is the cause of millions of deaths worldwide. To date, there is no clinically established treatment or specific vaccine for COVID-19. Nevertheless, according to the suggestions of different prudent infectious disease specialists, several countries have started using chloroquine combined with azithromycin as alternative drugs [15] . Several non-pharmaceutical interventions such as physical distancing, wearing face masks in public places, home quarantine, isolation and countrywide lockdown policies have been promoted to curb the spread of the virus [27] .

Cameroon, a country in sub-Saharan Africa, had its first confirmed case on March 05, 2020. It was a 58-year-old French citizen, arrived by Yaoundé airport on February 24, 2020 [22] . With an aim to quell the spread of the virus, the government undertook several policies such as closing the borders, confinement measures, prohibiting large gatherings of people, closures of all kinds of educational institutions (from kindergarten to university), wearing face masks in public places and so on. In addition, Cameroonian government implemented contact-tracing strategy together with mass-media sanitisation campaigns. Contact-tracing consists of search, or identify all individuals who have had a contact with a confirmed case and force them to be quarantined if they are COVID-19 positive [24] . The scarcity of medical resources in maximum sub-Saharan African countries compelled authorities to request asymptomatic cases of COVID-19 to be self-quarantined in their home while taking medication (conventional or not) following proper health guidelines [16] . The main problem encountered by health workers is the stigmatization of COVID-19 patients. Indeed, like the cases of HIV/AIDS in Africa, stigma pushes confirmed cases or those who have had direct contact with a confirmed case of COVID-19 to run away from hospitals, and in this way they are contributing to the rapid spread of the virus in the community [23, 37] .

Different types of mathematical models have been playing notable role in predicting the transmission dynamics of infectious diseases and different effective control measures can be designed to limit the community transmission. Several models have already been proposed to predict the evolution of COVID-19 and to study the impact of control measures imposed by different governments [1, 26] . Although these models have similarities, each of them presents specifications related to the evolution of COVID-19 and various control measures deployed by several respective governments. A model has been considered in [1] , describing the interactions among the bats and unknown hosts, human population and population of virus in reservoir (seafood market). Atangana-Baleanu derivative concept has been implemented successfully. However, confined or quarantined class and isolated class have not been considered in this study. In a recent study, Nabi [26] has proposed a new Susceptible-Exposed-Symptomatic Infectious-Asymptomatic Infectious-Quarantined-Hospitalized-Recovered-Dead (SEI D I U QHRD) compartmental mathematical model and calibrated model parameters to project the future dynamics of COVID-19 for various COVID-19 hotspots. Although he considered a quarantine compartment, he did not consider the possibility of those quarantined people to contract the virus. Although this consideration goes well with certain developed countries (France, Italy, Russia and Spain), certain realities [27] of COVID-19 pandemic in several developing countries and maximum African countries have not been taken into account in any model. In case of Cameroon, several crucial factors such as violation of containment measures and negligence towards confinement measures must be taken into account to understand the transmission dynamics of COVID-19. Importantly, numerous daily wage-earners are compelled to go out to fetch food for their family members contravening all confinement orders. In this way, they are getting infected and spreading the viruses to their family members as well as neighbors [27] . It is therefore prudent to envisage future outbreak dynamics of COVID-19 in developing countries with imperfect confinement. Moreover, some people flee quarantine in hospitals to return to their families due to the stigmatization of COVID-19 patients. Henceforth, considering perfect quarantine would be highly debatable if we want to model and forecast the transmission dynamics of COVID-19 in developing countries.

Fractional calculus or non-integer order calculus, a branch of mathematical analysis, contains the theory of fractional-order derivatives operators. The history of non-integer order calculus is more than 300 years old and modern fractional calculus is a rapidly evolving field both in theory, analysis and applications with a view to handling complicated real-world problems. A plethora of concepts regarding fractional derivatives have already been introduced and applied by researchers in various branches of Science and Engineering [13, 17, 32] . Atangana-Baleanu (AB), Caputo-Fabrizio (CF), and Caputo derivatives are the most commonly used derivatives in solving real-world problems. Grunwald-Letnikov fractional operators have been applied comprehensively in the field of image processing. AB, CF, and Caputo fractional derivatives have different Kernel properties. Caputo derivative is defined with power-law type kernel (non-local but singular), Caputo-Fabrizio with exponentially decaying type kernel (non-singular) and AB derivative in the Caputo sense has Mittag-Leffler type kernel [36] . A new generalised Caputo type fractional derivative has been introduced in [30] . Asymptotic stability of generalised Caputo FDEs has been introduced by Baleanu in [5] . Integer-order derivative operators have lucid geometric and physical interpretations, which clearly simplify their applications in tackling real-world challenges. Recently, various techniques have been proposed by researchers to describe the solution of non-linear fractional differential equations. A new technique for solving non-linear Volterra integro-differential equations in Atangana-Baleanu derivative sense has been introduced in [12] . Another technique with the uses of Genocchi polynomials has been proposed in [35] . A new technique to solve multi-variable order differential equations in AB sense has been introduced by Ganji et al. [11] .

In this paper, a new compartmental COVID-19 model has been studied rigorously with the help of Caputo fractional derivatives. The advantage of applying Caputo fractional derivatives to solve the proposed COVID-19 model is the dynamics of the model can be observed more deeply using the real-time Cameroon data. In this framework, real-time data can be compared with the model outputs in a precise manner. Fractional derivatives can be a suitable alternative to integer derivatives in studying complex dynamics as there are always some limitations of integer derivatives described in [2] . Podlubny has given the geometric and physical significance of Riemann-Liouville fractional integral in [33] . As compared to the Riemann-Liouville fractional derivative, Caputo fractional derivative has the identical physical interpretation [33] . In a physical model, fractional derivatives manage systems with memory that ensure that the evolving system state is reliant on its past states also.

The aim of this work is to forecast the probable time and size of the epidemic peaks of the novel coronavirus outbreak in Cameroon by studying a realistic compartmental model using the robust concept of Caputo fractional derivative. A compartmental mathematical model of COVID-19 progression dynamics has been proposed incorporating the efficacy of confinement measures and imperfect quarantine. Estimation of parameters has been performed by using real-time data, followed by a projection of the evolution of the disease. Global sensitivity analysis is applied to determine the influential mechanisms in the model that drive the transmission dynamics of the disease. Afterwards, the integer model is transformed into a fractional model in Caputo sense and the existence and uniqueness of solutions have been presented. Using the Adams-Bashforth-Moulton scheme, a numerical scheme for the fractional model has been constructed. Several numerical simulations are performed to compare the result obtained by the integer-order model with the fractional model outputs.

The entire paper is organized as follows. Model formulation and basic properties are presented in Section 2. Section 3 is devoted to model calibration using real data of reported cases of COVID-19 in Cameroon, global sensitivity analysis of the proposed model, and forecasting of the disease future dynamics. In section 4, the compartmental model is transformed into a fractional model in the Caputo sense, following by the existence and uniqueness of solutions, construction of a numerical scheme and numerical simulations in which the coefficient of fractional order vary. The paper ends with some insightful observations and fundamental findings.

A compartmental differential equation has been proposed to describe the transmission dynamics of the COVID-19. The spread starts with an introduction of at least one infected human into a susceptible population. The human population has been categorized into seven compartments according to their epidemiological status, i.e {S(t), C(t), E(t), A(t), Q(t), H(t), R(t)} which represent the number of susceptible individuals, confined individuals, infected individuals in incubation period, asymptomatic infectious individuals (undetected but infectious) including those who fled quarantine, quarantined individuals or confirmed infected individuals, hospitalised and recovered cases. So, the total human population, denoted by N , at any time can be represented by N (t) = S(t) + C(t) + E(t) + A(t) + Q(t) + H(t) + R(t). The formulation of the model are subject to some assumptions listed bellow.

• Migration is not considered. Indeed, concerning migration, the initial time of our model is taken when the disease is already inside the region or country and that international migrations are prohibited [26] .

• Vital dynamic (birth and natural deaths) is ignored. The idea is to observe what may happen in the short-term. Indeed, since for example to have a birth we need around 9 months on an average, to have a natural death rate which is the inverse of the average life expectancy i.e. 1/55 per year in some developing country like Cameroon and no one would like the disease to stay so long.

• Confinement is not perfect. Indeed, in some developing countries, because of insufficient financial income in households, people do not respect the confinement measures imposed by the government, and are forced to go out to work to provide for their families.

• The immunity is perfect. Indeed we assumed, without confirmed proof of lost of immunity, that recovered individuals develop certain immunity against the disease [9, 29] .

• Hospitalized patients (in respiratory support) can not spread the disease as they remain under closed supervision. The flow diagram of the proposed model is depicted on Figure 1 , where susceptible individuals decrease either via confinement at a rate c, or via infection by a direct contact with an infectious individual (A or Q).

The transmission rate from unquarantined asymptomatic carriers (A) (respectively quarantined symptomatic carriers (Q)) is β (respectively ηβ). The modification and dimensionless parameter 0 ≤ η < 1 accounts for the assumed reduction in transmissibility of quarantined symptomatic carriers relative to unquarantined asymptomatic carriers. Since the confinement is not perfect, confined individuals move to the latent class at a rate (1 − )λ where 0 ≤ ≤ 1 is the efficacy of confinement. This is in contrast with a variety of SEIR-style models recently employed in, e.g. [9, 29] , where the authors consider that confined individuals are not exposed to an infection. Indeed, in developing country such as Cameroon, since the confinement is not perfect, it is practical that an infected individual can infect a confined person by returning back home. This phenomenon was also observed at the beginning of the epidemic in a developed country such as Italy, where old confined individuals had been infected by their grandchildren who are asymptomatic carrier [16] . After the incubation period of 1/γ, a proportion q of infected individuals move to the asymptomatic class A and the other (1−q) to the quarantine class. It important to note that asymptomatic class A include all infectious individual who are not in quarantine and those who has fled the quarantine. The rate of transition from unquarantined asymptomatic infectious to quarantine infectious, unquarantined asymptomatic infectious to hospitalised class, unquarantined asymptomatic infectious to recovered class and hospitalised to recovered class are r 1 σ 1 , r 2 σ 1 , (1 − r 1 − r 2 )σ 1 and σ 3 respectively. r 3 σ 2 , r 4 σ 2 and (1 − r 3 − r 4 )σ 2 are transition rates of quarantined asymptomatic infectious to unquarantined asymptomatic infectious, quarantined asymptomatic Parameter Description β Infectious contact rate η

Infectiousness factor for quarantined infected carrier p Transition rate from confined class to unconfined class c Transition rate from unconfined class to confinement class Confinement efficacy γ

Transition rate from exposed class to infectious class q Fraction of exposed that become quarantined carriers σ 1 Transition rate from unquaratined class to quarantined class and recovered class σ 2

Transition rate from quaratined class to unquarantined class, hospitalised and recovered σ 3 Transition rate from hospitalised class to recovered class r 1

Fraction of unquaratined infectious that become quarantined infectious r 2

Fraction of unquaratined infectious that become hospitalised r 3 Fraction of quarantined infectious that become unquarantined infectious r 4 Fraction of quarantined infectious that become hospitalised δ a Disease induced death rate, unquarantined infectious δ q Disease induced death rate, quarantined infectious δ h Disease induced death rate, hospitalised infectious to hospitalised class, and quarantined asymptomatic infectious to recovered class respectively. δ a , δ q and δ h are the disease-induced death rate of unquarantined asymptomatic infectious, quarantined asymptomatic infectious and hospitalised patients respectively. The above assumptions led to the following nonlinear system of ordinary differential equations,

The parameters are described in Table 1 .

We set x = (S, C, E, A, Q, H, R) the vector of state variable, Let f : R 7 → R 7 the the right hand side of system (2), which is a continuously differentiable function on R 7 . According to [39, Theorem III.10 .VI], for any initial condition in Ω, a unique solution of (2) exists, at least locally, and remains in Ω for its maximal interval of existence [39, Theorem III. 10 .XVI]. Hence model (2) is biologically well-defined.

Here, we prove that all state variables of model (2) are non-negative for all time, i.e, solutions of the model (2) with positive initial data remain positive for all time t > 0. The following result can be obtained.

is a solution of (2) with positive initial conditions. Let us consider E(t) for t ≥ 0. It follows from the third equation of system (2) that

Since E(0) ≥ 0, it follows that E(t) ≥ 0 for t ≥ 0. We proceed with the same for A(t), Q(t), H(t) and R(t). Now, it remains to prove that S(t) and C(t) are also positive. Assume the contrary, and lett such that S(t) = 0 and C(t) ≥ 0. From the first equation of (2), it follows that dS(t) dt | t=t = pC(t), which means that S(t) < 0 for t ∈ (t − ζ,t) with ζ a small positive constant. This leads to a contradiction. Thus S(t) ≥ 0 for t ≥ 0. We proceed with the same to prove that C(t) ≥ 0. (2) with positive initial conditions are bounded by the total population N 0 .

where N 0 is equal to the total population. It follows that for all t ≥ 0, we have S(t) ≤ N 0 ,

In what follows, we study the model (2) in the following set

which is positively-invariant and attracting region for the model (2).

Now, we define the following manifold, in which any point is a disease-free equilibrium (DFE) of the model (2)

Note that W ⊂ D. In the following, we will work with the largest disease-free equilibrium point denoted by x 0 (see [14] ). Assuming S + C = N 0 , it follows that

Using notations in [38] , matrices F and V for the new infection terms and the remaining transfer terms are, respectively, given by

Then, the control reproduction ratio is defined, following [38] , as the spectral radius of the next generation matrix, F V −1 :

where ρ(·) represents the spectral radius operator.

The formula for control reproduction number has been formulated. Indeed, the insightful epidemic threshold, R c calculates the average number of new secondary COVID-19 cases generated by a COVID-19 positive individual in a population where a certain fraction of susceptible people are confined. The control of COVID-19 pandemic is governed by the application of some control measures which contribute to bring down R c under unity [38] . Hence, we claim the following result which is a direct consequence of the next generation operator method [38, Theorem 2] .

Lemma 3. If R c < 1, the disease-free equilibrium x 0 is locally asymptotically stable and unstable if R c > 1.

Proof. Note the two last equations of the system (2) are uncoupled to the remaining equations of the system. Since the total population, N 0 , is constant, we have S + C = N 0 − (E + A + Q + H + R). So, the local stability of the COVID-19 model (2) can be studied through the remaining system of state variables (E, A, Q). Thus, we obtain that the Jacobian matrix associated to these variables is given by

The eigenvalues of J are the roots of the following characteristic polynomial:

Note that a 2 is always positive. a 0 and a 1 are positive provide that R c < 1. Thus, it follows that J has all its eigenvalues with negative real parts. Hence, the disease-free equilibrium, x 0 , is locally asymptotically stable whenever R c < 1.

Remark 1. Lemma 3 implies that if R c < 1, then a sufficiently small flow of infected individuals will not generate an outbreak of COVID-19, whereas for R c > 1, epidemic curve reaches a peak by growing exponentially and then decreases to zero as t → ∞.

The better control of the COVID-19 can be established by the fact that the DFE x 0 is globally asymptotically stable (GAS). In this context, we claim the following result.

Theorem 1. if R c < 1, then the manifold, W, of disease-free equilibrium points of the model (2) is GAS in D.

Proof. Assume that R c < 1. Let X = (S, C, R) and Y = (E, A, Q) . Since

we deduce that

Thus from [38, Theorem 2] , it follows that s(F − V ) < 0 if and only if R c < 1 where s(M ) is the stability modulus of the matrix M . Therefore, the trajectories of the auxiliary system whose right-hand side (6) converges to zero whenever R c < 1. Since all the state variables are non-negative and F − V is a Metzler matrix, by the comparison theorem [18] , it follows that lim shows that W is an attractive manifold. Moreover, W is locally asymptotically stable when R 0 < 1. We conclude that in system (2), the manifold W is globally asymptotically stable when R c < 1.

In the absence of confinement measures, i.e. = 0, R c converges to the basic reproduction number, R 0 , given by

Using (5), it follows that

In the next section, numerical simulations have been performed rigorously. Parameters considered in our simulations are either estimated from real data or related to the COVID-19 outbreak data in Cameroon.

As in many countries in sub-Saharan Africa, the actual situation of the COVID-19 in Cameroon and its economic impact are mixed. Indeed, after the appearance of the first case confirmed on March 07, 2020 until the introduction of control measures by the government of Cameroon (closure of public places, in particular schools, universities, drinking places and other places of entertainment from 6 p.m. and so on), the statistics advanced by the Cameroonian health authorities do not reflect, according to certain non-governmental organizations, the actual situation of the COVID-19 epidemic in Cameroon. Since no nongovernmental organization is empowered to carry out tests, the only statistics we use here to calibrate our model are those communicated by the Cameroonian ministry of public health during her daily press briefings [25] . Also, these control measures decreed by the government have really impacted certain economics activity sectors. Indeed, the fact of the closing of the drinking places at other places of distractions, without forgetting the prohibition of the gatherings of more than 50 persons, a lot of people found themselves in technical unemployment, and of this fact was forced to leave for go to affect activities that may allow them to support themselves and their families.

The number of COVID-19 positive cases (NPC) per day, the number of recovered cases (NRC) per day and the number of death cases (NDC) per day in Cameroon, have been taken in the summarized table online on the evolution site of COVID-19 [34] , from early March to early July 2020. Time=0 corresponds to March 7, 2020 and time=124 stands for the day July 08, 2020. As of July 10, 2020, the total death toll is 359, the total number of recovered cases is 11,525 and the total number of positive cases is 14,916 in Cameroon [34] .

The model (2) calibration has been performed using a newly developed optimization algorithm based on trust-region-reflective (TRR) algorithm, which can be regarded as an evolution of Levenberg-Marquardt algorithm [26] . This robust optimization procedure can be used effectively for solving nonlinear least-squares problems. This algorithm has been implemented using the lsqcurvefit function, which is available in the Optimization Toolbox in MATLAB. Necessary model parameters have been estimated using this optimization technique. Daily infected cases data have been collected from a trusted data repository, which is available online [34] . A 7-day moving average of the daily reported cases has been used for our model calibration due to moderate volatile nature of real data Fig. 2 . It has been observed that the number of daily testing in Cameroon has been really inconsistent. With an aim to capture the real outbreak scenario, the 7-day moving average has been used in this regard. According to recent statistics, the total population of Cameroon is found to be 26,548,238 [34] . For initial conditions, we take S(0) = 26533750, C(0) = E(0) = A(0) = H(0) = R(0) = 0 and Q(0) = 1. The solutions have been obtained after resolving the following optimization problem numerically.

where ψ = {β, σ 1 , σ 2 , σ 3 , r 1 , r 2 , r 3 , r 4 , δ a , δ q , δ h , η, p, c, , γ, q}, is displayed in Table 2 Parameter Probable range base value TRR output References β 0. Table 2 : Calibration of the model parameters using a newly developed optimization algorithm based on well-known Trust-region-reflective algorithm. Partial Rank Correlation Coefficient (PRCC) method has been carried out to quantify the most dominant mechanisms, which are significantly responsible for the transmission dynamics of COVID-19 disease. PRCC is a global sensitivity analysis method, which is used to quantify the relationship between model response function (outputs) and model parameters (sampled by Latin Hypercube Sampling method) in an outbreak setting [28] . PRCC values are ranged between -1 and 1. A negative PRCC value indicates a negative correlation between the model output and respective input parameter, whereas a positive PRCC value depicts a positive correlation between response function and corresponding model parameter. Again, similar result has been found, when the confined infectious individual class (Q) has been used as the model response function Fig. 7 . In this case, the PRCC indexes are estimated to be 0.885 and −0.69 accordingly. Importantly, the public health implication of this scenario is that, the transmission dynamics of COVID-19 can be controlled effectively in the community by curtailing the direct infectious contact rate and strengthening the confinement measures. Direct transmission rate can be limited by promoting several non-pharmaceutical measures such as maintaining physical distance carefully, wearing face-masks in public places and following health guidelines properly. However, it is often really challenging for the government of a developing country to implement strict confinement policies [27] .

Thanks to the memory effect which represents an advantage of the fractional derivative compared to the ordinary derivative, the theory and the application of fractional calculus have been widely used to model dynamic processes in the fields of science, engineering and many others [3, 4] . Before presenting the fractional model with Caputo derivative, we give the definition of fractional derivative in Caputo sense, and some useful results (see [20] ).

Here, we remind a definition and some results about Caputo (C) derivative.

Definition 1. The Caputo non-integer order derivative of g ∈ C k −1 is presented as

Lemma 4. If 0 < β < 1 and m is an integer (non-negative), then ∃ the +ve constants C β,1 and C β,2 only dependent on β, s.t

and (m + 2) β+1 − 2(m + 1) β+1 + m β+1 ≤ C β,2 (m + 1) β−1 . where C is a +ve constant independent of m & h.

By Replacing the integer derivative order of the ODE COVID-19 model (2) by the fractional-order model in Caputo sense, we obtain the following fractional model

In this portion of the study, we will provide the existence of the unique solution for the proposed fractional-order model under the Caputo fractional operator by applying fixedpoint results. In this setting, B is the Banach space of real-valued continuous functions defined on an interval I with the associated norm

, with E(I) denotes the Banach space of real-valued continuous functions on I and the associated sup norm. For convenience the proposed system (11) can be rewritten in the equivalent form given below.

By applying the Caputo fractional integral operator, the above system (12) , reduces to the following integral equation of Volterra type with the Caputo fractional integral of order 0 < α < 1.

Now we prove that the kernels G 1 ; G 2 ; G 3 ; G 4 ; G 5 ; G 6 and G 7 fulfil the Lipschitz condition and contraction under some assumptions. In the following theorem, we have proved for G 1 and can proceeds for the rest in a similar patron.

Theorem 2. Let us consider the following inequality 0 ≤ π I ηp h αk 4 + ξ < 1.

The kernel G 1 satisfies the Lipschitz condition as well as contraction if the above inequality is satisfied.

Proof. For S and S 1 we proceed as below:

Since A(t) and Q(t) are bounded functions, i.e, A(t) ≤ k 1 and Q(t) ≤ k 2 , by the property of norm function, above inequality can be written as

.

Here V 1 = {β(k 1 + ηk 2 ) + c} which implies that

Hence, for G 1 the Lipschitz condition is obtained and if an additionally 0 ≤ β(k 1 + ηk 2 ) + c < 1, we obtain a contraction. The Lipschitz condition can be easily verified for the rest of the cases and given as follows:

Recursively, the expressions in (13) can be written as

The difference between successive terms of system (12) in recursive form is given below:

With below initial conditions S 0 (t) = S(0), C 0 (t) = C(0), E 0 (t) = E(0), A 0 (t) = A(0), Q 0 (t) = Q(0), H 0 (t) = H(0), and R 0 (t) = R(0).

Taking norm of the first equation of (17), we obtain

Using Lipschitz condition (14) we obtained

Thus, we have

Similarly, for the rest of equations in system (12) we obtained

From above we can write that

Now, we claim the following result which guaranteed the uniqueness of solution of model (12) .

Theorem 3. The proposed fractional epidemic model (12) has a unique solution for t ∈ [0, T ] if the the following inequality holds

Proof. As we have shown that the kernels condition given in (14) holds. So, by considering the Eqs. (20) and (29) , and by applying the recursive technique we obtained the succeeding results as below:

Therefore, the above mentioned sequences exist and satisfy φ 1n (24) and employing the triangle inequality for any k , we obtain

where 1 Γ(α) b α V i < 1 by assertion and T i = 1 Γ(α) V i b α n , i = 1, 2, ..., 7.

Therefore, S n , C n , E n , A n , Q n , H n , R n are regarded as the Cauchy sequence in G(J). Hence, they are uniformly convergent as described in [20] . Applying the limit theory on Eq. (16) when n → ∞ indicates that the limit of these sequences is the unique solution of the model (12) . Ultimately, the existence of a unique solution for (12) has been achieved.

In view of the fact that most of the fractional differential equations (FDEs) do not have exact analytic solutions, so approximation and numerical techniques must be employed. Numerous analytical and numerical methods have been put up to solve the FDEs. For numerical solutions of the system (11) one can apply the generalized Adams-Bashforth-Moulton method. To provide the exact solution by means of this algorithm, we take into account the following nonlinear fractional differential equation [19] :

This equation is equivalent to Volterra integral equation:

Diethelm and Neville [7] employed the predictor-correctors scheme based on the Adams-Bashforth-Moulton algorithm [6] . So we use again this same technique to find the solution of the projected model (11) .

For α ∈ [0, 1], 0 ≤ t ≤ T and setting h = T /N and t n = nh, for n = 0, 1, 2, ..., N ∈ Z + , the solution of the projected model is

where

and

Theorem 4. The numerical method (28)-(29) is conditionally stable.

Proof. LetS 0 ,S j (j = 0, ..., n + 1) andS P n+1 (n = 0, ..., N − 1) be perturbations of S 0 , S j and S P n+1 , respectively. Then, the following perturbation equations are obtained by using Eqs. (20) and (29) 

Using the Lipschitz condition, we obtain

where ζ 0 = max 0≤n≤N {|S 0 | + h α M a n,0 Γ(α+2) |S 0 |}. Also, from Eq. (3.18) in [20] we derive

where γ 0 = max{ζ 0 + h α M a n+1,n+1 Γ(α+2) η 0 }. C α,2 is a positive constant only depends on α (Lemma 4) and h is assumed to be small enough. Applying Lemma 2 concludes |S n+1 | ≤ Cγ 0 . which completes the proof.

We replace the integer order model to the fractional order model to study the true Cameroonian data. The advantage to use Caputo fractional derivative to study the current true data is to predict the epidemic peaks more clearly at different fractional order values α. Parameter values use here are found in Table 2 , and for which the basic reproduction number, R 0 is equal to 3.41. For initials conditions, we take S(0) = 26548238, C(0) = E(0) = A(0) = H(0) = R(0) = 0, Q(0) = 1, and a fixed time step size of h = 2 −8 . The value of the fractional-order derivatives are α = 1, 0.9, 0.8, 0.7, 0.6, 0.5.

Firstly, figures have been analysed at integer order (denoted by red line) for the given model and then we analyse the nature of the compartmental model at different fractional values to give more peak predictions for the real-time data. It has been observed that onset of the epidemic peaks is delayed when the value of order α decreases.

In Figures 8 and 9 , one can see that when α −→ 1 the solutions of our fractional model (11) converge to the solutions of the integer model (2) . From above simulations, it can be conveyed that there are so much uncertainties of epidemic peaks for the given real-time data range. Caputo derivative works well to study these biological phenomenons. Robust forecasting results of COVID-19 outbreak can be achieved by applying our model structure as imperfect quarantine and unsound confinement which are harsh realities in developing countries have been taken into account meticulously.

A compartmental mathematical model has been formulated and studied to predict the evolution of the novel coronavirus disease in Cameroon. After presenting the model with integer derivatives, the basic reproductive number R 0 of the model has been computed, which measures the number of secondary COVID-19 positive cases caused by a single infectious person in a completely susceptible Cameroonian population. The dynamics of the disease is impacted by the value of this epidemic threshold. Indeed, if R 0 < 1, then a sufficiently small influx of infected individuals will not generate an outbreak of the COVID-19, and if R 0 > 1, then a sufficiently small influx of infected individuals will generate an outbreak of the COVID-19. Using real data communicated by the Cameroonian government, our model has been calibrated using a newly developed optimization algorithm based on wellknown trust-region-reflective (TRR) algorithm. With the result of this model calibration, the probable peak date is estimated to be on August 01, 2020 with approximately 1073 (95% CI : 714 − 1654) daily new confirmed cases, and which corresponds to approximately ∼ 20, 100 (95% CI : 17, 343 − 24, 584) cumulative infected cases. As of July 11, according to the recent statistics communicated by the health authorities of Cameroon, nearly 14,916 confirmed cases have already been identified. This difference between the "real data" communicated by the government and the predictions of our model can be explained by the fact that the fear that this disease creates within populations pushes the wealthy individuals towards private hospitals which not only give false test results but overbidding the treatment of this disease [8] .

Global sensitivity analysis has been performed implementing partial rank correlation coefficient (PRCC) analysis to quantify the most crucial parameters, which significantly impact the progression dynamics of COVID-19. It has been unearthed in our analysis that the direct infectious contact rate (β) and the efficacy of confinement measures ( ) are the most influential parameters in controlling the outbreak of COVID-19. This ignites the public health implication that by reducing physical contact between people through social distancing, wearing face-masks and confinement measures, the disease outbreak can be controlled successfully.

Eventually, the concept of Caputo derivatives has been deployed to formulate the fractional model and the existence and uniqueness of the solutions have also been presented. Using generalized Adams-Bashforth-Moulton method, the numerical scheme has been constructed for the fractional model. Numerical simulations enlighten that for α = 1 the solutions of our fractional model converge to the solutions of the integer model, and the decrease of the fractional-order parameter (0 < α < 1) leads to delay the onset of the epidemic peaks.

On a comprehensive model of the novel coronavirus (covid-19) under mittag-leffler derivative

Fractal-fractional differentiation and integration: connecting fractal calculus and fractional calculus to predict complex system

Fractional discretization: the african's tortoise walk

Approximate solution of tuberculosis disease population dynamics model

Chaos analysis and asymptotic stability of generalized caputo fractional differential equations

An algorithm for the numerical solution of differential equations of fractional order

Analysis of fractional differential equations

Covid-19: la clinique Marie O et le scandale hors de prix

To mask or not to mask: Modeling the potential for face mask use by the general public to curtail the covid-19 pandemic

Impact of non-pharmaceutical interventions (npis) to reduce covid19 mortality and healthcare demand

A new approach for solving multi variable orders differential equations with mittag-leffler kernel

A new approach for solving nonlinear volterra integro-differential equations with mittag-leffler kernel

A new study of unreported cases of 2019-ncov epidemic outbreaks

Modeling the transmission dynamics of the covid-19 pandemic in south africa

Hydroxychloroquine and azithromycin as a treatment of covid-19: results of an open-label non-randomized clinical trial

Mgr Kleda reçoit de nombreux soutiens pour ses recherches contre le coronavirus

Theory and applications of fractional differential equations

On comparison systems for ordinary differential equations

On the fractional adams method

The finite difference methods for fractional ordinary differential equations. Numerical Functional Analysis and Optimization

Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (sars-cov-2)

Coronavirus : Un premier cas confirmé au cameroun

Stigma of people with hiv/aids in sub-saharan africa: a literature review

Covid-19 in cameroon: a crucial equation to resolve. The Lancet Infectious Diseases

Forecasting covid-19 pandemic: A data-driven analysis

Has countrywide lockdown worked as a feasible measure in bending the covid-19 curve in developing countries? medRxiv

Sensitivity analysis of chronic hepatitis c virus infection with immune response and cell proliferation

Mathematical assessment of the impact of non-pharmaceutical interventions on curtailing the 2019 novel coronavirus

Numerical simulation of initial value problems with generalized caputo-type fractional derivatives

The sars, mers and novel coronavirus (covid-19) epidemics, the newest and biggest global health threats: what lessons have we learned?

Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications

Geometric and physical interpretation of fractional integration and fractional differentiation

Operational matrix for atangana-baleanu derivative based on genocchi polynomials for solving fdes

Sir epidemic model with mittag-leffler fractional derivative

The health stigma and discrimination framework: a global, crosscutting framework to inform research, intervention development, and policy on health-related stigmas

Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission

Ordinary Differential Equations

The authors thank the Editor and the anonymous reviewer for their comments which permitted us to improve the manuscript.

This work does not have any conflict of interest. No funding is received for this study.

 Table 2