key: cord-1039131-399oxm32
authors: Das, Parthasakha; Nadim, Sk Shahid; Das, Samhita; Das, Pritha
title: Dynamics of COVID-19 transmission with comorbidity: a data driven modelling based approach
date: 2021-03-08
journal: Nonlinear Dyn
DOI: 10.1007/s11071-021-06324-3
sha: 9579102d486cc3a8ba6be1fcebc87b3292a9200a
doc_id: 1039131
cord_uid: 399oxm32

An outbreak of the COVID-19 pandemic is a major public health disease as well as a challenging task to people with comorbidity worldwide. According to a report, comorbidity enhances the risk factors with complications of COVID-19. Here, we propose and explore a mathematical framework to study the transmission dynamics of COVID-19 with comorbidity. Within this framework, the model is calibrated by using new daily confirmed COVID-19 cases in India. The qualitative properties of the model and the stability of feasible equilibrium are studied. The model experiences the scenario of backward bifurcation by parameter regime accounting for progress in susceptibility to acquire infection by comorbidity individuals. The endemic equilibrium is asymptotically stable if recruitment of comorbidity becomes higher without acquiring the infection. Moreover, a larger backward bifurcation regime indicates the possibility of more infection in susceptible individuals. A dynamics in the mean fluctuation of the force of infection is investigated with different parameter regimes. A significant correlation is established between the force of infection and corresponding Shannon entropy under the same parameters, which provides evidence that infection reaches a significant proportion of the susceptible.

COVID-19 was announced by the World health organization (WHO) as a major health hazard at the end of 2019. This was first seen in Wuhan province, China [1] . After that this infectious disease has been spread out in the whole world. To defeat the epidemic, the different countries of the world were in a lockdown scenario, the also successively unlock process due to socioeconomic and political impacts [2] [3] [4] [5] [6] . Moreover, the mortality and morbidity rate has varied across countries of the worlds [7] [8] [9] . COVID-19 is often present in human being as common symptoms, for example, cold-like illness including fever, muscle pain, fatigue, loss of taste and smell, sore throat, and dry cough [10] [11] [12] [13] .

COVID-19 can be transmitted from people to people through respiratory droplets from an infected human being or direct contact with contaminated surfaces or objects [14, 15] . As all vaccinations are in trial positions for the prevention of COVID-19. Various studies have pointed out the stage of comorbidity coinfection (like lung disease, heart disease, diabetes, etc) among the infected people. According to the Centers for Disease Control and Prevention, all people with comorbidity or medical condition have more risk to become infected than healthy or normal human being [16] [17] [18] . From the survey of clinical reports, this is further suggested that the confirmed COVID-19 patients are most having comorbidity co-infection which includes kidney disease, obesity those having basal metabolic rate (BMR) more than thirty, type-2 diabetes, COPD (chronic obstructive pulmonary disease) [19] [20] [21] . Moreover, sometimes this confirmed COVID-19 with comorbidity were admitted to the intensive care unit (ICU).

New daily cases of COVID-19 follow a very highly volatile pattern, i.e, deceasing-increasing arbitrarily [22, 23] . To implement fruitful public health measures in a scheduled time within a certain resource according to geographical state [24] , it is utmost essential to study the diffusion or force of infection among the wider population. In order to explain this observation, entropy helps us to determine the heterogeneity of a daily number of cases and time of highest diffusion or force of infection [25] . Various methods of measuring entropy, like weighted entropy [26] , maximum entropy [27] , evolutionary Entropy [28] , structural Entropy [29] are employed to forecast and study the transmission dynamics of pandemic like, COVID-19 [30, 31] . Here, the concept of Shannon entropy is indeed used and seemed to be appropriate to determine the spread of an epidemic as there is a significant analogy based on Boltzmann's classical thermodynamical paradigm [32] .

Various mathematical models have been formulated and studied the corresponding dynamics [33, 34] . Lie et al. proposed a data-driven COVID-19 model with distributed delay [35] . Global dynamics of a COVID-19 SEIR epidemic model has been investigated by Khyar et al. [36] . Further, different authors have suggested control strategies for disease transmission of COVID-19 [37] [38] [39] . The main aim of this paper is to study the qualitative behavior of COVID-19 transmission dynamics by bifurcation theory. Various qualitative dynamics have been studied in different mathematical model [40] [41] [42] . For this purpose, data-driven modeling is a powerful tool to investigate dynamical behav-ior in infectious disease [43] . Various mathematical models have been devolved as well as investigated the transmission dynamic of COVID-19 [44] [45] [46] [47] [48] . Indeed, considering comorbidity in modelling of COVID-19 is of noteworthy interest in ongoing studies and also has not been yet investigated according to literature.

