• Testing multispecies coalescent simulators with summary statistics

      Baños Cervantes, Hector Daniel; Allman, Elizabeth; Rhodes, John; Goddard, Scott; McIntyre, Julie; Barry, Ron (2018-12)
      The Multispecies coalescent model (MSC) is increasingly used in phylogenetics to describe the formation of gene trees (depicting the direct ancestral relationships of sampled lineages) within species trees (depicting the branching of species from their common ancestor). A number of MSC simulators have been implemented, and these are often used to test inference methods built on the model. However, it is not clear from the literature that these simulators are always adequately tested. In this project, we formulated tools for testing these simulators and use them to show that of four well-known coalescent simulators, Mesquite, Hybrid-Lambda, SimPhy, and Phybase, only SimPhy performs correctly according to these tests.
    • The linear algebra of interpolation with finite applications giving computational methods for multivariate polynomials

      Olmsted, Coert D.; Gislason, Gary A.; Lambert, J. P.; Lando, C. A.; Olson, J. V.; Piacenca, R. J. (1988)
      Linear representation and the duality of the biorthonormality relationship express the linear algebra of interpolation by way of the evaluation mapping. In the finite case the standard bases relate the maps to Gramian matrices. Five equivalent conditions on these objects are found which characterize the solution of the interpolation problem. This algebra succinctly describes the solution space of ordinary linear initial value problems. Multivariate polynomial spaces and multidimensional node sets are described by multi-index sets. Geometric considerations of normalization and dimensionality lead to cardinal bases for Lagrange interpolation on regular node sets. More general Hermite functional sets can also be solved by generalized Newton methods using geometry and multi-indices. Extended to countably infinite spaces, the method calls upon theorems of modern analysis.
    • Toward an optimal solver for the obstacle problem

      Heldman, Max; Bueler, Ed; Maxwell, David; Rhodes, John (2018-04)
      An optimal algorithm for solving a problem with m degrees of freedom is one that computes a solution in O (m) time. In this paper, we discuss a class of optimal algorithms for the numerical solution of PDEs called multigrid methods. We go on to examine numerical solvers for the obstacle problem, a constrained PDE, with the goal of demonstrating optimality. We discuss two known algorithms, the so-called reduced space method (RSP) [BM03] and the multigrid-based projected full-approximation scheme (PFAS) [BC83]. We compare the performance of PFAS and RSP on a few example problems, finding numerical evidence of optimality or near-optimality for PFAS.
    • The treatment of missing data on placement tools for predicting success in college algebra at the University of Alaska

      Crawford, Alyssa (2014-05)
      This project investigated the statistical significance of baccalaureate student placement tools such as tests scores and completion of a developmental course on predicting success in a college level algebra course at the University of Alaska (UA). Students included in the study had attempted Math 107 at UA for the first time between fiscal years 2007 and 2012. The student placement information had a high percentage of missing data. A simulation study was conducted to choose the best missing data method between complete case deletion, and multiple imputation for the student data. After the missing data methods were applied, a logistic regression with fitted with explanatory variables consisting of tests scores, developmental course grade, age (category) of scores and grade, and interactions. The relevant tests were SAT math, ACT math, AccuPlacer college level math, and the relevant developmental course was Devm /Math 105. The response variable was success in passing Math 107 with grade of C or above on the first attempt. The simulation study showed that under a high percentage of missing data and correlation, multiple imputation implemented by the R package Multivariate Imputation by Chained Equations (MICE) produced the least biased estimators and better confidence interval coverage compared to complete cases deletion when data are missing at random (MAR) and missing not at random (MNAR). Results from multiple imputation method on the student data showed that Devm /Math 105 grade was a significant predictor of passing Math 107. The age of Devm /Math 105, age of tests, and test scores were not significant predictors of student success in Math 107. Future studies may consider modeling with ALEKS scores, and high school math course information.
    • Tsunami runup in U and V shaped bays

      Garayshin, Viacheslav Valer'evich; Гарайшин, Вячеслав Валерьевич; Rybkin, Alexei; Rhodes, John; Nicolsky, Dmitry (2013-08)
      Tsunami runup can be effectively modeled using the shallow water wave equations. In 1958 Carrier and Greenspan in their work "Water waves of finite amplitude on a sloping beach" used this system to model tsunami runup on a uniformly sloping plane beach. They linearized this problem using a hodograph type transformation and obtained the Klein-Gordon equation which could be explicitly solved by using the Fourier-Bessel transform. In 2011, Efim Pelinovsky and Ira Didenkulova in their work "Runup of Tsunami Waves in U-Shaped Bays" used a similar hodograph type transformation and linearized the tsunami problem for a sloping bay with parabolic cross-section. They solved the linear system by using the D'Alembert formula. This method was generalized to sloping bays with cross-sections parameterized by power functions. However, an explicit solution was obtained only for the case of a bay with a quadratic cross-section. In this paper we will show that the Klein-Gordon equation can be solved by a spectral method for any inclined bathymetry with power function for any positive power. The result can be used to estimate tsunami runup in such bays with minimal numerical computations. This fact is very important because in many cases our numerical model can be substituted for fullscale numerical models which are computationally expensive, and time consuming, and not feasible to investigate tsunami behavior in the Alaskan coastal zone, due to the low population density in this area
    • Vertex arboricity of triangle-free graphs

      Warren, Samantha; Gimbel, John; Faudree, Jill; Allman, Elizabeth (2016-05)
      The vertex arboricity of a graph is the minimum number of colors needed to color the vertices so that the subgraph induced by each color class is a forest. In other words, the vertex arboricity of a graph is the fewest number of colors required in order to color a graph such that every cycle has at least two colors. Although not standard, we will refer to vertex arboricity simply as arboricity. In this paper, we discuss properties of chromatic number and k-defective chromatic number and how those properties relate to the arboricity of trianglefree graphs. In particular, we find bounds on the minimum order of a graph having arboricity three. Equivalently, we consider the largest possible vertex arboricity of triangle-free graphs of fixed order.