Analysis of transmission dynamics of COVID-19 via closed-form solutions of a susceptible-infectious-quarantined-diseased model with a quarantine-adjusted incidence function.

We analyze the disease control and prevention strategies in a susceptible-infectious-quarantined-diseased (SIQD) model with a quarantine-adjusted incidence function. We have established the closed-form solutions for all the variables of SIQD model with a quarantine-adjusted incidence function provided β ≠ γ + α by utilizing the classical techniques of solving ordinary differential equations (ODEs). The epidemic peak and time required to attain this peak are provided in closed form. We have provided closed-form expressions for force of infection and rate at which susceptible becomes infected. The management of epidemic perceptive using control and prevention strategies is explained as well. The epidemic starts when ρ 0 > 1, the peak of epidemic appears when number of infected attains peak value when ρ 0 = 1 , and the disease dies out ρ 0 < 1. We have provided the comparison of estimated and actual epidemic peak of COVID-19 in Pakistan. The forecast of epidemic peak for the United states, Brazil, India, and the Syrian Arab Republic is given as well.

INTRODUCTION

The novel coronavirus (COVID‐19) emerged in the early December 2019 in the city of Wuhan located in Hubei Province of China, and within a month, it spread outside China. On January 31, 2020, the World Health Organization (WHO) declared COVID‐19 as the global emergency. The containment strategies of COVID‐19 adopted by several countries include hand hygiene, social distancing, isolation, quarantine, and lockdown. It is crucial for the mathematicians, statisticians, and scientists to develop appropriate models to uncover the essential aspects of COVID‐19 spread and its transmission dynamics. The information about the total number of confirmed, infected, recovered, and mortality cases is essential for the public health decision makers in pandemic preparedness.

Many mathematical models have been developed. The mathematical models in epidemiology are governed by a system of first‐order ODEs. The classical susceptible‐infected‐recovered (SIR) model developed by Kermack‐McKendrick is widely used by researchers in several studies, and many extensions of model are proposed. The susceptible‐infectious‐recovered‐deceased (SIRD) model is utilized for forecasting of COVID‐19 by Mansoor et al for Pakistan, Chatterjee et al for India, and Al‐Raeei for the United States, Russia, China, and the Syrian Arab Republic. Several models are utilized to estimate and predict the COVID‐19 which include SIR model, susceptible‐exposed‐infectious‐removed (SEIR), , susceptible‐asymptomatic‐infected‐removed (SAIR), and extensions of these models.

The dynamic properties of COVID‐19 models are analyzed by utilizing different methods. Yadav et al utilized machine learning methods to analyze COVID‐19 epidemic. Ali et al applied optimal control theory to help public health decision makers to implement different control programs for the containment of this deadly virus. A small number of papers have utilized classical techniques to establish the closed‐form solutions of epidemiology models. , The stability analysis of the dynamical systems, numerical methods, and data drive approach are also utilized by several researchers to investigate transmission dynamics of COVID 19.

Hethcote et al studied the effects of quarantine in six endemic models for infectious diseases. The susceptible‐infectious‐quarantined‐susceptible (SIQS) and susceptible‐infectious‐quarantined‐removed (SIQR) models were analyzed by taking simple mass action incidence, the standard incidence, and the quarantine‐adjusted incidence. In Chinviriyasit S and Chinviriyasit W, global stability of an SIQ epidemic model with simple mass action incidence, constant immigration, and natural death is investigated. The SIQ model is extended to susceptible‐infectious‐quarantined‐diseased (SIQD) model by adding a death compartment, and a quarantine‐adjusted incidence function is considered. To the best of our knowledge, SIQD model with a quarantine‐adjusted incidence function has not been considered in literature to analyze the transmission dynamics of COVID‐19. Our newly developed SIQD model with a quarantine‐adjusted incidence function is analyzed by utilizing the closed‐form solutions. It is worthy to mention here that we focus on the infectious group rather than the exposed or asymptotic to understand the epidemic peak, force of infection (FOI), and rate at which susceptible becomes infected by utilizing the closed‐form solutions.

