Matching Items (6)
Filtering by
- All Subjects: Parameter estimation
- Creators: Kuang, Yang
Description
Predicting resistant prostate cancer is critical for lowering medical costs and improving the quality of life of advanced prostate cancer patients. I formulate, compare, and analyze two mathematical models that aim to forecast future levels of prostate-specific antigen (PSA). I accomplish these tasks by employing clinical data of locally advanced prostate cancer patients undergoing androgen deprivation therapy (ADT). I demonstrate that the inverse problem of parameter estimation might be too complicated and simply relying on data fitting can give incorrect conclusions, since there is a large error in parameter values estimated and parameters might be unidentifiable. I provide confidence intervals to give estimate forecasts using data assimilation via an ensemble Kalman Filter. Using the ensemble Kalman Filter, I perform dual estimation of parameters and state variables to test the prediction accuracy of the models. Finally, I present a novel model with time delay and a delay-dependent parameter. I provide a geometric stability result to study the behavior of this model and show that the inclusion of time delay may improve the accuracy of predictions. Also, I demonstrate with clinical data that the inclusion of the delay-dependent parameter facilitates the identification and estimation of parameters.
ContributorsBaez, Javier (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric (Committee member) / Crook, Sharon (Committee member) / Gardner, Carl (Committee member) / Nagy, John (Committee member) / Arizona State University (Publisher)
Created2017
Description
In recent decades, marine ecologists have conducted extensive field work and experiments to understand the interactions between bacteria and bacteriophage (phage) in marine environments. This dissertation provides a detailed rigorous framework for gaining deeper insight into these interactions. Specific features of the dissertation include the design of a new deterministic Lotka-Volterra model with n + 1 bacteria, n
+ 1 phage, with explicit nutrient, where the jth phage strain infects the first j bacterial strains, a perfectly nested infection network (NIN). This system is subject to trade-off conditions on the life-history traits of both bacteria and phage given in an earlier study Jover et al. (2013). Sufficient conditions are provided to show that a bacteria-phage community of arbitrary size with NIN can arise through the succession of permanent subcommunities, by the successive addition of one new population. Using uniform persistence theory, this entire community is shown to be permanent (uniformly persistent), meaning that all populations ultimately survive.
It is shown that a modified version of the original NIN Lotka-Volterra model with implicit nutrient considered by Jover et al. (2013) is permanent. A new one-to-one infection network (OIN) is also considered where each bacterium is infected by only one phage, and that phage infects only that bacterium. This model does not use the trade-offs on phage infection range, and bacterium resistance to phage. The OIN model is shown to be permanent, and using Lyapunov function theory, coupled with LaSalle’s Invariance Principle, the unique coexistence equilibrium associated with the NIN is globally asymptotically stable provided that the inter- and intra-specific bacterial competition coefficients are equal across all bacteria.
Finally, the OIN model is extended to a “Kill the Winner” (KtW) Lotka-Volterra model
of marine communities consisting of bacteria, phage, and zooplankton. The zooplankton
acts as a super bacteriophage, which infects all bacteria. This model is shown to be permanent.
+ 1 phage, with explicit nutrient, where the jth phage strain infects the first j bacterial strains, a perfectly nested infection network (NIN). This system is subject to trade-off conditions on the life-history traits of both bacteria and phage given in an earlier study Jover et al. (2013). Sufficient conditions are provided to show that a bacteria-phage community of arbitrary size with NIN can arise through the succession of permanent subcommunities, by the successive addition of one new population. Using uniform persistence theory, this entire community is shown to be permanent (uniformly persistent), meaning that all populations ultimately survive.
It is shown that a modified version of the original NIN Lotka-Volterra model with implicit nutrient considered by Jover et al. (2013) is permanent. A new one-to-one infection network (OIN) is also considered where each bacterium is infected by only one phage, and that phage infects only that bacterium. This model does not use the trade-offs on phage infection range, and bacterium resistance to phage. The OIN model is shown to be permanent, and using Lyapunov function theory, coupled with LaSalle’s Invariance Principle, the unique coexistence equilibrium associated with the NIN is globally asymptotically stable provided that the inter- and intra-specific bacterial competition coefficients are equal across all bacteria.
Finally, the OIN model is extended to a “Kill the Winner” (KtW) Lotka-Volterra model
of marine communities consisting of bacteria, phage, and zooplankton. The zooplankton
acts as a super bacteriophage, which infects all bacteria. This model is shown to be permanent.
ContributorsKorytowski, Daniel (Author) / Smith, Hal (Thesis advisor) / Gumel, Abba (Committee member) / Kuang, Yang (Committee member) / Gardner, Carl (Committee member) / Thieme, Horst (Committee member) / Arizona State University (Publisher)
Created2016
Description
Using a simple $SI$ infection model, I uncover the
overall dynamics of the system and how they depend on the incidence
function. I consider both an epidemic and endemic perspective of the
model, but in both cases, three classes of incidence
functions are identified.
In the epidemic form,
power incidences, where the infective portion $I^p$ has $p\in(0,1)$,
cause unconditional host extinction,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction. The case of non-extinction in upper
density-dependent
incidences extends to the case where a latent period is included.
Using data from experiments with rhanavirus and salamanders,
maximum likelihood estimates are applied to the data.
With these estimates,
I generate the corrected Akaike information criteria, which
reward a low likelihood and punish the use of more parameters.
This generates the Akaike weight, which is used to fit
parameters to the data, and determine which incidence functions
fit the data the best.
From an endemic perspective, I observe
that power incidences cause initial condition dependent host extinction for
some parameter constellations and global stability for others,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction.
The dynamics when the incidence function is homogeneous are deeply explored.
I expand the endemic considerations in the homogeneous case
by adding a predator into the model.
Using persistence theory, I show the conditions for the persistence of each of the
predator, prey, and parasite species. Potential dynamics of the system include parasite mediated
persistence of the predator, survival of the ecosystem at high initial predator levels and
ecosystem collapse at low initial predator levels, persistence of all three species, and much more.
overall dynamics of the system and how they depend on the incidence
function. I consider both an epidemic and endemic perspective of the
model, but in both cases, three classes of incidence
functions are identified.
In the epidemic form,
power incidences, where the infective portion $I^p$ has $p\in(0,1)$,
cause unconditional host extinction,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction. The case of non-extinction in upper
density-dependent
incidences extends to the case where a latent period is included.
Using data from experiments with rhanavirus and salamanders,
maximum likelihood estimates are applied to the data.
With these estimates,
I generate the corrected Akaike information criteria, which
reward a low likelihood and punish the use of more parameters.
This generates the Akaike weight, which is used to fit
parameters to the data, and determine which incidence functions
fit the data the best.
From an endemic perspective, I observe
that power incidences cause initial condition dependent host extinction for
some parameter constellations and global stability for others,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction.
The dynamics when the incidence function is homogeneous are deeply explored.
I expand the endemic considerations in the homogeneous case
by adding a predator into the model.
Using persistence theory, I show the conditions for the persistence of each of the
predator, prey, and parasite species. Potential dynamics of the system include parasite mediated
persistence of the predator, survival of the ecosystem at high initial predator levels and
ecosystem collapse at low initial predator levels, persistence of all three species, and much more.
ContributorsFarrell, Alexander E. (Author) / Thieme, Horst R (Thesis advisor) / Smith, Hal (Committee member) / Kuang, Yang (Committee member) / Tang, Wenbo (Committee member) / Collins, James (Committee member) / Arizona State University (Publisher)
Created2017
Description
A description of numerical and analytical work pertaining to models that describe the growth and progression of glioblastoma multiforme (GBM), an aggressive form of primary brain cancer. Two reaction-diffusion models are used: the Fisher-Kolmogorov-Petrovsky-Piskunov equation and a 2-population model that divides the tumor into actively proliferating and quiescent (or necrotic) cells. The numerical portion of this work (chapter 2) focuses on simulating GBM expansion in patients undergoing treatment for recurrence of tumor following initial surgery. The models are simulated on 3-dimensional brain geometries derived from magnetic resonance imaging (MRI) scans provided by the Barrow Neurological Institute. The study consists of 17 clinical time intervals across 10 patients that have been followed in detail, each of whom shows significant progression of tumor over a period of 1 to 3 months on sequential follow up scans. A Taguchi sampling design is implemented to estimate the variability of the predicted tumors to using 144 different choices of model parameters. In 9 cases, model parameters can be identified such that the simulated tumor contains at least 40 percent of the volume of the observed tumor. In the analytical portion of the paper (chapters 3 and 4), a positively invariant region for our 2-population model is identified. Then, a rigorous derivation of the critical patch size associated with the model is performed. The critical patch (KISS) size is the minimum habitat size needed for a population to survive in a region. Habitats larger than the critical patch size allow a population to persist, while smaller habitats lead to extinction. The critical patch size of the 2-population model is consistent with that of the Fisher-Kolmogorov-Petrovsky-Piskunov equation, one of the first reaction-diffusion models proposed for GBM. The critical patch size may indicate that GBM tumors have a minimum size depending on the location in the brain. A theoretical relationship between the size of a GBM tumor at steady-state and its maximum cell density is also derived, which has potential applications for patient-specific parameter estimation based on magnetic resonance imaging data.
ContributorsHarris, Duane C. (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric J. (Thesis advisor) / Preul, Mark C. (Committee member) / Crook, Sharon (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2023
Description
Efforts to treat prostate cancer have seen an uptick, as the world’s most commoncancer in men continues to have increasing global incidence. Clinically, metastatic
prostate cancer is most commonly treated with hormonal therapy. The idea behind
hormonal therapy is to reduce androgen production, which prostate cancer cells
require for growth. Recently, the exploration of the synergistic effects of the drugs
used in hormonal therapy has begun. The aim was to build off of these recent
advancements and further refine the synergistic drug model. The advancements I
implement come by addressing biological shortcomings and improving the model’s
internal mechanistic structure. The drug families being modeled, anti-androgens,
and gonadotropin-releasing hormone analogs, interact with androgen production in a
way that is not completely understood in the scientific community. Thus the models
representing the drugs show progress through their ability to capture their effect
on serum androgen. Prostate-specific antigen is the primary biomarker for prostate
cancer and is generally how population models on the subject are validated. Fitting
the model to clinical data and comparing it to other clinical models through the
ability to fit and forecast prostate-specific antigen and serum androgen is how this
improved model achieves validation. The improved model results further suggest that
the drugs’ dynamics should be considered in adaptive therapy for prostate cancer.
prostate cancer is most commonly treated with hormonal therapy. The idea behind
hormonal therapy is to reduce androgen production, which prostate cancer cells
require for growth. Recently, the exploration of the synergistic effects of the drugs
used in hormonal therapy has begun. The aim was to build off of these recent
advancements and further refine the synergistic drug model. The advancements I
implement come by addressing biological shortcomings and improving the model’s
internal mechanistic structure. The drug families being modeled, anti-androgens,
and gonadotropin-releasing hormone analogs, interact with androgen production in a
way that is not completely understood in the scientific community. Thus the models
representing the drugs show progress through their ability to capture their effect
on serum androgen. Prostate-specific antigen is the primary biomarker for prostate
cancer and is generally how population models on the subject are validated. Fitting
the model to clinical data and comparing it to other clinical models through the
ability to fit and forecast prostate-specific antigen and serum androgen is how this
improved model achieves validation. The improved model results further suggest that
the drugs’ dynamics should be considered in adaptive therapy for prostate cancer.
ContributorsReckell, Trevor (Author) / Kostelich, Eric (Thesis advisor) / Kuang, Yang (Committee member) / Mahalov, Alex (Committee member) / Arizona State University (Publisher)
Created2020
Description
Cancer is a disease involving abnormal growth of cells. Its growth dynamics is perplexing. Mathematical modeling is a way to shed light on this progress and its medical treatments. This dissertation is to study cancer invasion in time and space using a mathematical approach. Chapter 1 presents a detailed review of literature on cancer modeling.
Chapter 2 focuses sorely on time where the escape of a generic cancer out of immune control is described by stochastic delayed differential equations (SDDEs). Without time delay and noise, this system demonstrates bistability. The effects of response time of the immune system and stochasticity in the tumor proliferation rate are studied by including delay and noise in the model. Stability, persistence and extinction of the tumor are analyzed. The result shows that both time delay and noise can induce the transition from low tumor burden equilibrium to high tumor equilibrium. The aforementioned work has been published (Han et al., 2019b).
In Chapter 3, Glioblastoma multiforme (GBM) is studied using a partial differential equation (PDE) model. GBM is an aggressive brain cancer with a grim prognosis. A mathematical model of GBM growth with explicit motility, birth, and death processes is proposed. A novel method is developed to approximate key characteristics of the wave profile, which can be compared with MRI data. Several test cases of MRI data of GBM patients are used to yield personalized parameterizations of the model. The aforementioned work has been published (Han et al., 2019a).
Chapter 4 presents an innovative way of forecasting spatial cancer invasion. Most mathematical models, including the ones described in previous chapters, are formulated based on strong assumptions, which are hard, if not impossible, to verify due to complexity of biological processes and lack of quality data. Instead, a nonparametric forecasting method using Gaussian processes is proposed. By exploiting the local nature of the spatio-temporal process, sparse (in terms of time) data is sufficient for forecasting. Desirable properties of Gaussian processes facilitate selection of the size of the local neighborhood and computationally efficient propagation of uncertainty. The method is tested on synthetic data and demonstrates promising results.
Chapter 2 focuses sorely on time where the escape of a generic cancer out of immune control is described by stochastic delayed differential equations (SDDEs). Without time delay and noise, this system demonstrates bistability. The effects of response time of the immune system and stochasticity in the tumor proliferation rate are studied by including delay and noise in the model. Stability, persistence and extinction of the tumor are analyzed. The result shows that both time delay and noise can induce the transition from low tumor burden equilibrium to high tumor equilibrium. The aforementioned work has been published (Han et al., 2019b).
In Chapter 3, Glioblastoma multiforme (GBM) is studied using a partial differential equation (PDE) model. GBM is an aggressive brain cancer with a grim prognosis. A mathematical model of GBM growth with explicit motility, birth, and death processes is proposed. A novel method is developed to approximate key characteristics of the wave profile, which can be compared with MRI data. Several test cases of MRI data of GBM patients are used to yield personalized parameterizations of the model. The aforementioned work has been published (Han et al., 2019a).
Chapter 4 presents an innovative way of forecasting spatial cancer invasion. Most mathematical models, including the ones described in previous chapters, are formulated based on strong assumptions, which are hard, if not impossible, to verify due to complexity of biological processes and lack of quality data. Instead, a nonparametric forecasting method using Gaussian processes is proposed. By exploiting the local nature of the spatio-temporal process, sparse (in terms of time) data is sufficient for forecasting. Desirable properties of Gaussian processes facilitate selection of the size of the local neighborhood and computationally efficient propagation of uncertainty. The method is tested on synthetic data and demonstrates promising results.
ContributorsHan, Lifeng (Author) / Kuang, Yang (Thesis advisor) / Fricks, John (Thesis advisor) / Kostelich, Eric (Committee member) / Baer, Steve (Committee member) / Gumel, Abba (Committee member) / Arizona State University (Publisher)
Created2020