In this paper, we discuss a selection of tools from dynamical systems and order statistics, which are most often utilized separately, and combine them into an algorithm to estimate the parameters of mathematical models for infectious diseases in the case of small sample sizes and left censoring, which is relevant in the case of rapidly evolving infectious diseases and remote populations. The proposed method relies on the analogy between survival functions and the dynamics of the susceptible compartment in SIR-type models, which are both monotone decreasing in time and are both determined by a dual variable: the hazard function in survival prediction and the number of infected people in SIR-type models. We illustrate the methodology in the case of a continuous model in the presence of noisy measurements with different distributions (Normal, Poisson, Negative Binomial) and in a discrete model, reminiscent of the Ricker map, which admits chaotic dynamics. This estimation procedure shows stable results in experiments based on a popular benchmark dataset for SIR-type models and small samples. This manuscript illustrates how classical theoretical statistical methods and dynamical systems can be merged in interesting ways to study problems ranging from more fundamental small sample situations to more complex infectious disease and survival models, with the potential that this tools can be applied in the presence of a large number of covariates and different types of censored data.
- Award ID(s):
- 1920945
- PAR ID:
- 10292219
- Date Published:
- Journal Name:
- PeerJ Physical Chemistry
- Volume:
- 2
- ISSN:
- 2689-7733
- Page Range / eLocation ID:
- e14
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
We explore the effects of cross-diffusion dynamics in epidemiological models. Using reaction–diffusion models of infectious disease, we explicitly consider situations where an individual in a category will move according to the concentration of individuals in other categories. Namely, we model susceptible individuals moving away from infected and infectious individuals. Here, we show that including these cross-diffusion dynamics results in a delay in the onset of an epidemic and an increase in the total number of infectious individuals. This representation provides more realistic spatiotemporal dynamics of the disease classes in an Erlang SEIR model and allows us to study how spatial mobility, due to social behavior, can affect the spread of an epidemic. We found that tailored control measures, such as targeted testing, contact tracing, and isolation of infected individuals, can be more effective in mitigating the spread of infectious diseases while minimizing the negative impact on society and the economy.more » « less
-
more » « less
Abstract Background An understanding of epidemiological dynamics, once confined to mathematical epidemiologists and applied mathematicians, can be disseminated to a non-mathematical community of health care professionals and applied biologists through simple-to-use simulation applications. We used Numerus Model Builder RAMP
Ⓡ (Runtime Alterable Model Platform) technology, to construct deterministic and stochastic versions of compartmental SIR (Susceptible, Infectious, Recovered with immunity) models as simple-to-use, freely available, epidemic simulation application programs.Results We take the reader through simulations used to demonstrate the following concepts: 1) disease prevalence curves of unmitigated outbreaks have a single peak and result in epidemics that ‘burn’ through the population to become extinguished when the proportion of the susceptible population drops below a critical level; 2) if immunity in recovered individuals wanes sufficiently fast then the disease persists indefinitely as an endemic state, with possible dampening oscillations following the initial outbreak phase; 3) the steepness and initial peak of the prevalence curve are influenced by the basic reproductive value
R 0, which must exceed 1 for an epidemic to occur; 4) the probability that a single infectious individual in a closed population (i.e. no migration) gives rise to an epidemic increases with the value ofR 0>1; 5) behavior that adaptively decreases the contact rate among individuals with increasing prevalence has major effects on the prevalence curve including dramatic flattening of the prevalence curve along with the generation of multiple prevalence peaks; 6) the impacts of treatment are complicated to model because they effect multiple processes including transmission, recovery and mortality; 7) the impacts of vaccination policies, constrained by a fixed number of vaccination regimens and by the rate and timing of delivery, are crucially important to maximizing the ability of vaccination programs to reduce mortality.Conclusion Our presentation makes transparent the key assumptions underlying SIR epidemic models. Our RAMP simulators are meant to augment rather than replace classroom material when teaching epidemiological dynamics. They are sufficiently versatile to be used by students to address a range of research questions for term papers and even dissertations.
-
We study the effect of population mobility on the transmission dynamics of infectious diseases by considering a susceptible-exposed-infectious-recovered (SEIR) epidemic model with graph Laplacian diffusion, that is, on a weighted network. First, we establish the existence and uniqueness of solutions to the SEIR model defined on a weighed graph. Then by constructing Liapunov functions, we show that the disease-free equilibrium is globally asymptotically stable if the basic reproduction number is less than unity and the endemic equilibrium is globally asymptotically stable if the basic reproduction number is greater than unity. Finally, we apply our generalized weighed graph to Watts–Strogatz network and carry out numerical simulations, which demonstrate that degrees of nodes determine peak numbers of the infectious population as well as the time to reach these peaks. It also indicates that the network has an impact on the transient dynamical behaviour of the epidemic transmission.more » « less
-
The spread of an infectious disease depends on intrinsic properties of the disease as well as the connectivity and actions of the population. This study investigates the dynamics of an SIR type model which accounts for human tendency to avoid infection while also maintaining preexisting, interpersonal relationships. Specifically, we use a network model in which individuals probabilistically deactivate connections to infected individuals and later reconnect to the same individuals upon recovery. To analyze this network model, a mean field approximation consisting of a system of fourteen ordinary differential equations for the number of nodes and edges is developed. This system of equations is closed using a moment closure approximation for the number of triple links. By analyzing the differential equations, it is shown that, in addition to force of infection and recovery rate, the probability of deactivating edges and the average node degree of the underlying network determine if an epidemic occurs.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.7717/peerj-pchem.14