The subsequent parts of the paper are as follows: In Sect. 2, we propose and explore a mathematical model of COVID-19 with comorbidity. The qualitative dynamics is investigated at feasible equilibrium by performing bifurcation analysis in Sect. 3. In Sect. 4, the model is calibrated and sensitivity analysis is also performed. Moreover, the transmission dynamics of mean fluctuation in force of infection is studied by measuring the corresponding Shannon entropy. Finally, in Sect. 5, we discuss and conclude the results from our proposed study.

We develop here a SCEAIHR-type model by introducing threatened or comorbidity susceptible and hospitalized individuals in the dynamics of COVID-19 (SARS-Cov-2). The model consists of seven compartments, namely, susceptible individuals (S), comorbidity susceptible individuals (C), exposed individuals (E), asymptomatic individual (A), infected or infectious individuals with clinical symptoms (I), diagnosed and hospitalized individuals (H) and recovered individuals with no more infections (R). The total number of individuals is N = S + C + E + A + I + H + R. Our main aim is to study the dynamics of COVID-19 with comorbidity co-infection. Based on a biological viewpoint, we formulate a nonlinear mathematical model to investigate the COVID-19 or SARS-Cov-2 pandemic in a certain time window: 

The following initial values are considered in model (1) as follows:

Here, it is assumed that t ≥ t 0 , where t 0 indicates the initial date of the outbreak of COVID-19 for the model (1) . Susceptible individuals are generated by the recruitment of persons by birth or immigration to the community at a constant rate of s . Population decreases after infection and can be transmitted at a rate β s I N , δ a β s A N and δ 2 β s H N by direct interaction between susceptible individuals with symptomatic infected, asymptomatic infected and hospitalized individuals respectively and at a rate ρβ s I N , ρβ s δ a A N and ρβ s δ h H N by direct interaction between co-morbid susceptible individuals with symptomatic infected, asymptomatic individuals and hospitalized individuals respectively. Also, susceptible individuals often develop co-morbidity at a rate of φ s . It is assumed that there is a recovery rate for asymptomatic, symptomatic, and hospitalized individuals, such as γ a , γ i , and γ h respectively. The model incorporates certain demographic impacts by estimat-ing the disease-induced mortality rate η of hospitalized persons and natural death rate μ of each of the seven subpopulations. After the disease incubation period, 1 α e , a fraction ξ of the exposed individuals are symptomatically infected and the remaining fraction (1 − ξ) become asymptomatically infected. The exposed population also decreased at a rate of μ due to natural death. The pictorial diagram of our model is presented in Fig. 1 . The description of model parameters is provided in Table 1 .

In this section, we study some basic properties of the SCEAIHR model (1) such as positivity and boundedness, basic reproduction number, stability analysis at biologically feasible equilibrium points with nonnegative initial conditions (2) , the solutions of SCEAIHR model (1) remain positive through out the region in R 7 + for all t > 0.

Proof To establish the positivity of the system (1), we take any solution emerging from non-negative region R 7 + which remain positive for all t > 0. In order to do this, we prove that the points of vector field on the each hyper plane are bounded by non-negative region R 7 + . In the system (1), we see that d S dt | S=0 = S ≥ 0, 

which is further a positive invariant set to the SCEAIHR model (1) . Therefore, all solution trajectories initiating from interior of R 7 + always remain within the domain . this assures that the growth of all individuals can't be unbounded or exponential for time time window.

The basic reproduction number is a threshold value and characterizes the infectious disease. This is denoted by R 0 and defined by the number of secondarily infected individuals induced by single infected individuals as infectious during its incubation period within susceptible individuals. R 0 is a dimensionless number and quantifies the expectation of decreasing or increasing the disease outbreak. The system (1) has infection free equilibrium 0 , given by

, 0, 0, 0, 0, 0, 0). The compartment E, A, I, H in system (1) are explicitly associated with disease outbreak. Thus we get matrixF,Ṽ for new infection and transition part respectively, given bỹ

The infection-free equilibrium 0 is locally asymptotically stable if R 0 < 1 and unstable R 0 > 1.

