Sensitivity and elasticity analysis of novel corona virus transmission model: A mathematical approach.

The deadly corona virus continues to pound the globe mercilessly compelling mathematical models and computational simulations which might prove effective tools to enable global efforts to estimate key transmission parameters involved in the system. We propose a mathematical model using a set of non-linear differential equations to account for the spread of the COVID-19 infection with special compartment class isolation or quarantine and estimate the model parameters by fitting the model with reported data of the ongoing pandemic situation in India. The basic reproduction number is defined and local stability analysis is carried out at each equilibrium point in terms of the reproduction number R 0 . The model is fitted mathematically and makes the data India specific. Additionally, we examined sensitivity analysis of the model. These outcomes recommend how to control the spread of corona, keeping in mind contact and recovery rate. Also we have investigated the elasticity of the basic reproduction number as a measure of control parameters of the dynamical system. Numerical simulations were also done to show that the proposed model is valid for the type and spread of the outbreak which happened in India.

Introduction

Mathematical models relating to infectious disease transmission dynamics have assumed a degree of pervasiveness. The epic COVID-19 which started in recent times in Wuhan province of China has unleashed its tentacles across the globe. The World Health Organization declared it a pandemic. In the absence of a vaccine, social distancing has emerged as the most widely endorsed strategy for mitigation and control [1]. India took exceptional measures and suspended all traveler visas as most of con-solidified COVID cases were reported from different nations. The prime minister of India announced a 14 h voluntary public curfew (‘Janata Curfew’) on March 22, 2020 as a precautionary measure to combat COVID-19. The Govt. of India followed it up with lockdowns on March 23, 2020 to keep the fast approaching dangers from the infection in 75 locations across the nation over including significant urban communities where COVID-19 contamination was endemic [2]. Moreover, on March 24, 2020 the Govt. of India mandated country-wide lockdown for 21 days [3] and the lockdown was extended to June 07, 2020. Later, India loosened up lock down restrictions with open period. The MHA issued different guidelines for the last 5 lockdown lifting periods and on October 27, 2020; the Ministry of Home Affairs issued guidelines for activities permitted in Unlock 6.0 [4].

In the context of all that has been mentioned, the mathematical model proposed in the paper comes into its own is a very useful and important tool for the analysis of infectious diseases because this analysis enables prediction of future outbreaks while also detailing strategies to control the epidemic. Mathematical modeling based on a system of differential nonlinear equations may give a sensible tool to explain the elements of COVID-19 transmission. The global problem of outbreak of COVID-19 has attracted the interest of researchers in different disciplines. Okhuese et al. [5] developed an SIR model which showed how to control the infection and Batista et al. [6] developed a logistic growth regression model which is used for the estimation of the final size of the corona virus epidemic. In Ref. [7] Nagaraj et al. discussed invasion of novel COVID-19 in India, to create some awareness about the nature of the pandemic while Aritra Ghosh et al. [8] discussed how well or inadequately India was equipped to deal with an increasing number of COVID-19 cases, in terms of coping with the current adverse effects on the economy, livelihood and environment during the COVID-19 lockdown period along with various approaches undertaken to overcome this pandemic. Lin et al. [9] proposed a model for COVID-19, which effectively catches the time line of the present outbreak while T. Wu et al. [10] developed SEIR model, to clarify the transmission dynamics and projected national and global spread of disease. T. Chen et al. [11] suggested a model for reproducing stage-based transmissibility of a novel corona virus and in Ref. [12] Tang et al. proposed a model that would combine the clinical development of the disease, the epidemiological status of the patient and the requisite extent of intervention. Surbhi Sharma et al. [13] discussed current treatment protocol for COVID-19 in India and Aayushi Kundu et al. [14] discussed some factors such as dietary habits, vaccination, climate conditions which could be the explanation for the contrasting impact of COVID-19 in India and other developed nations. Some authors estimated the value of reproduction number, which helped to predict the outbreak of the disease. In Ref. [15] Read et al. calculated basic reproduction number using an assumption of poission-distributed daily time increments. Some authors [16,17] estimated the mean reproduction number for COVID-19 in the early phase of outbreak, which is slightly higher than that for SARS-CoV. Still other works where the basic reproduction number is estimated for different countries can be found in Refs. [[18], [19], [20], [21], [22], [23], [24], [25], [26]].

