Fractional model for the spread of COVID-19 subject to government intervention and public perception.

COVID-19 pandemic has impacted people all across the world. As a result, there has been a collective effort to monitor, predict, and control the spread of this disease. Among this effort is the development of mathematical models that could capture accurately the available data and simulate closely the futuristic scenarios. In this paper, a fractional-order memory-dependent model for simulating the spread of COVID-19 is proposed. In this model, the impact of governmental interventions and public perception are incorporated as part of the nonlinear time-varying transmission rate. In addition, an algorithm for approximating the optimal values of the fractional order and strength of governmental interventions is provided. This approach makes our model suitable for capturing the given data set and consequently reliable for future predictions. The model simulation is performed using the two-step generalized exponential time-differencing method and tested for data from Mainland China, Italy, Saudi Arabia and Brazil. The simulation results demonstrate that the fractional order model calibrates to the data better than its integer order counterpart. This observation is further endorsed by the calculated error metrics.

Introduction

The novel coronavirus pandemic, COVID-19, has impacted all aspects of our life and the disease continues to threaten global economy, public health, human interactions, and safety. Since the outbreak of the disease, there has been a collective effort to monitor, predict, and control its spread. For this purpose, different biological and non-biological techniques, tools, interventions, analysis, and experiments have been employed to understand the dynamics of disease spread.

Epidemiological mathematical models continue to play a pivotal role in understanding, capturing and predicting the spread of infectious diseases. In particular, first-order models such as SIR, SEIR, and many others, are intensively explored to mimic the collected data and make effective predictions. In the recent months, several such models have been developed to understand and predict the transmission of COVID-19 pandemic. See for example [1], [2].

Fractional-order models have been recognized as a powerful mathematical tool to study anomalous behaviors observed in many physical processes with prominent memory and hereditary properties. In general, such properties might not be satisfactorily captured by the integer-order counterpart models. This motivated researchers to consider fractional-order compartmental models for better understanding of dynamics of infectious diseases.

Several fractional-order epidemiological models have been proposed to study the outbreak and spread of different types of epidemics. They include models for Dengue fever [3], [4], [5], Influenza [6], Tuberculosis [7], HIV transmission [8], [9], and Malaria [10]. These models have been analyzed mathematically, validated epidemiologically, calibrated empirically and simulated numerically.

Since the outbreak of COVID-19 pandemic, a variety of fractional-order models for the spread of this disease are formulated [11], [12], [13], [14] These models differ mainly in the considered compartments, mitigation effects, multi-fractional orders, and parameter forms.

Fractional-order epidemiological models are often obtained from their integer-order counterparts by replacing the first-order derivative with a fractional one. However, this leads to a mismatch between the dimensions if the same parameters are used. This balance of dimensionality has been highlighted in some models but without much epidemiological interpretation of the new parameters. A common interpretation is that the parameters are observed in inhomogeneous time scale [15] corresponding to the fractional power of time unit. These parameters in the fractional model are commonly approximated by taking the fractional power of the classical parameters or by introducing an auxiliary parameter [16], [17], [18], [19], [20].

In this article, we propose and simulate a system of fractional order differential equations for the dynamics of COVID-19 pandemic. The proposed model uses the nonlinear time varying transmission rate given in [1], [21] which incorporates the effects of governmental interventions and public perception. In addition, an algorithm is developed so that for a given data set, the optimal fractional order of the model together with empirical values imitating the strength of governmental interventions are approximated. Consequently, reliable predictions for the future scenarios could be made.

The proposed model and algorithms are applied to the data from Mainland China, Saudi Arabia, Italy and Brazil. The simulation results suggest that the fractional order model has better agreement with the data than its integer order counterpart. This outcome is also supported by the error metrics considered in this paper. In addition, prediction trajectories for Saudi Arabia (a case with several changes in governmental interventions) and Brazil (a case with no governmental interventions) are included.

Since the focus is to provide adequate study and simulation of the fractional order model, we adopt the two-step generalized time differencing GETD2 scheme developed by Garrappa and Popolizio [22] to obtain the numerical results. However, to achieve a cost-effective implementation of this scheme, we use the partial fraction decomposition of the global Padé approximants of the two-parameter Mittag-Leffler function [23].

The rest of this paper is arranged as follows. Next, the memory-dependent SEIR model is introduced. Section 3 is devoted to the discussion of reproduction number of the disease and local stability analysis. The computational algorithm is described in Section 4. The simulation results are presented and discussed in Section 5.

