BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Elizabeth Gross

Elizabeth Gross appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.AGmath.COmath.ACmath.STMathematical SoftwareMethodologyMolecular NetworksNeurons and CognitionStatistics TheoryComputationGraphicsmath.NAPopulations and EvolutionQuantitative MethodsSymbolic Computation

Catalog footprint

What is connected

16works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

Broken Bracelets and Kostant's Partition Function

Inspired by the work of Amdeberhan, Can, and Moll on broken necklaces, we define a broken bracelet as a linear arrangement of marked and unmarked vertices and introduce a generalization called $n$-stars, which is a collection of $n$ broken bracelets whose final (unmarked) vertices are identified. Through these combinatorial objects, we provide a new framework for the study of Kostant's partition function, which counts the number of ways to express a vector as a nonnegative integer linear combination of the positive roots of a Lie algebra. Our main result establishes that (up to reflection) the number of broken bracelets with a fixed number of unmarked vertices with nonconsecutive marked vertices gives an upper bound for the value of Kostant's partition function for multiples of the highest root of a Lie algebra of type $A$. We connect this work to multiplex juggling sequences, as studied by Benedetti, Hanusa, Harris, Morales, and Simpson, by providing a correspondence to an equivalence relation on $n$-stars.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Binomial ideals of domino tilings

In this paper, we consider the set of all domino tilings of a cubiculated region. The primary question we explore is: How can we move from one tiling to another? Tiling spaces can be viewed as spaces of subgraphs of a fixed graph with a fixed degree sequence. Moves to connect such spaces have been explored in algebraic statistics. Thus, we approach this question from an applied algebra viewpoint, making new connections between domino tilings, algebraic statistics, and toric algebra. Using results from toric ideals of graphs, we are able to describe moves that connect the tiling space of a given cubiculated region of any dimension. This is done by studying binomials that arise from two distinct domino tilings of the same region. Additionally, we introduce tiling ideals and flip ideals and use these ideals to restate what it means for a tiling space to be flip connected. Finally, we show that if $R$ is a $2$-dimensional simply connected cubiculated region, any binomial arising from two distinct tilings of $R$ can be written in terms of quadratic binomials. As a corollary to our main result, we obtain an alternative proof to the fact that the set of domino tilings of a $2$-dimensional simply connected region is connected by flips.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The steady-state degree and mixed volume of a chemical reaction network

The steady-state degree of a chemical reaction network is the number of complex steady-states, which is a measure of the algebraic complexity of solving the steady-state system. In general, the steady-state degree may be difficult to compute. Here, we give an upper bound to the steady-state degree of a reaction network by utilizing the underlying polyhedral geometry associated with the corresponding polynomial system. We focus on three case studies of infinite families of networks, each generated by joining smaller networks to create larger ones. For each family, we give a formula for the steady-state degree and the mixed volume of the corresponding polynomial system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neural ideals and stimulus space visualization

A neural code $\mathcal{C}$ is a collection of binary vectors of a given length n that record the co-firing patterns of a set of neurons. Our focus is on neural codes arising from place cells, neurons that respond to geographic stimulus. In this setting, the stimulus space can be visualized as subset of $\mathbb{R}^2$ covered by a collection $\mathcal{U}$ of convex sets such that the arrangement $\mathcal{U}$ forms an Euler diagram for $\mathcal{C}$. There are some methods to determine whether such a convex realization $\mathcal{U}$ exists; however, these methods do not describe how to draw a realization. In this work, we look at the problem of algorithmically drawing Euler diagrams for neural codes using two polynomial ideals: the neural ideal, a pseudo-monomial ideal; and the neural toric ideal, a binomial ideal. In particular, we study how these objects are related to the theory of piercings in information visualization, and we show how minimal generating sets of the ideals reveal whether or not a code is $0$, $1$, or $2$-inductively pierced.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical algebraic geometry for model selection and its application to the life sciences

Researchers working with mathematical models are often confronted by the related problems of parameter estimation, model validation, and model selection. These are all optimization problems, well-known to be challenging due to non-linearity, non-convexity and multiple local optima. Furthermore, the challenges are compounded when only partial data is available. Here, we consider polynomial models (e.g., mass-action chemical reaction networks at steady state) and describe a framework for their analysis based on optimization using numerical algebraic geometry. Specifically, we use probability-one polynomial homotopy continuation methods to compute all critical points of the objective function, then filter to recover the global optima. Our approach exploits the geometric structures relating models and data, and we demonstrate its utility on examples from cell signaling, synthetic biology, and epidemiology.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phylogenetic trees

We introduce the package PhylogeneticTrees for Macaulay2 which allows users to compute phylogenetic invariants for group-based tree models. We provide some background information on phylogenetic algebraic geometry and show how the package PhylogeneticTrees can be used to calculate a generating set for a phylogenetic ideal as well as a lower bound for its dimension. Finally, we show how methods within the package can be used to compute a generating set for the join of any two ideals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

What makes a neural code convex?

