Matching Items (38)
Filtering by
- Genre: Doctoral Dissertation
Description
Parallel Monte Carlo applications require the pseudorandom numbers used on each processor to be independent in a probabilistic sense. The TestU01 software package is the standard testing suite for detecting stream dependence and other properties that make certain pseudorandom generators ineffective in parallel (as well as serial) settings. TestU01 employs two basic schemes for testing parallel generated streams. The first applies serial tests to the individual streams and then tests the resulting P-values for uniformity. The second turns all the parallel generated streams into one long vector and then applies serial tests to the resulting concatenated stream. Various forms of stream dependence can be missed by each approach because neither one fully addresses the multivariate nature of the accumulated data when generators are run in parallel. This dissertation identifies these potential faults in the parallel testing methodologies of TestU01 and investigates two different methods to better detect inter-stream dependencies: correlation motivated multivariate tests and vector time series based tests. These methods have been implemented in an extension to TestU01 built in C++ and the unique aspects of this extension are discussed. A variety of different generation scenarios are then examined using the TestU01 suite in concert with the extension. This enhanced software package is found to better detect certain forms of inter-stream dependencies than the original TestU01 suites of tests.
ContributorsIsmay, Chester (Author) / Eubank, Randall (Thesis advisor) / Young, Dennis (Committee member) / Kao, Ming-Hung (Committee member) / Lanchier, Nicolas (Committee member) / Reiser, Mark R. (Committee member) / Arizona State University (Publisher)
Created2013
Description
Mathematical modeling of infectious diseases can help public health officials to make decisions related to the mitigation of epidemic outbreaks. However, over or under estimations of the morbidity of any infectious disease can be problematic. Therefore, public health officials can always make use of better models to study the potential implication of their decisions and strategies prior to their implementation. Previous work focuses on the mechanisms underlying the different epidemic waves observed in Mexico during the novel swine origin influenza H1N1 pandemic of 2009 and showed extensions of classical models in epidemiology by adding temporal variations in different parameters that are likely to change during the time course of an epidemic, such as, the influence of media, social distancing, school closures, and how vaccination policies may affect different aspects of the dynamics of an epidemic. This current work further examines the influence of different factors considering the randomness of events by adding stochastic processes to meta-population models. I present three different approaches to compare different stochastic methods by considering discrete and continuous time. For the continuous time stochastic modeling approach I consider the continuous-time Markov chain process using forward Kolmogorov equations, for the discrete time stochastic modeling I consider stochastic differential equations using Wiener's increment and Poisson point increments, and also I consider the discrete-time Markov chain process. These first two stochastic modeling approaches will be presented in a one city and two city epidemic models using, as a base, our deterministic model. The last one will be discussed briefly on a one city SIS and SIR-type model.
ContributorsCruz-Aponte, Maytee (Author) / Wirkus, Stephen A. (Thesis advisor) / Castillo-Chavez, Carlos (Thesis advisor) / Camacho, Erika T. (Committee member) / Kang, Yun (Committee member) / Arizona State University (Publisher)
Created2014
Description
Rabies disease remains enzootic among raccoons, skunks, foxes and bats in the United States. It is of primary concern for public-health agencies to control spatial spread of rabies in wildlife and its potential spillover infection of domestic animals and humans. Rabies is invariably fatal in wildlife if untreated, with a non-negligible incubation period. Understanding how this latency affects spatial spread of rabies in wildlife is the concern of chapter 2 and 3. Chapter 1 deals with the background of mathematical models for rabies and lists main objectives. In chapter 2, a reaction-diffusion susceptible-exposed-infected (SEI) model and a delayed diffusive susceptible-infected (SI) model are constructed to describe the same epidemic process -- rabies spread in foxes. For the delayed diffusive model a non-local infection term with delay is resulted from modeling the dispersal during incubation stage. Comparison is made regarding minimum traveling wave speeds of the two models, which are verified using numerical experiments. In chapter 3, starting with two Kermack and McKendrick's models where infectivity, death rate and diffusion rate of infected individuals can depend on the age of infection, the asymptotic speed of spread $c^\ast$ for the cumulated force of infection can be analyzed. For the special case of fixed incubation period, the asymptotic speed of spread is governed by the same integral equation for both models. Although explicit solutions for $c^\ast$ are difficult to obtain, assuming that diffusion coefficient of incubating animals is small, $c^\ast$ can be estimated in terms of model parameter values. Chapter 4 considers the implementation of realistic landscape in simulation of rabies spread in skunks and bats in northeast Texas. The Finite Element Method (FEM) is adopted because the irregular shapes of realistic landscape naturally lead to unstructured grids in the spatial domain. This implementation leads to a more accurate description of skunk rabies cases distributions.
ContributorsLiu, Hao (Author) / Kuang, Yang (Thesis advisor) / Jackiewicz, Zdzislaw (Committee member) / Lanchier, Nicolas (Committee member) / Smith, Hal (Committee member) / Thieme, Horst (Committee member) / Arizona State University (Publisher)
Created2013
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
There has been important progress in understanding ecological dynamics through the development of the theory of ecological stoichiometry. This fast growing theory provides new constraints and mechanisms that can be formulated into mathematical models. Stoichiometric models incorporate the effects of both food quantity and food quality into a single framework that produce rich dynamics. While the effects of nutrient deficiency on consumer growth are well understood, recent discoveries in ecological stoichiometry suggest that consumer dynamics are not only affected by insufficient food nutrient content (low phosphorus (P): carbon (C) ratio) but also by excess food nutrient content (high P:C). This phenomenon, known as the stoichiometric knife edge, in which animal growth is reduced not only by food with low P content but also by food with high P content, needs to be incorporated into mathematical models. Here we present Lotka-Volterra type models to investigate the growth response of Daphnia to algae of varying P:C ratios. Using a nonsmooth system of two ordinary differential equations (ODEs), we formulate the first model to incorporate the phenomenon of the stoichiometric knife edge. We then extend this stoichiometric model by mechanistically deriving and tracking free P in the environment. This resulting full knife edge model is a nonsmooth system of three ODEs. Bifurcation analysis and numerical simulations of the full model, that explicitly tracks phosphorus, leads to quantitatively different predictions than previous models that neglect to track free nutrients. The full model shows that the grazer population is sensitive to excess nutrient concentrations as a dynamical free nutrient pool induces extreme grazer population density changes. These modeling efforts provide insight on the effects of excess nutrient content on grazer dynamics and deepen our understanding of the effects of stoichiometry on the mechanisms governing population dynamics and the interactions between trophic levels.
ContributorsPeace, Angela (Author) / Kuang, Yang (Thesis advisor) / Elser, James J (Committee member) / Baer, Steven (Committee member) / Tang, Wenbo (Committee member) / Kang, Yun (Committee member) / Arizona State University (Publisher)
Created2014
Description
In 1968, phycologist M.R. Droop published his famous discovery on the functional relationship between growth rate and internal nutrient status of algae in chemostat culture. The simple notion that growth is directly dependent on intracellular nutrient concentration is useful for understanding the dynamics in many ecological systems. The cell quota in particular lends itself to ecological stoichiometry, which is a powerful framework for mathematical ecology. Three models are developed based on the cell quota principal in order to demonstrate its applications beyond chemostat culture.
First, a data-driven model is derived for neutral lipid synthesis in green microalgae with respect to nitrogen limitation. This model synthesizes several established frameworks in phycology and ecological stoichiometry. The model demonstrates how the cell quota is a useful abstraction for understanding the metabolic shift to neutral lipid production that is observed in certain oleaginous species.
Next a producer-grazer model is developed based on the cell quota model and nutrient recycling. The model incorporates a novel feedback loop to account for animal toxicity due to accumulation of nitrogen waste. The model exhibits rich, complex dynamics which leave several open mathematical questions.
Lastly, disease dynamics in vivo are in many ways analogous to those of an ecosystem, giving natural extensions of the cell quota concept to disease modeling. Prostate cancer can be modeled within this framework, with androgen the limiting nutrient and the prostate and cancer cells as competing species. Here the cell quota model provides a useful abstraction for the dependence of cellular proliferation and apoptosis on androgen and the androgen receptor. Androgen ablation therapy is often used for patients in biochemical recurrence or late-stage disease progression and is in general initially effective. However, for many patients the cancer eventually develops resistance months to years after treatment begins. Understanding how and predicting when hormone therapy facilitates evolution of resistant phenotypes has immediate implications for treatment. Cell quota models for prostate cancer can be useful tools for this purpose and motivate applications to other diseases.
First, a data-driven model is derived for neutral lipid synthesis in green microalgae with respect to nitrogen limitation. This model synthesizes several established frameworks in phycology and ecological stoichiometry. The model demonstrates how the cell quota is a useful abstraction for understanding the metabolic shift to neutral lipid production that is observed in certain oleaginous species.
Next a producer-grazer model is developed based on the cell quota model and nutrient recycling. The model incorporates a novel feedback loop to account for animal toxicity due to accumulation of nitrogen waste. The model exhibits rich, complex dynamics which leave several open mathematical questions.
Lastly, disease dynamics in vivo are in many ways analogous to those of an ecosystem, giving natural extensions of the cell quota concept to disease modeling. Prostate cancer can be modeled within this framework, with androgen the limiting nutrient and the prostate and cancer cells as competing species. Here the cell quota model provides a useful abstraction for the dependence of cellular proliferation and apoptosis on androgen and the androgen receptor. Androgen ablation therapy is often used for patients in biochemical recurrence or late-stage disease progression and is in general initially effective. However, for many patients the cancer eventually develops resistance months to years after treatment begins. Understanding how and predicting when hormone therapy facilitates evolution of resistant phenotypes has immediate implications for treatment. Cell quota models for prostate cancer can be useful tools for this purpose and motivate applications to other diseases.
ContributorsPacker, Aaron (Author) / Kuang, Yang (Thesis advisor) / Nagy, John (Committee member) / Smith, Hal (Committee member) / Kostelich, Eric (Committee member) / Kang, Yun (Committee member) / Arizona State University (Publisher)
Created2014
Description
Pre-Exposure Prophylaxis (PrEP) is any medical or public health procedure used before exposure to the disease causing agent, its purpose is to prevent, rather than treat or cure a disease. Most commonly, PrEP refers to an experimental HIV-prevention strategy that would use antiretrovirals to protect HIV-negative people from HIV infection. A deterministic mathematical model of HIV transmission is developed to evaluate the public-health impact of oral PrEP interventions, and to compare PrEP effectiveness with respect to different evaluation methods. The effects of demographic, behavioral, and epidemic parameters on the PrEP impact are studied in a multivariate sensitivity analysis. Most of the published models on HIV intervention impact assume that the number of individuals joining the sexually active population per year is constant or proportional to the total population. In the second part of this study, three models are presented and analyzed to study the PrEP intervention, with constant, linear, and logistic recruitment rates. How different demographic assumptions can affect the evaluation of PrEP is studied. When provided with data, often least square fitting or similar approaches can be used to determine a single set of approximated parameter values that make the model fit the data best. However, least square fitting only provides point estimates and does not provide information on how strongly the data supports these particular estimates. Therefore, in the third part of this study, Bayesian parameter estimation is applied on fitting ODE model to the related HIV data. Starting with a set of prior distributions for the parameters as initial guess, Bayes' formula can be applied to obtain a set of posterior distributions for the parameters which makes the model fit the observed data best. Evaluating the posterior distribution often requires the integration of high-dimensional functions, which is usually difficult to calculate numerically. Therefore, the Markov chain Monte Carlo (MCMC) method is used to approximate the posterior distribution.
ContributorsZhao, Yuqin (Author) / Kuang, Yang (Thesis advisor) / Taylor, Jesse (Committee member) / Armbruster, Dieter (Committee member) / Tang, Wenbo (Committee member) / Kang, Yun (Committee member) / Arizona State University (Publisher)
Created2014
Description
By the von Neumann min-max theorem, a two person zero sum game with finitely many pure strategies has a unique value for each player (summing to zero) and each player has a non-empty set of optimal mixed strategies. If the payoffs are independent, identically distributed (iid) uniform (0,1) random variables, then with probability one, both players have unique optimal mixed strategies utilizing the same number of pure strategies with positive probability (Jonasson 2004). The pure strategies with positive probability in the unique optimal mixed strategies are called saddle squares. In 1957, Goldman evaluated the probability of a saddle point (a 1 by 1 saddle square), which was rediscovered by many authors including Thorp (1979). Thorp gave two proofs of the probability of a saddle point, one using combinatorics and one using a beta integral. In 1965, Falk and Thrall investigated the integrals required for the probabilities of a 2 by 2 saddle square for 2 × n and m × 2 games with iid uniform (0,1) payoffs, but they were not able to evaluate the integrals. This dissertation generalizes Thorp's beta integral proof of Goldman's probability of a saddle point, establishing an integral formula for the probability that a m × n game with iid uniform (0,1) payoffs has a k by k saddle square (k ≤ m,n). Additionally, the probabilities of a 2 by 2 and a 3 by 3 saddle square for a 3 × 3 game with iid uniform(0,1) payoffs are found. For these, the 14 integrals observed by Falk and Thrall are dissected into 38 disjoint domains, and the integrals are evaluated using the basic properties of the dilogarithm function. The final results for the probabilities of a 2 by 2 and a 3 by 3 saddle square in a 3 × 3 game are linear combinations of 1, π2, and ln(2) with rational coefficients.
ContributorsManley, Michael (Author) / Kadell, Kevin W. J. (Thesis advisor) / Kao, Ming-Hung (Committee member) / Lanchier, Nicolas (Committee member) / Lohr, Sharon (Committee member) / Reiser, Mark R. (Committee member) / Arizona State University (Publisher)
Created2011
Description
In this dissertation the potential impact of some social, cultural and economic factors on
Ebola Virus Disease (EVD) dynamics and control are studied. In Chapter two, the inability
to detect and isolate a large fraction of EVD-infected individuals before symptoms onset is
addressed. A mathematical model, calibrated with data from the 2014 West African outbreak,
is used to show the dynamics of EVD control under various quarantine and isolation
effectiveness regimes. It is shown that in order to make a difference it must reach a high
proportion of the infected population. The effect of EVD-dead bodies has been incorporated
in the quarantine effectiveness. In Chapter four, the potential impact of differential
risk is assessed. A two-patch model without explicitly incorporate quarantine is used to
assess the impact of mobility on communities at risk of EVD. It is shown that the
overall EVD burden may lessen when mobility in this artificial high-low risk society is allowed.
The cost that individuals in the low-risk patch must pay, as measured by secondary
cases is highlighted. In Chapter five a model explicitly incorporating patch-specific quarantine
levels is used to show that quarantine a large enough proportion of the population
under effective isolation leads to a measurable reduction of secondary cases in the presence
of mobility. It is shown that sharing limited resources can improve the effectiveness of
EVD effective control in the two-patch high-low risk system. Identifying the conditions
under which the low-risk community would be willing to accept the increases in EVD risk,
needed to reduce the total number of secondary cases in a community composed of two
patches with highly differentiated risks has not been addressed. In summary, this dissertation
looks at EVD dynamics within an idealized highly polarized world where resources
are primarily in the hands of a low-risk community – a community of lower density, higher
levels of education and reasonable health services – that shares a “border” with a high-risk
community that lacks minimal resources to survive an EVD outbreak.
Ebola Virus Disease (EVD) dynamics and control are studied. In Chapter two, the inability
to detect and isolate a large fraction of EVD-infected individuals before symptoms onset is
addressed. A mathematical model, calibrated with data from the 2014 West African outbreak,
is used to show the dynamics of EVD control under various quarantine and isolation
effectiveness regimes. It is shown that in order to make a difference it must reach a high
proportion of the infected population. The effect of EVD-dead bodies has been incorporated
in the quarantine effectiveness. In Chapter four, the potential impact of differential
risk is assessed. A two-patch model without explicitly incorporate quarantine is used to
assess the impact of mobility on communities at risk of EVD. It is shown that the
overall EVD burden may lessen when mobility in this artificial high-low risk society is allowed.
The cost that individuals in the low-risk patch must pay, as measured by secondary
cases is highlighted. In Chapter five a model explicitly incorporating patch-specific quarantine
levels is used to show that quarantine a large enough proportion of the population
under effective isolation leads to a measurable reduction of secondary cases in the presence
of mobility. It is shown that sharing limited resources can improve the effectiveness of
EVD effective control in the two-patch high-low risk system. Identifying the conditions
under which the low-risk community would be willing to accept the increases in EVD risk,
needed to reduce the total number of secondary cases in a community composed of two
patches with highly differentiated risks has not been addressed. In summary, this dissertation
looks at EVD dynamics within an idealized highly polarized world where resources
are primarily in the hands of a low-risk community – a community of lower density, higher
levels of education and reasonable health services – that shares a “border” with a high-risk
community that lacks minimal resources to survive an EVD outbreak.
ContributorsEspinoza Cortes, Baltazar (Author) / Castillo-Chavez, Carlos (Thesis advisor) / Kang, Yun (Committee member) / Safan, Muntaser (Committee member) / Arizona State University (Publisher)
Created2018
Description
This dissertation explores the impact of environmental dependent risk on disease dynamics within a Lagrangian modeling perspective; where the identity (defined by place of residency) of individuals is preserved throughout the epidemic process. In Chapter Three, the impact of individuals who refuse to be vaccinated is explored. MMR vaccination and birth rate data from the State of California are used to determine the impact of the anti-vaccine movement on the dynamics of growth of the anti-vaccine sub-population. Dissertation results suggest that under realistic California social dynamics scenarios, it is not possible to revert the influence of anti-vaccine
contagion. In Chapter Four, the dynamics of Zika virus are explored in two highly distinct idealized environments defined by a parameter that models highly distinctive levels of risk, the result of vector and host density and vector control measures. The underlying assumption is that these two communities are intimately connected due to economics with the impact of various patterns of mobility being incorporated via
the use of residency times. In short, a highly heterogeneous community is defined by its risk of acquiring a Zika infection within one of two "spaces," one lacking access to health services or effective vector control policies (lack of resources or ignored due to high levels of crime, or poverty, or both). Low risk regions are defined as those with access to solid health facilities and where vector control measures are implemented routinely. It was found that the better connected these communities are, the existence of communities where mobility between risk regions is not hampered, lower the overall, two patch Zika prevalence. Chapter Five focuses on the dynamics of tuberculosis (TB), a communicable disease, also on an idealized high-low risk set up. The impact of mobility within these two highly distinct TB-risk environments on the dynamics and control of this disease is systematically explored. It is found that collaboration and mobility, under some circumstances, can reduce the overall TB burden.
contagion. In Chapter Four, the dynamics of Zika virus are explored in two highly distinct idealized environments defined by a parameter that models highly distinctive levels of risk, the result of vector and host density and vector control measures. The underlying assumption is that these two communities are intimately connected due to economics with the impact of various patterns of mobility being incorporated via
the use of residency times. In short, a highly heterogeneous community is defined by its risk of acquiring a Zika infection within one of two "spaces," one lacking access to health services or effective vector control policies (lack of resources or ignored due to high levels of crime, or poverty, or both). Low risk regions are defined as those with access to solid health facilities and where vector control measures are implemented routinely. It was found that the better connected these communities are, the existence of communities where mobility between risk regions is not hampered, lower the overall, two patch Zika prevalence. Chapter Five focuses on the dynamics of tuberculosis (TB), a communicable disease, also on an idealized high-low risk set up. The impact of mobility within these two highly distinct TB-risk environments on the dynamics and control of this disease is systematically explored. It is found that collaboration and mobility, under some circumstances, can reduce the overall TB burden.
ContributorsMoreno Martínez, Victor Manuel (Author) / Castillo-Chavez, Carlos (Thesis advisor) / Kang, Yun (Committee member) / Mubayi, Anuj (Committee member) / Arizona State University (Publisher)
Created2018