Proof The variational matrix J of the system (1) is given by

For our convenience, we consider

The matrix J 0 at 0 has unique and repeated eigenvalues as −(ρφ s + μ) and −μ respectively. The rest of eigenvalues are the roots of following equation as follow:

which can be expressed as follow:

Denote

Consequently, n 11 (0) + n 12 (0) + n 13 (0) = n 1 (0) = R 0 < 1, which implies |n 1 (λ)| ≤ 1. So, eigenvalues corresponding characteristic equation n 1 (λ) = 1 have negative real parts for R 0 < 1. Thus, for R 0 < 1, all eigenvalues are negative. According to Routh-Hurwitz stability criterion, infection free equilibrium 0 is locally asymptotically stable for R 0 < 1. Again, if we consider R 0 > 1, i.e, n 1 (0) > 1, then lim λ→∞ n 1 (λ) = 0; then, there exits λ * 1 > 0 so that n 1 (λ) = 1. This means that there exits non-negative eigenvalue λ * 1 > 0 in matrix J 0 . Hence, infection free equilibrium 0 is unstable for R 0 > 1.

Here, we study the existence of endemic equilibrium. Let * = (S * , C * , E * , A * , I * , H * , R * ) be any endemic equilibrium of system (1). Moreover, The force of infection is

By solving the system (1) at equilibrium state, we

Substituting the all expressions into (5), we get the non-zero equilibrium of (1) which satisfies quadratic equation in ζ * as follow:

where,

From (6), we can obtain endemic equilibrium of (1).

Here, the coefficient A is always positive. The sign of coefficient C depends on lesser or greater value of R 0 than unity. It can be noted that multiple endemic equilibrium can be coexist for R 0 < 1. Thus, we can summarize as follows:

The system (1) has 1. two endemic equilibrium if P 0 > 0 and P 1 < 0 and P 2 1 − 4P 2 P 1 > 0.

2. at least one endemic equilibrium if P 0 < 0 iff R 0 > 1. 3. a unique endemic equilibrium if (P 0 = 0 and P 1 < 0) or P 2 1 − 4P 2 P 1 = 0. 4. no endemic equilibrium if P 0 > 0 iff R 0 < 1 and P 1 > 0.

From the above theorem, it can be noted that backward bifurcation is possible due to the coexistence of stable infection-free and endemic equilibrium for R 0 < 1.

For this purpose, we find a critical value R * 0 of R 0 for unique endemic equilibrium as follow:

Hence, backward bifurcation can occur with specific restriction R * 0 < R 0 < 1. Now, we turn on the analysis for backward bifurcation.

In order to investigate backward bifurcation at endemic equilibrium * (S * , C * , E * , A * , I * , H * , R * ), we apply centre manifold theorem by considering bifurcation parameter β s = β * s , corresponded to R 0 = 1. From variational matrix J in (3), we again consider few entries as follows:

Now, the variational matrix of system (1) at β s = β * s has right eigenvector corresponding to zero eigenvalue,

Similarly, at at β s = β * s , the variational matrix of system (1) has left eigenvector corresponding to zero eigenvalue, given by w = [w 1 , w 2 , w 3 , w 4 , w 5 , w 6 , w 7 ], where

We introduce a few notations for SCEAIHR model system as follows:

. . , 7. Now, we compute second order partial of f i at infection-free equilibrium 0 and obtain

,

,

The remaining partial derivatives at 0 are zero. Now we evaluate the coefficient a and b of well-established Theorem 4.1 in Castillo-Chavez et al. [49] as follows:

Now, we substitute all above values to find the coefficient a and b at threshold β * = β * s , we get

Here, the value of b is always positive, the system (1) undergoes backward bifurcation

where,

Theorem 3 The SCEAIHR model system (1) exhibits a backward bifurcation at R 0 = 1, if the parametric restriction in (9) holds.

Moreover, as comorbidity individual becomes always infected, backward bifurcation occurs in SCEAIHR model. Now for a case, it can be noted that comorbidity individuals do not get infection due to maintaining social distancing, wearing mask and proper sanitation, i.e, ρ = 0. The bifurcation coefficient a is as follow :

As, a < 0, b > 0, by Theorem 4.1 in Castillo-Chavez et al., a transcritical bifurcation occurs at R 0 = 1 for ρ = 0 and endemic equilibrium and the endemic equilibrium * is locally asymptotically stable for R 0 > 1. Here, a backward bifurcation do not exist for ρ = 0.