The human to human contact is the likely reason for episodes of COVID-19. Therefore, isolation of the infected human overall can reduce the risk of future COVID-19 spread. In order to do this, we divided the total population into five compartments, classifying them as susceptible, exposed, infected, isolated and recovered from the disease.

The remainder of our study is organized as follows. In section 2, we have proposed a mathematical model specific to COVID-19. In section 3, defines the basic reproduction numberwhile in section 4 discusses the local stability of two equilibrium points in terms of. The numerical simulations figure in section 5 along with a discussion of the sensitivity and elasticity analysis of the proposed model. The conclusions are given in section 6.

Mathematical formulation of the model

In this work we proposed a deterministic model for COVID-19 using a set of ordinary non-linear differential equations. The population has been divided into five compartments that is and. Many authors have carried out studies on COVID-19 with different compartments, but in our model we take into account asymptomatic cases, implying cases involving those who are unaware of infections or the limited capacity for testing. For Individuals without any symptoms, it is very difficult to detect the presence of the virus. Hence asymptomatic transmission is the most challenging one in COVID-19. Alongside this we consider another compartment i.e., isolated population which represents home quarantine and the hospitalized.

This model assumes a completely susceptible population with homogeneous mixing. The disease has an incubation period and after that the exposed individuals becomes infectious at the rate. The exposed individuals are removed and added to the infectious at the rate. Generally the infected individuals are isolated when the susceptible individuals are not infected. To factor this we added a compartment representing the quarantine/isolated. The isolated or quarantine individuals are removed from and added to the compartment at the rate. The recovered individuals from and are added to the compartment at a rate respectively. represents the natural death rate and the disease death rate. We assumed the total population to be constant.

From the all above assumptions we consider the following mathematical model for COVID-19 as follows

We can omit without generality the last equation because first four equations are independent of and the initial populations are defined in the following manners

The basic reproduction number

It measures the disease spread in the population and it is defined as the average number of cases produced if one infected person or individual introduced in the population. Here we defined this by using the next generation matrix approach [27]. System (1) has unique disease free equilibrium i.e., corona virus free equilibrium.

The following matrices and represents nonlinear terms of new infection and out flow terms respectively which are given by

Therefore, the linearized matrices at free corona virus equilibrium point is given by

The basic reproduction number is given by the dominant eigen value of i.e., .

Local stability analysis

The Jacobian matrix of system (1) at the free corona virus equilibrium is given by

The characteristic equation of (6) is given bywhere

One eigen value is negative i.e., and remaining three eigen values are determined by the equation.

It is easy to see that then which is equivalent to and to prove the free corona virus equilibrium is locally asymptotically stable we have to show that has only negative roots or roots with negative real parts. For this, we have to use Routh-Hurwitz criteria of three dimensions. The conditions are and. Already we have and clearly. The result gives inequality which gives.Since and on the other hand.

Therefore, .

Hence, by Routh-Hurwitz criteria the free corona virus equilibrium point is locally asymptotically stable if.

The unique corona virus equilibrium point of system (1) exists if which is given by

Clearly if and also when.

The Jacobian matrix for system (1) at is given by

The characteristic equation corresponding to (9) is given by

One Eigen value is and we have to show that remaining three Eigen values of the following equation are negative by Routh-Hurwitz criteria.Where

It can be verified that if. Clearly and by similar calculation we can show that and when. Therefore, by Routh-Hurwitz criteria the corona virus equilibrium point is locally asymptotically stable if.

Numerical simulation

The numerical simulations are performed for the proposed model and the outcomes are contrasted and the genuine information pertaining to India which were acquired from numerous reports distributed by the WHO and World meter [28,29].

We classified the outcomes into two sections-lockdown period and unlock period and our model fits additionally well with genuine information of day by day affirmed cases in previously mentioned two classifications, as shown in Fig. 1 and Fig. 2 . We accepted that the complete populace in the model was fixed. The starting point of simulation is first day of lockdown i.e., March 25 up to June 07 while the unlock period began from June 08 and went up to September 08. We assume the following parameter values to perform the numerical simulations in lockdown period.

Fig. 1

Flow chart of SEQIR model.

Fig. 2

Confirmed cases in India (lockdown period). Dots represent the real data and the line is corresponding to the real data which was obtained by solving system (1) numerically.

To fit the real time data to our model (1) in unlock period we changed two parameters , and the initial infected population is 8442.