Memory-dependent COVID-19 Model

Following [1], we assume the total population at any given time to be distributed among the four compartments: susceptible, exposed, infectious, and recovered, labeled by and respectively. That is,The assumed dynamics among these populations is depicted in Fig. 1 and can be describes as follows:

Fig. 1

Transfer diagram between the compartments.

Susceptible individuals enter the exposed compartment at the rate where is a nonlinear transmission rate.

Exposed individuals progress to the infectious compartment at the rate while the infectious individuals recover from the disease at a rate .

The population is assumed to have a natural death rate which is given by where denotes the life-expectancy.

The birthrate of the population is assumed to be equal to the death rate.

Furthermore, we consider a variable to represent the public perception of the risk of the pandemic. This quantity increases with the number of infections at a rate and diminishes over time at a rate . The parameter is the proportion of severe cases.

The values and can be interpreted as the mean latent period, mean infectious period, and mean duration of public reaction, respectively. These quantities together with and all have units of .

Based on the above assumptions, we propose the following model for the dynamics of COVID-19 pandemic:subject to the non-negative initial conditionsThe operator denotes the Caputo fractional derivative of order defined byThis derivative has been used successfully in the literature for compartmental models of different infectious diseases [3], [4], [6]. In addition, the Caputo derivative allows to specify the initial data in a natural way.

In the special case the model in (1) reduces to a modification of COVID-19 model introduced in [1].

As observed in the right-hand side of (1), the th power of the parameters is used as suggested by Diethelm [3] and Angstmann et al. [18]. This is necessary to ensure the consistency in units since the left-hand side has the units of (time).

For the transmission rate in (1), as in [1], we use the formulawhere is the baseline transmission rate, is the strength of governmental intervention. This quantity reflects the effect of government’s policies on the transmission of the disease. These policies include, closure of schools, stay home measures, mandatory use of face masks and physical distancing among others. On the other hand, is the intensity of individual reaction. It quantifies adherence of the population to public health measures such as the use of face masks, social distancing and proper hand-washing.

In the integer order SEIR models, the parameters and describe the transfer rates between the different compartments and are measured per unit time, i.e. . These parameters are the intuitive standard epidemiological rates such as, recovery, death, incubation, and transmission rates, which can be identified empirically or statistically. As for the fractional order model, these transfer rates retain similar epidemiological meanings. However, they are to be measured on time scale which might not be feasible. Raising these parameters to the power is an approach to approximate the corresponding transfer rates measured on time scale . An advantage of this approach is that the model reduces to the classical counterpart when .

Basic reproduction number and equilibria stability

The basic reproduction number, is the expected number of secondary infections that a single infectious individual will generate on average within a susceptible population. This number is arguably the most important quantity in determining the nature of an epidemic. In practice, a disease spreads whenever the reproduction number exceeds 1 (). On the other hand, an outbreak with a reproduction number below 1 will gradually disappear.

An equilibrium point is defined as follows.

[24] A point is an equilibrium point of the Caputo fractional dynamical system if and only if .

Using this definition, we see thatis an equilibrium point of (1). This point is called a disease-free equilibrium.

By applying the next-generation operator approach (see [25], [26]), the reproduction number is obtained as the spectral radius of the next generation matrix wherewith which is the transmission rate evaluated at . Thus, of (1) is given by

Next, we analyze the local asymptotic stability of the disease-free equilibrium.

For this purpose, we state the following result which follows by linearizing the systemabout the equilibrium point (as in Ahmed et al. [27]) and then applying Theorem 2 of Matignon [28].

Consider the following fractional-order system: If all eigenvalues of the Jacobian matrix satisfy then is locally asymptotically stable.

Based on Lemma 3.2, we have the following stability result for system (1).

The disease-free equilibria of the model (1) is locally asymptotically stable if .

The Jacobian of the right hand side of (1) at is given byThe eigenvalues of the matrix are given by

Clearly are negative. Moreover, if it is not hard to see thatand thus is also negative. Consequently, if then condition (3) is satisfied and Lemma 3.2 implies that is locally asymptotically stable. □

Computational algorithm

Here, we describe briefly the numerical algorithm used to simulate model (1). For this purpose, we consider the matrix form of the system,where