Neural codes allow the brain to represent, process, and store information about the world. Combinatorial codes, comprised of binary patterns of neural activity, encode information via the collective behavior of populations of neurons. A code is called convex if its codewords correspond to regions defined by an arrangement of convex open sets in Euclidean space. Convex codes have been observed experimentally in many brain areas, including sensory cortices and the hippocampus, where neurons exhibit convex receptive fields. What makes a neural code convex? That is, how can we tell from the intrinsic structure of a code if there exists a corresponding arrangement of convex open sets? In this work, we provide a complete characterization of local obstructions to convexity. This motivates us to define max intersection-complete codes, a family guaranteed to have no local obstructions. We then show how our characterization enables one to use free resolutions of Stanley-Reisner ideals in order to detect violations of convexity. Taken together, these results provide a significant advance in understanding the intrinsic combinatorial properties of convex codes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Algebraic Systems Biology: A Case Study for the Wnt Pathway

Steady state analysis of dynamical systems for biological networks give rise to algebraic varieties in high-dimensional spaces whose study is of interest in their own right. We demonstrate this for the shuttle model of the Wnt signaling pathway. Here the variety is described by a polynomial system in 19 unknowns and 36 parameters. Current methods from computational algebraic geometry and combinatorics are applied to analyze this model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Maximum Likelihood Threshold of a Graph

The maximum likelihood threshold of a graph is the smallest number of data points that guarantees that maximum likelihood estimates exist almost surely in the Gaussian graphical model associated to the graph. We show that this graph parameter is connected to the theory of combinatorial rigidity. In particular, if the edge set of a graph $G$ is an independent set in the $n-1$-dimensional generic rigidity matroid, then the maximum likelihood threshold of $G$ is less than or equal to $n$. This connection allows us to prove many results about the maximum likelihood threshold.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Goodness-of-fit for log-linear network models: Dynamic Markov bases using hypergraphs

Social networks and other large sparse data sets pose significant challenges for statistical inference, as many standard statistical methods for testing model fit are not applicable in such settings. Algebraic statistics offers a theoretically justified approach to goodness-of-fit testing that relies on the theory of Markov bases and is intimately connected with the geometry of the model as described by its fibers. Most current practices require the computation of the entire basis, which is infeasible in many practical settings. We present a dynamic approach to explore the fiber of a model, which bypasses this issue, and is based on the combinatorics of hypergraphs arising from the toric algebra structure of log-linear models. We demonstrate the approach on the Holland-Leinhardt $p_1$ model for random directed graphs that allows for reciprocated edges.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Maximum likelihood geometry in the presence of data zeros

Given a statistical model, the maximum likelihood degree is the number of complex solutions to the likelihood equations for generic data. We consider discrete algebraic statistical models and study the solutions to the likelihood equations when the data contain zeros and are no longer generic. Focusing on sampling and model zeros, we show that, in these cases, the solutions to the likelihood equations are contained in a previously studied variety, the likelihood correspondence. The number of these solutions give a lower bound on the ML degree, and the problem of finding critical points to the likelihood function can be partitioned into smaller and computationally easier problems involving sampling and model zeros. We use this technique to compute a lower bound on the ML degree for $2 \times 2 \times 2 \times 2$ tensors of border rank $\leq 2$ and $3 \times n$ tables of rank $\leq 2$ for $n=11, 12, 13, 14$, the first four values of $n$ for which the ML degree was previously unknown.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bertini for Macaulay2

Numerical algebraic geometry is the field of computational mathematics concerning the numerical solution of polynomial systems of equations. Bertini, a popular software package for computational applications of this field, includes implementations of a variety of algorithms based on polynomial homotopy continuation. The Macaulay2 package Bertini.m2 provides an interface to Bertini, making it possible to access the core run modes of Bertini in Macaulay2. With these run modes, users can find approximate solutions to zero-dimensional systems and positive-dimensional systems, test numerically whether a point lies on a variety, sample numerically from a variety, and perform parameter homotopy runs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Combinatorial degree bound for toric ideals of hypergraphs

Associated to any hypergraph is a toric ideal encoding the algebraic relations among its edges. We study these ideals and the combinatorics of their minimal generators, and derive general degree bounds for both uniform and non-uniform hypergraphs in terms of balanced hypergraph bicolorings, separators, and splitting sets. In turn, this provides complexity bounds for algebraic statistical models associated to hypergraphs. As two main applications, we recover a well-known complexity result for Markov bases of arbitrary 3-way tables, and we show that the defining ideal of the tangential variety is generated by quadratics and cubics in cumulant coordinates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

PHCpack in Macaulay2

The Macaulay2 package PHCpack.m2 provides an interface to PHCpack, a general-purpose polynomial system solver that uses homotopy continuation. The main method is a numerical blackbox solver which is implemented for all Laurent systems. The package also provides a fast mixed volume computation, the ability to filter solutions, homotopy path tracking, and a numerical irreducible decomposition method. As the size of many problems in applied algebraic geometry often surpasses the capabilities of symbolic software, this package will be of interest to those working on problems involving large polynomial systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Maximum likelihood degree of variance component models

Most statistical software packages implement numerical strategies for computation of maximum likelihood estimates in random effects models. Little is known, however, about the algebraic complexity of this problem. For the one-way layout with random effects and unbalanced group sizes, we give formulas for the algebraic degree of the likelihood equations as well as the equations for restricted maximum likelihood estimation. In particular, the latter approach is shown to be algebraically less complex. The formulas are obtained by studying a univariate rational equation whose solutions correspond to the solutions of the likelihood equations. Applying techniques from computational algebra, we also show that balanced two-way layouts with or without interaction have likelihood equations of degree four. Our work suggests that algebraic methods allow one to reliably find global optima of likelihood functions of linear mixed models with a small number of variance components.

0 citations0 reviewsNo code linkedOpen access