The detailed outline of paper is as follows. The SIQD model with quarantine‐adjusted incidence function is formulated in Section 2. The basic reproduction number is derived by utilizing the next generation operator approach by van den Driessche and Watmough. The classical techniques of solving ODEs are utilized to establish closed‐form solutions for all the variables of SIQD model with a quarantine‐adjusted incidence function. In Section 3, we have provided the closed‐form expressions for epidemic peak and maximum number of infected individuals. The closed‐form expressions for the FOI and rate at which susceptibles become infected are given in closed form. In Section 4, we have provided different control and prevention policies to slow down pace of epidemic. In Section 5, the model is utilized to study the COVID‐19 epidemic in Pakistan. Moreover, we have forecasted epidemic peak for USA, Brazil, India, and the Syrian Arab Republic. The concluding remarks are presented in Section 6.

THE SIQD MODEL WITH QUARANTINE‐ADJUSTED INCIDENCE FUNCTION

We consider a model of COVID‐19 transmission in a constant population of size N which consists of four compartments susceptible S(t), infected I(t), quarantined Q(t), deceased Δ(t), and total population . In literature, several types of incidence functions are considered to study the dynamics of infectious diseases; see, for example, Bailey, Jacquez et al, and Watson. Hethcote et al studied the SIQS and SIQR models by taking simple mass action incidence, the standard incidence, and the quarantine‐adjusted incidence. We follow the idea of Hethcote et al and consider a quarantine‐adjusted incidence function. The actively mixing population which is the pool of people with whom a susceptible can meet for our SIQD model is . The expression for the quarantine‐adjusted incidence function is , and it is obtained by replacing the denominator N in the standard incidence by actively mixed population . The FOI is βI/S + I. The model can be expressed as a system of following four ODEs: where dot above a variable denotes time differentiation and initial conditions are , and . The parameters are defined as follows: β > 0 is contact rate of susceptible individuals with infected individuals, γ > 0 is isolation rate of infected individuals, α > 0 is the virus‐induced mortality rate of infected class, and μ > 0 is the virus‐induced mortality rate of quarantined class.

The basic reproduction number for the SIQD model with quarantine‐adjusted incidence function

Next, we compute the basic reproduction number ρ 0 for this model by utilizing the next generation operator approach by van den Driessche and Watmough. It is an essential dimensionless quantity which provides the threshold in the analysis of a disease to predict outbreak and to evaluate its control strategies. Let represent the disease free equilibrium (DFE). The associated F and V matrices of SIQD model with quarantine‐adjusted incidence functions (1)–(4) for the next generation operator approach are given as and

The the basic reproduction number ρ 0 is the only eigenvalue of FV −1 at the DFE and is given by

It is important to mention here that the epidemic starts if . Setting in Equation (2), we have provided βS 0/S 0 + I 0 − γ − α > 0. This condition in terms of the basic reproduction number is ρ 0 > 1 + I 0/S 0. The epidemic starts when ρ 0 > 1 at initial time t 0, the peak of epidemic appears when I(t) attains peak value I at time t when , and the disease dies out ρ 0 < 1 at time t .

Closed‐form solutions for SIQD model with quarantine‐adjusted incidence function

In this section, we establish the closed‐form solutions of system of Equations (1)–(4). It is important to mention here that Equations (1) and (2) form a coupled system of nonlinear ODEs in terms of variables S and I. Equations (3) and (4) are linear ODEs for variables Q and Δ, respectively. We can solve Equations (3) and(4) if variable I is known. First, we solve Equations (1) and (2) for variables S and I. Then we use value of I in ODEs (3) and (4) to find S(t) and Δ(t).

Introduce a new variable then

With the aid of Equations (1) and (2), Equation (10) takes following form: Equations (1) and (2) can be expressed as a system of linear ODEs in terms of variable Z(t) where variable Z satisfies a linear ODE given by (11). It is straightforward to derive Z(t) from the linear ODE (11) and is given by where and . We substitute the expression for the variable Z(t) from (14) in ODE (12), and resulting the ODE can be expressed in variable separable as follows:

After some simplifications and using initial condition , Equation (15) yields

Substituting expression for Z(t) from (14) and S(t) from (16) in Equation (9), the closed‐form solution for I(t) is given by

It is worthy to mention here that, alternatively, we can establish closed‐form solution for the variable I(t) given in (17) by solving ODE (13) subject to initial condition . We substitute the expression for the variable I(t) from (17) in ODE (3), and we arrive at a linear ODE which finally yields and finally, Δ is given by

Alternatively, same expression for Δ(t) can be derived by solving ODE (4). The closed‐form solutions for all variables of model can be summarized as follows: and . *

The equilibrium point for the model