To solve this system numerically we use the GETD2 scheme developed in [22]. It has been shown that this scheme when applied to (4), has a rate of convergence. In addition, it has elegant stability properties that makes it robust for stiff systems without need to impose severe restrictions on the step size.where . To initialize the scheme we take and compute asFor efficient implementation of the GETD2 scheme, we use a technique based on the partial fraction decomposition of the global Padé approximant of type (7,2) developed in [23] for the Mittag-Leffler function.

Model simulation

In this section, we perform numerical simulations of the fractional order COVID-19 model (1) using the GETD2 scheme (5). We consider different case studies in which we calibrate the model to the data on confirmed COVID-19 infections in Mainland China, Italy, Saudi Arabia and Brazil. In addition, we present some predictions regarding possible spike and/or decay in the number of infections based on the intensity of behavioral response of population. For these purposes, we begin with a discussion about the parameters and their values.

The range of values of the parameters of the COVID-19 model (1) are listed in Table 1 . For our simulations, we use the following values.On the other hand, as outlined in Algorithm 1 , the optimal values of the fractional order and governmental intervention strength will be identified using the given data sets.

Table 1

Model parameters with possible range of values together with references.

Model Parameters	Range of values	References
Baseline Transmission rate, β_0	0.5 - 1.5	[29]
Governmental intervention strength, ν	[0,1)	[1], [21]
Intensity of behavioral response, κ	[0,∞)	[1], [21]
Mean latent period, σ^−1	2 - 14 days	[30]
Mean infectious period, γ^−1	5 days	[1]
Proportion of severe cases, ϕ	0.2	[1]
Mean duration of public reaction, ω^−1	11.2 days	[1], [21]

Algorithm 1

Approximating the optimal pair

Based on the timeline of governmental interventions, the numerical values assumed by the parameter may vary on different time intervals. For the implementation of Algorithm 1, feasible values of this parameter should reflect the intensity of the interventions. For example, periods with strict regulations should have no smaller values of than those with relaxed interventions. As such, at the initial phase of the pandemic and before any governmental intervention, we set to 0. On the other hand, based on the studies in [29], [31], substantial values of can be considered when accounting for the effect of mandatory use of face mask and physical distancing in public places.

The intensity was estimated in the study of Influenza pandemic by He et al. [21] as . However, by studying and matching different data sets of COVID-19 infections, we observed that to better reflect the human behavioral response, piecewise adjustments to this value are to be considered to reflect the interventions timeline. In particular, we use at the early stage of the outbreak since it is assumed that adequate awareness and precautions are not yet in place. For approximating and in Algorithm 1, we use the value for later time intervals. Nonetheless, when most of the governmental lockdown measures are being lifted, it is noted that the intensity of human behavioral response plays a major role in the decline of the number of infections. Thus, once the optimal values of and are approximated, we fine-tune in order to achieve better matching.

In all our case studies, we implement the GETD2 scheme with a step-size . Also, the root mean square error (RMSE) is used as the error metric in the implementation of Algorithm 1.

In all our simulations, the values of the RMSE between the approximated infected population and the actual data is provided for optimal and . Moreover, to give a perspective of the corresponding relative error, we also present the values of the mean absolute percentage error (MAPE).

Case study: outbreak in Mainland China

In Table 2 , we consider the reported data of cumulative COVID-19 infections in Mainland China from January 19, 2020 to March 6, 2020, as provided in [32]. To simulate the outbreak of the pandemic, we use thetogether with those given in (7). Furthermore, the initial populations are

Table 2

Cumulative data of number of infections in Mainland China between January 19, 2020 and March 6, 2020 as reported in [32].

January
19	20	21	22	23	24	25	26	27	28
198	291	440	571	830	1287	1975	2744	4515	5974
29	30	31
7711	9692	11791
February
1	2	3	4	5	6	7	8	9	10
14380	17205	20438	24324	28018	31161	34546	37198	40171	42638
11	12	13	14	15	16	17	18	19	20
44653	46472	48467	49970	51091	53139	55027	56776	57593	58482
21	22	23	24	25	26	27	28	29
58879	59527	59741	60249	60655	61088	61415	61842	62415
March
1	2	3	4	5	6
62617	62742	62861	63000	63143	63242

As per the lockdown records of Chinese Administrative Division, we consider the governmental intervention in two phases. The First phase corresponds to the period Jan 20 – Feb 02 and the second one to the period Feb 03 – March 06.