Sensitivity analysis and elasticity

In this we performed local sensitivity analysis by numerical simulations. The analysis examines the change in one input parameter value and then the change in output values. In this we fix all parameter values and we vary the particular parameter for several values. In this paper, we are investigating the sensitivity with respect to two parameters. The parameter values were. Now we fix and vary values 0.0001, 0.0005 and 0.0009. The corresponding graph is plotted in Fig. 3 . From this we observe that once the transmission rate or contact rate increases, the number of infected individuals in the population also increases.

Fig. 3

Confirmed cases in India (unlock period). Dots represent the real data and the line is corresponding to the real data which was obtained by solving system (1) numerically.

Here the value remains constant while values are varied 0.001, 0.005 and 0.009. The corresponding graph is given in Fig. 4 . From Fig. 4, we observe that when the recovery rate increases, the number of infected individuals in the population decreases.

Fig. 4

Sensitivity of the model (1) with respect to the transmission rate.

Elasticity of

The definition of sensitivity is local because the sensitivity is computed while all parameters, including parameter, are kept at their estimated values. However, this method does no longer absolutely explore the input space; in view that it does not take into account the simultaneous variation of input parameters. Another disadvantage is that it depends strongly on the magnitude of and the quantity. In this respect, a much greater useful concept is elasticity. The Elasticity of quantity with respect to the parameter is given by, which means the percentage change in with respect to the percentage change in the parameter. Now we are computing the elasticity of with respect to the parameters.

It is easy to show that the elasticity of with respect to the parameter is 1, i.e.,.

Computing the elasticity of reproduction number with respect to the is

This gives the result that 1% of increase in will produce 0.6% decrease in. From these results the elasticity's suggest that the recovery rate is more effect on, if we compare with transmission rate.

From these results we conclude that our model supports real data for the parameter values (Table 1 ) and the virus spread depend on transmission rate or contact rate and recovery rate from infection.

Table 1

Parameter values corresponding to the model (1) in India, for which.

Parameter	Value	Remark
A	1380 million	[30]
β	0.0009	Estimated
α	0.0037	Estimated
γ	0.0055	Estimated
q	0.00052	Estimated
θ	14 days	From MOHFW data
μ	0.0028	Estimated
S(0)	1380	[28]
I(0)	121	[28]

Conclusions and discussions

Recently in last decades as sensor-based diagnostic technologies are becoming very popular for importance in various interdisciplinary research field like as health monitoring, environmental monitoring etc., so in this pandemic situation we come to terms with the spread of COVID-19 infection where a wide scope of mathematical models have been called into play for inner dynamic analysis. In our present work, we analyze the observed data of COVID- 19 cases in India using the SEIQR mathematical model with quarantine or isolated compartment. In COVID-19 disease we can observe that human to human contact is the expected reason for outbreaks. Hence, isolation of the infected human generally can decrease the danger of future COVID-19 spread. For the proposed model we examined local stability at equilibrium points and these results defined in terms of. In this examination, we fitted our model based on data that was specific to India. These outcomes are introduced in Figs. 2 and 3. Further we discussed the sensitivity analysis of the model (1) with respect to the parameters transmission rate and recovery rate (Figs. 4 and 5 ). From the outcomes we can say that the best factor on spreading corona virus is contact rate and recovery rate. Additionally we examined the elasticity of basic reproduction number with respect to the transmission rate and recovery rate. These results indicate that the recovery rate highly impacts basic reproduction number. From the results, we observed that if the rate of transmission of disease is expanded, the disease will spread rapidly, and to control COVID-19 spread, we need to focus on contact rate. To control contact rate, we carefully follow physical distancing, wearing cover or face shield and so forth, in any case no control is possible. Among infected people, without isolation, tracking contacts is an important factor in fighting the epidemic. It is helpful to have a moderate analysis in the follow-up and isolation, as quickly as can be expected under the circumstances, so that those who may have been infected with the virus have not transmitted the disease. All the information utilized in this work has been obtained from authentic sources. By and large the numerical models have certain constraints, while our mathematical model performs rather well. Our model is based on data available on public platforms until September 8, 2020. Future models can include greater granularity as more data become available, and dynamics of the COVID-19 virus becomes better known.

Fig. 5

Sensitivity of the model (1) with respect to the recovery rate.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.