The model attains its equilibrium point when ρ 0 < 1, and thus, β < γ + α. We take limit t → ∞ for closed‐form solutions of all variables in Equation (20) and arrive at the following equilibrium point for the model provided β < γ + α:

The endemic will end as I(t) → 0 with S(t) approaching some positive value if t → ∞.

THE CLOSED‐FORM EXPRESSIONS OF SOME IMPORTANT INDICATORS OF EPIDEMIC

In this section, we have provided the closed‐form expressions for the epidemic peak and maximum number of infected individuals. The closed‐form expressions for the FOI and rate at which susceptibles become infected are also given in analytical form.

The epidemic peak and maximum number of infected individuals

It is worthy to mention here that we can find analytical expression for time t where peak of infected curve occurs and it is termed as the epidemic peak. Differentiate I(t) given in (20) with respect to t; we get

Setting , we get and t  > 0 as β/γ + α > 1 and S 0 > I 0. We check second‐order derivative at this peak value of t, and it is given by

which is negative as β/γ + α > 1. The maximum number of infected cases I is reported for this value of t and is given by

Equations (23) and (25) can expressed in terms of the basic reproduction number as follows: and

The FOI and rate at which susceptibles become infected

The FOI is the rate at which the susceptibles are infected. We utilize the closed‐form solutions for S(t) and I(t) from (20) to find following closed‐form expressions for the FOI:

The expression for FOI given in (28) can also be expressed in terms of the basic reproduction number ρ 0 and is given by

The rate at which susceptibles become infected is βSI/S + I, which is Ω × S(t), and is given as

This can be further simplified as follows:

The rate at which susceptibles † become infected given (31) can be expressed as

THE SIQD WITH A QUARANTINE‐ADJUSTED INCIDENCE FUNCTION FRAMEWORK TO MANAGE EPIDEMIC

In this section, we use the specific values of parameters and closed‐form solutions of all variables provided in (20) to explore the results derived in Sections 2 and 3. The important indicators of understanding the transmission dynamic of disease which are peak time of epidemic (23), epidemic peak (25), FOI (28), and rate at which susceptibles become infected (31) are studied with the aid of simulations. We also provide an effective framework to manage the epidemic.

It is important to provide the range of parameters of SIQD model with a quarantine‐adjusted incidence function in the context of the existing literature. The literature on the estimates of the basic reproduction number for the COVID‐19 suggests that it varies from 1.4 to 5.7. According to WHO, it can take values between 1.4 and 2.5. Zhao et al provided the preliminary estimates of the basic reproduction number of novel coronavirus between 3.6 and 4.0 and between 2.24 and 3.58. In China, ρ 0 was 5.6015 before sealing Wuhan city, and it was reduced to 3.4094 as of February 25, 2020. Musa et al estimated ρ 0 as 2.37 in Africa.

The estimated range of the parameters β, γ, and α for the SIQD model with a quarantine‐adjusted incidence function is presented in Table 1. The reader is refereed to see previous studies , , , , , , and references therein for further details. ‡

TABLE 1

The range of parameters for the SIQD model with a quarantine‐adjusted incidence function taken from clinical studies and modeling

Parameter	Range for parameters (day−1)	References
β	0.4–1.5	Previous studies 23 , 24 , 25
γ	0.2–0.8	Previous studies 10 , 25 , 28
α	0.001–0.1	Previous studies 25 , 26 , 27
μ	0.001	Assumed

The following comparative static analysis of the basic reproduction number, reveals that it has a direct relationship with infection rate β and inverse relationship with rate of isolation γ. The epidemic control strategies suggest to lower the value of the reproduction number and reduce it to less than one, that is, ρ 0 < 1. This can be achieved in different ways. When government enforces strict control measures in the form of lockdown (social distancing), then infection rate β declines. The height of the epidemic peak I reduces, the time span t to attain this height increases, and ρ 0 decreases. This slows down the transmission dynamics of the disease. And when community follows proper prevention measures (quarantine, self isolation, mask, disinfection of surfaces, and hands hygiene) then γ increases, and thus, it lowers the basic reproduction number. We can reduce the basic reproduction number either by reducing β or by increasing γ. One can also use appropriate combination of reducing β and increasing γ to lower the value of ρ 0.

The comparison of control strategies by government (lockdown) and prevention strategies (quarantine) is presented using simulations.

Simulations to study the effect of change of parameters on the epidemic peak and transmission of disease