Applying Algorithm 1, we obtain as an approximation of the optimal fractional order. The corresponding values of governmental intervention strength are given in Table 3 . For comparison purposes, we also provide the optimal values of for the integer order problem, . Furthermore, the values of RMSE and MAPE in approximating the infected population are also appended.

Table 3

Midland China, the piecewise values of and for and together with the error metrics in approximating the infected population.

Period	ν for q=1	ν for q=0.93	κ
Jan 20 – Feb 02	0	0	0
Feb 03 – Mar 06	0.945	0.875	1117.3
q	RMSE	MAPE
0.93	281.1	51.735
1	399.616	37.589

The error metrics in Table 3 suggest that the fractional order model performs better in terms of the RMSE. However, the MAPE seems to support the integer order model. This blended outcome can be understood in view of Fig. 2 . As can be observed, the absolute errors for the fractional order model are smaller throughout most of the period and thus smaller RMSE. On the other hand, although the relative errors are smaller in general for the fractional model, they attained some large values towards the end of the period, which in turn impacted the MAPE.

Fig. 2

Mainland China, error in matching the daily reported infections. Absolute error (left), relative error (right).

Figs. 3 and 4 include the outcomes of our simulations for the infection population and the parameters and . Fig. 5 provides an insight on the effect of governmental interventions. As observed, the change in governmental intervention strength on Feb 3rd yields an instantaneous change in the rate of change of exposed population. This phenomenon is expected, because governmental interventions such as inter-city lockdown, curfew, school closure among others, are expected to mitigate the contacts between susceptible and infectious individuals. However, some increase in the number of infections is still observed. This accounts for the individuals who were incubating the virus before implementation of the governmental interventions.

Fig. 3

Mainland China, matching of the reported daily infections (left). The values of and the corresponding transmission rate (right).

Fig. 4

Mainland China, matching of cumulative number of infections.

Fig. 5

Mainland China, Infected and Exposed population. = 0.93 (left), = 1 (right).

Our numerical experiments also demonstrate the robustness of the fractional-order model as a generalization of the integer-order model. It is observed that the fractional order model gives extra degree of freedom, via the order of differentiation for data fitting.

Case study: outbreak in Saudi Arabia

Since the outbreak of COVID-19 in Saudi Arabia, the Kingdom has witnessed several phases of governmental interventions aimed at mitigating spread of the diseases. Due to seasonality, special breaks and festivals, there have been a number of relaxations and tightening of these interventions.

To simulate the outbreak of the Pandemic within the Kingdom, we only adopt those interventions that were implemented at a national scale, as provided in https://www.garda.com/crisis24/news-alerts/334921/saudi-arabia-curfew-hours-modified-for-ramadan-april-21-update-29. We consider the parameters in (7) together withThe initial conditions areBased on the data provided in https://ourworldindata.org/coronavirus-source-data, our simulation is performed for the period of March 03 to July 17 excluding the days 4th, 5th and 7th of March, where there was no reported data.

Using Algorithm 1, we obtained an approximation of the optimal pair, with . For comparison purposes, we present in Table 4 the values corresponding to both and . In addition, the values of the RMSE and MAPE in approximating the infected population for these values of are presented. Brief remarks on the governmental actions and the values of intensity are also provided.

Table 4

Saudi Arabia, timelines of nationwide governmental interventions and the corresponding approximated values of and together with the RMSE and MAPE in approximating the infected population.

Date	Remark	ν for q=1	ν for q=0.87	κ
03-Mar	Discovery of the first cases in KSA	0	0	0
16-Mar	Closure of schools	0.1	0.1	0
23-Mar	Implementation of nationwide 7pm–6am curfew	0.4	0.25	111.73
06-Apr	Nationwide 3pm - 6am curfew and suspension of intercity travels	0.6	0.425	558.65
26-Apr	Relaxation of Nation curfew (5pm–9am) and intercity travels during Ramadan	0.5	0.4	558.65
23-May	Implementation of nationwide 24 h curfew	0.8478	0.8478	1117.3
28-May	Curfew 3pm–6am	0.45	0.5	1117.3
31-May	Relaxation of curfew (8pm–6am) and inter-region travels	0.4	0.3	4469.2
21-Jun	Return to normalcy with mandatory face masking and physical distancing measures	0.4	0.3	14524.9
q	RMSE	MAPE
0.87	425.917	7.907
1	603.735	10.632