Theorem 4 The SCEAIHR model system (1) experiences a transcritical bifurcation at R 0 = 1 with parametric restriction in (9) holds, i.e, ρ = 0.

In this section, we perform numerical simulation and its biological interpretation to complement the analytical findings. In the previous section, local stability for infection-free equilibrium and the existence of endemic equilibrium are studied. Moreover, backward bifurcation and transcritical bifurcation are observed. In order to validate the analytical findings, we estimate parameter values of the SCEAIHR model (1). (1) is fitted for new hospitalized COVID cases in India on a regular basis. Due to the high infectivity, the notified patients become hospitalized quickly and it is also easy to confirm the hospitalized cases from the recorded results. We have enlisted the main model parameters, which are estimated from the data in Table 1 . By adopting the model to the new regular cases reported, six parameters of the model have been estimated, such as (a) transmission rate (β s ), (b) modification factor for asymptomatic (δ a ) and symptomatic (δ h ), (c) rate of co-morbidity development by susceptible individuals (φ s ), (d) modification factor for comorbidity development (ρ) and (e) average hospitalized rate of infected individuals (θ ). The data will also estimate some of the initial conditions of the model (1). In MATLAB, the nonlinear least-square solver fmincon has been used during the defined time span to fit the simulated new daily data of COVID-19 recorded by India. The 95% confidence area is generated by using the Delayed Rejection Adaptive Metropolis (DRAM) algorithm. A description of this model fitting technique is given in [54] . Tables 1 and 2 provide the parameters and initial conditions respectively, which are estimated by the above technique. The fitting of India's daily new hospitalized COVID-19 cases is seen in Fig. 2 . The basic reproduction number is estimated at R 0 = 1.3607 Table 1 and parameters that are estimated from the model. In order to recognize the most significant parameters to infected individual, we next perform sensitivity analysis followed by uncertainty analysis.

We apply PRCC(Partial rank correlation coefficient) method for sensitivity analysis and LHS (Latin hypercube sampling) method for uncertainty analysis. A detailed technique can obtained in Marino et. al, [55] . We assign uniform distribution the parameters namely, β s , δ a , δ h , ρ, θ , φ s η and ξ with 95% confidence domain. We consider baseline values as parameter values. The outputs are given in Table 3 with bar diagram in Fig. 3a . From Fig. 3a , it can be seen that β s , ρ and ξ are influential parameter to infected individual. Moreover, we study the sensitivity indices of basic reproduction number R 0 . Here, only φ s is negatively correlated, and β s , δ a , ρ, and ξ are positively correlated to R 0 , given in Fig. 3b . As basic reproduction number (R 0 ) quantifies the expectation of decreasing or increasing the epidemic evolution, we next study the effects of parameter variation on R 0 .

We further study the effects of parameter variations on basic reproduction number Fig. 4 . We observe that only the increasing value of ρ can change R 0 < 1 to R 0 > 1 in Fig. 4a. From Fig. 4b , the increasing value simultaneously changes the value of R 0 . Further, it can be also seen that the value of R 0 always remains less than one within the parametric plane of (φ s , β s ) ∈ (0, 0.1] × (0, 2].

On investigation of the existence of endemic equilibrium, we derive the sub-threshold range of bistable equilibrium in the SCEAIHR model (1). Thus backward bifurcation can be seen within the domain R * 0 < R 0 < 1 (R * 0 = 0.4617) with ρ ∈ [0.01, 0.8], given in Fig. 5a , where two positive endemic equilibrium coexist, i.e, asymptotically stable equilibrium (blue line) and unstable equilibrium (red line). From an epidemiological viewpoint, the scenario of backward bifurcation indicates that the model possesses two endemic equilibrium. However, infection-free equilibrium can exist only for R 0 < R * 0 . Indeed, this restriction is sufficient for the eradication of infection from the system. Moreover, it sometimes depends on the initial size of the sub-population of model.