We consider the following parameters as an example to explore different control strategies mentioned as above: N is 10 million, , , and initial conditions are The graphical analysis of all variables S(t), I(t), Q(t), and Δ(t) using closed‐from solutions (20) is given in Figure 1. §

FIGURE 1

The graphs of S(t), I(t), Q(t), and Δ(t) for SIQD model with a quarantine‐adjusted incidence function [Colour figure can be viewed at wileyonlinelibrary.com]

Figure 2A represents the effect of change of β while keeping other parameters as fixed. The effect of change of γ while keeping other parameters as fixed is given in Figure 2B. A closer look at Figure 2A,B shows that when β reduces (lockdown) or γ increases (quarantine and prevention), the epidemic peak reduces, and time span to attain this peak increases. This is an effective strategy to control spread of disease.

FIGURE 2

Effect of change of parameters on epidemic peak I and time span t [Colour figure can be viewed at wileyonlinelibrary.com]

One can also use appropriate combination of reducing β and increasing γ to lower the value of ρ 0. We have set initially and , and this is shown in Figure 3 in green solid line. The value of ; the epidemic peak is on Day 55 with 4 178 373 number of infected individuals. Then we changed parameters by a fixed proportion according to following three scenarios:

FIGURE 3

The effect of change of β and γ in a fixed proportion on epidemic peak I and time span t [Colour figure can be viewed at wileyonlinelibrary.com]

Increase β to ϕ β and reduce γ to γ/ϕ. This is shown in Figure 3 in blue dashed line. The value of ; the epidemic peak is on Day 32 with 6 171 341 number of infected individuals.

Reduce β to β/ϕ and increase γ to ϕ γ. This is shown in Figure 3 in red dotted line. The value of ; the epidemic peak is on Day 144 with 1 642 341 number of infected individuals. Scenario 2 is the best scenario to reduce epidemic peak and increase time span to attain the epidemic peak.

Increase γ to ϕ γ with no change in β. This is shown in Figure 3 in cyan long dash line. The value of ; the epidemic peak is on Day 65 with 3 054 914 number of infected individuals.

Scenario 2 or 3 is the best scenario to reduce epidemic peak and increase time span to attain the epidemic peak. We can conclude that increasing γ (quarantine and prevention) strategy is an effective way to slow down transmission dynamics than the strategy of reducing β (lockdown). The prolonged lockdown results in social and economic loss. One should use appropriate combination of reducing β (lockdown) and increasing γ (quarantine and prevention) to lower the value of ρ 0.

Another important strategy is to set the value of basic reproduction number at some reasonable fixed level close to the threshold value where transmission dynamics is slow. And then change parameters β and γ to reduce epidemic peak and increase the time span to attain this epidemic peak. This scenario is presented in Figure 4. This is a good strategy if the health care facilities are adequate to treat only I ∗ infected for time period t ∗ and is an effective way to develop herd immunity. To prevent second epidemic wave, it is necessary to slow down transmission dynamics instead of trying to quickly flatten the curve by strict lockdown measures. One should use appropriate combination of reducing β and increasing γ to maintain ρ 0 which will eventually help to attain the DFE.

FIGURE 4

Effect of change of β and γ for fixed value of [Colour figure can be viewed at wileyonlinelibrary.com]

Effect of change of parameters on the FOI and rate of infection

Next, we analyze the effect of change of parameters on the FOI and rate of infection. Figure 5A,B represents the effect of change of β while keeping other parameters as fixed. The FOI and rate of infection slow down as β reduces. This means that lockdown works efficiently for a short period in reducing force and rate of infection but at the cost of social and economics loss. But this is not an effective strategy in long run.

FIGURE 5

Effect of change of parameters on force of infection and rate of infection [Colour figure can be viewed at wileyonlinelibrary.com]

The effect of change of γ while keeping other parameters as fixed on the FOI and rate of infection is given in Figure 5C,D. A closer look at Figure 5C shows that when γ increases (quarantine and prevention), the time span to arrive at a certain point increases, but intensity of force remains same. The rate at which susceptible moves to infected class is reduced with an increase in γ.

The control and prevention policies include lockdown by government, quarantine, self isolation, social distancing, mask, disinfection of surfaces, and hands hygiene and are helpful in containment of corona virus. The effective control and prevention policies help in slow down pace of epidemic. The slow transmission dynamics will be helpful for a country to improve health care facilities by increasing capacity of beds, intensive care units, ventilators, and hiring more health care staff to deal with large number of patients.

