Matching Items (210)
Description
In this research, I surveyed existing methods of characterizing Epilepsy from Electroencephalogram (EEG) data, including the Random Forest algorithm, which was claimed by many researchers to be the most effective at detecting epileptic seizures [7]. I observed that although many papers claimed a detection of >99% using Random Forest, it was not specified “when” the detection was declared within the 23.6 second interval of the seizure event. In this research, I created a time-series procedure to detect the seizure as early as possible within the 23.6 second epileptic seizure window and found that the detection is effective (> 92%) as early as the first few seconds of the epileptic episode. I intend to use this research as a stepping stone towards my upcoming Masters thesis research where I plan to expand the time-series detection mechanism to the pre-ictal stage, which will require a different dataset.
ContributorsBou-Ghazale, Carine (Author) / Lai, Ying-Cheng (Thesis director) / Berisha, Visar (Committee member) / Barrett, The Honors College (Contributor) / Electrical Engineering Program (Contributor)
Created2022-05
Description
Since the 20th century, Arizona has undergone shifts in agricultural practices, driven by urban expansion and crop irrigation regulations. These changes present environmental challenges, altering atmospheric processes and influencing climate dynamics. Given the potential threats of climate change and drought on water availability for agriculture, further modifications in the agricultural landscape are expected. To understand these land use changes and their impact on carbon dynamics, our study quantified aboveground carbon storage in both cultivated and
abandoned agricultural fields. To accomplish this, we employed Python and various geospatial libraries in Jupyter Notebook files, for thorough dataset assembly and visual, quantitative analysis. We focused on nine counties known for high cultivation levels, primarily located in the lower latitudes of Arizona. Our analysis investigated carbon dynamics across not only abandoned and actively cultivated croplands but also neighboring uncultivated land, for which we estimated the extent. Additionally, we compared these trends with those observed in developed land areas.
The findings revealed a hierarchy in aboveground carbon storage, with currently cultivated lands having the lowest levels, followed by abandoned croplands and uncultivated wilderness. However, wilderness areas exhibited significant variation in carbon storage by county compared to cultivated and abandoned lands. Developed lands ranked highest in aboveground carbon storage, with the median value being the highest. Despite county-wide variations, abandoned croplands generally contained more carbon than currently cultivated areas, with adjacent wilderness lands containing even more than both. This trend suggests that cultivating croplands in the region reduces aboveground carbon stores, while abandonment allows for some replenishment, though only to a limited extent. Enhancing carbon stores in Arizona can be achieved through active restoration efforts on abandoned cropland. By promoting native plant regeneration and boosting aboveground carbon levels, these measures are crucial for improving carbon sequestration. We strongly advocate for implementing this step to facilitate the regrowth of native plants and enhance overall carbon storage in the region.
ContributorsGoodwin, Emily (Author) / Eikenberry, Steffen (Thesis director) / Kuang, Yang (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor)
Created2024-05
Description
Glioblastoma Multiforme is a prevalent and aggressive brain tumor. It has an average 5-year survival rate of 6% and average survival time of 14 months. Using patient-specific MRI data from the Barrow Neurological Institute, this thesis investigates the impact of parameter manipulation on reaction-diffusion models for predicting and simulating glioblastoma growth. The study aims to explore key factors influencing tumor morphology and to contribute to enhancing prediction techniques for treatment.
ContributorsShayegan, Tara (Author) / Kostelich, Eric (Thesis director) / Kuang, Yang (Committee member) / Barrett, The Honors College (Contributor) / School of Human Evolution & Social Change (Contributor)
Created2024-05
Description
Predicting nonlinear dynamical systems has been a long-standing challenge in science. This field is currently witnessing a revolution with the advent of machine learning methods. Concurrently, the analysis of dynamics in various nonlinear complex systems continues to be crucial. Guided by these directions, I conduct the following studies. Predicting critical transitions and transient states in nonlinear dynamics is a complex problem. I developed a solution called parameter-aware reservoir computing, which uses machine learning to track how system dynamics change with a driving parameter. I show that the transition point can be accurately predicted while trained in a sustained functioning regime before the transition. Notably, it can also predict if the system will enter a transient state, the distribution of transient lifetimes, and their average before a final collapse, which are crucial for management. I introduce a machine-learning-based digital twin for monitoring and predicting the evolution of externally driven nonlinear dynamical systems, where reservoir computing is exploited. Extensive tests on various models, encompassing optics, ecology, and climate, verify the approach’s effectiveness. The digital twins can extrapolate unknown system dynamics, continually forecast and monitor under non-stationary external driving, infer hidden variables, adapt to different driving waveforms, and extrapolate bifurcation behaviors across varying system sizes. Integrating engineered gene circuits into host cells poses a significant challenge in synthetic biology due to circuit-host interactions, such as growth feedback. I conducted systematic studies on hundreds of circuit structures exhibiting various functionalities, and identified a comprehensive categorization of growth-induced failures. I discerned three dynamical mechanisms behind these circuit failures. Moreover, my comprehensive computations reveal a scaling law between the circuit robustness and the intensity of growth feedback. A class of circuits with optimal robustness is also identified. Chimera states, a phenomenon of symmetry-breaking in oscillator networks, traditionally have transient lifetimes that grow exponentially with system size. However, my research on high-dimensional oscillators leads to the discovery of ’short-lived’ chimera states. Their lifetime increases logarithmically with system size and decreases logarithmically with random perturbations, indicating a unique fragility. To understand these states, I use a transverse stability analysis supported by simulations.
ContributorsKong, Lingwei (Author) / Lai, Ying-Cheng (Thesis advisor) / Tian, Xiaojun (Committee member) / Papandreou-Suppappola, Antonia (Committee member) / Alkhateeb, Ahmed (Committee member) / Arizona State University (Publisher)
Created2023
Description
A notable challenge when assembling synthetic gene circuits is that modularity often fails to function as intended. A crucial underlying reason for this modularity failure is the existence of competition for shared and limited gene expression resources. By designing a synthetic cascading bistable switches (Syn-CBS) circuit in a single strain with two coupled self-activation modules to achieve successive cell fate transitions, nonlinear resource competition within synthetic gene circuits is unveiled. However, in vivo it can be seen that the transition path was redirected with the activation of one switch always prevailing over that of the other, contradictory to coactivation theoretically expected. This behavior is a result of resource competition between genes and follows a ‘winner-takes-all’ rule, where the winner is determined by the relative connection strength between the two modules. Despite investigation demonstrating that resource competition between gene modules can significantly alter circuit deterministic behaviors, how resource competition contributes to gene expression noise and how this noise can be controlled is still an open issue of fundamental importance in systems biology and biological physics. By utilizing a two-gene circuit, the effects of resource competition on protein expression noise levels can be closely studied. A surprising double-edged role is discovered: the competition for these resources decreases noise while the constraint on resource availability adds its own term of noise into the system, denoted “resource competitive” noise. Noise reduction effects are then studied using orthogonal resources. Results indicate that orthogonal resources are a good strategy for eliminating the contribution of resource competition to gene expression noise.
Noise propagation through a cascading circuit has been considered without resource competition. It has been noted that the noise from upstream genes can be transmitted downstream. However, resource competition’s effects on this cascading noise have yet to be studied. When studied, it is found that resource competition can induce stochastic state
switching and perturb noise propagation. Orthogonal resources can remove some of the resource competitive behavior and allow for a system with less noise.
ContributorsGoetz, Hanah Elizabeth (Author) / Tian, Xiaojun (Thesis advisor) / Wang, Xiao (Committee member) / Lai, Ying-Cheng (Committee member) / Arizona State University (Publisher)
Created2022
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
The most advanced social insects, the eusocial insects, form often large societies in which there is reproductive division of labor, queens and workers, have overlapping generations, and cooperative brood care where daughter workers remain in the nest with their queen mother and care for their siblings. The eusocial insects are composed of representative species of bees and wasps, and all species of ants and termites. Much is known about their organizational structure, but remains to be discovered.
The success of social insects is dependent upon cooperative behavior and adaptive strategies shaped by natural selection that respond to internal or external conditions. The objective of my research was to investigate specific mechanisms that have helped shaped the structure of division of labor observed in social insect colonies, including age polyethism and nutrition, and phenomena known to increase colony survival such as egg cannibalism. I developed various Ordinary Differential Equation (ODE) models in which I applied dynamical, bifurcation, and sensitivity analysis to carefully study and visualize biological outcomes in social organisms to answer questions regarding the conditions under which a colony can survive. First, I investigated how the population and evolutionary dynamics of egg cannibalism and division of labor can promote colony survival. I then introduced a model of social conflict behavior to study the inclusion of different response functions that explore the benefits of cannibalistic behavior and how it contributes to age polyethism, the change in behavior of workers as they age, and its biological relevance. Finally, I introduced a model to investigate the importance of pollen nutritional status in a honeybee colony, how it affects population growth and influences division of labor within the worker caste. My results first reveal that both cannibalism and division of labor are adaptive strategies that increase the size of the worker population, and therefore, the persistence of the colony. I show the importance of food collection, consumption, and processing rates to promote good colony nutrition leading to the coexistence of brood and adult workers. Lastly, I show how taking into account seasonality for pollen collection improves the prediction of long term consequences.
The success of social insects is dependent upon cooperative behavior and adaptive strategies shaped by natural selection that respond to internal or external conditions. The objective of my research was to investigate specific mechanisms that have helped shaped the structure of division of labor observed in social insect colonies, including age polyethism and nutrition, and phenomena known to increase colony survival such as egg cannibalism. I developed various Ordinary Differential Equation (ODE) models in which I applied dynamical, bifurcation, and sensitivity analysis to carefully study and visualize biological outcomes in social organisms to answer questions regarding the conditions under which a colony can survive. First, I investigated how the population and evolutionary dynamics of egg cannibalism and division of labor can promote colony survival. I then introduced a model of social conflict behavior to study the inclusion of different response functions that explore the benefits of cannibalistic behavior and how it contributes to age polyethism, the change in behavior of workers as they age, and its biological relevance. Finally, I introduced a model to investigate the importance of pollen nutritional status in a honeybee colony, how it affects population growth and influences division of labor within the worker caste. My results first reveal that both cannibalism and division of labor are adaptive strategies that increase the size of the worker population, and therefore, the persistence of the colony. I show the importance of food collection, consumption, and processing rates to promote good colony nutrition leading to the coexistence of brood and adult workers. Lastly, I show how taking into account seasonality for pollen collection improves the prediction of long term consequences.
ContributorsRodríguez Messan, Marisabel (Author) / Kang, Yun (Thesis advisor) / Castillo-Chavez, Carlos (Thesis advisor) / Kuang, Yang (Committee member) / Page Jr., Robert E (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2018
Description
The role of climate change, as measured in terms of changes in the climatology of geophysical variables (such as temperature and rainfall), on the global distribution and burden of vector-borne diseases (VBDs) remains a subject of considerable debate. This dissertation attempts to contribute to this debate via the use of mathematical (compartmental) modeling and statistical data analysis. In particular, the objective is to find suitable values and/or ranges of the climate variables considered (typically temperature and rainfall) for maximum vector abundance and consequently, maximum transmission intensity of the disease(s) they cause.
Motivated by the fact that understanding the dynamics of disease vector is crucial to understanding the transmission and control of the VBDs they cause, a novel weather-driven deterministic model for the population biology of the mosquito is formulated and rigorously analyzed. Numerical simulations, using relevant weather and entomological data for Anopheles mosquito (the vector for malaria), show that maximum mosquito abundance occurs when temperature and rainfall values lie in the range [20-25]C and [105-115] mm, respectively.
The Anopheles mosquito ecology model is extended to incorporate human dynamics. The resulting weather-driven malaria transmission model, which includes many of the key aspects of malaria (such as disease transmission by asymptomatically-infectious humans, and enhanced malaria immunity due to repeated exposure), was rigorously analyzed. The model which also incorporates the effect of diurnal temperature range (DTR) on malaria transmission dynamics shows that increasing DTR shifts the peak temperature value for malaria transmission from 29C (when DTR is 0C) to about 25C (when DTR is 15C).
Finally, the malaria model is adapted and used to study the transmission dynamics of chikungunya, dengue and Zika, three diseases co-circulating in the Americas caused by the same vector (Aedes aegypti). The resulting model, which is fitted using data from Mexico, is used to assess a few hypotheses (such as those associated with the possible impact the newly-released dengue vaccine will have on Zika) and the impact of variability in climate variables on the dynamics of the three diseases. Suitable temperature and rainfall ranges for the maximum transmission intensity of the three diseases are obtained.
Motivated by the fact that understanding the dynamics of disease vector is crucial to understanding the transmission and control of the VBDs they cause, a novel weather-driven deterministic model for the population biology of the mosquito is formulated and rigorously analyzed. Numerical simulations, using relevant weather and entomological data for Anopheles mosquito (the vector for malaria), show that maximum mosquito abundance occurs when temperature and rainfall values lie in the range [20-25]C and [105-115] mm, respectively.
The Anopheles mosquito ecology model is extended to incorporate human dynamics. The resulting weather-driven malaria transmission model, which includes many of the key aspects of malaria (such as disease transmission by asymptomatically-infectious humans, and enhanced malaria immunity due to repeated exposure), was rigorously analyzed. The model which also incorporates the effect of diurnal temperature range (DTR) on malaria transmission dynamics shows that increasing DTR shifts the peak temperature value for malaria transmission from 29C (when DTR is 0C) to about 25C (when DTR is 15C).
Finally, the malaria model is adapted and used to study the transmission dynamics of chikungunya, dengue and Zika, three diseases co-circulating in the Americas caused by the same vector (Aedes aegypti). The resulting model, which is fitted using data from Mexico, is used to assess a few hypotheses (such as those associated with the possible impact the newly-released dengue vaccine will have on Zika) and the impact of variability in climate variables on the dynamics of the three diseases. Suitable temperature and rainfall ranges for the maximum transmission intensity of the three diseases are obtained.
ContributorsOkuneye, Kamaldeen O (Author) / Gumel, Abba B (Thesis advisor) / Kuang, Yang (Committee member) / Smith, Hal (Committee member) / Thieme, Horst (Committee member) / Nagy, John (Committee member) / Arizona State University (Publisher)
Created2018
Description
Rabies is an infectious viral disease. It is usually fatal if a victim reaches the rabid stage, which starts after the appearance of disease symptoms. The disease virus attacks the central nervous system, and then it migrates from peripheral nerves to the spinal cord and brain. At the time when the rabies virus reaches the brain, the incubation period is over and the symptoms of clinical disease appear on the victim. From the brain, the virus travels via nerves to the salivary glands and saliva.
A mathematical model is developed for the spread of rabies in a spatially distributed fox population to model the spread of the rabies epizootic through middle Europe that occurred in the second half of the 20th century. The model considers both territorial and wandering rabid foxes and includes a latent period for the infection. Since the model assumes these two kinds of rabid foxes, it is a system of both partial differential and integral equations (with integration
over space and, occasionally, also over time). To study the spreading speeds of the rabies epidemic, the model is reduced to a scalar Volterra-Hammerstein integral equation, and space-time Laplace transform of the integral equation is used to derive implicit formulas for the spreading speed. The spreading speeds are discussed and implicit formulas are given for latent periods of fixed length, exponentially distributed length, Gamma distributed length, and log-normally distributed length. A number of analytic and numerical results are shown pertaining to the spreading speeds.
Further, a numerical algorithm is described for the simulation
of the spread of rabies in a spatially distributed fox population on a bounded domain with Dirichlet boundary conditions. I propose the following methods for the numerical approximation of solutions. The partial differential and integral equations are discretized in the space variable by central differences of second order and by
the composite trapezoidal rule. Next, the ordinary or delay differential equations that are obtained this way are discretized in time by explicit
continuous Runge-Kutta methods of fourth order for ordinary and delay differential systems. My particular interest
is in how the partition of rabid foxes into
territorial and diffusing rabid foxes influences
the spreading speed, a question that can be answered by purely analytic means only for small basic reproduction numbers. I will restrict the numerical analysis
to latent periods of fixed length and to exponentially
distributed latent periods.
The results of the numerical calculations
are compared for latent periods
of fixed and exponentially distributed length
and for various proportions of territorial
and wandering rabid foxes.
The speeds of spread observed in the
simulations are compared
to spreading speeds obtained by numerically solving the analytic formulas
and to observed speeds of epizootic frontlines
in the European rabies outbreak 1940 to 1980.
A mathematical model is developed for the spread of rabies in a spatially distributed fox population to model the spread of the rabies epizootic through middle Europe that occurred in the second half of the 20th century. The model considers both territorial and wandering rabid foxes and includes a latent period for the infection. Since the model assumes these two kinds of rabid foxes, it is a system of both partial differential and integral equations (with integration
over space and, occasionally, also over time). To study the spreading speeds of the rabies epidemic, the model is reduced to a scalar Volterra-Hammerstein integral equation, and space-time Laplace transform of the integral equation is used to derive implicit formulas for the spreading speed. The spreading speeds are discussed and implicit formulas are given for latent periods of fixed length, exponentially distributed length, Gamma distributed length, and log-normally distributed length. A number of analytic and numerical results are shown pertaining to the spreading speeds.
Further, a numerical algorithm is described for the simulation
of the spread of rabies in a spatially distributed fox population on a bounded domain with Dirichlet boundary conditions. I propose the following methods for the numerical approximation of solutions. The partial differential and integral equations are discretized in the space variable by central differences of second order and by
the composite trapezoidal rule. Next, the ordinary or delay differential equations that are obtained this way are discretized in time by explicit
continuous Runge-Kutta methods of fourth order for ordinary and delay differential systems. My particular interest
is in how the partition of rabid foxes into
territorial and diffusing rabid foxes influences
the spreading speed, a question that can be answered by purely analytic means only for small basic reproduction numbers. I will restrict the numerical analysis
to latent periods of fixed length and to exponentially
distributed latent periods.
The results of the numerical calculations
are compared for latent periods
of fixed and exponentially distributed length
and for various proportions of territorial
and wandering rabid foxes.
The speeds of spread observed in the
simulations are compared
to spreading speeds obtained by numerically solving the analytic formulas
and to observed speeds of epizootic frontlines
in the European rabies outbreak 1940 to 1980.
ContributorsAlanazi, Khalaf Matar (Author) / Thieme, Horst R. (Thesis advisor) / Jackiewicz, Zdzislaw (Committee member) / Baer, Steven (Committee member) / Gardner, Carl (Committee member) / Kuang, Yang (Committee member) / Smith, Hal (Committee member) / Arizona State University (Publisher)
Created2018
Description
Ideas from coding theory are employed to theoretically demonstrate the engineering of mutation-tolerant genes, genes that can sustain up to some arbitrarily chosen number of mutations and still express the originally intended protein. Attention is restricted to tolerating substitution mutations. Future advances in genomic engineering will make possible the ability to synthesize entire genomes from scratch. This presents an opportunity to embed desirable capabilities like mutation-tolerance, which will be useful in preventing cell deaths in organisms intended for research or industrial applications in highly mutagenic environments. In the extreme case, mutation-tolerant genes (mutols) can make organisms resistant to retroviral infections.
An algebraic representation of the nucleotide bases is developed. This algebraic representation makes it possible to convert nucleotide sequences into algebraic sequences, apply mathematical ideas and convert results back into nucleotide terms. Using the algebra developed, a mapping is found from the naturally-occurring codons to an alternative set of codons which makes genes constructed from them mutation-tolerant, provided no more than one substitution mutation occurs per codon. The ideas discussed naturally extend to finding codons that can tolerate t arbitrarily chosen number of mutations per codon. Finally, random substitution events are simulated in both a wild-type green fluorescent protein (GFP) gene and its mutol variant and the amino acid sequence expressed from each post-mutation is compared with the amino acid sequence pre-mutation.
This work assumes the existence of synthetic protein-assembling entities that function like tRNAs but can read k nucleotides at a time, with k greater than or equal to 5. The realization of this assumption is presented as a challenge to the research community.
An algebraic representation of the nucleotide bases is developed. This algebraic representation makes it possible to convert nucleotide sequences into algebraic sequences, apply mathematical ideas and convert results back into nucleotide terms. Using the algebra developed, a mapping is found from the naturally-occurring codons to an alternative set of codons which makes genes constructed from them mutation-tolerant, provided no more than one substitution mutation occurs per codon. The ideas discussed naturally extend to finding codons that can tolerate t arbitrarily chosen number of mutations per codon. Finally, random substitution events are simulated in both a wild-type green fluorescent protein (GFP) gene and its mutol variant and the amino acid sequence expressed from each post-mutation is compared with the amino acid sequence pre-mutation.
This work assumes the existence of synthetic protein-assembling entities that function like tRNAs but can read k nucleotides at a time, with k greater than or equal to 5. The realization of this assumption is presented as a challenge to the research community.
ContributorsAmpofo, Prince Kwame (Author) / Tian, Xiaojun (Thesis advisor) / Kiani, Samira (Committee member) / Kuang, Yang (Committee member) / Arizona State University (Publisher)
Created2019