
Citation: Xuejuan Lu, Shaokai Wang, Shengqiang Liu, Jia Li. An SEI infection model incorporating media impact[J]. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1317-1335. doi: 10.3934/mbe.2017068
[1] | Rina Su, Chunrui Zhang . The generation mechanism of Turing-pattern in a Tree-grass competition model with cross diffusion and time delay. Mathematical Biosciences and Engineering, 2022, 19(12): 12073-12103. doi: 10.3934/mbe.2022562 |
[2] | Yue Xing, Weihua Jiang, Xun Cao . Multi-stable and spatiotemporal staggered patterns in a predator-prey model with predator-taxis and delay. Mathematical Biosciences and Engineering, 2023, 20(10): 18413-18444. doi: 10.3934/mbe.2023818 |
[3] | Mingzhu Qu, Chunrui Zhang, Xingjian Wang . Analysis of dynamic properties on forest restoration-population pressure model. Mathematical Biosciences and Engineering, 2020, 17(4): 3567-3581. doi: 10.3934/mbe.2020201 |
[4] | Fang Liu, Yanfei Du . Spatiotemporal dynamics of a diffusive predator-prey model with delay and Allee effect in predator. Mathematical Biosciences and Engineering, 2023, 20(11): 19372-19400. doi: 10.3934/mbe.2023857 |
[5] | Ranjit Kumar Upadhyay, Swati Mishra . Population dynamic consequences of fearful prey in a spatiotemporal predator-prey system. Mathematical Biosciences and Engineering, 2019, 16(1): 338-372. doi: 10.3934/mbe.2019017 |
[6] | Sourav Kumar Sasmal, Jeet Banerjee, Yasuhiro Takeuchi . Dynamics and spatio-temporal patterns in a prey–predator system with aposematic prey. Mathematical Biosciences and Engineering, 2019, 16(5): 3864-3884. doi: 10.3934/mbe.2019191 |
[7] | Swadesh Pal, Malay Banerjee, Vitaly Volpert . Spatio-temporal Bazykin’s model with space-time nonlocality. Mathematical Biosciences and Engineering, 2020, 17(5): 4801-4824. doi: 10.3934/mbe.2020262 |
[8] | Meiling Zhu, Huijun Xu . Dynamics of a delayed reaction-diffusion predator-prey model with the effect of the toxins. Mathematical Biosciences and Engineering, 2023, 20(4): 6894-6911. doi: 10.3934/mbe.2023297 |
[9] | Tingting Ma, Xinzhu Meng . Global analysis and Hopf-bifurcation in a cross-diffusion prey-predator system with fear effect and predator cannibalism. Mathematical Biosciences and Engineering, 2022, 19(6): 6040-6071. doi: 10.3934/mbe.2022282 |
[10] | Xiaomei Bao, Canrong Tian . Turing patterns in a networked vegetation model. Mathematical Biosciences and Engineering, 2024, 21(11): 7601-7620. doi: 10.3934/mbe.2024334 |
Respiratory syncytial virus (RSV), although the leading cause of hospitalization of infants in developed countries, is still awaiting development of a safe and effective vaccine [1]. Infants typically contract RSV within the first two years of life [2], although there has been a recent suppression of RSV due to COVID-19 mitigation measures [3,4]. This has resulted in earlier and larger epidemics of RSV in many countries this year as COVID mitigation measures have been relaxed [5,6]. RSV can also cause serious illness and death in elderly adults [7]. While the severity of RSV in these two groups is thought to be in part due to a limited immune response [8], other factors might also play a role. A recent study by Felt et al. examined the role of defective viral genomes (DVGs) [9], finding that they might be an indicator of disease severity in RSV.
DVGs, are a feature of a number of viral infections besides RSV [10], including influenza [11,12,13], Sindbis virus [14], vesicular stomatitis virus [15,16], poliovirus [17,18,19], measles [20], SARS-CoV-2 [21], and parvovirus [22]. Defective genomes are thought to have a number of effects on the course of an infection. They interfere with production of fully competent infectious virus when they coinfect cells [11,23]. Interestingly, the within-cell interference can lead to persistence of the infection within the host [23,24,25]. More recently, studies have determined that DVGs also alter the immune response during an infection [9,26,27,28], with evidence suggesting that the presence of DVGs stimulates the innate immune response. Some DVGs are also being investigated as potential vaccines since they have the ability to also stimulate the adaptive immune response [29]. DVGs are also thought to play a role in the evolution of viruses providing increased genetic diversity during infections in a single host [30]. Thus there is a need to develop a better understanding of how the appearance of DVGs alters viral infections.
Mathematical modeling can help us develop an understanding of the role of DVGs in viral infections. There have already been a number of studies that examined models incorporating defective interfering particles (DIPs). DIPs occur when DVGs are packaged and exported out of the cell. DIPs can infect cells, but cannot produce new progeny virions on their own. They can, however, produce new DIPs when in a cell that is infected with both the DIP and a fully competent virus, since the full virus provides the missing components needed for replication, virion formation, and release [31,32]. For example, Liao et al. used a model to determine conditions under which current DIP quantification techniques resulted in correct estimates of the number of DIPs [11]. Shirogane et al. [19] examined how the DIPs competed with wild-type virus in a competitive coinfection model. Other within host models examined the possibility of using defective particles to treat viral infections [33,34,35]. More recent modeling efforts combined intracellular models with within host models in order to understand not only the effect of DIPs at the within host level, but to understand how DVGs might interfere with normal viral replication within the cell [36,37,38].
A limitation of all of these studies is that the models were developed based on assumptions of the effect of DIPs and the models were then used to examine possible outcomes of these assumptions. Unfortunately, many of the known effects of DVGs, such as their effect on the immune response [9,26,27,28,29], were not incorporated into these models, perhaps missing an important process that could alter the interplay between DVGs and fully competent virions. A recent study by Felt et al. [9] could allow for assessment of how DVGs alter the time course of a viral infection — information that can then be incorporated into future modeling efforts in order to produce more accurate mathematical models.
Here, we use mathematical models and the Felt at el. study to examine how DVGs alter the time course of viral load. Without directly incorporating DVGs, we fit a simple mathematical model to experimental data in order to identify differences in viral kinetics parameters under different infection conditions, such as infections with no DVGs and infections where DVGs appear. This type of approach has been used to quantify differences in viral kinetics parameters between wild-type and drug-resistant strains of influenza [39,40], quantify differences in aging animals with RSV [41,42] or SARS-CoV-2 [43], and examine the effect of an antiviral [44]. We find that the appearance of DVGs primarily changes the infection rate and the clearance rate of infections, but only when the DVGs appear late in the infection or are prolonged.
We use a basic model of viral infection that does not include an eclipse phase and assumes that cells produce and release virus immediately upon infection [45],
dTdt=−βTVdIdt=βTV−δIdVdt=pI−cV. | (2.1) |
In this model, T, I, and V are the variables of the model and β, δ, c, and p are the parameters. T represents the number of target cells infected with virus V at a rate β. Once infected, the cells enter the infectious phase I, in which they produce virus at a rate p and die at rate δ. Virus is cleared at rate c. While this model excludes the eclipse phase — where the cell is infected, but not yet producing virus — inclusion of this phase would add an additional parameter that is known to not be identifiable from viral titer data alone [46].
The experimental data used in this study was digitized from Felt et al. [9], which describes a challenge study of RSV infections in healthy adults. In this study, data was collected from nasal swabs taken from healthy adults challenged with RSV. 61 experimentally infected healthy non-smoking adults, aged 18-50 years old, were inoculated experimentally via intranasal drops. Nasal washes and blood samples were taken immediately before inoculation, daily for the first 10 days post infection, and again on days 14 and 28. For all patients, quantitative polymerase chain reaction with reverse transcription methods were used to estimate viral loads. Patients were separated into different groups based on the dynamics of DVGs during the infection and median viral titer curves were used to compare the different groups. One comparison was between patients who had no DVGs detected, those who had DVGs appear early (before 3, 5, or 6 days post infection), and those who had DVGs appear late (after 3, 5, or 6 days post infection). Another comparison was between patients who had no DVGs detected, those who had transient DVG detection (not detected after 6 dpi), and those who had prolonged DVG detection (detected at least once before 6 dpi and detected again after 6 dpi).
Data was extracted from figures using Webplotdigitizer (Webplotdigitizer.com). To fit the data and calculate the parameter estimates for data sets from this study, we used Python's built-in optimizing function scipy.optimize.minimize and minimized the sum of squared residuals (SSR) using the Nelder-Mead algorithm. In addition to estimating the free parameters of the model (β, p, δ, c), we use the estimates to calculate the basic reproductive number, R0=βpcδ, which is the average number of secondary infections caused by a single infected cell, and the infecting time tinf=√2pβ [47], which is the average time between release of virus from one cell and infection of the next cell.
Bootstrapping was used to estimate the distributions of the parameter estimates [48]. We shuffled the residuals at each time point using the numpy shuffle function, which uses a Fisher-Yates shuffle. We then repeated the minimization process to produce another set of parameter estimates. After repeating this process 1000 times, we calculated the 95% confidence interval for each set of parameters.
To compare the parameter distributions, we performed Mann-Whitney U tests using Python's built-in mannwhitneyu function to check if parameter distributions were distinct. The Mann-Whitney test is used since parameter distributions cannot all be assumed to be normal. p values less than 0.05 are considered statistically significant.
In the study by Felt et al. [9], infections are compared using two main features: the time at which DVGs are first detected, and the duration of DVG detection. In the first comparison, participants are divided into three groups: those where DVGs are not detected, those that have DVGs detected early and those that have DVGs detected late. Cutoffs to define early versus late start were set at 3, 5, or 6 days post-infection. The second comparison also consists of three groups: those who do not shed DVGs, those who shed DVGs transiently, and those who shed DVGs for a prolonged period of time. Prolonged DVG was defined as detection of DVGs on two or more days after 6 days post-infection. We perform our analysis using the same groupings.
We show the viral load curves of Figure 5C in the Felt et al. manuscript and our best fit curves along with distributions for the parameter estimates in Figure 1. These viral load curves represent the mean viral load of people with no DVGs detected (A), early DVGs (B), and late DVGs (C) with the early/late cutoff set at 5 days post-infection. The best fit parameter values are given in Table 1. Distributions for the parameters of the viral kinetics model, β, p, c, and δ, as well as distributions for the calculated values R0 and tinf, for all three groups are given in Figure 1 (bottom two rows). p values from the Mann-Whitney test are also presented in Table 1.
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.93×10−6 | 1.77×109 | 26.9 | 181 | 1.07 | 0.0200 | 1.82 |
95% CI | (0.569–221)×10−6 | 3.76×104–1.01×1015 | 0.864–401 | 0.614–96.1 | 1.02–6.82 | 1.96×10−5–0.506 | |
Early | 4.38×10−5 | 5.06×106 | 7.62 | 23.1 | 1.26 | 0.0950 | 2.09 |
95% CI | (0.106–42.5)×10−5 | 2.44×104–5.78×1014 | 1.28–2066 | 1.64–472 | 1.01–3.52 | 3.88×10−5–0.442 | |
Late | 2.21×10−8 | 4.95×1012 | 226 | 476 | 1.01 | 0.00428 | 3.26 |
95% CI | (2.10–186)×10−8 | (0.000309–12.6)×1012 | 2.34–408 | 50.6–1350 | 1.02–2.35 | 0.00279–0.0703 | |
p values | |||||||
Early/DVG- | 0.389 | 0.356 | 0.247 | 0.076 | 0.405 | 0.431 | |
Early/Late | 0.0005 | 0.098 | 0.365 | 0.008 | 0.440 | 0.269 | |
DVG-/Late | 0.001 | 0.478 | 0.141 | 5.82×10−5 | 0.436 | 0.475 |
The model is able to fit the viral load curves quite well. The R0 values are fairly low, just over the threshold needed to initiate an infection (R0=1), which is smaller than other estimates found in literature [42,49], although those RSV infections were in animals other than humans. The infecting times range from about 6 minutes to 2 hours, which are more in line with previous estimates [42,49,50]. When comparing parameter estimates between the 3 groups, we note that parameter differences show up primarily in the infection rate, β, and the clearance rate c. When DVGs appear late, the infection rate appears to be lower and the clearance rate appears to be higher, than when DVGs appear early or not at all.
Since the cut-off of 5 days for early vs. late is arbitrary, two other cut-off times were tested. An earlier cut-off of 3 days is shown in Figure 2. In this case, the mean viral load in the early DVG group stayed below the threshold of detection, so we show only the viral time courses and best fit curves for no DVGs (A) and late DVGs (B). The best fit parameter values along with 95% confidence intervals are given in Table 2. Distributions for the parameters of the viral kinetics model, β, p, c, and δ, as well as distributions for the calculated values R0 and tinf, for all three groups are given in Figure 2 (bottom two rows). p values from the Mann-Whitney test are also presented in Table 2.
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 3.10×10−5 | 3.66×107 | 2.41 | 251 | 1.88 | 0.0420 | 0.755 |
95% CI | (0.0948–6.86)×10−5 | 2.56×106–1.53×1016 | 1.53–558 | 16.5–3880 | 1.00–2.48 | 9.41×10−6–0.117 | |
Late | 3.27×10−7 | 5.48×109 | 9.01 | 151 | 1.32 | 0.0334 | 0.821 |
95% CI | (0.453–12.3)×10−6 | (0.262–1280)×109 | 3.47–159 | 31.2–365 | 1.03–1.86 | 0.00653–0.0838 | |
p values | |||||||
DVG-/Late | 0.015 | 0.483 | 0.228 | 0.521 | 0.167 | 0.281 |
The mean viral titer curve for the early DVG group remained below threshold for all time points, making it impossible to fit the model to this data. However, there were three patients in the early DVG group with the earliest cut-off of 3 days. This means that the early DVG patients were shedding very low amounts throughout the infection. For the remaining two groups, we see slightly higher basic reproduction numbers, although they are still fairly close to the threshold of 1. The infecting times are 45 minutes to 1 hour, which are similar to previous estimates [42,49,50]. We can only compare the no DVG and late DVG groups in this case, and we note that the only parameter that is statistically significantly different between the two groups is the infection rate.
A later cut-off date of 6 days was also tested. This is shown in Figure 3. The best fit parameter values along with 95% confidence intervals are given in Table 3. Distributions of the parameter estimates of the viral kinetics model for all three groups are given in Figure 3 (bottom two rows). p values from the Mann-Whitney test are also presented in Table 3.
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.20×10−5 | 1.82×107 | 3.25 | 78.0 | 1.58 | 0.0707 | 1.47 |
95% CI | (0.0667–24.8)×10−5 | (0.00348–100)×107 | 0.771–1840 | 0.542–92.3 | 1.02–6.57 | 8.26×10−6–0.521 | |
Early | 5.00×10−6 | 2.19×109 | 31.7 | 326 | 1.06 | 0.0135 | 2.06 |
95% CI | (0.508–180)×10−6 | (0.000128–572)×109 | 1.58–1760 | 2.70–869 | 1.01–2.77 | 4.10×10−6–0.296 | |
Late | 1.68×10−8 | 1.31×1013 | 69.0 | 3060 | 1.04 | 0.00301 | 2.26 |
95% CI | (0.864–91.3)×10−8 | 9.15×108–3.92×1013 | 2.67–476 | 57.7–2080 | 1.01–2.24 | 0.00257–0.0562 | |
p values | |||||||
Early/DVG- | 0.441 | 0.488 | 0.514 | 0.005 | 0.426 | 0.476 | |
Early/Late | 0.005 | 0.238 | 0.365 | 0.024 | 0.468 | 0.446 | |
DVG-/Late | 7.44×10−5 | 0.216 | 0.197 | 2.17×10−5 | 0.503 | 0.478 |
As with the two previous assessments of the effect of time of appearance of DVGs, we see a statistically significant difference in the infection rates of early or no DVGs and late DVGs. We also see a statistically significant difference in the clearance rate between early or no DVGs and late DVGs. In this case, there is also a statistically significant difference between the clearance rates of early and no DVGs. So even when the cut-off day marking the difference between early and late DVGs is shifted, there is a consistent difference in infection rates between early or no DVGs and late DVGs. The differences in clearance rate appear to be more dependent on the chosen cut-off point.
We show the viral load curves of Figure 6C in the Felt et al. manuscript and our best fit curves along with distributions for the parameter estimates in Figure 4. These viral load curves represent the mean viral load of people shedding no DVGs (A), transient DVGs (B), and prolonged DVGs (C). The best fit parameter values are given in Table 4. Distributions for the parameters of the viral kinetics model, β, p, c, and δ as well as distributions for the calculated values R0 and tinf, for all three groups are given in Figure 4 (bottom two rows). p values from the Mann-Whitney test are presented in Table 4.
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.55×10−7 | 1.89×1012 | 387 | 1240 | 1.01 | 0.00200 | 2.21 |
95% CI | (0.377–2400)×10−7 | 4.25×104–1.47×1017 | 1.61–486 | 0.853–3910 | 1.00–3.98 | 2.76×10−6–0.455 | |
Transient | 1.01×10−5 | 1.78×108 | 11.8 | 131 | 1.16 | 0.0334 | 1.53 |
95% CI | (0.0654–20.2)×10−5 | 5.50×104–1.08×1016 | 1.17–1817 | 0.961–139 | 1.01–4.82 | 7.22×10−6–0.461 | |
Prolonged | 1.02×10−8 | 2.18×1013 | 61.1 | 3460 | 1.05 | 0.00300 | 2.83 |
95% CI | (0.496–24.9)×10−8 | 5.35×108–7.91×1014 | 18.1–2600 | 3.56–1330 | 1.01–1.96 | 0.00106–0.0908 | |
p values | |||||||
Transient/DVG- | 0.371 | 0.411 | 0.510 | 0.129 | 0.469 | 0.443 | |
Transient/Prolonged | 1.30×10−5 | 0.206 | 0.249 | 0.001 | 0.493 | 0.519 | |
DVG-/Prolonged | 2.37×10−5 | 0.428 | 0.227 | 0.034 | 0.343 | 0.382 |
The model, again, fits the viral load curves well. R0 values are again just over the threshold for infection. The infecting times range from 3 to 48 minutes. The distribution for β estimates for prolonged DVGs is clearly lower than the β distributions for no DVGs or transient DVGs. This is reflected in the statistical tests, which show p values less than 0.05 for no DVG/prolonged DVGs and transient DVGs/prolonged DVGs comparisons of β estimates. We also see similar statistically significant differences in the estimates for the clearance rate c when comparing the prolonged DVGs to the other two groups.
In this manuscript, we fit a mathematical model of viral infections to data examining the appearance of DVGs in RSV infections. We found that the infection rate and the clearance rate are the parameters most affected by the presence of DVGs. Specifically, the infection rate is lower when DVGs appear later and are prolonged during the infection, while the clearance rate is higher in these two cases. In the case of the clearance rate, a possible mechanism is fairly straightforward. When DVGs occur late and are prolonged, they are essentially displacing fully competent virions, which then decrease more quickly than in infections that do not have DVGs or where DVGs are not present later in the infection. The link between later and longer DVGs and lower infection rate is less clear. Defective particles are known to interfere directly with intracellular replication of fully competent virus [51], but they can also interfere with infection dynamics by stimulating the immune response [52]. Since DVGs are known to stimulate the innate immune response in particular [9,26,27,28], this could be an indirect mechanism that reduces the infection rate since the innate immune response typically rises early and rapidly during infections [53]. It is also possible that a lower infection rate leads to a later appearance of DVGs. Since replication is slower in the case of low infection rates, replication of DVGs would presumably also be slower causing them to be detected later in the infection.
While our study could not quantify any differences in the viral time courses of early and no DVGs, the original study by Felt et al. noted that there was a marked difference in the clinical severity of the infection in these two cases. Specifically, they found that the appearance of DVGs early in the infection reduced the severity of the illness. Other recent studies have examined the correlation between viral load and clinical severity, with mixed results. Some studies have found no link between single viral load measures and severity [54,55] while others noted a correlation [56,57,58,59]. More detailed challenge studies note that there appears to be a time delay between rise of viral load and rise of symptoms [60,61], suggesting that there might be some intermediary steps linking viral load and symptoms.
There has been interest in using DVGs to potentially treat viral infections [51,52]. This treatment strategy has been tested in a number of animal models [62,63,64,65,66,67,68,69]. Our results indicate that the timing of treatment initiation could be an important factor since the timing of the appearance of DVGs is associated with different changes in viral kinetics parameters. Duration of production of the DVGs could also be an important factor in the efficacy of treatment since prolonged persistence of DVGs was associated with a lower infection rate and faster clearance of virus. While the idea is that a single injection of DVGs is sufficient because the DVGs should continue to replicate and will persist, this might not always happen since there were some patients in the Felt et al. study who had transient DVGs [9]. Transient DVGs showed no statistically significant difference in infection rate or clearance rate from no DVG patients in this study, suggesting that DVGs need to persist for some period of time in order to affect the time course of the infection.
This study used a very simple model of viral infection that neglects many details of the infection. Our understanding of how DVGs alter the time course of infections would be helped by fitting models that incorporate DVGs [11,33,34,36,37,38] to clinical data. In order to properly parameterize such a model, more data needs to be collected. In addition to viral time courses, we also need the time course of DVGs. A model that explicitly incorporates DVGs could provide insight into how late appearance of DVGs leads to lower infection rates, higher viral loads and more rapid clearance. The explicit inclusion of DVGs will also allow for estimation of parameters related to DVGs such as the DVG production rate, DVG infection rate, and DVG clearance rate.
Despite the simplicity of the model, we might not have found a global minimum and there was still difficulty in identifying all the parameters. It is known that for this simple model δ and c are not both uniquely identifiable from viral titer data alone [46]. We note that some of the confidence intervals are fairly broad, so there is room for improvement on the identifiability of the parameters. In order to improve parameter estimates measurement of other quantities will be needed. For example, measuring both infectious virus and RNA copies can narrow the confidence intervals around parameters estimates [70], while measurement of infected cells in addition to viral measurements can help separate the effects of δ and c [71].
In spite of these limitations, our application of a mathematical model to challenge infections of RSV found a robust link between the appearance and duration of DVGs and the infection rate. Further studies to refine both the model and the parameter estimates will lead to a better understanding of causal link between the appearance of DVGs and the infection rate.
The authors declare there is no conflict of interest.
[1] | [ M. Becker,J. Joseph, AIDS and behavioral change to reduce risk: A review, Amer. J. Public Health, 78 (1988): 394-410. |
[2] | [ R. J. Blendon,J. M. Benson,C. M. DesRoches,E. Raleigh,K. Taylor-Clark, The public's response to severe acute respiratory syndormw in Toronto and the United States, Clinic Infectious Diseases, 38 (2004): 925-931. |
[3] | [ S. M. Blower,A. R. McLean, Prophylactic vaccines, risk behavior change, and the probability of eradicating HIV in San Francisco, Science, 265 (1994): 1451-1454. |
[4] | [ F. Brauer and C. Castillo-Chavez, Mathematical Models in Population Biology and Epidemics, Springer-Verlag, New York, 2001. |
[5] | [ S. Busenberg and K. Cooke, Vertically Transmitted Diseases, Springer-Verlag, New York, 1993. |
[6] | [ L. Cai,X. Li, Analysis of a SEIV epidemic model with a nonlinear incidence rate, Applied Mathematical Modelling, 33 (2009): 2919-2926. |
[7] | [ V. Capasso,G. Serio, A generalisation of the Kermack-McKendrick deterministic epidemic model, Math. Biosci., 42 (1978): 43-61. |
[8] | [ R. Chunara,J. R. Andrews,J. S. Brownstein, Social and news media enable estimation of epidemiological patterns early in the 2010 Haitian cholera outbreak, Am. J. Trop. Med. Hyg., 86 (2012): 39-45. |
[9] | [ S. Collinson and J. M. Heffernan, Modelling the effects of media during an influenza epidemic, BMC Public Health, 14 (2014), 376. |
[10] | [ J. Cui,Y. Sun,H. Zhu, The impact of media on the control of infectious diseases, Journal of Dynamics and Differential Equations, 20 (2008): 31-53. |
[11] | [ J. Cui,X. Tao,H. Zhu, An SIS infection model incorporating media coverage, Rocky Mountain Journal of Mathematics, 38 (2008): 1323-1334. |
[12] | [ O. Diekmann and J. Heesterbeek, Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation, Wiley, New York, 2000. |
[13] | [ K. A. Fanti,E. Vanman,C. C. Henrich,M. N. Avraamide, Desensitization to media violence over a short period of time, Aggressive Behavior, 35 (2009): 179-187. |
[14] | [ K. Frost,E. Frank,E. Maibach, Relative risk in the news media: A quantification of misrepresentation, American Journal of Public Health, 87 (1997): 842-845. |
[15] | [ H. Hethcote, The mathematics of infectious diseases, SIAM Revi., 42 (2000): 599-653. |
[16] | [ X. Lai,S. Liu,R. Lin, Rich dynamical behaviors for predator-prey model with weak Allee effect, Applicable Analysis, 89 (2010): 1271-1292. |
[17] | [ J. LaSalle, The Stability of Dynamical Systems, Reg. Conf. Ser. Appl. Math., SIAM, Philadelphia, 1976. |
[18] | [ J. LaSalle and S. Lefschetz, Stability by Lyapunov's Direct Method, Academic Press, New York, 1961. |
[19] | [ G. M. Leung, T. H. Lam, L. M. Ho, S. Y. Ho, B. H. Chang, I. O. Wong and A. J. Hedley, The impac of community sphychological response on outbreak control for severe acute respiratory syndrome in Hong Kong, J. Epid. Comm. Health 2003, p995. |
[20] | [ S. Levin, T. Hallam and L. Gross, Applied Mathematical Ecology, Lecture Notes in Biomathematics, vol. 18, Springer, Berlin, Germany, 1989. |
[21] | [ B. Li, S. Liu, J. Cui and J. Li, A simple predator-prey population model with rich dynamics, App. Sci., 6 (2016), 151. |
[22] | [ R. Lin,S. Liu,X. Lai, Bifurcations of a Predator-prey System with Weak Allee effects, J. Korean Math. Soc., 50 (2013): 695-713. |
[23] | [ R. Liu,J. Wu,H. Zhu, Media/psychological impact on multiple outbreaks of emerging infectious diseases, Computational and Mathematical Methods in Medicine, 8 (2007): 153-164. |
[24] | [ W. Liu,H. Hethcote,S. Levin, Dynamical behaviour of epidemiological models with nonlinear incidence rates, J. Math. Biol., 25 (1987): 359-380. |
[25] | [ W. Liu,S. Levin,Y. Iwasa, Influence of nonlinear incidence rates upon the behaviour of SIRS epidemiological models, J. Math. Biol., 23 (1986): 187-204. |
[26] | [ Y. Liu,J. Cui, The impact of media coverage on the dynamics of infectious disease, International Journal of Biomathematics, 1 (2008): 65-74. |
[27] | [ A. K. Misra,A. Sharma,J. Li, A mathematical model for control of vector borne diseases through media campaigns, Discret. Contin. Dyn. S., 18 (2013): 1909-1927. |
[28] | [ A. K. Misra,A. Sharma,J. B. Shukla, Modeling and analysis of effects of awareness programs by media on the spread of infectious diseases, Math. Comput. Model, 53 (2011): 1221-1228. |
[29] | [ A. K. Misra,A. Sharma,V. Singh, Effect of awareness programs in controlling the prevalence of an epidemic with time delay, J. Biol. Sys., 19 (2011): 389-402. |
[30] | [ A. Mummert and H. Weiss, Get the news out loudly and quickly: The Influence of the media on limiting emerging infectious disease outbreaks, PloS One, 8 (2013), e71692. |
[31] | [ J. Murray, Mathematical Biology, Springer-Verlag, Berlin, 1989. |
[32] | [ C. Sun,W. Yang,J. Arino,K. Khan, Effect of media-induced social distancing on disease transmission in a two patch setting, Math. Biosci., 230 (2011): 87-95. |
[33] | [ J. Tchuenche, N. Dube, C. Bhunu, R. Smith and C. Bauch, The impact of media coverage on the transmission dynamics of human inuenza, BMC Public Health, 11 (2011), S5. |
[34] | [ P. van den Driessche,J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of diseases transmission, Math. Biosci., 180 (2002): 29-48. |
[35] | [ A. A. F. Wahlberg,L. Sjoberg, Risk perception and the media, Journal of Risk Research, 3 (2001): 31-50. |
[36] | [ Y. Xiao, S. Tang and J. Wu, Media impact switching surface during an infectious disease outbreak, Scientific Reports, 5 (2015), 7838. |
[37] | [ Y. Xiao,T. Zhao,S. Tang, Dynamics of an infectious diseases with media/psychology induced non-smooth incidence, Math. Biosci. and Eng., 10 (2013): 445-461. |
[38] | [ M. E. Young, G. R. Norman and K. R. Humphreys, Medicine in the popular press: The Influence of the media on perceptions of disease, PloS One, (3) (2008), e3552. |
[39] | [ X. Zhao, Uniform persistence and periodic coexistence states in infinite-dimensional periodic semiflow with applications, Can. Appl. Math. Q., 3 (1995): 473-495. |
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.93×10−6 | 1.77×109 | 26.9 | 181 | 1.07 | 0.0200 | 1.82 |
95% CI | (0.569–221)×10−6 | 3.76×104–1.01×1015 | 0.864–401 | 0.614–96.1 | 1.02–6.82 | 1.96×10−5–0.506 | |
Early | 4.38×10−5 | 5.06×106 | 7.62 | 23.1 | 1.26 | 0.0950 | 2.09 |
95% CI | (0.106–42.5)×10−5 | 2.44×104–5.78×1014 | 1.28–2066 | 1.64–472 | 1.01–3.52 | 3.88×10−5–0.442 | |
Late | 2.21×10−8 | 4.95×1012 | 226 | 476 | 1.01 | 0.00428 | 3.26 |
95% CI | (2.10–186)×10−8 | (0.000309–12.6)×1012 | 2.34–408 | 50.6–1350 | 1.02–2.35 | 0.00279–0.0703 | |
p values | |||||||
Early/DVG- | 0.389 | 0.356 | 0.247 | 0.076 | 0.405 | 0.431 | |
Early/Late | 0.0005 | 0.098 | 0.365 | 0.008 | 0.440 | 0.269 | |
DVG-/Late | 0.001 | 0.478 | 0.141 | 5.82×10−5 | 0.436 | 0.475 |
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 3.10×10−5 | 3.66×107 | 2.41 | 251 | 1.88 | 0.0420 | 0.755 |
95% CI | (0.0948–6.86)×10−5 | 2.56×106–1.53×1016 | 1.53–558 | 16.5–3880 | 1.00–2.48 | 9.41×10−6–0.117 | |
Late | 3.27×10−7 | 5.48×109 | 9.01 | 151 | 1.32 | 0.0334 | 0.821 |
95% CI | (0.453–12.3)×10−6 | (0.262–1280)×109 | 3.47–159 | 31.2–365 | 1.03–1.86 | 0.00653–0.0838 | |
p values | |||||||
DVG-/Late | 0.015 | 0.483 | 0.228 | 0.521 | 0.167 | 0.281 |
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.20×10−5 | 1.82×107 | 3.25 | 78.0 | 1.58 | 0.0707 | 1.47 |
95% CI | (0.0667–24.8)×10−5 | (0.00348–100)×107 | 0.771–1840 | 0.542–92.3 | 1.02–6.57 | 8.26×10−6–0.521 | |
Early | 5.00×10−6 | 2.19×109 | 31.7 | 326 | 1.06 | 0.0135 | 2.06 |
95% CI | (0.508–180)×10−6 | (0.000128–572)×109 | 1.58–1760 | 2.70–869 | 1.01–2.77 | 4.10×10−6–0.296 | |
Late | 1.68×10−8 | 1.31×1013 | 69.0 | 3060 | 1.04 | 0.00301 | 2.26 |
95% CI | (0.864–91.3)×10−8 | 9.15×108–3.92×1013 | 2.67–476 | 57.7–2080 | 1.01–2.24 | 0.00257–0.0562 | |
p values | |||||||
Early/DVG- | 0.441 | 0.488 | 0.514 | 0.005 | 0.426 | 0.476 | |
Early/Late | 0.005 | 0.238 | 0.365 | 0.024 | 0.468 | 0.446 | |
DVG-/Late | 7.44×10−5 | 0.216 | 0.197 | 2.17×10−5 | 0.503 | 0.478 |
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.55×10−7 | 1.89×1012 | 387 | 1240 | 1.01 | 0.00200 | 2.21 |
95% CI | (0.377–2400)×10−7 | 4.25×104–1.47×1017 | 1.61–486 | 0.853–3910 | 1.00–3.98 | 2.76×10−6–0.455 | |
Transient | 1.01×10−5 | 1.78×108 | 11.8 | 131 | 1.16 | 0.0334 | 1.53 |
95% CI | (0.0654–20.2)×10−5 | 5.50×104–1.08×1016 | 1.17–1817 | 0.961–139 | 1.01–4.82 | 7.22×10−6–0.461 | |
Prolonged | 1.02×10−8 | 2.18×1013 | 61.1 | 3460 | 1.05 | 0.00300 | 2.83 |
95% CI | (0.496–24.9)×10−8 | 5.35×108–7.91×1014 | 18.1–2600 | 3.56–1330 | 1.01–1.96 | 0.00106–0.0908 | |
p values | |||||||
Transient/DVG- | 0.371 | 0.411 | 0.510 | 0.129 | 0.469 | 0.443 | |
Transient/Prolonged | 1.30×10−5 | 0.206 | 0.249 | 0.001 | 0.493 | 0.519 | |
DVG-/Prolonged | 2.37×10−5 | 0.428 | 0.227 | 0.034 | 0.343 | 0.382 |
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.93×10−6 | 1.77×109 | 26.9 | 181 | 1.07 | 0.0200 | 1.82 |
95% CI | (0.569–221)×10−6 | 3.76×104–1.01×1015 | 0.864–401 | 0.614–96.1 | 1.02–6.82 | 1.96×10−5–0.506 | |
Early | 4.38×10−5 | 5.06×106 | 7.62 | 23.1 | 1.26 | 0.0950 | 2.09 |
95% CI | (0.106–42.5)×10−5 | 2.44×104–5.78×1014 | 1.28–2066 | 1.64–472 | 1.01–3.52 | 3.88×10−5–0.442 | |
Late | 2.21×10−8 | 4.95×1012 | 226 | 476 | 1.01 | 0.00428 | 3.26 |
95% CI | (2.10–186)×10−8 | (0.000309–12.6)×1012 | 2.34–408 | 50.6–1350 | 1.02–2.35 | 0.00279–0.0703 | |
p values | |||||||
Early/DVG- | 0.389 | 0.356 | 0.247 | 0.076 | 0.405 | 0.431 | |
Early/Late | 0.0005 | 0.098 | 0.365 | 0.008 | 0.440 | 0.269 | |
DVG-/Late | 0.001 | 0.478 | 0.141 | 5.82×10−5 | 0.436 | 0.475 |
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 3.10×10−5 | 3.66×107 | 2.41 | 251 | 1.88 | 0.0420 | 0.755 |
95% CI | (0.0948–6.86)×10−5 | 2.56×106–1.53×1016 | 1.53–558 | 16.5–3880 | 1.00–2.48 | 9.41×10−6–0.117 | |
Late | 3.27×10−7 | 5.48×109 | 9.01 | 151 | 1.32 | 0.0334 | 0.821 |
95% CI | (0.453–12.3)×10−6 | (0.262–1280)×109 | 3.47–159 | 31.2–365 | 1.03–1.86 | 0.00653–0.0838 | |
p values | |||||||
DVG-/Late | 0.015 | 0.483 | 0.228 | 0.521 | 0.167 | 0.281 |
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.20×10−5 | 1.82×107 | 3.25 | 78.0 | 1.58 | 0.0707 | 1.47 |
95% CI | (0.0667–24.8)×10−5 | (0.00348–100)×107 | 0.771–1840 | 0.542–92.3 | 1.02–6.57 | 8.26×10−6–0.521 | |
Early | 5.00×10−6 | 2.19×109 | 31.7 | 326 | 1.06 | 0.0135 | 2.06 |
95% CI | (0.508–180)×10−6 | (0.000128–572)×109 | 1.58–1760 | 2.70–869 | 1.01–2.77 | 4.10×10−6–0.296 | |
Late | 1.68×10−8 | 1.31×1013 | 69.0 | 3060 | 1.04 | 0.00301 | 2.26 |
95% CI | (0.864–91.3)×10−8 | 9.15×108–3.92×1013 | 2.67–476 | 57.7–2080 | 1.01–2.24 | 0.00257–0.0562 | |
p values | |||||||
Early/DVG- | 0.441 | 0.488 | 0.514 | 0.005 | 0.426 | 0.476 | |
Early/Late | 0.005 | 0.238 | 0.365 | 0.024 | 0.468 | 0.446 | |
DVG-/Late | 7.44×10−5 | 0.216 | 0.197 | 2.17×10−5 | 0.503 | 0.478 |
Group | β (day)−1⋅(copies/ml)−1 | p (copies/ml)⋅(day)−1 | δ (day)−1 | c (day)−1 | R0 | tinf (day) | SSR |
DVG- | 2.55×10−7 | 1.89×1012 | 387 | 1240 | 1.01 | 0.00200 | 2.21 |
95% CI | (0.377–2400)×10−7 | 4.25×104–1.47×1017 | 1.61–486 | 0.853–3910 | 1.00–3.98 | 2.76×10−6–0.455 | |
Transient | 1.01×10−5 | 1.78×108 | 11.8 | 131 | 1.16 | 0.0334 | 1.53 |
95% CI | (0.0654–20.2)×10−5 | 5.50×104–1.08×1016 | 1.17–1817 | 0.961–139 | 1.01–4.82 | 7.22×10−6–0.461 | |
Prolonged | 1.02×10−8 | 2.18×1013 | 61.1 | 3460 | 1.05 | 0.00300 | 2.83 |
95% CI | (0.496–24.9)×10−8 | 5.35×108–7.91×1014 | 18.1–2600 | 3.56–1330 | 1.01–1.96 | 0.00106–0.0908 | |
p values | |||||||
Transient/DVG- | 0.371 | 0.411 | 0.510 | 0.129 | 0.469 | 0.443 | |
Transient/Prolonged | 1.30×10−5 | 0.206 | 0.249 | 0.001 | 0.493 | 0.519 | |
DVG-/Prolonged | 2.37×10−5 | 0.428 | 0.227 | 0.034 | 0.343 | 0.382 |