Lockdown policies implemented by government in most of countries have negative impacts in terms of social and economic activities. The saving lives strategy has resulted in social and economic losses: decline in earnings, depreciation in gross domestic product (GDP), high unemployment rates, stagnant revenue, increased costs, and loss of social capital. We conclude that quarantine and prevention strategy measures are more efficient to slow down transmission dynamics than strict lockdown. The prolonged lockdown results in social and economic loss. One should use appropriate combinations of lockdown and prevention measures to flatten the curve of epidemic. The partial lockdown with effective quarantine and prevention measures is best way to handle the epidemic.

THE SIMULATIONS OF SIQD WITH A QUARANTINE‐ADJUSTED INCIDENCE FUNCTION BY UTILIZING LEAST SQUARE METHOD

In this section, we illustrate the simulation of the SIQD model with a quarantine‐adjusted incidence function by utilizing the closed‐form solutions of all variables provided in (20) for the new coronavirus disease. We had concluded in Section 4 that reducing infection rate β (lockdown) and increasing isolation rate γ (quarantine and prevention) in a fixed proportion are best ways to handle the epidemic. This can be done by reducing β to β/ϕ and increasing γ to ϕ γ for some fixed proportion ϕ. First, we analyze the estimated and actual epidemic peak for Pakistan. Then we apply the model for the United States, Brazil, India, and the Syrian Arab Republic. First, we fit the collected data of the new coronavirus disease to find the infection rate β, isolation rate γ, virus‐induced mortality rate of infected individuals α, and scaling factor ϕ for the model. We assume that the virus‐induced mortality rate of quarantined class μ is 10−7 days −1 for the pandemic. We apply the least square method for the fitting of the collected cases with the the aid of closed‐form solutions of all variables provided in (20). It is worthy to mention here that we have used the least square method with a general function form instead of a linear function. In Figure 6, we illustrate the flowchart of the fitting.

FIGURE 6

Flowchart of the fitting

The SIQD with a quarantine‐adjusted incidence function to analyze the estimated and actual epidemic peak for Pakistan

We study the epidemic peak of Pakistan which has already arrived on July 1, 2020. We use the collected cases of the new coronavirus disease up to date August 2, 2020, to fit the parameters of the new coronavirus disease. , The start date of pandemic in Pakistan is February 26, 2020, and we take this as t 0. The initial values of the susceptible, infected, quarantine, and death cases are , , . The values of parameters estimated by utilizing the closed‐form solution and least square method on actual data of Pakistan are given in Table 2.

TABLE 2

The parameter values for the SIQD model with quarantine‐adjusted incidence function to estimate COVID‐19 transmission dynamics for Pakistan

β	γ	α	μ	ϕ	N
(day −1)	(day −1)	(day −1)	(day −1)
71.0248	0.70909	0.0213	10−7	100	220 892 331

The basic reproduction number for Pakistan for set of parameters given in Table 2 is . A closer look at Figure 7 reveals that the infection cases increase up to the Day 126, and total number of infected cases are 108 193. Now, we compare these estimated values with actual data of coronavirus cases in Pakistan. The actual peak occurred on Day 127, and number of infected at this peak value are 108 642. We conclude that our model provided close estimates to the actual data.

FIGURE 7

The graph of active infected cases I(t) for Pakistan [Colour figure can be viewed at wileyonlinelibrary.com]

The SIQD with a quarantine‐adjusted incidence function to forecast epidemic peak for USA, Brazil, India, and the Syrian Arab Republic

First, we fit the collected cases of the pandemic for the United States where we use the collected cases of the new coronavirus disease up to date August 2, 2020, to fit the parameters of the new coronavirus disease. , The first case of the pandemic in the United States appeared on January 13, 2020, and from this day, we use the initial values where the initial value of the infection cases is and the initial values of the quarantine and death . The parameters of model for the United States are given in Table 3. A similar procedure is carried out for the collected cases of the pandemic for Brazil, India, and the Syrian Arab republic, , and collected cases of the new coronavirus disease up to date August 2, 2020 are considered to fit the parameters of the new coronavirus disease. The initial conditions, parameter values, and scaling factor for the SIQD model with a quarantine‐adjusted incidence function to forecast COVID‐19 transmission dynamics for the United States, Brazil, India, and the Syrian Arab Republic are summarized in Table 3. The graphs for the active cases of infected individuals for United States, Brazil, India, and the Syrian Arab republic are presented in Figure 8. The peak time t (days), number of infected at peak I , and basic reproduction number are provided in Table 4.

