
We propose a mathematical model based on a system of differential equations, which incorporates the impact of the chronic health conditions of the host population, to investigate the transmission dynamics of COVID-19. The model divides the total population into two groups, depending on whether they have underlying conditions, and describes the disease transmission both within and between the groups. As an application of this model, we perform a case study for Hamilton County, the fourth-most populous county in the US state of Tennessee and a region with high prevalence of chronic conditions. Our data fitting and simulation results quantify the high risk of COVID-19 for the population group with underlying health conditions. The findings suggest that weakening the disease transmission route between the exposed and susceptible individuals, including the reduction of the between-group contact, would be an effective approach to protect the most vulnerable people in this population group.
Citation: Chayu Yang, Jin Wang. COVID-19 and underlying health conditions: A modeling investigation[J]. Mathematical Biosciences and Engineering, 2021, 18(4): 3790-3812. doi: 10.3934/mbe.2021191
[1] | Honglei Wang, Wenliang Zeng, Xiaoling Huang, Zhaoyang Liu, Yanjing Sun, Lin Zhang . MTTLm6A: A multi-task transfer learning approach for base-resolution mRNA m6A site prediction based on an improved transformer. Mathematical Biosciences and Engineering, 2024, 21(1): 272-299. doi: 10.3934/mbe.2024013 |
[2] | Shanzheng Wang, Xinhui Xie, Chao Li, Jun Jia, Changhong Chen . Integrative network analysis of N6 methylation-related genes reveal potential therapeutic targets for spinal cord injury. Mathematical Biosciences and Engineering, 2021, 18(6): 8174-8187. doi: 10.3934/mbe.2021405 |
[3] | Pingping Sun, Yongbing Chen, Bo Liu, Yanxin Gao, Ye Han, Fei He, Jinchao Ji . DeepMRMP: A new predictor for multiple types of RNA modification sites using deep learning. Mathematical Biosciences and Engineering, 2019, 16(6): 6231-6241. doi: 10.3934/mbe.2019310 |
[4] | Yong Zhu, Zhipeng Jiang, Xiaohui Mo, Bo Zhang, Abdullah Al-Dhelaan, Fahad Al-Dhelaan . A study on the design methodology of TAC3 for edge computing. Mathematical Biosciences and Engineering, 2020, 17(5): 4406-4421. doi: 10.3934/mbe.2020243 |
[5] | Atefeh Afsar, Filipe Martins, Bruno M. P. M. Oliveira, Alberto A. Pinto . A fit of CD4+ T cell immune response to an infection by lymphocytic choriomeningitis virus. Mathematical Biosciences and Engineering, 2019, 16(6): 7009-7021. doi: 10.3934/mbe.2019352 |
[6] | Tamás Tekeli, Attila Dénes, Gergely Röst . Adaptive group testing in a compartmental model of COVID-19*. Mathematical Biosciences and Engineering, 2022, 19(11): 11018-11033. doi: 10.3934/mbe.2022513 |
[7] | Wenli Cheng, Jiajia Jiao . An adversarially consensus model of augmented unlabeled data for cardiac image segmentation (CAU+). Mathematical Biosciences and Engineering, 2023, 20(8): 13521-13541. doi: 10.3934/mbe.2023603 |
[8] | Tongmeng Jiang, Pan Jin, Guoxiu Huang, Shi-Cheng Li . The function of guanylate binding protein 3 (GBP3) in human cancers by pan-cancer bioinformatics. Mathematical Biosciences and Engineering, 2023, 20(5): 9511-9529. doi: 10.3934/mbe.2023418 |
[9] | Xin Yu, Jun Liu, Ruiwen Xie, Mengling Chang, Bichun Xu, Yangqing Zhu, Yuancai Xie, Shengli Yang . Construction of a prognostic model for lung squamous cell carcinoma based on seven N6-methylandenosine-related autophagy genes. Mathematical Biosciences and Engineering, 2021, 18(5): 6709-6723. doi: 10.3934/mbe.2021333 |
[10] | Tahir Rasheed, Faran Nabeel, Muhammad Bilal, Yuping Zhao, Muhammad Adeel, Hafiz. M. N. Iqbal . Aqueous monitoring of toxic mercury through a rhodamine-based fluorescent sensor. Mathematical Biosciences and Engineering, 2019, 16(4): 1861-1873. doi: 10.3934/mbe.2019090 |
We propose a mathematical model based on a system of differential equations, which incorporates the impact of the chronic health conditions of the host population, to investigate the transmission dynamics of COVID-19. The model divides the total population into two groups, depending on whether they have underlying conditions, and describes the disease transmission both within and between the groups. As an application of this model, we perform a case study for Hamilton County, the fourth-most populous county in the US state of Tennessee and a region with high prevalence of chronic conditions. Our data fitting and simulation results quantify the high risk of COVID-19 for the population group with underlying health conditions. The findings suggest that weakening the disease transmission route between the exposed and susceptible individuals, including the reduction of the between-group contact, would be an effective approach to protect the most vulnerable people in this population group.
The constituent members in a system mainly found in nature can be interacting with each other through cooperation and competition. Demonstrations for such systems involve biological species, countries, businesses, and many more. It's very much intriguing to investigate in a comprehensive manner numerous social as well as biological interactions existent in dissimilar species/entities utilizing mathematical modeling. The predation and the competition species are the most famous interactions among all such types of interactions. Importantly, Lotka [1] and Volterra [2] in the 1920s have announced individually the classic equations portraying population dynamics. Such illustrious equations are notably described as predator-prey (PP) equations or Lotka-Volterra (LV) equations. In this structure, PP/LV model represents the most influential model for interacting populations. The interplay between prey and predator together with additional factors has been a prominent topic in mathematical ecology for a long period. Arneodo et al. [3] have established in 1980 that a generalized Lotka-Volterra biological system (GLVBS) would depict chaos phenomena in an ecosystem for some explicitly selected system parameters and initial conditions. Additionally, Samardzija and Greller [4] demonstrated in 1988 that GLVBS would procure chaotic reign from the stabled state via rising fractal torus. LV model was initially developed as a biological concept, yet it is utilized in enormous diversified branches for research [5,6,7,8]. Synchronization essentially is a methodology of having different chaotic systems (non-identical or identical) following exactly a similar trajectory, i.e., the dynamical attributes of the slave system are locked finally into the master system. Specifically, synchronization and control have a wide spectrum for applications in engineering and science, namely, secure communication [9], encryption [10,11], ecological model [12], robotics [13], neural network [14], etc. Recently, numerous types of secure communication approaches have been explored [15,16,17,18] such as chaos modulation [18,19,20,21], chaos shift keying [22,23] and chaos masking [9,17,20,24]. In chaos communication schemes, the typical key idea for transmitting a message through chaotic/hyperchaotic models is that a message signal is nested in the transmitter system/model which originates a chaotic/ disturbed signal. Afterwards, this disturbed signal has been emitted to the receiver through a universal channel. The message signal would finally be recovered by the receiver. A chaotic model has been intrinsically employed both as receiver and transmitter. Consequently, this area of chaotic synchronization & control has sought remarkable considerations among differential research fields.
Most prominently, synchronization theory has been in existence for over 30 years due to the phenomenal research of Pecora and Carroll [25] established in 1990 using drive-response/master-slave/leader-follower configuration. Consequently, many authors and researchers have started introducing and studying numerous control and synchronization methods [9,26,27,28,29,30,31,32,33,34,35,36] etc. to achieve stabilized chaotic systems for possessing stability. In [37], researchers discussed optimal synchronization issues in similar GLVBSs via optimal control methodology. In [38,39], the researchers studied the adaptive control method (ACM) to synchronize chaotic GLVBSs. Also, researchers [40] introduced a combination difference anti-synchronization scheme in similar chaotic GLVBSs via ACM. In addition, authors [41] investigated a combination synchronization scheme to control chaos existing in GLVBSs using active control strategy (ACS). Bai and Lonngren [42] first proposed ACS in 1997 for synchronizing and controlling chaos found in nonlinear dynamical systems. Furthermore, compound synchronization using ACS was first advocated by Sun et al. [43] in 2013. In [44], authors discussed compound difference anti-synchronization scheme in four chaotic systems out of which two chaotic systems are considered as GLVBSs using ACS and ACM along with applications in secure communications of chaos masking type in 2019. Some further research works [45,46] based on ACS have been reported in this direction. The considered chaotic GLVBS offers a generalization that allows higher-order biological terms. As a result, it may be of interest in cases where biological systems experience cataclysmic changes. Unfortunately, some species will be under competitive pressure in the coming years and decades. This work may be comprised as a step toward preserving as many currently living species as possible by using the proposed synchronization approach which is based on master-slave configuration and Lyapunov stability analysis.
In consideration of the aforementioned discussions and observations, our primary focus here is to develop a systematic approach for investigating compound difference anti-synchronization (CDAS) approach in 4 similar chaotic GLVBSs via ACS. The considered ACS is a very efficient yet theoretically rigorous approach for controlling chaos found in GLVBSs. Additionally, in view of widely known Lyapunov stability analysis (LSA) [47], we discuss actively designed biological control law & convergence for synchronization errors to attain CDAS synchronized states.
The major attributes for our proposed research in the present manuscript are:
● The proposed CDAS methodology considers four chaotic GLVBSs.
● It outlines a robust CDAS approach based active controller to achieve compound difference anti-synchronization in discussed GLVBSs & conducts oscillation in synchronization errors along with extremely fast convergence.
● The construction of the active control inputs has been executed in a much simplified fashion utilizing LSA & master-salve/ drive-response configuration.
● The proposed CDAS approach in four identical chaotic GLVBSs of integer order utilizing ACS has not yet been analyzed up to now. This depicts the novelty of our proposed research work.
This manuscript is outlined as follows: Section 2 presents the problem formulation of the CDAS scheme. Section 3 designs comprehensively the CDAS scheme using ACS. Section 4 consists of a few structural characteristics of considered GLVBS on which CDAS is investigated. Furthermore, the proper active controllers having nonlinear terms are designed to achieve the proposed CDAS strategy. Moreover, in view of Lyapunov's stability analysis (LSA), we have examined comprehensively the biological controlling laws for achieving global asymptotical stability of the error dynamics for the discussed model. In Section 5, numerical simulations through MATLAB are performed for the illustration of the efficacy and superiority of the given scheme. Lastly, we also have presented some conclusions and the future prospects of the discussed research work in Section 6.
We here formulate a methodology to examine compound difference anti-synchronization (CDAS) scheme viewing master-slave framework in four chaotic systems which would be utilized in the coming up sections.
Let the scaling master system be
˙wm1= f1(wm1), | (2.1) |
and the base second master systems be
˙wm2= f2(wm2), | (2.2) |
˙wm3= f3(wm3). | (2.3) |
Corresponding to the aforementioned master systems, let the slave system be
˙ws4= f4(ws4)+U(wm1,wm2,wm3,ws4), | (2.4) |
where wm1=(wm11,wm12,...,wm1n)T∈Rn, wm2=(wm21,wm22,...,wm2n)T∈Rn, wm3=(wm31,wm32,...,wm3n)T∈Rn, ws4=(ws41,ws42,...,ws4n)T∈Rn are the state variables of the respective chaotic systems (2.1)–(2.4), f1,f2,f3,f4:Rn→Rn are four continuous vector functions, U=(U1,U2,...,Un)T:Rn×Rn×Rn×Rn→Rn are appropriately constructed active controllers.
Compound difference anti-synchronization error (CDAS) is defined as
E=Sws4+Pwm1(Rwm3−Qwm2), |
where P=diag(p1,p2,.....,pn),Q=diag(q1,q2,.....,qn),R=diag(r1,r2,.....,rn),S=diag(s1,s2,.....,sn) and S≠0.
Definition: The master chaotic systems (2.1)–(2.3) are said to achieve CDAS with slave chaotic system (2.4) if
limt→∞‖E(t)‖=limt→∞‖Sws4(t)+Pwm1(t)(Rwm3(t)−Qwm2(t))‖=0. |
We now present our proposed CDAS approach in three master systems (2.1)–(2.3) and one slave system (2.4). We next construct the controllers based on CDAS approach by
Ui= ηisi−(f4)i−KiEisi, | (3.1) |
where ηi=pi(f1)i(riwm3i−qiwm2i)+piwm1i(ri(f3)i−qi(f2)i), for i=1,2,...,n.
Theorem: The systems (2.1)–(2.4) will attain the investigated CDAS approach globally and asymptotically if the active control functions are constructed in accordance with (3.1).
Proof. Considering the error as
Ei= siws4i+piwm1i(riwm3i−qiwm2i),fori=1,2,3,.....,n. |
Error dynamical system takes the form
˙Ei= si˙ws4i+pi˙wm1i(riwm3i−qiwm2i)+piwm1i(ri˙wm3i−qi˙wm2i)= si((f4)i+Ui)+pi(f1)i(riwm3i−qiwm2i)+piwm1i(ri(f3)i−qi(f2)i)= si((f4)i+Ui)+ηi, |
where ηi=pi(f1)i(riwm3i−qiwm2i)+piwm1i(ri(f3)i−qi(f2)i), i=1,2,3,....,n. This implies that
˙Ei= si((f4)i−ηisi−(f4)i−KiEisi)+ηi= −KiEi | (3.2) |
The classic Lyapunov function V(E(t)) is described by
V(E(t))= 12ETE= 12ΣE2i |
Differentiation of V(E(t)) gives
˙V(E(t))=ΣEi˙Ei |
Using Eq (3.2), one finds that
˙V(E(t))=ΣEi(−KiEi)= −ΣKiE2i). | (3.3) |
An appropriate selection of (K1,K1,.......,Kn) makes ˙V(E(t)) of eq (3.3), a negative definite. Consequently, by LSA [47], we obtain
limt→∞Ei(t)=0,(i=1,2,3). |
Hence, the master systems (2.1)–(2.3) and slave system (2.4) have attained desired CDAS strategy.
We now describe GLVBS as the scaling master system:
{˙wm11=wm11−wm11wm12+b3w2m11−b1w2m11wm13,˙wm12=−wm12+wm11wm12,˙wm13=b2wm13+b1w2m11wm13, | (4.1) |
where (wm11,wm12,wm13)T∈R3 is state vector of (4.1). Also, wm11 represents the prey population and wm12, wm13 denote the predator populations. For parameters b1=2.9851, b2=3, b3=2 and initial conditions (27.5,23.1,11.4), scaling master GLVBS displays chaotic/disturbed behaviour as depicted in Figure 1(a).
The base master systems are the identical chaotic GLVBSs prescribed respectively as:
{˙wm21=wm21−wm21wm22+b3w2m21−b1w2m21wm23,˙wm22=−wm22+wm21wm22,˙wm23=b2wm23+b1w2m21wm23, | (4.2) |
where (wm21,wm22,wm23)T∈R3 is state vector of (4.2). For parameter values b1=2.9851, b2=3, b3=2, this base master GLVBS shows chaotic/disturbed behaviour for initial conditions (1.2,1.2,1.2) as displayed in Figure 1(b).
{˙wm31=wm31−wm31wm32+b3w2m31−b1w2m31wm33,˙wm32=−wm32+wm31wm32,˙wm33=b2wm33+b1w2m31wm33, | (4.3) |
where (wm31,wm32,wm33)T∈R3 is state vector of (4.3). For parameters b1=2.9851, b2=3, b3=2, this second base master GLVBS displays chaotic/disturbed behaviour for initial conditions (2.9,12.8,20.3) as shown in Figure 1(c).
The slave system, represented by similar GLVBS, is presented by
{˙ws41=ws41−ws41ws42+b3w2s41−b1w2s41ws43+U1,˙ws42=−ws42+ws41ws42+U2,˙ws43=b2ws43+b1w2s41ws43+U3, | (4.4) |
where (ws41,ws42,ws43)T∈R3 is state vector of (4.4). For parameter values, b1=2.9851, b2=3, b3=2 and initial conditions (5.1,7.4,20.8), the slave GLVBS exhibits chaotic/disturbed behaviour as mentioned in Figure 1(d).
Moreover, the detailed theoretical study for (4.1)–(4.4) can be found in [4]. Further, U1, U2 and U3 are controllers to be determined.
Next, the CDAS technique has been discussed for synchronizing the states of chaotic GLVBS. Also, LSA-based ACS is explored & the necessary stability criterion is established.
Here, we assume P=diag(p1,p2,p3), Q=diag(q1,q2,q3), R=diag(r1,r2,r3), S=diag(s1,s2,s3). The scaling factors pi,qi,ri,si for i=1,2,3 are selected as required and can assume the same or different values.
The error functions (E1,E2,E3) are defined as:
{E1=s1ws41+p1wm11(r1wm31−q1wm21),E2=s2ws42+p2wm12(r2wm32−q2wm22),E3=s3ws43+p3wm13(r3wm33−q3wm23). | (4.5) |
The major objective of the given work is the designing of active control functions Ui,(i=1,2,3) ensuring that the error functions represented in (4.5) must satisfy
limt→∞Ei(t)=0for(i=1,2,3). |
Therefore, subsequent error dynamics become
{˙E1=s1˙ws41+p1˙wm11(r1wm31−q1wm21)+p1wm11(r1˙wm31−q1˙wm21),˙E2=s2˙ws42+p2˙wm12(r2wm32−q2wm22)+p2wm12(r2˙wm32−q2˙wm22),˙E3=s3˙ws43+p3˙wm13(r3wm33−q3wm23)+p3wm13(r3˙wm33−q3˙wm23). | (4.6) |
Using (4.1), (4.2), (4.3), and (4.5) in (4.6), the error dynamics simplifies to
{˙E1=s1(ws41−ws41ws42+b3w2s41−b1w2s41ws43+U1)+p1(wm11−wm11wm12+b3w2m11−b1w2m11wm13)(r1wm31−q1wm21)+p1wm11(r1(wm31−wm31wm32+b3w2m31−b1w2m31wm33)−q1(wm21−wm21wm22+b3w2m21−b1w2m21wm23),˙E2=s2(−ws42+ws41ws42+U2)+p2(−wm12+wm11wm12)(r2wm32−q2wm22)+p2wm12(r2(−wm32+wm31wm32)−q2(−wm22+wm21wm22)),˙E3=s3(b2ws43+b1w2s41ws43+U3)+p3(b2wm13+b1w2m11wm13)(r3wm33−q3wm23)+p3wm13(r3(b2wm33+b1w2m31wm33)−q3(b2wm23+b1w2m21wm23)). | (4.7) |
Let us now choose the active controllers:
U1= η1s1−(f4)1−K1E1s1, | (4.8) |
where η1=p1(f1)1(r1wm31−q1wm21)+p1wm11(r1(f3)1−q1(f2)1), as described in (3.1).
U2= η2s2−(f4)2−K2E2s2, | (4.9) |
where η2=p2(f1)2(r2wm32−q2wm22)+p2wm12(r2(f3)2−q2(f2)2).
U3= η3s3−(f4)3−K3E3s3, | (4.10) |
where η3=p3(f1)3(r3wm33−q3wm23)+p3wm13(r3(f3)3−q3(f2)3) and K1>0,K2>0,K3>0 are gaining constants.
By substituting the controllers (4.8), (4.9) and (4.10) in (4.7), we obtain
{˙E1=−K1E1,˙E2=−K2E2,˙E3=−K3E3. | (4.11) |
Lyapunov function V(E(t)) is now described by
V(E(t))= 12[E21+E22+E23]. | (4.12) |
Obviously, the Lyapunov function V(E(t)) is +ve definite in R3. Therefore, the derivative of V(E(t)) as given in (4.12) can be formulated as:
˙V(E(t))= E1˙E1+E2˙E2+E3˙E3. | (4.13) |
Using (4.11) in (4.13), one finds that
˙V(E(t))= −K1E21−K2E22−K3E23<0, |
which displays that ˙V(E(t)) is -ve definite.
In view of LSA [47], we, therefore, understand that CDAS error dynamics is globally as well as asymptotically stable, i.e., CDAS error E(t)→0 asymptotically for t→∞ to each initial value E(0)∈R3.
This section conducts a few simulation results for illustrating the efficacy of the investigated CDAS scheme in identical chaotic GLVBSs using ACS. We use 4th order Runge-Kutta algorithm for solving the considered ordinary differential equations. Initial conditions for three master systems (4.1)–(4.3) and slave system (4.4) are (27.5,23.1,11.4), (1.2,1.2,1.2), (2.9,12.8,20.3) and (14.5,3.4,10.1) respectively. We attain the CDAS technique among three masters (4.1)–(4.3) and corresponding one slave system (4.4) by taking pi=qi=ri=si=1, which implies that the slave system would be entirely anti-synchronized with the compound of three master models for i=1,2,3. In addition, the control gains (K1,K2,K3) are taken as 2. Also, Figure 2(a)–(c) indicates the CDAS synchronized trajectories of three master (4.1)–(4.3) & one slave system (4.4) respectively. Moreover, synchronization error functions (E1,E2,E3)=(51.85,275.36,238.54) approach 0 as t tends to infinity which is exhibited via Figure 2(d). Hence, the proposed CDAS strategy in three masters and one slave models/systems has been demonstrated computationally.
In this work, the investigated CDAS approach in similar four chaotic GLVBSs using ACS has been analyzed. Lyapunov's stability analysis has been used to construct proper active nonlinear controllers. The considered error system, on the evolution of time, converges to zero globally & asymptotically via our appropriately designed simple active controllers. Additionally, numerical simulations via MATLAB suggest that the newly described nonlinear control functions are immensely efficient in synchronizing the chaotic regime found in GLVBSs to fitting set points which exhibit the efficacy and supremacy of our proposed CDAS strategy. Exceptionally, both analytic theory and computational results are in complete agreement. Our proposed approach is simple yet analytically precise. The control and synchronization among the complex GLVBSs with the complex dynamical network would be an open research problem. Also, in this direction, we may extend the considered CDAS technique on chaotic systems that interfered with model uncertainties as well as external disturbances.
The authors gratefully acknowledge Qassim University, represented by the Deanship of Scientific Research, on the financial support for this research under the number 10163-qec-2020-1-3-I during the academic year 1441 AH/2020 AD.
The authors declare there is no conflict of interest.
[1] | Centers for Disease Control and Prevention, Coronavirus (COVID-19). Available from: https://www.cdc.gov/coronavirus/2019-ncov. |
[2] | World Health Organization, Coronavirus disease (COVID-19) pandemic. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019. |
[3] |
A. Clark, M. Jit, C. Warren-Gash, B. Guthrie, H. H. X. Wang, S. W. Mercer, et al., Global, regional, and national estimates of the population at increased risk of severe COVID-19 due to underlying health conditions in 2020: A modelling study, Lancet Glob. Health, 8 (2020), e1003–e1017. doi: 10.1016/S2214-109X(20)30264-3
![]() |
[4] |
E. K. Stokes, L. D. Zambrano, K. N. Anderson, E. P. Marder, K. M. Raz, S. E. B. Felix, et al., Coronavirus disease 2019 case surveillance – United States, January 22–May 30, 2020, Morb. Mortal. Wkly. Rep., 69 (2020), 759-765. doi: 10.15585/mmwr.mm6924e2
![]() |
[5] |
H. Razzaghi, Y. Wang, H. Lu, K. E. Marshall, N. F. Dowling, G. Paz-Bailey, et al., Estimated county-level prevalence of selected underlying medical conditions associated with increased risk for severe COVID-19 illness – United States, 2018, Morb. Mortal. Wkly. Rep., 69 (2020), 945–950. doi: 10.15585/mmwr.mm6929a1
![]() |
[6] |
Z. J. Cheng, J. Shan, 2019 Novel coronavirus: where we are and what we know, Infection, 48 (2020), 155–163. doi: 10.1007/s15010-020-01401-y
![]() |
[7] | A. Sahin, A. Erdogan, P. M. Agaoglu, Y. Dineri, A. Cakirci, M. Senel, et al., 2019 novel coronavirus (COVID-19) outbreak: A review of the current literature, Eurasian J. Med. Oncol., 4 (2020), 1–7. |
[8] |
K. Leung, J. T. Wu, D. Liu, G. M. Leung, First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment, Lancet, 395 (2020), 1382–1393. doi: 10.1016/S0140-6736(20)30746-7
![]() |
[9] |
R. Li, S. Pei, B. Chen, Y. Song, T. Zhang, W. Yang, et al., Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2), Science, 368 (2020), 489–493. doi: 10.1126/science.abb3221
![]() |
[10] | J. M. Read, J. R. E. Bridgen, D. A. T. Cummings, A. Ho, C. P. Jewell, Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions, MedRxiv, 2020. Available from: https://doi.org/10.1101/2020.01.23.20018549. |
[11] |
B. Tang, X. Wang, Q. Li, N. L. Bragazzi, S. Tang, Y. Xiao, et al., Estimation of the Transmission Risk of 2019-nCoV and Its Implication for Public Health Interventions, J. Clin. Med., 9 (2020), 462. doi: 10.3390/jcm9020462
![]() |
[12] |
J. Wang, Mathematical models for COVID-19: applications, limitations, and potentials, J. Public Health Emerg., 4 (2020), 9. doi: 10.21037/jphe-2020-05
![]() |
[13] |
J. T. Wu, K. Leung, G. M. Leung, Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study, Lancet, 395 (2020), 689–697. doi: 10.1016/S0140-6736(20)30260-9
![]() |
[14] |
C. Yang, J. Wang, A mathematical model for the novel coronavirus epidemic in Wuhan, China, Math. Bios. Eng., 17 (2020), 2708–2724. doi: 10.3934/mbe.2020148
![]() |
[15] |
C. Yang, J. Wang, Modeling the transmission of COVID-19 in the US – A case study, Infect. Disease Modelling, 6 (2021), 195–211. doi: 10.1016/j.idm.2020.12.006
![]() |
[16] |
C. Yang, J. Wang, Transmission rates and environmental reservoirs for COVID-19: A modeling study, J. Biol. Dyn., 15 (2021), 86–108. doi: 10.1080/17513758.2020.1869844
![]() |
[17] |
H. Zhong, W. Wang, Mathematical analysis for COVID-19 resurgence in the contaminated environment, Math. Bios. Eng., 17 (2020), 6909–6927. doi: 10.3934/mbe.2020357
![]() |
[18] | Data USA: Chattanooga, TN. Available from: https://datausa.io/profile/geo/chattanooga-tn/. |
[19] | Hamilton County Health Department. Available from: https://health.hamiltontn.org/default.aspx. |
[20] | COVID-19 in Hamilton County, TN. Available from: https://sites.google.com/view/hamiltoncounty-tn-covid19. |
[21] |
C. Yang, X. Wang, D. Gao, J. Wang, Impact of awareness programs on cholera dynamics: Two modeling approaches, Bull. Math. Biol., 79 (2017), 2109–2131. doi: 10.1007/s11538-017-0322-1
![]() |
[22] |
J. F. W. Chan, S. Yuan, K. H. Kok, K. K. W. To, H. Chu, J. Yang, et al., A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster, Lancet, 395 (2020), 514–523. doi: 10.1016/S0140-6736(20)30154-9
![]() |
[23] |
A. Kimball, K. M. Hatfield, M. Arons, A. James, J. Taylor, K. Spicer, et al., Asymptomatic and presymptomatic SARS-CoV-2 infections in residents of a long-term care skilled nursing facility – King County, Washington, March 2020, Morb. Mortal. Wkly. Rep., 69 (2020), 377–381. doi: 10.15585/mmwr.mm6913e1
![]() |
[24] |
W. E. Wei, Z. Li, C. J. Chiew, S. E. Yong, M. P. Toh, V. J. Lee, Presymptomatic transmission of SARS-CoV-2 – Singapore, January 23 - March 16, 2020, Morb. Mortal. Wkly. Rep., 69 (2020), 411–415. doi: 10.15585/mmwr.mm6914e1
![]() |
[25] |
M. A. Johansson, T. M. Quandelacy, S. Kada, P. V. Prasad, M. Steele, J. T. Brooks, et al., SARS-CoV-2 transmission from people without COVID-19 symptoms, JAMA Network Open, 4 (2021), e2035057. doi: 10.1001/jamanetworkopen.2020.35057
![]() |
[26] | J. A. Spencer, D. P. Shutt, S. K. Moser, H. Clegg, H. J. Wearing, H. Mukundan, et al., Epidemiological parameter review and comparative dynamics of influenza, respiratory syncytial virus, rhinovirus, human coronvirus, and adenovirus, MedRxiv, 2020. Available from: http://dx.doi.org/10.1101/2020.02.04.20020404. |
[27] | World Health Organization, Coronavirus disease (COVID-19) situation reports. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports. |
[28] |
Z. S. Almalki, M. F. Khan, S. Almazrou, A. S. Alanazi, M. S. Iqbal, A. Alqahtani, et al., Clinical characteristics and outcomes among COVID-19 hospitalized patients with chronic conditions: A retrospective single-center study, J. Multidiscip. Healthc., 13 (2020), 1089–1097. doi: 10.2147/JMDH.S273918
![]() |
[29] |
L. Ellwein, H. Tran, C. Zapata, V. Novak, M. Olufsen, Sensitivity analysis and model assessment: Mathematical models for arterial blood flow and blood pressure, J. Cardiovasc. Eng., 8 (2008), 94–108. doi: 10.1007/s10558-007-9047-3
![]() |
[30] | COVID-19 Vaccinations in Hamilton County, TN. Available from: https://datastudio.google.com/reporting/1d5c1095-e7a1-4a05-80db-e81e33f323fb/page/mxd0B?s=pj_BWXEzAU8. |
[31] |
H. J. Wearing, P. Rohani, M. J. Keeling, Appropriate models for the management of infectious diseases, PLoS Med., 2 (2005), e174. doi: 10.1371/journal.pmed.0020174
![]() |
[32] |
P. van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Bios., 180 (2002), 29–48. doi: 10.1016/S0025-5564(02)00108-6
![]() |
[33] | J. P. LaSalle, The Stability of Dynamical Systems, CBMS-NSF Regional Conference Series in Applied Mathematics, Philadelphia, 1976. |
1. | Muhammad Zubair Mehboob, Arslan Hamid, Jeevotham Senthil Kumar, Xia Lei, Comprehensive characterization of pathogenic missense CTRP6 variants and their association with cancer, 2025, 25, 1471-2407, 10.1186/s12885-025-13685-0 |