-
Evaluating COVID-19 Surveillance Testing Strategies at Colorado School of Mines: A Stochastic Modeling Approach
Authors:
Laura Albrecht,
Karin Leiderman,
Suzanne Sindi,
Douglas Nychka
Abstract:
This study introduces a stochastic model of COVID-19 transmission tailored to the Colorado School of Mines campus and evaluates surveillance testing strategies within a university context. Enhancing the conventional SEIR framework with stochastic transitions, our model accounts for the unique characteristics of disease spread in a residential college, including specific states for testing, quarant…
▽ More
This study introduces a stochastic model of COVID-19 transmission tailored to the Colorado School of Mines campus and evaluates surveillance testing strategies within a university context. Enhancing the conventional SEIR framework with stochastic transitions, our model accounts for the unique characteristics of disease spread in a residential college, including specific states for testing, quarantine, and isolation. Employing an approximate Bayesian computation (ABC) method for parameter estimation, we navigate the complexities inherent in stochastic models, enabling an accurate fit of the model with the campus case data. We then studied our model under the estimated parameters to evaluate the efficacy of different testing policies that could be implemented on a university campus. This framework not only advances understanding of COVID-19 dynamics on the Mines campus but serves as a blueprint for comparable settings, providing insights for informed strategies against infectious diseases.
△ Less
Submitted 5 April, 2024;
originally announced April 2024.
-
N-ACT: An Interpretable Deep Learning Model for Automatic Cell Type and Salient Gene Identification
Authors:
A. Ali Heydari,
Oscar A. Davalos,
Katrina K. Hoyer,
Suzanne S. Sindi
Abstract:
Single-cell RNA sequencing (scRNAseq) is rapidly advancing our understanding of cellular composition within complex tissues and organisms. A major limitation in most scRNAseq analysis pipelines is the reliance on manual annotations to determine cell identities, which are time consuming, subjective, and require expertise. Given the surge in cell sequencing, supervised methods-especially deep learni…
▽ More
Single-cell RNA sequencing (scRNAseq) is rapidly advancing our understanding of cellular composition within complex tissues and organisms. A major limitation in most scRNAseq analysis pipelines is the reliance on manual annotations to determine cell identities, which are time consuming, subjective, and require expertise. Given the surge in cell sequencing, supervised methods-especially deep learning models-have been developed for automatic cell type identification (ACTI), which achieve high accuracy and scalability. However, all existing deep learning frameworks for ACTI lack interpretability and are used as "black-box" models. We present N-ACT (Neural-Attention for Cell Type identification): the first-of-its-kind interpretable deep neural network for ACTI utilizing neural-attention to detect salient genes for use in cell-type identification. We compare N-ACT to conventional annotation methods on two previously manually annotated data sets, demonstrating that N-ACT accurately identifies marker genes and cell types in an unsupervised manner, while performing comparably on multiple data sets to current state-of-the-art model in traditional supervised ACTI.
△ Less
Submitted 8 May, 2022;
originally announced June 2022.
-
A Statistical Model of Word Rank Evolution
Authors:
Alex John Quijano,
Rick Dale,
Suzanne Sindi
Abstract:
The availability of large linguistic data sets enables data-driven approaches to study linguistic change. The Google Books corpus unigram frequency data set is used to investigate the word rank dynamics in eight languages. We observed the rank changes of the unigrams from 1900 to 2008 and compared it to a Wright-Fisher inspired model that we developed for our analysis. The model simulates a neutra…
▽ More
The availability of large linguistic data sets enables data-driven approaches to study linguistic change. The Google Books corpus unigram frequency data set is used to investigate the word rank dynamics in eight languages. We observed the rank changes of the unigrams from 1900 to 2008 and compared it to a Wright-Fisher inspired model that we developed for our analysis. The model simulates a neutral evolutionary process with the restriction of having no disappearing and added words. This work explains the mathematical framework of the model - written as a Markov Chain with multinomial transition probabilities - to show how frequencies of words change in time. From our observations in the data and our model, word rank stability shows two types of characteristics: (1) the increase/decrease in ranks are monotonic, or (2) the rank stays the same. Based on our model, high-ranked words tend to be more stable while low-ranked words tend to be more volatile. Some words change in ranks in two ways: (a) by an accumulation of small increasing/decreasing rank changes in time and (b) by shocks of increase/decrease in ranks. Most words in all of the languages we have looked at are rank stable, but not as stable as a neutral model would predict. The stopwords and Swadesh words are observed to be rank stable across eight languages indicating linguistic conformity in established languages. These signatures suggest unigram frequencies in all languages have changed in a manner inconsistent with a purely neutral evolutionary process.
△ Less
Submitted 14 February, 2022; v1 submitted 21 July, 2021;
originally announced July 2021.
-
A field guide to cultivating computational biology
Authors:
Anne E Carpenter,
Casey S Greene,
Piero Carnici,
Benilton S Carvalho,
Michiel de Hoon,
Stacey Finley,
Kim-Anh Le Cao,
Jerry SH Lee,
Luigi Marchionni,
Suzanne Sindi,
Fabian J Theis,
Gregory P Way,
Jean YH Yang,
Elana J Fertig
Abstract:
Biomedical research centers can empower basic discovery and novel therapeutic strategies by leveraging their large-scale datasets from experiments and patients. This data, together with new technologies to create and analyze it, has ushered in an era of data-driven discovery which requires moving beyond the traditional individual, single-discipline investigator research model. This interdisciplina…
▽ More
Biomedical research centers can empower basic discovery and novel therapeutic strategies by leveraging their large-scale datasets from experiments and patients. This data, together with new technologies to create and analyze it, has ushered in an era of data-driven discovery which requires moving beyond the traditional individual, single-discipline investigator research model. This interdisciplinary niche is where computational biology thrives. It has matured over the past three decades and made major contributions to scientific knowledge and human health, yet researchers in the field often languish in career advancement, publication, and grant review. We propose solutions for individual scientists, institutions, journal publishers, funding agencies, and educators.
△ Less
Submitted 22 April, 2021;
originally announced April 2021.
-
Slack Reactants: A State-Space Truncation Framework to Estimate Quantitative Behavior of the Chemical Master Equation
Authors:
Jinsu Kim,
Jason Dark,
German Enciso,
Suzanne Sindi
Abstract:
State space truncation methods are widely used to approximate solutions of the chemical master equation. While most methods of this kind focus on truncating the state space directly, in this work we propose modifying the underlying chemical reaction network by introducing \emph{slack reactants} that indirectly truncate the state space. More specifically, slack reactants introduce an expanded chemi…
▽ More
State space truncation methods are widely used to approximate solutions of the chemical master equation. While most methods of this kind focus on truncating the state space directly, in this work we propose modifying the underlying chemical reaction network by introducing \emph{slack reactants} that indirectly truncate the state space. More specifically, slack reactants introduce an expanded chemical reaction network and impose a truncation scheme based on user defined properties, such as mass-conservation. This network structure allows us to prove inheritance of special properties of the original model, such as irreducibility and complex balancing.
We use the network structure imposed by slack reactants to prove the convergence of the stationary distribution and first arrival times. We then provide examples comparing our method with the stationary finite state projection and finite buffer methods. Our slack reactant system appears to be more robust than some competing methods with respect to calculating first arrival times.
△ Less
Submitted 11 May, 2020;
originally announced May 2020.
-
Ensemble-based Topological Entropy Calculation (E-tec)
Authors:
Eric Roberts,
Suzanne Sindi,
Spencer Smith,
Kevin Mitchell
Abstract:
Topological entropy measures the number of distinguishable orbits in a dynamical system, thereby quantifying the complexity of chaotic dynamics. One approach to computing topological entropy in a two-dimensional space is to analyze the collective motion of an ensemble of system trajectories taking into account how trajectories "braid" around one another. In this spirit, we introduce the Ensemble-b…
▽ More
Topological entropy measures the number of distinguishable orbits in a dynamical system, thereby quantifying the complexity of chaotic dynamics. One approach to computing topological entropy in a two-dimensional space is to analyze the collective motion of an ensemble of system trajectories taking into account how trajectories "braid" around one another. In this spirit, we introduce the Ensemble-based Topological Entropy Calculation, or E-tec, a method to derive a lower-bound on topological entropy of two-dimensional systems by considering the evolution of a "rubber band" (piece-wise linear curve) wrapped around the data points and evolving with their trajectories. The topological entropy is bounded below by the exponential growth rate of this band. We use tools from computational geometry to track the evolution of the rubber band as data points strike and deform it. Because we maintain information about the configuration of trajectories with respect to one another, updating the band configuration is performed locally, which allows E-tec to be more computationally efficient than some competing methods. In this work, we validate and illustrate many features of E-tec on a chaotic lid-driven cavity flow. In particular, we demonstrate convergence of E-tec's approximation with respect to both the number of trajectories (ensemble size) and the duration of trajectories in time.
△ Less
Submitted 17 April, 2019; v1 submitted 30 August, 2018;
originally announced September 2018.
-
Topological Chaos in a Three-Dimensional Spherical Fluid Vortex
Authors:
Spencer A. Smith,
Joshua Arenson,
Eric Roberts,
Suzanne Sindi,
Kevin A. Mitchell
Abstract:
In chaotic deterministic systems, seemingly stochastic behavior is generated by relatively simple, though hidden, organizing rules and structures. Prominent among the tools used to characterize this complexity in 1D and 2D systems are techniques which exploit the topology of dynamically invariant structures. However, the path to extending many such topological techniques to three dimensions is fil…
▽ More
In chaotic deterministic systems, seemingly stochastic behavior is generated by relatively simple, though hidden, organizing rules and structures. Prominent among the tools used to characterize this complexity in 1D and 2D systems are techniques which exploit the topology of dynamically invariant structures. However, the path to extending many such topological techniques to three dimensions is filled with roadblocks that prevent their application to a wider variety of physical systems. Here, we overcome these roadblocks and successfully analyze a realistic model of 3D fluid advection, by extending the homotopic lobe dynamics (HLD) technique, previously developed for 2D area-preserving dynamics, to 3D volume-preserving dynamics. We start with numerically-generated finite-time chaotic-scattering data for particles entrained in a spherical fluid vortex, and use this data to build a symbolic representation of the dynamics. We then use this symbolic representation to explain and predict the self-similar fractal structure of the scattering data, to compute bounds on the topological entropy, a fundamental measure of mixing, and to discover two different mixing mechanisms, which stretch 2D material surfaces and 1D material curves in distinct ways.
△ Less
Submitted 15 March, 2017;
originally announced March 2017.