TABLE 3

The parameter values for the SIQD model with a quarantine‐adjusted incidence function to forecast COVID‐19 transmission dynamics for the United States, Brazil, India, and the Syrian Arab Republic

	USA	Brazil	India	The Syrian Arab Republic
β (day −1)	2.129	6.9	37.2711	17.6764
γ (day −1)	0.20469	0.68090	0.3719029	0.176009
α (day −1)	0.0203	0.0230	0.0189	0.02
μ (day −1)	10−7	10−7	10−7	10−7
ϕ	10	10	100	100
N	33 100 2647	212 559 409	1 380 004 385	17 500 657
t 0	Jan 13, 2020	Feb 26, 2020	Jan 30, 2020	March 22, 2020
I 0	1	1	1	1
Q 0	0	0	0	0
Δ0	0	0	0	0

FIGURE 8

The graphs of I(t) for SIQD model with a quarantine‐adjusted incidence function for USA, Brazil, India, and the syrian Arab republic [Colour figure can be viewed at wileyonlinelibrary.com]

TABLE 4

The values of t , I , and ρ 0 for the SIQD model with a quarantine‐adjusted incidence function to forecast COVID‐19 transmission dynamics for the United States, Brazil, India, and the Syrian Arab Republic

Country	t p (days)	I p	ρ 0
USA	261	3 586 865	1.0299
Brazil	214	774 449	1.0099
India	237	843 987	1.0017
Syrian	197	20 246	1.0031
ArabRepublic

In this section, we discussed the forecasting of the new coronavirus disease based on the closed‐form solutions (20) of the SIQD model with a quarantine‐adjusted incidence function for five countries, namely, Pakistan, the United States, Brazil, India, and the Syrian Arab Republic. We used the collected data of the infectious cases, the total cases, and the mortality cases up to date August 2 in the five countries to calculate the infection rate β, isolation rate γ, virus‐induced mortality rate of infected individuals α, and scaling factor ϕ for the model. We assume that the virus‐induced mortality rate of quarantined class μ is 10−7 days −1 for the pandemic. We used the least square method with the general function for this estimating. Also, we computed the peak of the infection function and the predict date of this peak for the five countries. We have compared the estimated and actual epidemic peak t and the number of infected individuals I at the peak for Pakistan. We conclude that our model provided close estimates to the actual data. Moreover, we have computed the basic reproduction number based on the fitting parameter for this model which lies in the range [1.0013–1.0299]. The smallest value of the basic reproduction number is for Pakistan, and the highest value is for the United States. It is important to mention here that highest number of the new coronavirus appeared in the United States up to date of writing this work.

CONCLUDING REMARKS

We have developed an SIQD model with a quarantine‐adjusted incidence function. The closed‐form solutions for all variables of model are established by utilizing the classical techniques of solving ODEs, and these hold provided . We have provided the closed‐form expressions for FOI, rate at which susceptible becomes infected, the epidemic peak and time required to attain this peak. The management of epidemic perceptive using control and prevention strategies is explained as well. The epidemic starts when ρ 0 > 1, the peak of epidemic appears when number of infected attains peak value when , and the disease dies out ρ 0 < 1. The epidemic control strategies suggest to lower the value of the reproduction number and reduce it to less than, that is, ρ 0 < 1. This can be achieved in different ways. The epidemic peak reduces, and time span to attain this peak increases by reducing infection rate β or increasing isolation rate of infected individuals γ. Another effective strategy to reduce epidemic peak and increase the time span to attain this epidemic peak is to set the value of basic reproduction number at some fixed level close to the threshold value and then change parameters β and γ in a fixed proportion. The effect of change of parameter on the FOI and rate of infection is also presented in detail. We have provided the forecasting of the new coronavirus disease based on the closed‐form solutions of SIQD model with a quarantine‐adjusted incidence function and the scaling factor for five countries, namely, Pakistan, the United States, Brazil, India, and the Syrian Arab Republic.

The effective control and prevention policies help in slow down pace of epidemic. The slow transmission dynamics will be helpful for a country to improve health care facilities by increasing capacity of beds, intensive care units, ventilators, and hiring more health care staff to deal with large number of patients. One can also use appropriate combination of reducing infection rate (lockdown) and increasing quarantine rate to lower the value of the reproduction number to its threshold value. This will help in epidemic end. This model can be extended by adding exposed or asymptotic compartment and will be considered in a future work.