
The spreading of COVID-19 has been considered a worldwide issue, and many global efforts have been suggested. Suggested control strategies to minimize the impact of the disease have effectively worked with computational simulations and mathematical models. Model critical transmissions and sensitivities are also key elements to study this pandemic more widely. This work reviews and discusses susceptible–exposed–infected–recovered (SEIR) model to predict the spreading of this disease. Accordingly, the basic reproduction number and its parameter elasticity are considered at the equilibrium points. Furthermore, the real data of confirmed cases in the Kurdistan region of Iraq are used in estimating model parameters and model validating. Computational model results provide some important model improvements and suggest control strategies. Firstly, the model population states have different model dynamics using the estimated parameters and the initial values. Another result is that almost all model states are sensitive to the model parameters at different levels. Interestingly, contact rate, transition rate from exposed class to the infected class and natural recovery rate are the most important controllable parameters to reduce the basic reproduction number R0, and they become the model critical parameters. More interestingly, computational results for the real data provide that the basic reproduction number in the Kurdistan Region was about 1.28, which is greater than unity. This means that the new coronavirus still has a high potential to spread among individuals, and it will require more interventions and new strategies to control this disease further.
Citation: Sarbaz H. A. Khoshnaw, Kawther Y. Abdulrahman, Arkan N. Mustafa. Identifying key critical model parameters in spreading of COVID-19 pandemic[J]. AIMS Bioengineering, 2022, 9(2): 163-177. doi: 10.3934/bioeng.2022012
[1] | Bashdar A. Salam, Sarbaz H. A. Khoshnaw, Abubakr M. Adarbar, Hedayat M. Sharifi, Azhi S. Mohammed . Model predictions and data fitting can effectively work in spreading COVID-19 pandemic. AIMS Bioengineering, 2022, 9(2): 197-212. doi: 10.3934/bioeng.2022014 |
[2] | Sarbaz H. A. Khoshnaw, Azhi Sabir Mohammed . Computational simulations of the effects of social distancing interventions on the COVID-19 pandemic. AIMS Bioengineering, 2022, 9(3): 239-251. doi: 10.3934/bioeng.2022016 |
[3] | Mehmet Yavuz, Waled Yavız Ahmed Haydar . A new mathematical modelling and parameter estimation of COVID-19: a case study in Iraq. AIMS Bioengineering, 2022, 9(4): 420-446. doi: 10.3934/bioeng.2022030 |
[4] | Fırat Evirgen, Fatma Özköse, Mehmet Yavuz, Necati Özdemir . Real data-based optimal control strategies for assessing the impact of the Omicron variant on heart attacks. AIMS Bioengineering, 2023, 10(3): 218-239. doi: 10.3934/bioeng.2023015 |
[5] | Ayub Ahmed, Bashdar Salam, Mahmud Mohammad, Ali Akgül, Sarbaz H. A. Khoshnaw . Analysis coronavirus disease (COVID-19) model using numerical approaches and logistic model. AIMS Bioengineering, 2020, 7(3): 130-146. doi: 10.3934/bioeng.2020013 |
[6] | Prasantha Bharathi Dhandapani, Dumitru Baleanu, Jayakumar Thippan, Vinoth Sivakumar . On stiff, fuzzy IRD-14 day average transmission model of COVID-19 pandemic disease. AIMS Bioengineering, 2020, 7(4): 208-223. doi: 10.3934/bioeng.2020018 |
[7] | Júlia Novaes Matias, Gyovanna Sorrentino dos Santos Campanari, Gabriela Achete de Souza, Vinícius Marinho Lima, Ricardo José Tofano, Claudia Rucco Penteado Detregiachi, Sandra M. Barbalho . Metabolic syndrome and COVID-19. AIMS Bioengineering, 2020, 7(4): 242-253. doi: 10.3934/bioeng.2020021 |
[8] | Kranthi Kumar Lella, Alphonse PJA . A literature review on COVID-19 disease diagnosis from respiratory sound data. AIMS Bioengineering, 2021, 8(2): 140-153. doi: 10.3934/bioeng.2021013 |
[9] | Mariem Jelassi, Kayode Oshinubi, Mustapha Rachdi, Jacques Demongeot . Epidemic dynamics on social interaction networks. AIMS Bioengineering, 2022, 9(4): 348-361. doi: 10.3934/bioeng.2022025 |
[10] | Honar J. Hamad, Sarbaz H. A. Khoshnaw, Muhammad Shahzad . Model analysis for an HIV infectious disease using elasticity and sensitivity techniques. AIMS Bioengineering, 2024, 11(3): 281-300. doi: 10.3934/bioeng.2024015 |
The spreading of COVID-19 has been considered a worldwide issue, and many global efforts have been suggested. Suggested control strategies to minimize the impact of the disease have effectively worked with computational simulations and mathematical models. Model critical transmissions and sensitivities are also key elements to study this pandemic more widely. This work reviews and discusses susceptible–exposed–infected–recovered (SEIR) model to predict the spreading of this disease. Accordingly, the basic reproduction number and its parameter elasticity are considered at the equilibrium points. Furthermore, the real data of confirmed cases in the Kurdistan region of Iraq are used in estimating model parameters and model validating. Computational model results provide some important model improvements and suggest control strategies. Firstly, the model population states have different model dynamics using the estimated parameters and the initial values. Another result is that almost all model states are sensitive to the model parameters at different levels. Interestingly, contact rate, transition rate from exposed class to the infected class and natural recovery rate are the most important controllable parameters to reduce the basic reproduction number R0, and they become the model critical parameters. More interestingly, computational results for the real data provide that the basic reproduction number in the Kurdistan Region was about 1.28, which is greater than unity. This means that the new coronavirus still has a high potential to spread among individuals, and it will require more interventions and new strategies to control this disease further.
Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The virus was first identified in Wuhan, China, in December 2019, and it quickly spread to other parts of mainland China and many other countries [1]–[6]. Due to the rapid spread of the virus with consequences worldwide, the World Health Organization (WHO) declared a public health emergency of international concern on January 30, 2020, and a pandemic on March 11, 2020 [7],[8]. The transmission of COVID-19 has become a global health threat. As of April 17, 2020, there was a total of 2,230,439 COVID-19 cases reported worldwide, with 150,810 deaths and 564,210 recovered cases [9]. According to the WHO's COVID-19 global report, over 21.2 million people had been infected with the virus as of August 16, 2020, with over 761,779 deaths [10]. Moreover, till October 14, 2021, more than 239,426,728 cases of infection and at least 4,878,373 deaths have been confirmed in Africa, Asia, Europe, North America, Oceania, South America and Antarctica. This has become one of the deadliest pandemics in history according to the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU) [11].
Studying and analyzing complex infectious disease models are difficult tasks and require some computational tools and mathematical approaches. In systems biology, mathematical presentations and computational simulations are effective tools and will help further understanding and future predictions about model dynamics. The models of COVID-19 are complicated, and they are required some mathematical tools to improve interventions and healthcare strategies. Recently, several studies about COVID-19 were suggested based on mathematical tools and computational simulations to understand this pandemic and its control approaches. In [12]–[14], the SEIR model was modified by several researchers to understand the impacts of social distancing, rapid testing, vaccination and contaminated objects in controlling the spread of COVID-19. One of the most important key elements in the spreading of COVID-19 is the basic reproduction number (also called basic reproductive ratio). This quantity, denoted by R0, represents the number of secondary infections induced by a single primary infection in a population where everyone is susceptible [15]. According to developed techniques of mathematical epidemiology models, R0 plays an important role for understanding epidemiological ideas and identifying key critical parameters for such models, and R0 may have different values for reported cases of coronavirus. According to a published study, the basic reproduction number was computed using incidence data in the city of Jakarta from March 3, 2020, to May 10, 2020. The median value of R0 was about 1.75 during the first interval. After the rigorous social distancing implemented, this value was about 1.22. In addition, the basic reproduction number has been calculated for different infectious disease models. Consequently, this quantity was calculated for a mathematical model of the CD8+ T-cell response to human T cell leukemia/lymphoma virus type I (HTLV-I) infection, given in [16].
According to published studies given in [14],[17],[18], the concept of sensitivity analysis has been presented to detect model sensitive parameters, including local sensitivities for each model state with respect to model parameters. Such sensitivities are computed by using three different methods: non-normalization, half-normalization and full-normalization. Accordingly, there was another method of sensitivity analysis given in [19]. They suggested the partial rank correlation coefficient (PRCC) technique to determine how the glioma-immune model output is affected by changes in specific parameter values.
Optimal control analysis is another effective method to identify the model critical parameters and to suggest control strategies. This technique has been used with the COVID-19 pandemic for deriving control policies. Recently, this technique has been used to find some control strategies for the COVID-19 pandemic in Indonesia, given in [20]. Then, there was a study about mathematical modeling with optimal control analysis for the spread of the COVID-19 pandemic with a case study of Brazil and cost-effectiveness; this was given in [21]. More recently, there was another approach that suggested a mathematical model with optimal control to highlight the negative impact of quarantine on diabetic people, with cost-effectiveness, presented in [22]. The reader can also find more relevant works about different approaches and mathematical models of the COVID-19 pandemic and its control strategies in [23]–[29].
Although there are many mathematical models that have been suggested for prediction of the COVID-19 pandemic, such models can be further improved and analyzed in terms of real data. Calculating the basic reproduction number with local sensitivities can improve the model outcomes further. An issue that has not been considered enough is elasticity of basic reproduction number with respect to model parameters. This helps us to identify the most critical model parameters in spreading this disease. Another issue of modeling such phenomena is validating the model results. Therefore, comparing the model predictions with real data can also indicate the accuracy and reliability of such models.
Estimating the model initial parameters for the suggested SEIR model is an important step forward in this study compared to the other existing studies about the confirmed cases of COVID-19 in the Kurdistan Region of Iraq. Therefore, we suggest the Metropolis–Hastings (MH) algorithm to estimate such model parameters, and then we use them in computational simulations. Such estimated parameters can also provide a good agreement between the model results and the daily infected cases. Another novelty here is the data collection of COVID-19 infected cases in the area that provides the future predictions and interventions. Accordingly, we have collected the daily confirmed cases of COVID-19 from January 10, 2021, to March 22, 2021. The collected real data here give us a great understanding about this disease in the area and suggested more interventions and strategies to control this pandemic. Interestingly, identifying the critical model parameters based on sensitivity and elasticity methods becomes another novelty for the suggested SEIR model.
This paper has some contributions. The main contribution is suggesting the simplified SEIR model for COVID-19 disease to predict the spreading of the disease in the Kurdistan Region of Iraq. We use this model to compare the model results and the daily infected cases more easily and correctly. Another contribution here is identifying the model critical parameters. Consequently, we suggest the local sensitivity methods to calculate the model sensitivities between the model variables and parameters. Additionally, the elasticity between the basic reproduction number, R0, and model parameters are also calculated and helps to identify the key critical parameters in spreading this disease.
The organization of this paper is as follows. In Section 2, we construct the SEIR model based on our assumptions and model transmissions. Using the next generation matrix, we derive the formula of the basic reproduction number, calculate the elasticity between R0 and model parameters and show the model local sensitivities based on the estimated values, as discussed in Section 3. We perform our computational model simulations of daily infected cases in Section 4. Finally, some conclusions and suggestions are given in Section 5.
Many infectious disease models have been suggested to understand the transmission rates of COVID-19. The SEIR model has been commonly applied in many infectious diseases. This model consists of four compartments, namely, susceptible S(t), exposed E(t), infected I(t) and recovered R(t) individuals. Several authors modified the SEIR model for studying the transmission of pathogens. Among them, in [30], a modified SEIR model was used to describe the 2009 influenza A (H1N1) pandemic within local regions in Japan. It was also used for the Ebola epidemic occurring in the West African nations of Guinea, Liberia and Sierra Leone [31]. In this paper, we use the SEIR model for COVID-19 disease that describes the spreading of COVID-19 in the Kurdistan Region of Iraq. Accordingly, we develop the previously published work given in [32], and we add a parameter that shows a rate of recovery of the individuals from the exposed class. The model network of SEIR shown in Figure 1 includes various compartments with their transmissions. The model compartments and transmission rates are defined in Table 1.
The chemical reactions of the model can be expressed as follows:
Furthermore, all parameters and initial states with their biological meanings are shown in Table 1. We use the concepts of standard chemical kinetics and the mass action law to derive a system of non–linear ordinary differential equations [33],[34]. Then, the system of nonlinear ordinary differential equations takes the following form:
with non-negative initial conditions S(0) = S0 > 0, E(0) = E0 > 0, I(0) = I0 > 0, R(0) = R0 ≥ 0. All parameters are defined as contact rates between individual categories.
Symbol | Biological definition |
S(0) | Initial susceptible individuals |
E(0) | Initial exposed individuals |
I(0) | Initial symptomatic infected individuals |
R(0) | Initial recovered individuals |
β | Infection rate |
Λ | Recruitment rate |
k | Transition rate of exposed individuals to the infected class |
α | Transition rate of exposed individuals to the recovered class |
γ | Natural recovery rate |
µ | Natural death rate |
Infectious disease models cannot be well understood only by biological techniques. Model analysis techniques and computational simulations can serve a great role in predicting and analyzing the model dynamics and identifying critical model parameters. Recently, some techniques of model analysis have been applied such as sensitivity and elasticity methods. They are useful tools to determine which variable or parameter is sensitive to a specific condition [12],[17],[18]. Calculating the basic reproduction numbers at the equilibrium points can also provide a great environment to discuss infectious disease models more widely and accurately [13],[14]. Here, in this study, we focus on three important techniques of model analysis: basic reproduction numbers, local sensitivity and elasticity methods. More details about the suggested techniques are discussed in the following sections.
In epidemiology, the next-generation matrix is used to derive the basic reproduction number for a compartmental model of the spread of infectious disease. This method was proposed by Diekmann et al. (1990) [35] and van den Driessche and Watmough (2002) [36]. One of the most important parameters for epidemiological models is the basic reproduction number [37]. Using Equation 2, we can easily derive R0 at an equilibrium point E0. The suggested model here includes two infected states, E and I, and two uninfected states, S and R. Let X be a vector of infected classes and Y be a vector of uninfected classes. They are shown below:
To find equilibrium points for the model, let us set the right-hand side of Equation 2 equal to zero. Then, the equilibrium points are obtained. The first equilibrium is the COVID-19-free equilibrium point, which is represented by
The second type of equilibrium is called the endemic equilibrium point. It is given below:
where
Using the next generation method, we have the vector functions
Then, the Jacobian matrices for M, G at E0 are given below:
Since W is a non-singular matrix, W−1 can be calculated:
Therefore, the basic reproduction number is the spectral radius of HW−1. Thus, the basic reproduction number for Equation 2 is
The concept of elasticity for infectious disease models has an effective role in determining the sensitivity between R0 and the model parameters. For calculating such coefficients, the following main formula is used:
where ω is the set of all parameters. Using Equation 5, we can calculate the elasticity index of R0 for each model parameter. They are given below:
Using the estimated parameters, we can compute the model elasticity; see Table 2.
Parameter | |
Λ | 1 |
k | 0.14039612 |
α | −0.403487536 |
γ | −0.9996225141 |
β | 1 |
µ | −1.0004274882 |
The results given in Table 2 show the positive and negative signs for the model elasticity. These give us the direct or indirect relations between R0 and the model parameters. For example, this coefficient is positive for the set of parameters {Λ, k, β}. This means that when these parameters, {Λ, k, β}, are increased, the value of R0 is also increased, and this virus spreads further. On the other hand, the signs of {γ, α, µ} are negative. This means that the value of the basic reproduction number can be reduced by increasing such parameters. Interestingly, β is the most significant positive parameter that can be utilized to reduce R0. For example, if the value of β is decreased by 10%, then the basic reproduction number will be reduced to R0 = 0.52886708. Furthermore, reducing the contact rate between individuals can impact the number of infected people from the disease.
The concept of local sensitivity analysis is another step forward for studying infection disease models and identifying the model critical elements. This is simply used for determining the sensitivity of each model species with respect to the model parameters. This method was recently considered in [14],[17],[18],[38]. Assume that there are p compartments (yi for i = 1,2,...,p) and q parameters (bj for j = 1,2,...,q) in an infectious disease model. The model balanced equations are represented by the following system of differential equations:
where
The full-normalization sensitivities are given as follows:
The half-normalization sensitivities are defined as follows:
The non-normalization sensitivities can be expressed as follows:
where
Computational techniques are effective tools and frequently required to analyze some complex infectious disease models. The estimated parameter values and initial populations are obtained from the real data in the Kurdistan Region of Iraq given in the Appendices. The model estimated parameters are α = 0.23, k = 0.34, γ = 0.08, β = 0.17, µ = 0.00003, Λ = 182.46, and the model initial variables are S0 = 6 × 106, E0 = 105.26, I0 = 1027, R0 = 1027. Using System Biology Toolbox (SBedit) in MATLAB, we determine the numerical approximate solutions of the model (Equation 2); see Figures 5 and 6. Numerical simulations are calculated in two-dimensional planes. All model computational results here are very interesting and provide some biological aspects. Firstly, it can be seen that Figure 5 shows the model dynamics of susceptible, exposed, asymptomatic infected and recovered people. We have noticed that the number of healthy people decreases gradually and becomes stable after 500 days, while the dynamics of recovered people increases and then becomes more flat after 500 days as well. Then, from 200 to 600 days, the numbers of exposed and infected individuals change dramatically. More interestingly, Figure 6 shows the comparison of total populations between the real data cases shown in Table 3 and the model equations given in Equation 2. For example, the total populations of infected people for the model equations and real data are given in Figure 6a. Then, Figure 6b shows the total populations of recovered individuals for model equations and real data. The results show good agreement between the real data and the model results. Finally, it is clear that our model results are fitted for the confirmed cases in Kurdistan. This means that the suggested model here can further be used to predict and give more control strategies.
Control strategies and preventions for the COVID-19 pandemic may not be well understood only by healthcare techniques and biological concepts. Mathematical models and computational simulations have been widely used for further improvements and identifying the model critical parameters. Therefore, identifying such parameters becomes an issue that can be further studied, and more techniques can be proposed. Sensitivity of model parameters helps international efforts to suggest more control strategies and preventions. Accordingly, the suggested SEIR model in this study plays an effective role in identifying the model critical elements using elasticity and sensitivity of R0. Thus, the model sensitivities were calculated using three different techniques: non normalizations, half normalizations, and full normalizations. The results showed that the most effective parameters for spreading this disease are contact rate among people, transition rate from exposed class to the infected class and natural recovery rate. Interestingly, identifying the model critical parameters here will help to suggest further interventions and control strategies in the Kurdistan Region of Iraq. The suggested SEIR model can be further developed and can have added more transmission rates between model states. For example, the vaccination compartment with its contact rates can be added to the model. This helps to study this disease more widely and accurately. Further techniques of model analysis such as optimal control will also be applied to the model to find more control strategies.
[1] |
Li Q, Guan X, Wu P, et al. (2020) Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. N Engl J Med 382: 1199-1207. https://doi.org/10.1056/NEJMoa2001316 ![]() |
[2] | Gorbalenya AE, Baker SC, Baric R, et al. Severe acute respiratory syndrome-related coronavirus: The species and its viruses–a statement of the Coronavirus Study Group (2020). https://doi.org/10.1101/2020.02.07.937862 |
[3] |
Chen N, Zhou M, Dong X, et al. (2020) Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study. The Lancet 395: 507-513. https://doi.org/10.1016/S0140-6736(20)30211-7 ![]() |
[4] |
Huang C, Wang Y, Li X, et al. (2020) Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The Lancet 395: 497-506. https://doi.org/10.1016/S0140-6736(20)30183-5 ![]() |
[5] |
Wang C, Horby PW, Hayden FG, et al. (2020) A novel coronavirus outbreak of global health concern. The Lancet 395: 470-473. https://doi.org/10.1016/S0140-6736(20)30185-9 ![]() |
[6] |
Holshue ML, DeBolt C, Lindquist S, et al. (2020) First case of 2019 novel coronavirus in the United States. N Engl J Med 382: 929-936. https://doi.org/10.1056/NEJMoa2001191 ![]() |
[7] | WHO G, Statement on the second meeting of the International Health Regulations.Emergency Committee regarding the outbreak of novel coronavirus (2019-nCoV). World Health Organization, 2020 (2005) . Available from: https://www.who.int/news/item/30-01-2020-statement-on-the-second-meeting-of-the-international-health-regulations-(2005)-emergency-committee-regarding-the-outbreak-of-novel-coronavirus-(2019-ncov) |
[8] | Organization, W.H.WHO Director-General's opening remarks at the media briefing on COVID-19, 2020. Available from: https://reliefweb.int/report/world/who-director-generals-opening-remarks-media-briefing-covid-19-20-november-2020?gclid=CjwKCAjwur-SBhB6EiwA5sKtjlyi9UUQaeFUMljHfMiVkWeSgVVARurASa9mFnEB5GmjY95h2o2dKhoCHw8QAvD_BwE |
[9] | Coronavirus Outbreak [Internet] COVID-19 CORONAVIRUS PANDEMIC, 2020. Available from: https://www.worldometers.info/coronavirus/ |
[10] | World Health Organization, Coronavirus disease (COVID-19) situation report 153, 2020. Available from: https://apps.who.int/iris/handle/10665/332556 |
[11] | COVID J, dashboard by the center for systems science and engineering (CSSE) at Johns Hopkins university (JHU), 2020. Available from: https://coronavirus.jhu.edu/map.html |
[12] |
Aldila D, Samiadji BM, Simorangkir GM, et al. (2021) Impact of early detection and vaccination strategy in COVID-19 eradication program in Jakarta, Indonesia. BMC Res Notes 14: 1-7. https://doi.org/10.1186/s13104-021-05540-9 ![]() |
[13] |
Mekonen KG, Habtemicheal TG, Balcha SF (2021) Modeling the effect of contaminated objects for the transmission dynamics of COVID-19 pandemic with self-protection behavior changes. Results Appl Math 9: 100134. https://doi.org/10.1016/j.rinam.2020.100134 ![]() |
[14] |
Aldila D, Khoshnaw SH, Safitri E, et al. (2020) A mathematical study on the spread of COVID-19 considering social distancing and rapid assessment: The case of Jakarta, Indonesia. Chaos Soliton Fract 139: 110042. https://doi.org/10.1016/j.chaos.2020.110042 ![]() |
[15] |
Perasso A (2018) An introduction to the basic reproduction number in mathematical epidemiology. ESAIM: ProcS 62: 123-138. https://doi.org/10.1051/proc/201862123 ![]() |
[16] |
Khajanchi S, Bera S, Roy TK (2021) Mathematical analysis of the global dynamics of a HTLV-I infection model, considering the role of cytotoxic T-lymphocytes. Math Comput Simulat 180: 354-378. https://doi.org/10.1016/j.matcom.2020.09.009 ![]() |
[17] |
Khoshnaw SH, Salih RH, Sulaimany S (2020) Mathematical modelling for coronavirus disease (COVID-19) in predicting future behaviours and sensitivity analysis. Math Model Nat Pheno 15: 33. https://doi.org/10.1051/mmnp/2020020 ![]() |
[18] |
Khoshnaw SH, Shahzad M, Ali M, et al. (2020) A quantitative and qualitative analysis of the COVID–19 pandemic model. Chaos Soliton Fract 138: 109932. https://doi.org/10.1016/j.chaos.2020.109932 ![]() |
[19] |
Khajanchi S, Nieto JJ (2021) Spatiotemporal dynamics of a glioma immune interaction model. Sci Rep 11: 22385. https://doi.org/10.1038/s41598-021-00985-1 ![]() |
[20] |
Aldila D, Ndii MZ, Samiadji BM (2020) Optimal control on COVID-19 eradication program in Indonesia under the effect of community awareness. Math Biosci Eng 17: 6355-6389. https://doi.org/10.3934/mbe.2020335 ![]() |
[21] |
Kouidere A, Kada D, Balatif O, et al. (2021) Optimal control approach of a mathematical modeling with multiple delays of the negative impact of delays in applying preventive precautions against the spread of the COVID-19 pandemic with a case study of Brazil and cost-effectiveness. Chaos Soliton Fract 142: 110438. https://doi.org/10.1016/j.chaos.2020.110438 ![]() |
[22] |
Kouidere A, Youssoufi LE, Ferjouchia H, et al. (2021) Optimal control of mathematical modeling of the spread of the COVID-19 pandemic with highlighting the negative impact of quarantine on diabetics people with cost-effectiveness. Chaos Soliton Fract 145: 110777. https://doi.org/10.1016/j.chaos.2021.110777 ![]() |
[23] |
Khajanchi S, Sarkar K, Banerjee S (2022) Modeling the dynamics of COVID-19 pandemic with implementation of intervention strategies. Eur Phys J Plus 137: 129. https://doi.org/10.1140/epjp/s13360-022-02347-w ![]() |
[24] | Mondal J, Khajanchi S Mathematical modeling and optimal intervention strategies of the COVID-19 outbreak (2022). https://doi.org/10.1007/s11071-022-07235-7 |
[25] |
Rai RK, Khajanchi S, Tiwari PK, et al. (2022) Impact of social media advertisements on the transmission dynamics of COVID-19 pandemic in India. J Appl Math Comput 68: 19-44. https://doi.org/10.1007/s12190-021-01507-y ![]() |
[26] |
Khajanchi S, Sarkar K, Mondal J, et al. (2021) Mathematical modeling of the COVID-19 pandemic with intervention strategies. Results Phys 25: 104285. https://doi.org/10.1016/j.rinp.2021.104285 ![]() |
[27] |
Khajanchi S, Sarkar K (2020) Forecasting the daily and cumulative number of cases for the COVID-19 pandemic in India. Chaos 30: 071101. https://doi.org/10.1063/5.0016240 ![]() |
[28] |
Samui P, Mondal J, Khajanchi S (2020) A mathematical model for COVID-19 transmission dynamics with a case study of India. Chaos Soliton Fract 140: 110173. https://doi.org/10.1016/j.chaos.2020.110173 ![]() |
[29] |
Sarkar K, Khajanchi S, Nieto JJ (2020) Modeling and forecasting the COVID-19 pandemic in India. Chaos Soliton Fract 139: 110049. https://doi.org/10.1016/j.chaos.2020.110049 ![]() |
[30] |
Saito MM, Imoto S, Yamaguchi R, et al. (2013) Extension and verification of the SEIR model on the 2009 influenza A (H1N1) pandemic in Japan. Math Biosci 246: 47-54. https://doi.org/10.1016/j.mbs.2013.08.009 ![]() |
[31] | Diaz P, Constantine P, Kalmbach K, et al. (2018) A modified SEIR model for the spread of Ebola in Western Africa and metrics for resource allocation. Appl Math Comput 324: 141-155. https://doi.org/10.1016/j.amc.2017.11.039 |
[32] |
Li MY, Graef JR, Wang L, et al. (1999) Global dynamics of a SEIR model with varying total population size. Math Biosci 160: 191-213. https://doi.org/10.1016/S0025-5564(99)00030-9 ![]() |
[33] |
Chellaboina V, Bhat SP, Haddad WM, et al. (2009) Modeling and analysis of mass-action kinetics. IEEE Contr Syst Mag 29: 60-78. https://doi.org/10.1109/MCS.2009.932926 ![]() |
[34] | Khoshnaw SH Dynamic analysis of a predator and prey model with some computational simulations (2017) arXiv: 1709.00053. https://doi.org/10.4172/2329-9533.1000137 |
[35] |
Diekmann O, Heesterbeek JAP, Metz JA (1990) On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J Math Biol 28: 365-382. https://link.springer.com/article/10.1007/BF00178324 ![]() |
[36] |
Van den Driessche P, Watmough J (2002) Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci 180: 29-48. https://doi.org/10.1016/S0025-5564(02)00108-6 ![]() |
[37] |
Heffernan JM, Smith RJ, Wahl LM (2005) Perspectives on the basic reproductive ratio. J R Soc Interface 2: 281-293. https://doi.org/10.1098/rsif.2005.0042 ![]() |
[38] |
Rahman B, Khoshnaw SH, Agaba GO, et al. (2021) How containment can effectively suppress the outbreak of COVID-19: a mathematical modeling. Axioms 10: 204. https://doi.org/10.3390/axioms10030204 ![]() |
[39] |
Biegler LT, Damiano JJ, Blau GE (1986) Nonlinear parameter estimation: A case study comparison. AIChE J 32: 29-45. https://doi.org/10.1002/aic.690320105 ![]() |
[40] |
Varah JM (1982) A spline least squares method for numerical parameter estimation in differential equations. SIAM J Sci Stat Comput 3: 28-46. https://doi.org/10.1137/0903003 ![]() |
[41] | Kurdistan Regional Government, GOV.KRD. Coronavirus(COVID-19). Available from: https://gov.krd/coronavirus-en/situation-update/ |
1. | Sreenitha Kasarapu, Rakibul Hassan, Houman Homayoun, Sai Manoj Pudukotai Dinakarrao, Scalable and Demography-Agnostic Confinement Strategies for COVID-19 Pandemic with Game Theory and Graph Algorithms, 2022, 2, 2673-8112, 767, 10.3390/covid2060058 | |
2. | Alaa Falih Mahdi, Hussein K. Asker, A Sensitivity Analysis of the PSITPS Epidemic Model’s Parameters for COVID-19, 2024, 11, 2518-0010, 79, 10.31642/JoKMC/2018/110210 |
Symbol | Biological definition |
S(0) | Initial susceptible individuals |
E(0) | Initial exposed individuals |
I(0) | Initial symptomatic infected individuals |
R(0) | Initial recovered individuals |
β | Infection rate |
Λ | Recruitment rate |
k | Transition rate of exposed individuals to the infected class |
α | Transition rate of exposed individuals to the recovered class |
γ | Natural recovery rate |
µ | Natural death rate |
Parameter | |
Λ | 1 |
k | 0.14039612 |
α | −0.403487536 |
γ | −0.9996225141 |
β | 1 |
µ | −1.0004274882 |
Symbol | Biological definition |
S(0) | Initial susceptible individuals |
E(0) | Initial exposed individuals |
I(0) | Initial symptomatic infected individuals |
R(0) | Initial recovered individuals |
β | Infection rate |
Λ | Recruitment rate |
k | Transition rate of exposed individuals to the infected class |
α | Transition rate of exposed individuals to the recovered class |
γ | Natural recovery rate |
µ | Natural death rate |
Parameter | |
Λ | 1 |
k | 0.14039612 |
α | −0.403487536 |
γ | −0.9996225141 |
β | 1 |
µ | −1.0004274882 |