Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
Nonauton. Dyn. Syst. 2021; 8:75–86 Research Article Open Access Swati Tyagi, Shaifu Gupta, Syed Abbas*, Krishna Pada Das, and Baazaoui Riadh Analysis of infectious disease transmission and prediction through SEIQR epidemic model https://doi.org/10.1515/msds-2020-0126 Received November 3, 2020; accepted March 25, 2021 Abstract: In literature, various mathematical models have been developed to have a better insight into the transmission dynamics and control the spread of infectious diseases. Aiming to explore more about various aspects of infectious diseases, in this work, we propose conceptual mathematical model through a SEIQR (Susceptible-Exposed-Infected-Quarantined-Recovered) mathematical model and its control measurement. We establish the positivity and boundedness of the solutions. We also compute the basic reproduction num- ber and investigate the stability of equilibria for its epidemiological relevance. To validate the model and estimate the parameters to predict the disease spread, we consider the special case for COVID-19 to study the real cases of infected cases from [2] for Russia and India. For better insight, in addition to mathematical model, a history based LSTM model is trained to learn temporal patterns in COVID-19 time series and predict future trends. In the end, the future predictions from mathematical model and the LSTM based model are compared to generate reliable results. Keywords: Infectious disease; Mathematical model; Stability analysis; Long Short Term Memory networks (LSTM); Parameter estimation MSC: 00A71; 92B20; 34D20 1 Introduction In history, various pandemics and epidemics have emerged out several times, which have demolished human- ity, sometimes resulting into end of civilizations and tremendous change in the course of history. Mathemati- cal modelling plays an important role in understanding the complexities of such infectious diseases and their control. The initial models of infectious diseases primarily focused on the theoretical disease research. More recently, there have been eorts to broaden the eld still further by incorporating many aspects of complex- ity of natural systems. Provided this rich history in development of fundamental theory, the applied aspect of the questions being asked, and the extensive data collected on many disease systems, the eld of infec- tious disease modelling represents one of the richest areas of research at the interface of pure theory and data. Mathematical modelling can be used to study the mechanisms underlying observed epidemiological patterns, assessing the eectiveness of control strategies, and predicting epidemiological trends. The investigation of dynamics of infectious disease transmission can be reached from various frameworks, such as relatively sim- ple curve tting techniques to standard SIR, SEIR etc. compartmental models to complex stochastic models using simulations. To model the outbreaks of infectious diseases, there are two major components, namely; Swati Tyagi: Department of Mathematics, Amity University, Punjab, 140306, India Shaifu Gupta: Department of Computer Science & Engineering, Indian Institute of Technology Jammu, Jammu, Jammu and Kashmir, 181221, India *Corresponding Author: Syed Abbas: School of Basic Sciences, Indian Institute of Technology Mandi, Mandi, 175005, H.P., India, E-mail: sabbas.iitk@gmail.com; abbas@iitmandi.ac.in Krishna Pada Das: Mahadevananda Mahavidyalaya Monirampore, P.O.-Barrackpore, Kolkata 700120, India Baazaoui Riadh: College of Sciences and Humanities, Prince Sattam bin Abdulaziz University, Saudi Arabia Open Access. © 2021 Swati Tyagi et al., published by De Gruyter. This work is licensed under the Creative Commons Attribution alone 4.0 License. 76 Ë Swati Tyagi et al. (i) the basic reproduction number (R , the average number of secondary infections produced by primary in- fection while coming in contact with the entire susceptible population) and (ii) the generation time (average time from the onset of symptom in primary infectious individual to symptom onset in secondary case). These two components jointly determine the likelihood and speed of the outbreaks of epidemic. Furthermore, all the modelling techniques largely depend on the availability of data for the corresponding infectious disease under investigation, in order to provide the robust estimate the parameters and their distribution. In case of emergence of a new infectious disease, rapid growth of literature on mathematical modelling also plays as a key factor. Mathematical modelling also helps in assessment of potential long-term impacts of such diseases and provides means to evaluate and predict the eects of possible interventions, even with the availability of limited data. Some typical examples are profound early models for HIV in 1980s, pandemic inuenza, Creutzfeldt–Jakob Disease, outbreak of severe acute respiratory syndrome (SARS) in 2003, Nipah virus etc. Very recently, the novel coronavirus, also known as COVID-19, 2019-nCOV or SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) has provided the world with an unparallel challenge and is crucially threatening the mankind [1, 2]. The rst case of COVID-19 was identied in the Wuhan city, Hubei Province of China on December 31, 2019, and since then at this stage, it has spread widely in 215 countries and territories around the whole world [1, 2] with the infections increasing exponentially. The coronavirus is a zoonotic disease, where primary host was animals, from there it transmitted to humans [3]. The infected people face u-like symptoms, which leads to critical illness, respiratory problems, primarily to those which have other diseases like diabetes, asthma, blood pressure, heart disease etc. [2]. On 30 January 2020, the epidemic of coronavirus disease 2019 (COVID-19) was declared as a Public Health Emergency of International Concern, the highest level in the World Health Organization’s (WHO) emergency response for infectious diseases. The number of cases was accelerating in China and subsequently all over the world; however, the publications of the COVID-19’s potential transmission and eectiveness of government interventions were limited. It is urgent to provide more scientic information for a better understanding of the novel coronavirus and further containment of the outbreak. Till date, more than 9.9 million people are in- fected with COVID-19 according to WHO and 0.49 million people have died due to the infection. As a preventive measure, various countries have implemented a strict lock-down and their administrations are encouraging and advising people to stay at home and maintain social distancing and proper hygiene. This pandemic has impacted immensely in a negative manner on world politics, education, economy, sociality, and various other global aspects [4, 5]. Knowledge of early spread dynamics of infection and guring out the possible control measures is critical to assess the extensive outspread. Besides medical and biological research, it becomes an emergency to for- mulate a mathematical model which eectively describes the development and transmission of disease. This will help to make signicant decisions and preventive measures based on the eective assumptions provided by the model. This paper is organized as follows. In Section 2, we propose a SEIQR mathematical model to predict the dynamics of infectious disease. In Section 3, the positivity and boundedness of the solutions are described. We also compute the basic reproduction number and investigate the stability of the equilibria. In Section 4, we perform the data analysis as a special case for COVID-19 in comparison with the model solution and also predict the disease dynamics. We train the proposed mathematical model using the data available on [2] to estimate the parameters, as a case for Russia and India. In Section 5, we present the deep learning model to be used for learning past trends and to predict future trends. The experimental results of this model are presented in Section 6. Section 7 presents a comparative evaluation of results of mathematical and temporal analysis based machine learning model. At the end we conclude our paper in Section 8. 2 Model formulation In this work we are investigating SEIQR model based on the fact that infected/symptomatic individuals are put into quarantine for treatment (hospital or home remedies), from which they move into recovered class. The Analysis of infectious disease transmission and prediction through SEIQR epidemic model Ë 77 SEIQR, an extensively used epidemic model, can reect the ows of people between ve states: susceptible (S), exposed (E), infectious (I), quarantine (Q) and recovered (R). The total host population is partitioned into susceptibles, exposed (in the latent period), infectious, quarantined and recovered, with the densities, respectively, denoted by S(t), E(t),I(t), Q(t) and R(t). The total population size at time t is denoted by N(t), with N = S + E + I + Q + R. The model is described by the following system of dierential equations: dS = A − β SE − β SI − μS 10 20 dt dE = β SE + β SI − E − μE (1) 10 20 0 dt dI = E − k I − α I − μI 0 0 0 dt dQ = k I − k Q − α Q − μQ 0 1 0 dt dR = k Q − μR. dt with non-negative initial conditions; the positive parameters μ is the rate of natural death; α is non- negative constant and denote the rate of disease caused death. The parameter denote the transfer rates between the exposed and the infectious. The constant k is the rate at which the infectious individuals are quarantined; the constant k is the rate at which the quarantine individuals recover; β , β are the rate of 1 10 20 the ecient contact in the latent and infected; A are constant recruitment of susceptibles respectively. The initial conditions are S(0) > 0; E(0) ≥ 0; I(0) ≥ 0; Q(0) ≥ 0 and R(0) ≥ 0. 3 Analysis of the system 3.1 Non-negativity and boundedness of solution Theorem 3.1. The solutions of the system (1) are non-negative and bounded. dS Proof. From (1), we have ≥ −μS. This implies dt S(t) ≥ S(0)exp(−μt) > 0. Similarly we have E(t) ≥ E(0)exp(−( + μ)t) ≥ 0, I(t) ≥ I(0)exp(−(k + α + μ)t) ≥ 0, 0 0 Q(t) ≥ Q(0)exp(−(k + α + μ)t) ≥ 0, 1 0 R(t) ≥ R(0)exp(−μt) ≥ 0. Next, let N = S + E + I + Q + R. Then we have dN + μN ≤ A − α I − α Q ≤ A. 0 0 dt Integrating both sides of the above inequality and using theory of dierential inequalities, we obtain −μt −μt 0 < N(S, E, Q, I, R) ≤ 1 − e + N S(0), E(0), Q(0), I(0), R(0) e and letting t ! ∞, we have 0 < N ≤ . Thus all solutions of (1) that initiate in fR \ 0g are conned in the n o 5 A region Ω = (S, E, Q, I, R) 2 R : N = + ϵ for any ϵ > 0 and for t ! ∞. This proves the theorem. μ 78 Ë Swati Tyagi et al. 3.2 Computation of Reproduction number Following the similar approaches as in [6], it can be veried that the system (1) always has the disease free 0 A equilibrium (DFE) P = (S , 0, 0, 0, 0) = ( , 0, 0, 0, 0) in the absence of infection. Furthermore, let R = 0 0 −1 −1 ρ(FV ) represents the spectral radius of the matrix FV , then R is given by β S −1 0 20 R = ρ(FV ) = . (2) (k + α + μ)( + μ − β S ) 0 0 0 10 The parameter R is referred to as the Basic reproduction number. The number R indicates the transmission 0 0 of infectious disease, for R > 1, the outbreak is expected to continue and expected to end for R < 1. 0 0 It can be observed from (2) that R signicantly depends upon the average rate of contact between sus- ceptible and infected, and transfer rate between susceptible and latent individuals. Thus, we can say lock- down/quarantining individuals plays a signicant role in controlling the spread of infectious disease. 3.3 Stability analysis Stability analysis of DFE P : Theorem 3.2. If β S < + μ, then the disease-free equilibrium point is locally asymptotically stable. 10 0 Proof. In order to check the stability of the equilibrium points, we calculate the eigenvalues using Jacobian matrix. The Jacobian at DFE is given by 0 1 0 0 −μ −β S −β S 0 0 10 20 B 0 0 C 0 β S − ( + μ) β S 0 0 B 10 0 20 C B C J = B 0 −(k + α + μ) 0 0 C DFE 0 0 0 B C @ A 0 0 k −(k + α + μ) 0 0 1 0 0 0 0 k −μ β S − −(k +α +2μ) 10 0 0 0 The characteristic roots of the Jacobian matrix J are −μ, −μ, −(k + α + μ) and ± DFE 1 0 0 2 0 (k +α +μ+β S − −μ) +4 β S 0 0 10 0 0 20 . Clearly all ve roots are negative .Thus the disease-free equilibrium point (P ) is locally asymptotically stable. Further to prove the global asymptotic stability of DFE, we use theorem by Castilo-Chavez et.al. [7]. Theorem 3.3. [8] The given mathematical model (1) can be represented in the form dX dY = F(X, Y), = G(X, Y), G(X, 0) = 0 (3) dt dt where X = S; Y = (E, I, Q) denote the classes of uninfected and infected individuals respectively. The DFE is represented by P = (X , 0) = ( , 0). 0 0 For global asymptotic stability of P , condition (H1) and (H2) must be satised. dX (H1) For = F(X , 0), X is globally asymptotically stable. 0 0 dt ˆ ˆ (H2) G(X, Y) = MY − G(X, Y), G(X, Y) ≥ 0. Here M = D G(X , 0) is M-matrix, that is all o-diagonal elements Y 0 are non-negative. If given system (1) satises (3), then P = (X , 0) is a globally asymptotically stable equilibrium of given 0 0 mathematical model. Theorem 3.4. The equilibrium P = (X , 0) of system (1) is globally asymptotically stable provided the condi- 0 0 tions in (3.3) are satised. Proof: Using Theorem 3.3 to our model system (1), we obtain F(X , 0) = A − μS, G(X, Y) = MY − G(X, Y) 0 Analysis of infectious disease transmission and prediction through SEIQR epidemic model Ë 79 where 2 3 2 3 −( + μ) 0 0 β SE + β SI 0 1 2 6 7 6 7 M = −(k + α + μ) 0 , G(X, Y) = 0 . 4 5 4 5 0 0 0 0 k −(k + α + μ) 0 0 1 0 We can observe that G(X, Y) ≥ 0 8 (X, Y). Furthermore, the matrix A is M-matrix. Thus the DFE P is globally asymptotically stable. Existence of Endemic equilibria and its stability: * * * * * In addition to DFE P , the model (1) admits a unique endemic equilibrium point P = (S , E , Q , I , R ), where + μ * 0 S = β + d β 10 1 20 Ad β R − 1 ( ) * 1 20 0 E = R ( + μ)(β + d β ) 0 0 10 1 20 * * I = d E * * Q = d d E 2 1 d d k * 1 2 1 * R = E d = k + α + μ 0 0 d = k + α + μ 1 0 From above expressions, we can observe that P is positive only when R > 1. Thus the endemic equilibrium point P exists whenever R > 1. Theorem 3.5. The endemic equilibrium point (EE) P (S , E , I , Q , R ) is locally asymptotically stable. * * * * * * Proof. The Jacobian matrix at EE is given by 0 1 * * * * −β E − β I − μ −β S −β S 0 0 10 20 10 20 B C * * * * β E + β I β S − ( + μ) β S 0 0 B 10 20 10 0 20 C B C J = B 0 −(k + α + μ) 0 0 C EE 0 0 0 B C @ A 0 0 k −(k + α + μ) 0 0 1 0 0 0 0 k −μ Solving the above for det(J − λI) = 0, we obtain the following characteristic equation EE 3 2 (−μ − λ)(−k − α − μ − λ)(λ + Aλ + Bλ + C) = 0, (4) 1 0 * * where A = −(a + d + h), B = ad + ah + dh − bc − eg, C = aeg + bch − adh − e, with a = −β E − β I − μ, 10 20 * * * * * * b = β E + β I , c = −β S , d = β S − ( + μ), e = , f = −β S , g = β S and h = −(k + α + μ). 10 20 10 10 0 0 20 20 0 0 Since roots of the characteristic equation (4) corresponding to Jacobian matrix evaluated at P are−μ < 0, 3 2 −(k + α + μ) < 0, and remaining roots satisfy the characteristic equation λ + Aλ + Bλ + C = 0. Applying 1 0 Routh-Hurwitz criterion, the endemic equilibrium P is locally asymptotically stable, provided A > 0, C > 0 and AB − C > 0. This completes the proof. 4 Parameter estimation, model validation and prediction In this section, we perform numerical simulations to support our analytical results to observe the dynamics of the proposed nonlinear system. In order to do so, we rst estimate the parameters by analyzing the present 80 Ë Swati Tyagi et al. ×10 Predicted cases 12 Actual cases 0 50 100 150 200 250 300 350 Day (after 1 may 2020) Figure 1: Fitting of model for prediction and validity for the COVID-19 outbreak in Russia state in the country. In this work ‘fmincon’ function of MATLAB is used as an optimization function to perform numerical simulations. It is used to optimize the unknown parameters of the model equations. The format of the fmincon function is as follows: [b, fmin,ag] = fmincon(@fun,b ,[],[],[],[],lb,ub,options); Here, ‘fun’ is the objective function. In this study, the objective function is N = S + E + I + Q + R; b refers to the initial state of parameters. The initial state of parameters are chosen randomly. The terms lb and ub represent the lower bound and upper bounds of the parameters that are chosen randomly. In the options set, the maximum number of iterations required for attaining global optimum are specied. Case study: Russia We rst check for the validity of model and predictions for COVID-19 in Russia. The model parameters have been estimated considering the total infectious cases from May 1 to September 30 as provided by [2] using minimization techniques. We have considered the initial population size as on 1 May as S(0) = 140000000; E(0) = 140000; I(0) = 106498; Q(0)=100000; and R(0) = 13220. The parameter values are given in the following Table: Parameter Value β 0.020 β 0.254 0.456 k 0.645 α 1.292 A 100000.649 μ 0.007 k 0.650 Cases Analysis of infectious disease transmission and prediction through SEIQR epidemic model Ë 81 ×10 Predicted cases Actual Cases 0 20 40 60 80 100 120 140 160 180 200 Day (after 27 June 2020) Figure 2: Fitting of model for prediction and validity for the COVID-19 outbreak in India In Figure 1, we provide the time-series analysis of total cases (infected, quarantined and recovered). It is clear from the gure that disease will possibly grow above 1300000 by end of November 2020 in Russia, if all the conditions remain same with a controllable number of infections. Case Study: India We further perform rigorous numerical simulation to analyze the scenario of COVID-19 epidemic in India. We rst estimate the model parameters using the total infectious cases from June 27 to September 30 as provided by [2] using minimization techniques. We have considered the initial population size as S(0) = 130000000; E(0) = 100000; I(0) = 508953; Q(0)=100000; and R(0) = 295881. The estimated parameter values are given in the following Table: Parameter Value β 0.020 β 0.001 0.023 k 0.650 α 0.020 A 100000.650 μ 0.001 k 0.650 In Figure 2, we have presented the best tted curve for the total number of COVID-19 and the same of real cases [2] of the best tted curve and the model predictions. It can be observed from the model predictions that the disease may continue upto even after end of November, if all the conditions remain same and the total infections may reach 11000000. Cases 82 Ë Swati Tyagi et al. 5 Predictions using Long short term memory (LSTM) model Estimation of future COVID-19 cases can be viewed as a time series analysis problem where past trends in the increase of COVID-19 cases are used to predict future trend. The time series to be used for prediction of COVID- t+1 19 in future can be formulated as given in Eq. 5. In Eq. 5, let ˆx denote the predicted number of COVID-19 cases at time instant t + 1. The number of cases are predicted as a function f (.) of the cases at previous time steps. In this Eq. p represents the past time lags under consideration for the study and the terms x represent actual cases at any time instant t. t+1 t t−1 t−p+1 ˆx = f (x , x , . . . , x ) (5) Multi-step ahead prediction includes predicting two or more values in future, In this study, they are generated using an iterative method [9] where the predicted value is repeatedly passed back to the model to generate t+2 the next step forecast. For example ˆx can be predicted as: t+2 t t t+1−p−1 x = f (ˆx + 1, x , . . . , x ) (6) t+M Using the iterative process given in Eq. 6, the multi-step ahead prediction ˆx can be obtained. An M step ahead forecast requires M one-step ahead forecasts. The term M is called prediction horizon. . The nature of the function f(.) in Eq. 5 can be linear or non-linear. There are dierent models in literature for learning this linear or non-linear function for learning patterns in data. Autoregressive integrated moving average (ARIMA) method uses a linear function for prediction [10]. However it is seen in dierent predictive analytic applications such as future resource workload prediction [11] that non-linear models perform better in comparison to linear models. This is due to the fact that most of the real life patterns are non-linear in nature. Neural networks overcome the limitation of ARIMA models and can learn non-linear relationships in past time steps and predicted values [12]. A neural network uses non-linear activation function in its layers which enables to learn non-linear relations. Dierent works in literature have also used hybrid of ARIMA and neural networks for better modelling of patterns in a time series [13, 14]. Recently, deep neural networks are well favoured over conventional neural networks as they have shown promising results. Among deep neural networks, recurrent neural network and its variation long short term memory model (LSTM) is very popular. LSTM has shown its ecacy in dealing with short-term and long-term dependency problems. It has been successfully applied in dierent domains for modelling sequential and temporal dependencies. For exam- ple, LSTM is used for estimating future resource usage prediction in cloud data centres in [15, 16], resource contention failure identication [17], trac ow prediction [18] and many other use cases. These applications from literature indicate that LSTM is playing signicant role for modelling temporal and sequential depen- dencies to predict COVID-19 patterns. The architecture of LSTM model is explained in next subsection. 5.1 Architecture of LSTM model In this section, we provide basic introduction to use LSTM for time series prediction. LSTM is a variant of recurrent neural network (RNN). RNNs are used for modelling sequential dependencies in data. A basic RNN consists of neurons which are arranged in the form of an array in dierent layers of the RNN. Each neuron in a given layer is connected to every neuron in the next layer. In these networks, the output from the previous time step is passed as input to the next time step which allows the network to retain memory and hence enables the model to learn sequential and temporal dependencies in data [17]. The prediction is computed as: t+1 t ˆx = W h + b (7) hy The terms W and b represent the weights and bias respectively between output layer and last RNN layer hy whose output is h . The network is trained using a mean squared error loss function. LSTM model is better than basic RNN as unlike RNN it does not have vanishing gradient problem. The time series prediction architecture of LSTM is similar to RNN except that a neuron of basic RNN model is Analysis of infectious disease transmission and prediction through SEIQR epidemic model Ë 83 replaced by a cell in LSTM. Figure 3 shows the architecture of basic LSTM cell. Each cell has three gates namely, input gate, output gate and the forget gate. The forget gate identies the information that is to be forgotten by the cell memory. The term ˜c denotes the vector of new candidate values. The role of the input gate is to decide the values that are to be updated and the output gate controls the ow of information out of the cell. Together these three gates control the ow of information in and out of cell memory enabling it learn long and short term dependencies [17]. The rectangular boxes in Figure 3 represent the layer of neurons with a desired activation function. Estimated COVID-19 cases can be predicted by analysing its temporal dependencies with other measur- able metrics as well. We explain this prediction model in the next subsection. 5.2 Univariate and multivariate modelling of COVID-19 cases Univariate prediction models analyse temporal patterns in the past history of the number of COVID-19 cases itself to predict future trends. However it is seen in real life scenarios that each decision involves analysing multiple variables together. This is due to the fact that studying multiple related variables together can aid in the understanding of variations in the trend of desired metric to be predicted [19]. Multivariate analysis in- volves analysing multiple variables together to analyse the relationships between them and their combined eect on the prediction of desired metric. It is shown in the context of prediction of resource usage and per- formance in cloud data centres [15, 20] and prediction of mortality rates [21] that better prediction capabilities in forecasting new values are achieved by analysing metrics which are related with the desired metric. This indicates that multivariate analysis provides better statistical power to the model and aids in achieving more realistic predictions for the desired metric. Therefore, we analyse both univariate as well as multivariate pre- diction models. Here the desired metric is COVID-19 cases. To build multivariate prediction model, temporal dependencies in deaths as well as recovered cases are also analysed in addition to the number of COVID-19 cases. In the next section, we analyse the results of LSTM based model for estimating COVID-19 cases for India and Russia. 6 Experimental analysis using LSTM model In this section, we present the results of the experimental studies of COVID-19 cases. Here, we rst analyse the nature of COVID-19 cases in the next subsection which is followed by its modelling and prediction by LSTM models in subsequent sections. The data of COVID-19 cases, deaths and recovered cases in India and Russia is obtained from [2]. t−1 tanh t tt f oo forget gate output gate t t ˜c i t−1 t t−1 t input gate tanh h x h x t−1 t t−1 t h x h x Figure 3: Structure of LSTM cell. 1 84 Ë Swati Tyagi et al. 6.1 Pre-processing We analysed the time varying trends of COVID-19 cases in India and Russia. It is seen that time series of COVID-19 cases is non-stationary in nature. This is due to the reason that the series is a monotonically in- creasing series and hence the statistical properties of the series are not constant. We also experimentally veried this fact by using augmented Dickey-Fuller test [22]. Since the series is non-stationary, we use Box- Cox transformation to pre-process the series. Many machine learning and deep learning models are based on the assumption of normality in a time series. Therefore, here Box-Cox transformation is used so that the transformed series is closer to the normal distribution where the mean, mode and median lie at the centre. This pre-processing helps LSTM models to learn the trends in the non-stationary time series better. 6.2 Model specications and results To train our model we have used univariate as well as multivariate data for model training and prediction. To train the LSTM based neural network, we have used tanh(·) activation function. It is seen empirically that for India, a two layered LSTM with 4 and 8 cells in individual layers performs the best. For Russia, a two layered LSTM network with 32 cells in each layer performs the best. The model is trained for the data from May 1 to September 30 for Russia and from June 27 to September 30 for India and then predictions are made using the test results. We have predicted the total COVID-19 cases 60 days ahead for both India and Russia. The obtained prediction results for India are given in Fig. 4(a)-4(b) and for Russia are given in Fig. 5(a)-5(b). 1e7 1e7 1.2 Actual Actual 1.75 Predicted Predicted 1.0 1.50 1.25 0.8 1.00 0.6 0.75 0.4 0.50 0.2 0.25 0.00 0.0 0 50 100 150 200 0 50 100 150 200 Day (after 1 May, 2020) Day (after 1 May, 2020) (a) For univariate data (b) For multivariate data Figure 4: Model training and Prediction analysis for India 7 Comparative discussion of results From the predictions and obtained results in Section 4 and 6, we can observe that the virus infection will increase rapidly with time resulting into a vast infected population and the situation may worsen, if no nec- essary steps are taken or vaccine is not invented. Thus it is high time to take necessary precautions and preven- tive measures. People should themselves need to be aware and follow strict precautionary measures. There should be minimum public mobility along with strict hygienic and social distancing norms. Implementing Government restrictions mathematically, we can say that the total cases may peak upto 1300000 for Russia Cases Cases Analysis of infectious disease transmission and prediction through SEIQR epidemic model Ë 85 1e6 1e6 Actual Actual 1.4 1.6 Predicted Predicted 1.2 1.4 1.2 1.0 1.0 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 50 100 150 200 0 50 100 150 200 Day (after 1 May, 2020) Day (after 1 May, 2020) (a) For univariate data (b) For multivariate data Figure 5: Model training and Prediction analysis for Russia and 11000000 for India by end of November. We can also observe that LSTM is giving more precise predic- tions. Since, Government is changing its restriction measures and policies from time to time, so the system parameters are also changing after that period of time. Therefore, in the current position we are not interested about infected population for long time of period, but through these estimates, we can say that we need to follow strict policies to reduce the infection. 8 Conclusion In this work, we have formulated and proposed a SEIQR epidemic model to study the dynamics of spread of infectious diseases. The proposed model is an attempt to provide a deep insight to analyze the disease spread and design possible control strategies. We have derived various analytical results, such as positivity and boundendess of the solutions, local stability of equilibria, basic reproduction number. There exist two equilibria, namely the disease-free equilibrium (DFE) and the endemic equilibrium (EE). We have estimated the parameters using the real world data as a special case for COVID-19. Our main aim is to predict mathematically the number of infected individual due COVID- 19 virus in Russia and India in this study. To full our aim, we have performed a detailed numerical simulation of our proposed model through minimization techniques and long short term memory networks. We have tried to establish an acceptable agreement between the data analysis and numerical solutions. Analysing the results of our proposed model, we can say that India may be in a big trouble due to COVID-19 virus in near future. Indian Government should take stricter measures and policies other than lock down, quarantine etc. Since the government policies are changing continuously to protect India, thus parameters are also changing with time. Thus, we can assume that if the Indian Government and all state government may take proper steps time to time and people also follow strict social distancing measures along with proper hygiene, yoga, and etc., then the infected number of population may dier from our predictions with time and India will recover in recent future from this pandemic. At the end, we can summarize and conclude that mathematical modelling is an ecient method to estimate the situation and dynamics of this global pandemic COVID-19, if the parameters are estimated properly. Due to limited data and short time span, we have made certain assumptions while framing model. But with the release of more epidemic data, the key parameters may undergo signicant changes on account of health, education, religion, per capita income, inuencing the pandemic spread among the masses. Cases Cases 86 Ë Swati Tyagi et al. AcknowledgementsL=: The authors are thankful to the reviewers for their constructive comments and suggestions, which have helped in improving the manuscript considerably. Conict of interest: On behalf of all authors, the corresponding author states that there is no conict of interest. Data Availability Statement: Data sharing is not applicable to this article as no datasets were generated or analysed during the current study. References [1] https://www.who.int/emergencies/diseases/novel-coronavirus-2019 [2] https://www.worldometers.info/coronavirus [3] Kamrujjaman, M., Mahmud, M. S., & Islam, M. S. (2020). Coronavirus outbreak and the mathematical growth map of Covid- 19. Annual Research & Review in Biology, 72-78. [4] Verdière, N., Lanza, V., Charrier, R., Provitolo, D., Dubos-Paillard, E., Bertelle, C., & Alaoui, A. (2014, June). Mathematical modeling of human behaviors during catastrophic events. [5] Ivorra, B., & Ramos, A. M. (2020). Application of the Be-CoDiS mathematical model to forecast the international spread of the 2019–20 Wuhan coronavirus outbreak. Research Gate Preprint, 9, 1-13. [6] Diekmann, O., Heesterbeek, J. A. P., & Roberts, M. G. (2010). The construction of next-generation matrices for compartmental epidemic models. Journal of the Royal Society Interface, 7(47), 873-885. [7] Li, M. Y., & Muldowney, J. S. (1996). A geometric approach to global-stability problems. SIAM Journal on Mathematical Anal- ysis, 27(4), 1070-1083. [8] Srivastav, A. K., & Ghosh, M. (2016). Modeling and analysis of the symptomatic and asymptomatic infections of swine flu with optimal control. Modeling Earth Systems and Environment, 2(4), 1-9. [9] Zhang, L., Zhou, W. D., Chang, P. C., Yang, J. W., & Li, F. Z. (2013). Iterated time series prediction with multiple support vector regression models. Neurocomputing, 99, 411-422. [10] Caglar, F., & Gokhale, A. (2014, June). iOverbook: intelligent resource-overbooking to support soft real-time applications in the cloud. In 2014 IEEE 7th International Conference on Cloud Computing (pp. 538-545). [11] Gupta, S., & Dinesh, D. A. (2017). Online adaptation models for resource usage prediction in cloud network. In 2017 Twenty- third National Conference on Communications (NCC) (pp. 1-6). [12] Gupta, S., Dileep, A. D., & Gonsalves, T. A. (2016). Fractional dierence based hybrid model for resource prediction in cloud network. In Proceedings of the Fifth International Conference on Network, Communication and Computing (pp. 93-97). [13] Zhang, G. P. (2003). Time series forecasting using a hybrid ARIMA and neural network model. Neurocomputing, 50, 159-175. [14] Chaâbane, N. (2014). A hybrid ARFIMA and neural network model for electricity price prediction. International journal of electrical power & energy systems, 55, 187-194. [15] Gupta, S., Dileep, A. D., & Gonsalves, T. A. (2018). A joint feature selection framework for multivariate resource usage pre- diction in cloud servers using stability and prediction performance. The Journal of Supercomputing, 74(11), 6033-6068. [16] Song, B., Yu, Y., Zhou, Y., Wang, Z., & Du, S. (2018). Host load prediction with long short-term memory in cloud computing. The Journal of Supercomputing, 74(12), 6554-6568. [17] Gupta, S., & Dileep, A. D. (2020). Relevance feedback based online learning model for resource bottleneck prediction in cloud servers. Neurocomputing. Vol. 402, 307-322. [18] Zhao, F., Zeng, G. Q., & Lu, K. D. (2019). EnLSTM-WPEO: short-term trac flow prediction by ensemble LSTM, NNCT weight integration, and population extremal optimization. IEEE Transactions on Vehicular Technology, 69(1), 101-113. [19] Chakraborty, K., Mehrotra, K., Mohan, C. K., & Ranka, S. (1992). Forecasting the behavior of multivariate time series using neural networks. Neural networks, 5(6), 961-970. [20] Hirwa J.S. & Cao J (2014) An ensemble multivariate model for resource performance prediction in the cloud. In: Network and Parallel Computing NPC 2014, Lecture Notes in Computer Science, Vol. 8707, pp 333–346. [21] Peña, D., & Sãnchez, I. (2007). Measuring the advantages of multivariate vs. univariate forecasts. Journal of Time Series Analysis, 28(6), 886-909. [22] Harris, R. I. (1992). Testing for unit roots using the augmented Dickey-Fuller test: Some issues relating to the size, power and the lag structure of the test. Economics letters, 38(4), 381-386.
Nonautonomous and Stochastic Dynamical Systems – de Gruyter
Published: Jan 1, 2021
Keywords: Infectious disease; Mathematical model; Stability analysis; Long Short Term Memory networks (LSTM); Parameter estimation; 00A71; 92B20; 34D20
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
Access the full text.
Sign up today, get DeepDyve free for 14 days.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.