On the other hand, if comorbidity individuals do not acquire infection, i.e, ρ = 0, the model (1) experiences transcritical with φ s ∈ [10 × 10 −6 , 2.7 × 10 −5 ] given in Fig. 5b . Here, only one endemic equilibrium (blue line) and infection-free equilibrium (red line) exist for R 0 > 1. It indicates that force of infection is not strong enough to spread out as all susceptible individuals take the precaution of disease spreading, like, wearing a mask, social distancing, maintaining sanitation. Moreover, it can be also noted that increasing the value of φ s becomes a challenging task for the eradication of infection. Numerically, it can be seen that a backward bifurcation regime increases gradually with increases of φ s with ρ ∈ [0.01, 0.8], given in Fig. 6 . Consequently, the value of R 0 is essential to be reduced to ensure the eradication of infection in such a way that R * 0 and R 0 is close to one. Now we study the variation of infection intensity. For this purpose, we consider average of force of infection for daily basis cases. (1). Moreover, in order to understand the combined effect under (φ s , β s ) ∈ (0, 0.1] × (0, 2] on < ζ * >, a matrix plot is given in Fig. 7c. In Fig. 7c , it can be mentioned that the average force of infection is dependent on value of φ s for any value of β s , defined in efferent color region of increased value of φ s ∈ (0, 0.1]. In order to quantitative measure the degree of disorder in force of infection. Shannon entropy [32] is employed for this purpose. Here, disorder means the non-uniform distribution of infection under distribution of susceptible.

Shannon entropy is generally defined as

N being length of event ζ * i and p(ζ * i ) is probability of event ζ * i , i.e, relative frequency of occurrence of non-recurrent event ζ * i . The trend of En(ζ * i ) on φ s ∈ (0, 0.1] and β s ∈ (0, 2] are shown in Fig. 8a , b. Figure 8a shows that deceasing trend of En(ζ * i ) with increase of φ s . On the contrary, En(ζ * i ) increases with β s in Fig. 8b. From Figs.7a, b and 8a , b, it can be noted that the trend of En(ζ * i ) is highly correlated with the fluctuation in force of infection ζ * . Moreover, we also plot in Fig. 8c to study the dependence of En(ζ * i ) on Table 1 (φ s , β s ) ∈ (0, 0.1]×(0, 2]. Further comparing Figs. 7c and 8, we observe similar pattern of increasing entropy between ζ * and En(ζ * ). This result reveals the additional evidence about the measure of dynamical disorder in force of infection in the model (1) for parameter regimes.

In this paper, we have proposed and studied a mathematical model for disease dynamics of COVID-19 accounting of impacts of comorbidities on the complication of COVID-19. The proposed model has been calibrated by using new daily cases of India. The qual-itative properties of the model have been studied and the basic reproduction number has been calculated by the method of the next-generation matrix. The model has asymptotically stable infection-free equilibrium for less than a unity of basic reproduction number. Using the center manifold theorem, the phenomena of backward bifurcation have been observed for increasing the modification factor of comorbidity development. This indicates that infection persists for less than a unity of basic reproduction number. The occurrence of backward bifurcation assures the rich dynamics of the model. From an epidemiological viewpoint, comorbidity individuals acquire more re-infection due to lack of surveillance and precautions like wearing masks, social distancing, proper sanitation, etc. Further, endemic equilibrium has become asymptotically stable if the modification factor for comorbidity is zero (ρ = 0), i,e. none of the comorbidity individual gets an infection and keeps maintained themselves with high surveillance. In these circumstances, more susceptible individuals have acquired infection and become exposed individuals as basic reproduction is more than unity, indeed, exposed individuals have a strong immunity to prevent COVID-19 infection in this context. Moreover, the increase of comorbidity development and modification factor simultaneously has increased the backward bifurcation regimes. Consequently, this can lead to a disaster situation.

Sensitivity analysis of the model reveals that transmission rate (β s ), modification factor of comorbidity development (ρ), comorbidity development (φ s ) have a significant influence on infected individuals and basic reproduction number (R 0 ). In order to control disease transmission, a policy might be implemented by enhancing surveillance with proper maintaining sanitation and keep safe from infected individuals.

On investigation of the pattern of transmission dynamics, the mean fluctuations of the force of infection have shown decreasing as well as increasing trend with comorbidity development (φ s ) and transmission rate (β s ). Consequently, the combined effect has shown an increasing trend to the average force of infection. From the dynamical viewpoint, the measure of disorder in force of infection, i.e, the Shannon entropy has provided evidence towards the higher entropy production indicating the more powerful force of infection, i.e, the infection reaches a significant proportion in susceptible individuals. This might lead to a devastating situation in society regarding disease transmission.