We observe that the values of presented in Table 4 reflect the intensity of the control measures in the sense that periods with more restrictive interventions have higher values of . From an epidemiology perspective, these interventions such as imposing a curfew, restricting the gatherings size, among others, reduce interactions and consequently mitigate the spread of the disease. As a result, the fitted values of tend to reflect the epidemiological reality since bigger values of will imply smaller values of the transmission rate in (2).

Fig. 6 includes plots of the daily infections obtained from the fractional model (1) and plots of the parameters and . It is observed that the polynomial behavior of the data at the early stage of the outbreak is better captured by the power-law corresponding to than by the exponential or almost exponential law, or . However, the final stage of simulation is better captured with . This Further illustrates the effectiveness of fractional order models for making predictions. The outcome of our simulation for the cumulative infections is presented in Fig. 7 . We note that yields a satisfactory matching.

Fig. 6

Saudi Arabia, matching of the reported daily infections (left). The values of and the corresponding transmission rate (right).

Fig. 7

Saudi Arabia, matching of cumulative number of infections.

Using the calibrated model (1), predictions for the spread of the disease subject to different values of are shown in Fig. 8 . We assume that the governmental intervention continues to be the same as recent period. It is shown that, a drop in the intensity below the value between 21-June and 17-July will yield a spike in the infection numbers. On the other hand, maintaining same intensity will help obtain a slow decline of the infection numbers. However, a notable decay in the number of infections could be achieved if there is a significant increase in the intensity of the behavioral response. This corresponds to improved compliance with proper hand-washing, face masking, physical distancing, among others.

Fig. 8

Saudi Arabia, matching of reported daily infections up to July 17 with and predictions subject to different human behavioral response for the period July 18–August 7.

Case study: outbreak in Italy

In the early outbreak of COVID-19 in Italy, the country was hit badly and a number of governmental interventions were initiated. To simulate this outbreak using (1), we considertogether with those provided in (7). The initial conditions areThe timeline of governmental interventions considered for this simulation are based on the report by Briscese et al. [33], up to 3-April, together with other reopening information available on https://www.bbc.com/news/world-europe-52435273. Our simulation covers the period of 31-January 2020 to 17-July 2020, based on the data provided in https://ourworldindata.org/coronavirus-source-data. However for the period February 01-21, 2020, no new infections were reported. Consequently, we excluded these dates from our simulation.

Algorithm 1 applied to the fractional order model (1) yields an approximation for the optimal and pair, with . The corresponding values of are provided in Table 5 . In addition, brief descriptions of the governmental actions together with the adopted values of are also provided. Similar to Saudi Arabia case, we observe that the fitted values of reflect the intensity of governmental actions and are consistent from an epidemiological point of view.

Table 5

Italy, timelines of nationwide governmental interventions and the corresponding approximated values of and together with the RMSE and MAPE in approximating the infected population.

Date	Remarks	ν for q=1	ν for q=0.88	κ
22-Feb		0	0	0
4-Mar	Closure of Schools	0.3	0.15	0
9-Mar	Lockdown Measures	0.4	0.2	111.73
11-Mar	Closure ofmost retail stores	0.6	0.375	558.65
22-Mar	Suspension of all non-essential economic activities	0.8	0.75	558.65
14-Apr	Re-opening of some retail activities	0.75	0.75	11173
4-May	Re-opening of salon, museums, bars among others	0.75	0.75	15642.2
1-Jun	Re-opening of more economic and sporting activities	0.75	0.75	20111.4
q	RMSE	MAPE
0.88	425.917	7.907
1	547.329	8.740

Plots of the daily and cumulative infections obtained from the model (1) are presented in Figs. 9 and 10 , respectively. We observe from these plots that the fractional order model has a better matching. This observation is also supported by the values of the error metrics included in Table 5.

Fig. 9

Italy, matching of the reported daily infections (left). The values of and the corresponding transmission rate (right).

Fig. 10

Italy, matching of cumulative number of infections.

Case study: outbreak in Brazil

Information from various sources indicate that there were no governmental control measures to mitigate the COVID-19 outbreak in Brazil. Some of the events in the early days of the outbreak can be found in an article by Reuters https://www.reuters.com/article/us-health-coronavirus-brazil-timeline/timeline-key-moments-in-bolsonaros-handling-of-covid-19-crisis-idINKBN2482UL. Thus, for simulating the outbreak we set . The simulation covers the period from 26th February to 17th July using the data provided in https://ourworldindata.org/coronavirus-source-data. We exclude the days February 27–29 and March 02–04, 08 and 10, since no infections were reported. The other parameters used in simulating (1) for the outbreak in Brazil aretogether with those provided in (7). The initial conditions areWith for the whole simulation time, implementation of Algorithm 1reduces to approximating only the optimal fractional order which is approximated as .

