Matching Items (14)
Description
In the field of infectious disease epidemiology, the assessment of model robustness outcomes plays a significant role in the identification, reformulation, and evaluation of preparedness strategies aimed at limiting the impact of catastrophic events (pandemics or the deliberate release of biological agents) or used in the management of disease prevention strategies, or employed in the identification and evaluation of control or mitigation measures. The research work in this dissertation focuses on: The comparison and assessment of the role of exponentially distributed waiting times versus the use of generalized non-exponential parametric distributed waiting times of infectious periods on the quantitative and qualitative outcomes generated by Susceptible-Infectious-Removed (SIR) models. Specifically, Gamma distributed infectious periods are considered in the three research projects developed following the applications found in (Bailey 1964, Anderson 1980, Wearing 2005, Feng 2007, Feng 2007, Yan 2008, lloyd 2009, Vergu 2010). i) The first project focuses on the influence of input model parameters, such as the transmission rate, mean and variance of Gamma distributed infectious periods, on disease prevalence, the peak epidemic size and its timing, final epidemic size, epidemic duration and basic reproduction number. Global uncertainty and sensitivity analyses are carried out using a deterministic Susceptible-Infectious-Recovered (SIR) model. The quantitative effect and qualitative relation between input model parameters and outcome variables are established using Latin Hypercube Sampling (LHS) and Partial rank correlation coefficient (PRCC) and Spearman rank correlation coefficient (RCC) sensitivity indices. We learnt that: For relatively low (R0 close to one) to high (mean of R0 equals 15) transmissibility, the variance of the Gamma distribution for the infectious period, input parameter of the deterministic age-of-infection SIR model, is key (statistically significant) on the predictability of the epidemiological variables such as the epidemic duration and the peak size and timing of the prevalence of infectious individuals and therefore, for the predictability these variables, it is preferable to utilize a nonlinear system of Volterra integral equations, rather than a nonlinear system of ordinary differential equations. The predictability of epidemiological variables such as the final epidemic size and the basic reproduction number are unaffected by (or independent of) the variance of the Gamma distribution for the infectious period and therefore for the choice on which type of nonlinear system for the description of the SIR model (VIE's or ODE's) is irrelevant. Although, for practical proposes, with the aim of lowering the complexity and number operations in the numerical methods, a nonlinear system of ordinary differential equations is preferred. The main contribution lies in the development of a model based decision-tool that helps determine when SIR models given in terms of Volterra integral equations are equivalent or better suited than SIR models that only consider exponentially distributed infectious periods. ii) The second project addresses the question of whether or not there is sufficient evidence to conclude that two empirical distributions for a single epidemiological outcome, one generated using a stochastic SIR model under exponentially distributed infectious periods and the other under the non-exponentially distributed infectious period, are statistically dissimilar. The stochastic formulations are modeled via a continuous time Markov chain model. The statistical hypothesis test is conducted using the non-parametric Kolmogorov-Smirnov test. We found evidence that shows that for low to moderate transmissibility, all empirical distribution pairs (generated from exponential and non-exponential distributions) for each of the epidemiological quantities considered are statistically dissimilar. The research in this project helps determine whether the weakening exponential distribution assumption must be considered in the estimation of probability of events defined from the empirical distribution of specific random variables. iii) The third project involves the assessment of the effect of exponentially distributed infectious periods on estimates of input parameter and the associated outcome variable predictions. Quantities unaffected by the use of exponentially distributed infectious period within low transmissibility scenarios include, the prevalence peak time, final epidemic size, epidemic duration and basic reproduction number and for high transmissibility scenarios only the prevalence peak time and final epidemic size. An application designed to determine from incidence data whether there is sufficient statistical evidence to conclude that the infectious period distribution should not be modeled by an exponential distribution is developed. A method for estimating explicitly specified non-exponential parametric probability density functions for the infectious period from epidemiological data is developed. The methodologies presented in this dissertation may be applicable to models where waiting times are used to model transitions between stages, a process that is common in the study of life-history dynamics of many ecological systems.
ContributorsMorales Butler, Emmanuel J (Author) / Castillo-Chavez, Carlos (Thesis advisor) / Aparicio, Juan P (Thesis advisor) / Camacho, Erika T (Committee member) / Kang, Yun (Committee member) / Arizona State University (Publisher)
Created2014
Description
Diseases have been part of human life for generations and evolve within the population, sometimes dying out while other times becoming endemic or the cause of recurrent outbreaks. The long term influence of a disease stems from different dynamics within or between pathogen-host, that have been analyzed and studied by many researchers using mathematical models. Co-infection with different pathogens is common, yet little is known about how infection with one pathogen affects the host's immunological response to another. Moreover, no work has been found in the literature that considers the variability of the host immune health or that examines a disease at the population level and its corresponding interconnectedness with the host immune system. Knowing that the spread of the disease in the population starts at the individual level, this thesis explores how variability in immune system response within an endemic environment affects an individual's vulnerability, and how prone it is to co-infections. Immunology-based models of Malaria and Tuberculosis (TB) are constructed by extending and modifying existing mathematical models in the literature. The two are then combined to give a single nine-variable model of co-infection with Malaria and TB. Because these models are difficult to gain any insight analytically due to the large number of parameters, a phenomenological model of co-infection is proposed with subsystems corresponding to the individual immunology-based model of a single infection. Within this phenomenological model, the variability of the host immune health is also incorporated through three different pathogen response curves using nonlinear bounded Michaelis-Menten functions that describe the level or state of immune system (healthy, moderate and severely compromised). The immunology-based models of Malaria and TB give numerical results that agree with the biological observations. The Malaria--TB co-infection model gives reasonable results and these suggest that the order in which the two diseases are introduced have an impact on the behavior of both. The subsystems of the phenomenological models that correspond to a single infection (either of Malaria or TB) mimic much of the observed behavior of the immunology-based counterpart and can demonstrate different behavior depending on the chosen pathogen response curve. In addition, varying some of the parameters and initial conditions in the phenomenological model yields a range of topologically different mathematical behaviors, which suggests that this behavior may be able to be observed in the immunology-based models as well. The phenomenological models clearly replicate the qualitative behavior of primary and secondary infection as well as co-infection. The mathematical solutions of the models correspond to the fundamental states described by immunologists: virgin state, immune state and tolerance state. The phenomenological model of co-infection also demonstrates a range of parameter values and initial conditions in which the introduction of a second disease causes both diseases to grow without bound even though those same parameters and initial conditions did not yield unbounded growth in the corresponding subsystems. This results applies to all three states of the host immune system. In terms of the immunology-based system, this would suggest the following: there may be parameter values and initial conditions in which a person can clear Malaria or TB (separately) from their system but in which the presence of both can result in the person dying of one of the diseases. Finally, this thesis studies links between epidemiology (population level) and immunology in an effort to assess the impact of pathogen's spread within the population on the immune response of individuals. Models of Malaria and TB are proposed that incorporate the immune system of the host into a mathematical model of an epidemic at the population level.
ContributorsSoho, Edmé L (Author) / Wirkus, Stephen (Thesis advisor) / Castillo-Chavez, Carlos (Thesis advisor) / Chowell-Puente, Gerardo (Committee member) / Arizona State University (Publisher)
Created2011
Description
Mathematical models are important tools for addressing problems that exceed experimental capabilities. In this work, I present ordinary and partial differential equation (ODE, PDE) models for two problems: Vicodin abuse and impact cratering.
The prescription opioid Vicodin is the nation's most widely prescribed pain reliever. The majority of Vicodin abusers are first introduced via prescription, distinguishing it from other drugs in which the most common path to abuse begins with experimentation. I develop and analyze two mathematical models of Vicodin use and abuse, considering only those patients with an initial Vicodin prescription. Through adjoint sensitivity analysis, I show that focusing efforts on prevention rather than treatment has greater success at reducing the total population of abusers. I prove that solutions to each model exist, are unique, and are non-negative. I also derive conditions for which these solutions are asymptotically stable.
Verification and Validation (V&V) are necessary processes to ensure accuracy of computational methods. Simulations are essential for addressing impact cratering problems, because these problems often exceed experimental capabilities. I show that the Free Lagrange (FLAG) hydrocode, developed and maintained by Los Alamos National Laboratory, can be used for impact cratering simulations by verifying FLAG against two analytical models of aluminum-on-aluminum impacts at different impact velocities and validating FLAG against a glass-into-water laboratory impact experiment. My verification results show good agreement with the theoretical maximum pressures, and my mesh resolution study shows that FLAG converges at resolutions low enough to reduce the required computation time from about 28 hours to about 25 minutes.
Asteroid 16 Psyche is the largest M-type (metallic) asteroid in the Main Asteroid Belt. Radar albedo data indicate Psyche's surface is rich in metallic content, but estimates for Psyche's composition vary widely. Psyche has two large impact structures in its Southern hemisphere, with estimated diameters from 50 km to 70 km and estimated depths up to 6.4 km. I use the FLAG hydrocode to model the formation of the largest of these impact structures. My results indicate an oblique angle of impact rather than a vertical impact. These results also support previous claims that Psyche is metallic and porous.
The prescription opioid Vicodin is the nation's most widely prescribed pain reliever. The majority of Vicodin abusers are first introduced via prescription, distinguishing it from other drugs in which the most common path to abuse begins with experimentation. I develop and analyze two mathematical models of Vicodin use and abuse, considering only those patients with an initial Vicodin prescription. Through adjoint sensitivity analysis, I show that focusing efforts on prevention rather than treatment has greater success at reducing the total population of abusers. I prove that solutions to each model exist, are unique, and are non-negative. I also derive conditions for which these solutions are asymptotically stable.
Verification and Validation (V&V) are necessary processes to ensure accuracy of computational methods. Simulations are essential for addressing impact cratering problems, because these problems often exceed experimental capabilities. I show that the Free Lagrange (FLAG) hydrocode, developed and maintained by Los Alamos National Laboratory, can be used for impact cratering simulations by verifying FLAG against two analytical models of aluminum-on-aluminum impacts at different impact velocities and validating FLAG against a glass-into-water laboratory impact experiment. My verification results show good agreement with the theoretical maximum pressures, and my mesh resolution study shows that FLAG converges at resolutions low enough to reduce the required computation time from about 28 hours to about 25 minutes.
Asteroid 16 Psyche is the largest M-type (metallic) asteroid in the Main Asteroid Belt. Radar albedo data indicate Psyche's surface is rich in metallic content, but estimates for Psyche's composition vary widely. Psyche has two large impact structures in its Southern hemisphere, with estimated diameters from 50 km to 70 km and estimated depths up to 6.4 km. I use the FLAG hydrocode to model the formation of the largest of these impact structures. My results indicate an oblique angle of impact rather than a vertical impact. These results also support previous claims that Psyche is metallic and porous.
ContributorsCaldwell, Wendy K (Author) / Wirkus, Stephen (Thesis advisor) / Asphaug, Erik (Committee member) / Camacho, Erika T (Committee member) / Crook, Sharon (Committee member) / Plesko, Catherine S (Committee member) / Smith, Hal (Committee member) / Arizona State University (Publisher)
Created2019
Description
Pharmacokinetics describes the movement and processing of a drug in a body, while Pharmacodynamics describes the drug's effect on a given subject. Pharmacokinetic/Pharmacodynamic(Pk/Pd) models have become a fundamental tool when predicting bacterial behavior and drug development. In November of 2009, Katsube et al. published their paper detailing their Pk/Pd model for the drug Doripenem and the bacteria P. aeruginosa. In their paper, they determined that there is a dependent relationship between the drug's effectiveness and the dosing strategy of the drug. Therefore, this thesis has applied optimal control in order to optimize the drug's effectiveness, while not burdening the subject with the side effects of the drug. Optimal Control is a mathematical tool used to balance two competing factors. As a result, it has become a useful tool used to make decisions involving complex behavior. By using Optimal Control, the model will maximize the drug's effect on the bacterial population of P. aeruginosa, while minimizing the drug concentration of Doripenem. In doing so, our research will enable doctors and clinicians to maximize a drug's effectiveness on the body, while minimizing side effects.
ContributorsSawkins, Bryan Thomas (Author) / Camacho, Erika (Thesis director) / Wirkus, Stephen (Committee member) / School of Mathematical and Natural Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2018-05
Description
The retina is the lining in the back of the eye responsible for vision. When light photons hits the retina, the photoreceptors within the retina respond by sending impulses to the optic nerve, which connects to the brain. If there is injury to the eye or heredity retinal problems, this part can become detached. Detachment leads to loss of nutrients, such as oxygen and glucose, to the cells in the eye and causes cell death. Sometimes the retina is able to be surgically reattached. If the photoreceptor cells have not died and the reattachment is successful, then these cells are able to regenerate their outer segments (OS) which are essential for their functionality and vitality. In this work we will explore how the regrowth of the photoreceptor cells in a healthy eye after retinal detachment can lead to a deeper understanding of how eye cells take up nutrients and regenerate. This work uses a mathematical model for a healthy eye in conjunction with data for photoreceptors' regrowth and decay. The parameters for the healthy eye model are estimated from the data and the ranges of these parameter values are centered +/- 10\% away from these values are used for sensitivity analysis. Using parameter estimation and sensitivity analysis we can better understand how certain processes represented by these parameters change within the model as a result of retinal detachment. Having a deeper understanding for any sort of photoreceptor death and growth can be used by the greater scientific community to help with these currently irreversible conditions that lead to blindness, such as retinal detachment. The analysis in this work shows that maximizing the carrying capacity of the trophic pool and the rate of RDCVF, as well as minimizing nutrient withdrawal of the rods and the cones from the trophic pool results in both the most regrowth and least cell death in retinal detachment.
ContributorsGoldman, Miriam Ayla (Author) / Camacho, Erikia (Thesis director) / Wirkus, Stephen (Committee member) / School of Mathematical and Natural Sciences (Contributor, Contributor) / Barrett, The Honors College (Contributor)
Created2017-05
Description
Infectious diseases are a leading cause of death worldwide. With the development of drugs, vaccines and antibiotics, it was believed that for the first time in human history diseases would no longer be a major cause of mortality. Newly emerging diseases, re-emerging diseases and the emergence of microorganisms resistant to existing treatment have forced us to re-evaluate our optimistic perspective. In this study, a simple mathematical framework for super-infection is considered in order to explore the transmission dynamics of drug-resistance. Through its theoretical analysis, we identify the conditions necessary for the coexistence between sensitive strains and drug-resistant strains. Farther, in order to investigate the effectiveness of control measures, the model is extended so as to include vaccination and treatment. The impact that these preventive and control measures may have on its disease dynamics is evaluated. Theoretical results being confirmed via numerical simulations. Our theoretical results on two-strain drug-resistance models are applied in the context of Malaria, antimalarial drugs, and the administration of a possible partially effective vaccine. The objective is to develop a monitoring epidemiological framework that help evaluate the impact of antimalarial drugs and partially-effective vaccine in reducing the disease burden at the population level. Optimal control theory is applied in the context of this framework in order to assess the impact of time dependent cost-effective treatment efforts. It is shown that cost-effective combinations of treatment efforts depend on the population size, cost of implementing treatment controls, and the parameters of the model. We use these results to identify optimal control strategies for several scenarios.
ContributorsUrdapilleta, Alicia (Author) / Castillo-Chavez, Carlos (Thesis advisor) / Wang, Xiaohong (Thesis advisor) / Wirkus, Stephen (Committee member) / Camacho, Erika (Committee member) / Arizona State University (Publisher)
Created2011
Description
Adderall remains to be one of the most commonly abused drug among college campuses. Although it is a prescription drug that is primarily used to treat attention deficit hyperactivity disorder (ADHD), it has become one of the many "study drugs" due to its usage among college students during stressful school times, such as exams, where increased concentration and energy levels are thought to improve work efficiency. However, Adderall is notable for having a high potential for abuse and a risk of psychological and physical side effects. In this paper, we conducted a mathematical analysis on an existing epidemiological model of Adderall abuse. We started by verifying the positivity of solutions using techniques from dynamical systems because this is a population model dealing with people. Then, we found and investigated different equilibrium solutions to analyze their stability using both analytical and graphical approaches. Finally, the results were tied back into the Adderall model and conclusions were drawn.
ContributorsKerseg, Cassidy (Author) / Wirkus, Stephen (Thesis director) / Brager, Danielle (Committee member) / School of Mathematical and Natural Sciences (Contributor, Contributor) / Barrett, The Honors College (Contributor)
Created2021-05
Description
This outlines a mathematical model created in MATLAB for the purposes of predicting nitrous oxide emissions from wastewater treatment plants with updated an updated understanding of AOB metabolic pathway.
ContributorsOverbey, Jorja (Author) / Hart, Steven (Thesis director) / Young, Michelle (Committee member) / Wirkus, Stephen (Committee member) / Barrett, The Honors College (Contributor) / School of Music, Dance and Theatre (Contributor) / School of Mathematical and Natural Sciences (Contributor)
Created2023-05
Description
Acetylcholinesterase (AChE) inhibition by chemical toxicants such as organophosphates, nerve agents, and carbamates can lead to a series of adverse health outcomes including seizures, coma, and death. An adverse outcome pathway (AOP) is a framework that describes a series of biologically measurable key events (KEs) leading from some molecular initiating event (MIE) to an adverse outcome (AO) of regulatory significance, all developed and hosted in the AOP Wiki. A quantitative AOP (qAOP) is a mathematical model that predicts how perturbations in the MIE affect KEs based on the key event relationships (KERs) that define the AOP. The purpose of this thesis was to expand upon the KERs that define the AOP for AChE inhibition leading to neurodegeneration in order to better understand the effects of AChE inhibitors and the risks they pose to ecosystems, wildlife, and human health. In order to reduce the resources and time spent for chemical toxicity testing, a qAOP was developed based on the available quantitative data and models that supported the AOP. A literature review for the collection of qualitative evidence and quantitative data in support of the AOP was performed resulting in further expansion of the relationships between key events (KERs) through construction of additional KER description pages. A model evaluation was performed by comparing the qAOP model predictions with experimental data, with a subsequent sensitivity analysis of unknown parameters. The qAOP model simulates the MIE through its fifth KE (KE 5) and KE 7. Model predictions compared to experimentally measured data either under- or overpredicting multiple KEs warranting additional refinement such as a formal parameter optimization. Overall, more data amenable to qAOP model development are needed. To aid qAOP model development, the presentation of data in the AOPWiki may be improved by presenting the quantitative data in the AOP Wiki in a tabular format and allowing for the hosting of mathematical models or raw data. With these recommendations in mind, and through continued AOP construction in the AOP Wiki, new qAOP models will be developed, ultimately supporting chemical risk assessment and the mitigation of effects upon exposed individuals and wildlife populations.
ContributorsSinitsyn, Dennis (Author) / Watanabe, Karen (Thesis advisor) / Vinas, Natalia (Committee member) / Wirkus, Stephen (Committee member) / Arizona State University (Publisher)
Created2023
Description
A leading crisis in the United States is the opioid use disorder (OUD) epidemic. Opioid overdose deaths have been increasing, with over 100,000 deaths due to overdose from April 2020 to April 2021. This dissertation presents two mathematical models to address illicit OUD (IOUD), treatment, and recovery within an epidemiological framework. In the first model, individuals remain in the recovery class unless they relapse. Due to the limited availability of specialty treatment facilities for individuals with OUD, a saturation treat- ment function was incorporated. The second model is an extension of the first, where a casual user class and its corresponding specialty treatment class were added. Using U.S. population data, the data was scaled to a population of 200,000 to find parameter estimates. While the first model used the heroin-only dataset, the second model used both the heroin and all-illicit opioids datasets. Backward bifurcation was found in the first IOUD model for realistic parameter values. Additionally, bistability was observed in the second IOUD model with the heroin-only dataset. This result implies that it would be beneficial to increase the availability of treatment. An alarming effect was discovered about the high overdose death rate: by 2038, the disease-free equilibrium would be the only stable equilibrium. This consequence is concerning because although the goal is for the epidemic to end, it would be preferable to end it through treatment rather than overdose. The IOUD model with a casual user class, its sensitivity results, and the comparison of parameters for both datasets, showed the importance of not overlooking the influence that casual users have in driving the all-illicit opioid epidemic. Casual users stay in the casual user class longer and are not going to treatment as quickly as the users of the heroin epidemic. Another result was that the users of the all-illicit opioids were going to the recovered class by means other than specialty treatment. However, the relapse rates for those individuals were much more significant than in the heroin-only epidemic. The results above from analyzing these models may inform health and policy officials, leading to more effective treatment options and prevention efforts.
ContributorsCole, Sandra (Author) / Wirkus, Stephen (Thesis advisor) / Gardner, Carl (Committee member) / Lanchier, Nicolas (Committee member) / Camacho, Erika (Committee member) / Fricks, John (Committee member) / Arizona State University (Publisher)
Created2022