Finally, our study assures that the increase of force of infection builds up with the enhanced entropy production attributed to highly disorder in epidemic evolution. The production of entropy are concepts that can be helpful to guide computational biologists to elucidate the associations between mean fluctuation and changes in parameters to develop policy in disease control and novel vaccination strategies. Thus, it would be useful to include in future publications how vaccination policy can be fruitful. 

A novel coronavirus outbreak of global health concern

Transmission dynamics of covid-19 in wuhan, china: effects of lockdown and medical resources

An investigation of transmission control measures during the first 50 days of the covid-19 epidemic in china

After the lockdown: simulating mobility, public health and economic recovery scenarios

Lockdown and unlock for covid-19 and its impact on residential mobility in india: an analysis of the covid-19 community mobility reports

A vulnerability index for the management of and response to the covid-19 epidemic in india: an ecological study

Modeling the effect of contaminated objects for the transmission dynamics of covid-19 pandemic with self protection behavior changes

The impact of covid-19 and strategies for mitigation and suppression in low-and middle-income countries

Impact of covid-19 and comorbidities on health and economics: Focus on developing countries and india

Rapid simulation of spatial epidemics: a spectral method

Return of the coronavirus: 2019-ncov

Clinical features of patients infected with 2019 novel coronavirus in wuhan, china

Potential impact of the covid-19 pandemic on hiv, tuberculosis, and malaria in low-income and middle-income countries: a modelling study

The continuing 2019-ncov epidemic threat of novel coronaviruses to global health-the latest 2019 novel coronavirus outbreak in wuhan, china

Pandemic potential of 2019-ncov

Prevalence of comorbidities and its effects in patients infected with sars-cov-2: a systematic review and meta-analysis

Comorbidity and its impact on 1590 patients with covid-19 in china: A nationwide analysis

Prevalence of covid-19-related risk factors and risk of severe influenza outcomes in cancer survivors: A matched cohort study using linked english electronic health records data

Diabetes and covid-19: evidence, current status and unanswered research questions

Impact of comorbid asthma on severity of coronavirus disease (covid-19)

A novel covid-19 and its effects on cardiovascular disease

Spatial dynamics of an epidemic model with nonlocal infection

Pattern transitions in spatial epidemics: mechanisms and emergent properties

Migration rate estimation in an epidemic network

Entropy-based pandemics forecasting

Exploring dynamical complexity in a time-delayed tumor-immune model

Maximum entropy method for estimating the reproduction number: An investigation for covid-19 in china and the united states

Evolutionary entropy determines invasion success in emergent epidemics

Structural entropy of daily number of covid-19 related fatalities

The dynamics of entropy in the covid-19 outbreaks

Entropy ratio and entropy concentration coefficient, with application to the covid-19 pandemic

A mathematical theory of communication

Prediction of bifurcations by varying critical parameters of covid-19

Seir modeling of the covid-19 and its dynamics

Covid-19: datadriven dynamics, statistical and distributed delay models, and observations

Global dynamics of a multi-strain seir epidemic model with general incidence rates: application to covid-19 pandemic

Dynamics and control of covid-19 pandemic with nonlinear incidence rates

Optimal treatment strategies for delayed cancer-immune system with multiple therapeutic approach

Effects of control measures on the dynamics of covid-19 and double-peak behavior in spain

Stochastic dynamics of michaelis-menten kinetics based tumor-immune interactions

An investigation on michaelis -menten kinetics based complex dynamics of tumorimmune interaction

An investigation on monod-haldane immune response based tumor-effector-interleukin-2 interactions with treatments

Modeling infectious disease dynamics in the complex landscape of global health

Modeling shield immunity to reduce covid-19 epidemic spread

Modelling the covid-19 epidemic and implementation of population-wide interventions in italy

Characterizing chaos and multifractality in noise-assisted tumor-immune interplay

Occurrence of backward bifurcation and prediction of disease transmission with imperfect lockdown: A case study on covid-19

Short-term predictions and prevention strategies for covid-2019: A model based study

Dynamical models of tuberculosis and their applications

COVID-19 coronavirus outbreak

Early transmission dynamics in wuhan, china, of novel coronavirus-infected pneumonia

Assessment of lockdown effect in some states and overall india: A predictive mathematical study on covid-19 outbreak

Mathematical analysis of a powerlaw form time dependent vector-borne disease transmission model

A methodology for performing global uncertainty and sensitivity analysis in systems biology

The authors have no conflict of interest concerning the publication of this manuscript.