Our examination of the infection data in Fig. 11 shows that there is a significant change in the reported numbers around day sixty. Under no governmental interventions, the parameter in our model that can be used to account for this dramatic change is the human behavioral response . For this purpose, we split our time interval into [0, 60] and [61, 134]. Then, we fit appropriate values of for both (the optimal value of ) and . The calculated values are shown in Table 6 together with the RMSE and MAPE in approximating the daily number of infections.

Fig. 11

Reported daily infections in Brazil.

Table 6

Brazil, the piecewise values of for and together with the RMSE and MAPE in approximating the infected population.

Period	κ for q=1	κ for q=0.70
Feb 26 – May 04	80445.6	26815.2
May 05 – July 17	27932.5	17876.8
q	RMSE	MAPE
0.70	6732.751	48.414
1	11783.101	145.99

Plots of the daily infections obtained from (1) are presented in Fig. 12 . As can be observed, the integer order model could not provide a reasonable matching of the data. On the other hand, the fractional order model with is able to capture well the initial spread of the disease. Plots of the cumulative infections are also presented in 13 .

Fig. 12

Brazil, matching of the reported daily infections in (left). The values of and the corresponding transmission rate (right).

Fig. 13

Brazil, Matching of cumulative number of infections.

Overall, we remark that modeling the COVID-19 outbreak in Brazil is not a trivial task due to the erratic nature of the reported data. However, using the extra degree of freedom offered by the fractional order model through the fractional order we are able to find a possible fit for daily infection data.

Using the calibrated model (1), predictions for the spread of the disease subject to different values of are shown in Fig. 14 . For the period from February 26 to July 17, we use the values given in Table 6 for . Then we predict possible infection numbers for the period from July 18 to August 16. The plots suggest that for a significant reduction in infections, a substantial increase in is required. This corresponds to improved compliance with proper hand-washing, face masking, physical distancing, among other public health measures.

Fig. 14

Brazil, matching of reported daily infections up to 17-July with and predictions subject different human behavioral response from 18-July to 16-August.

Discussion

The proposed model aims for improved fitting of the real data and thus more realistic predictions. It provides extra degrees of freedom through the fractional order and the transmission rate that includes the government intervention and intensity of human behavioral response.

From the outcomes of our numerical simulations, we note the following observations which can help in policy making for public health interventions.

The data matching reveals that strict governmental lockdown measures result in large values of the parameter and thus rapid decrease in the transmission rate. This is notably observed in the case of China.

The intensity of human behavioral response, such as adherence to masking, hand washing, physical distancing, among others, compliments the governmental interventions in mitigating the transmission of the disease. This is demonstrated by the predictions in Fig. 8, where we see that for fixed governmental interventions, the number of infections is inversely proportional to the behavioral response .

Complete absence of governmental interventions results in growing transmission of the disease even with strong human behavioral response. This scenario is demonstrated in the Brazil model, where we set to zero since there is no record of governmental interventions. As can be observed, the human behavioral response alone is not adequate to mitigate the spread of the disease. Rather, it complements the governmental interventions.

As a consequence of the above observations, gradual lifting of governmental measures is advisable. For example, interventions such as mandatory physical distancing and masking in public places should be in place. These less strict interventions can still be well complemented by the behavioral responses and consequently help flatten the curve. The prediction for Saudi Arabia in Fig. 8 is a typical reference.

Conclusion

In this paper, we proposed a nonlinear fractional order model to study the dynamics of COVID-19 pandemic. Using the data of the outbreak in Mainland China, Saudi Arabia, Italy, and Brazil, the effectiveness of the proposed model in simulating the spread of the disease is demonstrated. Predictions based on varying intensity of human behavioral responses are provided. By examining the mean square error of the predicted number of infections and the reported data, we found that the proposed fractional model yields better matching than the integer order model. From the perspective of relative error, the values of mean absolute percentage error confirm this conclusion for the simulation of outbreak in Saudi Arabia, Italy and Brazil. On the other hand, for China, these two error metrics have mixed outcomes.