From RNA-seq Time Series Data to Models of Regulatory Networks
|
|
- Jean Adams
- 5 years ago
- Views:
Transcription
1 From RNA-seq Time Series Data to Models of Regulatory Networks Konstantin Mischaikow Dept. of Mathematics, Rutgers Brown, Dec 2017
2 A Model Problem
3 Malaria Estimated number of malaria cases in 2010: between 219 and 550 million Estimated number of deaths due to malaria in 2010: 600,000 to 1,240,000 Malaria may have killed half of all the people that ever lived. And more people are now infected than at any point in history. There are up to half a billion cases every year, and about 2 million deaths - half of those are children in sub-saharan Africa. J. Whitfield, Nature, 2002 Malaria is of great public health concern, and seems likely to be the vector-borne disease most sensitive to long-term climate change. World Health Organization Resistance is now common against all classes of antimalarial drugs apart from artemisinins. Malaria strains found on the Cambodia Thailand border are resistant to combination therapies that include artemisinins, and may therefore be untreatable. World Health Organization
4 Malaria: P. falciparum Task: Understand the regulation on the genetic/biomolecular level with the goal of affecting the dynamics with drugs. All genes (5409) 48 hour cycle 1-2 minutes periodic$genes$(43)$ 1.5$ 0.0$ High Standard$devia0ons$from$ mean$expression$(z&score)$ A differential equation dx dt = f(x, ) is probably a reasonable model for the dynamics, but I do not have an analytic description of f or estimates of the parameters. &1.5$ Low 0 0$ 10$ 20$ 30$ 40$ 50$ 60$ 0me$in#vitro#(hours)$ Malaria is Sequenced Poorly annotated Walter Reed Army Inst. Research Duke University A proposed network
5 A Philosophical Interlude
6 The Lac Operon Ozbudak et al. Nature 2004 ODES are great modeling tools, but should be handled with care. Network Model 1 R T ẏ = y R T + R(x) 1 ẋ = y x x R(x) = 1+ R T n x x 0 ODE Model y = Data (G/8.1) 1.2 =... parameter values
7 What does it mean to solve an ODE? truth model Precise Not Accurate Not Rigorous Classical Qualita7ve Representa7on of Dynamics parameter Conley-Morse Chain Complex Dynamic Signature (Morse Graph) Not Precise Accurate Rigorous
8 Combinatorial Dynamics
9 State Transition Graph F : X! X p3 Linear time Algorithm! Simple decomposition of Dynamics: Recurrent p2 Nonrecurrent (gradient-like) p1 p0 Vertices: States Edges: Dynamics Don t know exact current state, so don t know exact next state Morse Graph of state transition graph Poset
10 What is observable? A X is an attractor if F(A) =A Birkhoff s Theorem implies that the Morse graph and the lattice of Attractors are equivalent. p3 Observable p3,p2,p1,p0 p2,p1,p0 p1,p0 Lower Sets O(M) Join Irreducible J _ (A) p2 P O S E T Computable p1 p0 p1 p0 ; Lattice of Attractors _ = [ of F : X! X ^ = maximal attractor in \ Morse Graph of F : X! X
11 Topology (differential equations are not defined on discrete sets)
12 Topology Let X be a compact metric space. Let R(X) denote the lattice of regular closed subsets of X. Let L be a finite bounded sublattice of R(X). G(L) denoted atoms of L phase space Infinite unbounded lattice Level of measurement Applicable scale for model smallest elements of L Dynamics Declare a bounded sublattice A L to be a lattice of forward invariant regions (attracting blocks). Use Birkhoff to define poset (P := J _ (A),<) For each p 2 P define a Morse tile M(p) :=cl(a \ pred(a)) Remark: I have purposefully ignored the relation between L and F : X! X
13 Example Phase space: X =[ 4, 4] R Atoms of lattice: G(L) ={[n, n + 1] n = 4,...,3} Lattice of attracting blocks: A = {[ 3, 1], [1, 3], [ 3, 1] [ [1, 3], [ 4, 4]} A Birkhoff P Let F 0 (x) = How does this relate to a differential equation dx dt = f(x)? f(x). Attracting blocks are regions of phase space that are forward invariant with time. Morse tiles M(p) F F Remark: This leads to a homology theory F
14 Switching Systems (an example of how to use these ideas) Choosing L and F : X!! X
15 Biological Model How do I want to interpret this information? What differential equation do I want to use? x i denotes amount of species i. Assume x i decays. Proposed model: dx dx i i dx i dt = ix i dt dt = ix ix i i + i (x) i j ) Parameters 1/node 3/edge 1 2 x 1 represses the production of x 2. dx 2 dt u 2,1 l 2,1 1 2 x 1 activates the production of x 2. x 1 2,1 For x 1 < 2,1 we ask about sign ( 2x 2 + u 2,1 ). For x 1 > 2,1 we ask about sign ( 2x 2 + l 2,1 ). j,i (x i)= ( u j,i `j,i if x i < j,i if x i > j,i Focus on sign of ix i + i,j (x j ) ix i + + i,j (x j)
16 Example 1 2 We care about sign of: 1 2,1 + 1,2 (x 2 ) 2 1, ,1 (x 1) Parameter space is a subset of (0, 1) 8 Fix z a regular parameter value. x 2 If x 1 < 2,1 and 2 1, ,1 (x 1) > 0 If x 1 < 2,1 and 2 1, ,1 (x 1) < 0 1,2 z is a regular parameter value if 2,1 Phase space: X =(0, 1) 2 x 1 0 < i 0 < `i,j <u i,j, 0 < i,k 6= j,k, and 0 6= i j,i + i (x)
17 Fix z a regular parameter value. Need to Construct State Transition Graph F z : X! X x 2 Vertices X corresponds to all rectangular domains and co-dimension 1 faces defined by thresholds. 1,2 Edges Faces pointing in map to their domain. 2,1 x 1 Domains map to their faces pointing out. If no outpointing faces domain maps to itself.
18 Example 1 2 Constructing state transition graph F z : X! X Fix z a regular parameter value. Check signs of i j,i + ± i,j (x j) Assume: l 1,2 < 1 2,1 <u 1,2 2 1,2 <l 2,1 <u 2,1 l 1,2 < 1 2,1 <u 1,2 l 2,1 < 2 1,2 <u 2,1 x 2 x 2 1,2 1,2 x 1 x 1 Morse Graph 2,1 FP{0,1} FC 2,1 Dynamics orders maxima and minima (M 1,M 2,m 1,m 2 )
19 Input: Regulatory Network DSGRN Database 1 2 Output: DSGRN database (7) FP(1,0) 1 2,1 <l 1,2 <u 1,2 l 2,1 <u 2,1 < 2 1,2 (8) FP(1,0) l 1,2 < 1 2,1 <u 1,2 l 2,1 <u 2,1 < 2 1,2 (9) FP(0,0) l 1,2 <u 1,2 < 1 2,1 l 2,1 <u 2,1 < 2 1,2 (4) FP(1,1) 1 2,1 <l 1,2 <u 1,2 l 2,1 < 2 1,2 <u 2,1 (5) FC l 1,2 < 1 2,1 <u 1,2 l 2,1 < 2 1,2 <u 2,1 (6) FP(0,0) l 1,2 <u 1,2 < 1 2,1 l 2,1 < 2 1,2 <u 2,1 (1) FP(1,1) 1 2,1 <l 1,2 <u 1,2 (2) FP(0,1) l 1,2 < 1 2,1 <u 1,2 (3) FP(0,1) l 1,2 <u 1,2 < 1 2,1 2 1,2 <l 2,1 <u 2,1 2 1,2 <l 2,1 <u 2,1 2 1,2 <l 2,1 <u 2,1 Parameter graph provides explicit partition of entire 8-D parameter space. Observe that we can query this database for local or global dynamics.
20 Back to Malaria Remark: there are a variety of statistical methods for generating possible regulatory networks from this type of time series data.
21 In vitro data (WRAIR+ Haase) Remarks about dynamics: 1. Gene expression is cyclic in nature. 2. We know relative times of expression of genes. Putative TF genes (456) High Standard deviations from mean expression (z-score) -1.5 Low Assumption: Expression of important functions must be robust to perturbations Time in vitro (hours)
22 Simple Test Cyclic feedback system: well understood using classical dynamical systems techniques. Experimental 7me series for associated genes Under the assump7on of monotone switches if parameter values are chosen such that there exists a stable periodic orbit, then the maxima in the network must occur in the order: (188,93,184, 395) (green, blue, cyan,red) Conclusion: This network does not generate observed dynamics
23 No mathema7cal theory DSGRN computa7on produces a parameter graph with approximately 45,000 nodes. Computa7on 7me on laptop approximately 1 second. Time series for associated genes SQL Query: A stable cycle involving oscilla7ons in all genes 96 parameter graph nodes with Morse graph that has a minimal node consis7ng of a Full Cycle (FC).
24 DSGRN Analysis (II): Max-Min Matching M m M m M m M M m M m m Have developed polynomial 7me algorithm that take paths in state transi7on graph and iden7fies sequences of possible maxima and minima. Tested all max-min sequences from state transi7on graphs from all 96 parameter graph nodes against 17,280 experimental pa`erns. No Match Conclusion: This network does not generate observed dynamics
25 DSGRN strategy I. Start with a proposed gene regulatory network. Extract from experimental data the poset indicating possible max-min orderings. Perform DSGRN computation to identify parameter node for which minimal Morse node is FC Reject parameter node if max-min sequences of FC are not linear extensions of poset. Compute fraction of parameter nodes that match experimental data. Assumption: Expression of important functions must be robust to perturbations. IF fraction is small, THEN reject regulatory network.
26 DSGRN strategy II. Start with an acceptable gene regulatory network. Create new regulatory network via random perturbations: Add/Remove edge(s) Add/Remove node I. Start with the newly proposed gene regulatory network. 90% of parameter nodes* induce minimal FC node Current favorite network 80% of parameter nodes* induce minimal FC which agrees with experimentally determined max-min ordering.
27 Thank-you for your Attention Rutgers S. Harker MSU T. Gedeon B. Cummings FAU W. Kalies VU Amsterdam R. Vandervorst Homology + Database Software chomp.rutgers.edu
28 Example (continued) X =[ 4, 4] R A (P,<) The homology Conley index of M(p) is (k, 0,...) (0, k, 0,...) 3 5 (k, 0,...) F F CH (p) :=H (A, pred(a); k) k a field Conley index can be used to guarantee existence of equilibria, periodic orbits, heteroclinic and homoclinic orbits, and chaotic dynamics. Theorem: (R. Franzosa) There exists a strictly upper triangular (with respect to <) boundary operator : M p2p CH (p)! M p2p CH (p) such that the induced homology is isomorphic to H (X).
TOWARDS AUTOMATED CHAOS VERIFICATION
TOWARDS AUTOMATED CHAOS VERIFICATION SARAH DAY CDSNS, Georgia Institute of Technology, Atlanta, GA 30332 OLIVER JUNGE Institute for Mathematics, University of Paderborn, 33100 Paderborn, Germany KONSTANTIN
More informationLattice Structures for Attractors I
Lattice Structures for Attractors I WILLIAM D. KALIES Florida Atlantic University 777 Glades Road Boca Raton, FL 33431, USA KONSTANTIN MISCHAIKOW Rutgers University 110 Frelinghusen Road Piscataway, NJ
More informationEFFICIENT COMPUTATION OF LYAPUNOV FUNCTIONS FOR MORSE DECOMPOSITIONS
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX EFFICIENT COMPUTATION OF LYAPUNOV FUNCTIONS FOR MORSE DECOMPOSITIONS ARNAUD GOULLET, SHAUN HARKER, KONSTANTIN
More informationComputing All Isolated Invariant Sets at a Finite Resolution
College of William and Mary W&M ScholarWorks Undergraduate Honors Theses Theses, Dissertations, & Master Projects 5-2016 Computing All Isolated Invariant Sets at a Finite Resolution Martin Salgado-Flores
More informationMath 4200, Problem set 3
Math, Problem set 3 Solutions September, 13 Problem 1. ẍ = ω x. Solution. Following the general theory of conservative systems with one degree of freedom let us define the kinetic energy T and potential
More informationLattice Structures for Attractors II
Lattice Structures for Attractors II WILLIAM D. KALIES Florida Atlantic University 777 Glades Road Boca Raton, FL 33431, USA KONSTANTIN MISCHAIKOW Rutgers University 110 Frelinghusen Road Piscataway, NJ
More informationMorse-Forman-Conley theory for combinatorial multivector fields
1 Discrete, Computational and lgebraic Topology Copenhagen, 12th November 2014 Morse-Forman-Conley theory for combinatorial multivector fields (research in progress) P Q L M N O G H I J K C D E F Marian
More informationAn Algorithmic Approach to Chain Recurrence
An Algorithmic Approach to Chain Recurrence W.D. Kalies, K. Mischaikow, and R.C.A.M. VanderVorst October 10, 2005 ABSTRACT. In this paper we give a new definition of the chain recurrent set of a continuous
More informationThe Existence of Chaos in the Lorenz System
The Existence of Chaos in the Lorenz System Sheldon E. Newhouse Mathematics Department Michigan State University E. Lansing, MI 48864 joint with M. Berz, K. Makino, A. Wittig Physics, MSU Y. Zou, Math,
More informationDynamic Systems and Applications 11 (2002) TOWARDS A RIGOROUS NUMERICAL STUDY OF THE KOT SCHAFFER MODEL
Dynamic Systems and Applications (2002) 87-97 TOWARDS A RIGOROUS NUMERICAL STUDY OF THE KOT SCHAFFER MODEL SARAH DAY Georgia Institute of Technology, Department of Mathematics, Atlanta, GA 30332 090, USA.
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationKey words. database, dynamical system, Conley index, Morse decomposition, Leslie population models, combinatorial dynamics, multiparameter system
DATABASES FOR THE GLOBAL DYNAMICS OF MULTIPARAMETER SYSTEMS ZIN ARAI, WILLIAM KALIES, HIROSHI KOKUBU, KONSTANTIN MISCHAIKOW, HIROE OKA, AND PAWE L PILARCZYK Abstract. A generally applicable, automatic
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationAlgebra Homework, Edition 2 9 September 2010
Algebra Homework, Edition 2 9 September 2010 Problem 6. (1) Let I and J be ideals of a commutative ring R with I + J = R. Prove that IJ = I J. (2) Let I, J, and K be ideals of a principal ideal domain.
More informationDynamical Systems and Ergodic Theory PhD Exam Spring Topics: Topological Dynamics Definitions and basic results about the following for maps and
Dynamical Systems and Ergodic Theory PhD Exam Spring 2012 Introduction: This is the first year for the Dynamical Systems and Ergodic Theory PhD exam. As such the material on the exam will conform fairly
More informationNonlinear Dynamical Systems Eighth Class
Nonlinear Dynamical Systems Eighth Class Alexandre Nolasco de Carvalho September 19, 2017 Now we exhibit a Morse decomposition for a dynamically gradient semigroup and use it to prove that a dynamically
More informationComputational homology in dynamical systems
1 Computational homology in dynamical systems Marian Mrozek Jagiellonian University, Kraków National University of Ireland, Galway 29 June - 10 July 2009 Outline 2 Dynamical systems Rigorous numerics of
More informationB553 Lecture 1: Calculus Review
B553 Lecture 1: Calculus Review Kris Hauser January 10, 2012 This course requires a familiarity with basic calculus, some multivariate calculus, linear algebra, and some basic notions of metric topology.
More informationComparative Network Analysis
Comparative Network Analysis BMI/CS 776 www.biostat.wisc.edu/bmi776/ Spring 2016 Anthony Gitter gitter@biostat.wisc.edu These slides, excluding third-party material, are licensed under CC BY-NC 4.0 by
More informationProblem 1. Let f n (x, y), n Z, be the sequence of rational functions in two variables x and y given by the initial conditions.
18.217 Problem Set (due Monday, December 03, 2018) Solve as many problems as you want. Turn in your favorite solutions. You can also solve and turn any other claims that were given in class without proofs,
More informationNonlinear Dynamical Systems Ninth Class
Nonlinear Dynamical Systems Ninth Class Alexandre Nolasco de Carvalho September 21, 2017 Lemma Let {T (t) : t 0} be a dynamically gradient semigroup in a metric space X, with a global attractor A and a
More informationThe Structure of Hyperbolic Sets
The Structure of Hyperbolic Sets p. 1/35 The Structure of Hyperbolic Sets Todd Fisher tfisher@math.umd.edu Department of Mathematics University of Maryland, College Park The Structure of Hyperbolic Sets
More informationKey words. database, dynamical system, Conley index, Morse decomposition, Leslie population models, combinatorial dynamics, multiparameter system
A DATABASE SCHEMA FOR THE ANALYSIS OF GLOBAL DYNAMICS OF MULTIPARAMETER SYSTEMS ZIN ARAI, WILLIAM KALIES, HIROSHI KOKUBU, KONSTANTIN MISCHAIKOW, HIROE OKA, AND PAWE L PILARCZYK Abstract. A generally applicable,
More informationThe Conley Index and Rigorous Numerics of Attracting Periodic Orbits
The Conley Index and Rigorous Numerics of Attracting Periodic Orbits Marian Mrozek Pawe l Pilarczyk Conference on Variational and Topological Methods in the Study of Nonlinear Phenomena (Pisa, 2000) 1
More informationEquivalence of Graded Module Braids and Interlocking Sequences
J Math Kyoto Univ (JMKYAZ) (), Equivalence of Graded Module Braids and Interlocking Sequences By Zin ARAI Abstract The category of totally ordered graded module braids and that of the exact interlocking
More informationKevin James. MTHSC 102 Section 4.3 Absolute Extreme Points
MTHSC 102 Section 4.3 Absolute Extreme Points Definition (Relative Extreme Points and Relative Extreme Values) Suppose that f(x) is a function defined on an interval I (possibly I = (, ). 1 We say that
More informationIntroduction Knot Theory Nonlinear Dynamics Topology in Chaos Open Questions Summary. Topology in Chaos
Introduction Knot Theory Nonlinear Dynamics Open Questions Summary A tangled tale about knot, link, template, and strange attractor Centre for Chaos & Complex Networks City University of Hong Kong Email:
More informationHyperbolic Dynamics. Todd Fisher. Department of Mathematics University of Maryland, College Park. Hyperbolic Dynamics p.
Hyperbolic Dynamics p. 1/36 Hyperbolic Dynamics Todd Fisher tfisher@math.umd.edu Department of Mathematics University of Maryland, College Park Hyperbolic Dynamics p. 2/36 What is a dynamical system? Phase
More informationMorse Theory for Lagrange Multipliers
Morse Theory for Lagrange Multipliers γ=0 grad γ grad f Guangbo Xu Princeton University and Steve Schecter North Carolina State University 1 2 Outline 3 (1) History of Mathematics 1950 2050. Chapter 1.
More information1 Introduction Definitons Markov... 2
Compact course notes Dynamic systems Fall 2011 Professor: Y. Kudryashov transcribed by: J. Lazovskis Independent University of Moscow December 23, 2011 Contents 1 Introduction 2 1.1 Definitons...............................................
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Centre for Research on Inner City Health St Michael s Hospital Toronto On leave from Department of Mathematics University of Manitoba Julien Arino@umanitoba.ca
More informationSolution to Homework #4 Roy Malka
1. Show that the map: Solution to Homework #4 Roy Malka F µ : x n+1 = x n + µ x 2 n x n R (1) undergoes a saddle-node bifurcation at (x, µ) = (0, 0); show that for µ < 0 it has no fixed points whereas
More informationWILLIAM D. KALIES CURRICULUM VITAE FEBRUARY 2018
WILLIAM D. KALIES CURRICULUM VITAE FEBRUARY 2018 Professor Phone: (561) 297-1107 Department of Mathematical Sciences FAX: (561) 297-2436 Florida Atlantic University E-mail: wkalies@fau.edu Boca Raton,
More informationGLOBEX Bioinformatics (Summer 2015) Genetic networks and gene expression data
GLOBEX Bioinformatics (Summer 2015) Genetic networks and gene expression data 1 Gene Networks Definition: A gene network is a set of molecular components, such as genes and proteins, and interactions between
More informationFOURTH ORDER CONSERVATIVE TWIST SYSTEMS: SIMPLE CLOSED CHARACTERISTICS
FOURTH ORDER CONSERVATIVE TWIST SYSTEMS: SIMPLE CLOSED CHARACTERISTICS J.B. VAN DEN BERG AND R.C.A.M. VANDERVORST ABSTRACT. On the energy manifolds of fourth order conservative systems closed characteristics
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationCOMPUTATION OF CONNECTION MATRICES USING THE SOFTWARE PACKAGE conley
COMPUTATION OF CONNECTION MATRICES USING THE SOFTWARE PACKAGE conley MOHAMED BARAKAT AND STANISLAUS MAIER-PAAPE Abstract. In this paper we demonstrate the power of the computer algebra package conley,
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More information1. Group Theory Permutations.
1.1. Permutations. 1. Group Theory Problem 1.1. Let G be a subgroup of S n of index 2. Show that G = A n. Problem 1.2. Find two elements of S 7 that have the same order but are not conjugate. Let π S 7
More informationMATH 614 Dynamical Systems and Chaos Lecture 2: Periodic points. Hyperbolicity.
MATH 614 Dynamical Systems and Chaos Lecture 2: Periodic points. Hyperbolicity. Orbit Let f : X X be a map defining a discrete dynamical system. We use notation f n for the n-th iteration of f defined
More informationWHAT IS A CHAOTIC ATTRACTOR?
WHAT IS A CHAOTIC ATTRACTOR? CLARK ROBINSON Abstract. Devaney gave a mathematical definition of the term chaos, which had earlier been introduced by Yorke. We discuss issues involved in choosing the properties
More informationMA102: Multivariable Calculus
MA102: Multivariable Calculus Rupam Barman and Shreemayee Bora Department of Mathematics IIT Guwahati Differentiability of f : U R n R m Definition: Let U R n be open. Then f : U R n R m is differentiable
More informationEquivalence relations
Equivalence relations R A A is an equivalence relation if R is 1. reflexive (a, a) R 2. symmetric, and (a, b) R (b, a) R 3. transitive. (a, b), (b, c) R (a, c) R Example: Let S be a relation on people
More information(Non-)Existence of periodic orbits in dynamical systems
(Non-)Existence of periodic orbits in dynamical systems Konstantin Athanassopoulos Department of Mathematics and Applied Mathematics University of Crete June 3, 2014 onstantin Athanassopoulos (Univ. of
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationPropp-Wilson Algorithm (and sampling the Ising model)
Propp-Wilson Algorithm (and sampling the Ising model) Danny Leshem, Nov 2009 References: Haggstrom, O. (2002) Finite Markov Chains and Algorithmic Applications, ch. 10-11 Propp, J. & Wilson, D. (1996)
More informationDYNAMICAL SYSTEMS PROBLEMS. asgor/ (1) Which of the following maps are topologically transitive (minimal,
DYNAMICAL SYSTEMS PROBLEMS http://www.math.uci.edu/ asgor/ (1) Which of the following maps are topologically transitive (minimal, topologically mixing)? identity map on a circle; irrational rotation of
More informationTEST CODE: PMB SYLLABUS
TEST CODE: PMB SYLLABUS Convergence and divergence of sequence and series; Cauchy sequence and completeness; Bolzano-Weierstrass theorem; continuity, uniform continuity, differentiability; directional
More informationfy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f))
1. Basic algebra of vector fields Let V be a finite dimensional vector space over R. Recall that V = {L : V R} is defined to be the set of all linear maps to R. V is isomorphic to V, but there is no canonical
More informationReading: Chapter 9.3. Carnegie Mellon
I II Lecture 3 Foundation of Data Flow Analysis Semi-lattice (set of values, meet operator) Transfer functions III Correctness, precision and convergence IV Meaning of Data Flow Solution Reading: Chapter
More informationHandlebody Decomposition of a Manifold
Handlebody Decomposition of a Manifold Mahuya Datta Statistics and Mathematics Unit Indian Statistical Institute, Kolkata mahuya@isical.ac.in January 12, 2012 contents Introduction What is a handlebody
More informationIsomorphism of free G-subflows (preliminary draft)
Isomorphism of free G-subflows (preliminary draft) John D. Clemens August 3, 21 Abstract We show that for G a countable group which is not locally finite, the isomorphism of free G-subflows is bi-reducible
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationSIMPLE MULTIVARIATE OPTIMIZATION
SIMPLE MULTIVARIATE OPTIMIZATION 1. DEFINITION OF LOCAL MAXIMA AND LOCAL MINIMA 1.1. Functions of variables. Let f(, x ) be defined on a region D in R containing the point (a, b). Then a: f(a, b) is a
More informationComputation of homoclinic and heteroclinic orbits for flows
Computation of homoclinic and heteroclinic orbits for flows Jean-Philippe Lessard BCAM BCAM Mini-symposium on Computational Math February 1st, 2011 Rigorous Computations Connecting Orbits Compute a set
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Department of Mathematics University of Manitoba Winnipeg Julien Arino@umanitoba.ca 19 May 2012 1 Introduction 2 Stochastic processes 3 The SIS model
More information- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs
LP-Duality ( Approximation Algorithms by V. Vazirani, Chapter 12) - Well-characterized problems, min-max relations, approximate certificates - LP problems in the standard form, primal and dual linear programs
More informationComplexity Theory VU , SS The Polynomial Hierarchy. Reinhard Pichler
Complexity Theory Complexity Theory VU 181.142, SS 2018 6. The Polynomial Hierarchy Reinhard Pichler Institut für Informationssysteme Arbeitsbereich DBAI Technische Universität Wien 15 May, 2018 Reinhard
More informationOutline. Complexity Theory EXACT TSP. The Class DP. Definition. Problem EXACT TSP. Complexity of EXACT TSP. Proposition VU 181.
Complexity Theory Complexity Theory Outline Complexity Theory VU 181.142, SS 2018 6. The Polynomial Hierarchy Reinhard Pichler Institut für Informationssysteme Arbeitsbereich DBAI Technische Universität
More informationSHIFT EQUIVALENCE AND A COMBINATORIAL-TOPOLOGICAL APPROACH TO DISCRETE-TIME DYNAMICAL SYSTEMS
SHIFT EQUIVALENCE AND A COMBINATORIAL-TOPOLOGICAL APPROACH TO DISCRETE-TIME DYNAMICAL SYSTEMS BY JUSTIN BUSH A dissertation submitted to the Graduate School New Brunswick Rutgers, The State University
More informationAnalysis-3 lecture schemes
Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationLorenz like flows. Maria José Pacifico. IM-UFRJ Rio de Janeiro - Brasil. Lorenz like flows p. 1
Lorenz like flows Maria José Pacifico pacifico@im.ufrj.br IM-UFRJ Rio de Janeiro - Brasil Lorenz like flows p. 1 Main goals The main goal is to explain the results (Galatolo-P) Theorem A. (decay of correlation
More informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
More informationApril 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.
April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set
More informationMTH 132 Solutions to Exam 2 Nov. 23rd 2015
Name: Section: Instructor: READ THE FOLLOWING INSTRUCTIONS. Do not open your exam until told to do so. No calculators, cell phones or any other electronic devices can be used on this exam. Clear your desk
More information10725/36725 Optimization Homework 2 Solutions
10725/36725 Optimization Homework 2 Solutions 1 Convexity (Kevin) 1.1 Sets Let A R n be a closed set with non-empty interior that has a supporting hyperplane at every point on its boundary. (a) Show that
More informationApplications. More Counting Problems. Complexity of Algorithms
Recurrences Applications More Counting Problems Complexity of Algorithms Part I Recurrences and Binomial Coefficients Paths in a Triangle P(0, 0) P(1, 0) P(1,1) P(2, 0) P(2,1) P(2, 2) P(3, 0) P(3,1) P(3,
More informationMATH 61-02: PRACTICE PROBLEMS FOR FINAL EXAM
MATH 61-02: PRACTICE PROBLEMS FOR FINAL EXAM (FP1) The exclusive or operation, denoted by and sometimes known as XOR, is defined so that P Q is true iff P is true or Q is true, but not both. Prove (through
More informationList of topics for the preliminary exam in algebra
List of topics for the preliminary exam in algebra 1 Basic concepts 1. Binary relations. Reflexive, symmetric/antisymmetryc, and transitive relations. Order and equivalence relations. Equivalence classes.
More informationg 2 (x) (1/3)M 1 = (1/3)(2/3)M.
COMPACTNESS If C R n is closed and bounded, then by B-W it is sequentially compact: any sequence of points in C has a subsequence converging to a point in C Conversely, any sequentially compact C R n is
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationImmerse Metric Space Homework
Immerse Metric Space Homework (Exercises -2). In R n, define d(x, y) = x y +... + x n y n. Show that d is a metric that induces the usual topology. Sketch the basis elements when n = 2. Solution: Steps
More informationDynamical modelling of systems of coupled oscillators
Dynamical modelling of systems of coupled oscillators Mathematical Neuroscience Network Training Workshop Edinburgh Peter Ashwin University of Exeter 22nd March 2009 Peter Ashwin (University of Exeter)
More informationPopulation Games and Evolutionary Dynamics
Population Games and Evolutionary Dynamics William H. Sandholm The MIT Press Cambridge, Massachusetts London, England in Brief Series Foreword Preface xvii xix 1 Introduction 1 1 Population Games 2 Population
More informationSingular perturbation analysis of an additive increase multiplicative decrease control algorithm under time-varying buffering delays.
Singular perturbation analysis of an additive increase multiplicative decrease control algorithm under time-varying buffering delays. V. Guffens 1 and G. Bastin 2 Intelligent Systems and Networks Research
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More informationSynchronous state transition graph
Heike Siebert, FU Berlin, Molecular Networks WS10/11 2-1 Synchronous state transition graph (0, 2) (1, 2) vertex set X (state space) edges (x,f(x)) every state has only one successor attractors are fixed
More informationProof: The coding of T (x) is the left shift of the coding of x. φ(t x) n = L if T n+1 (x) L
Lecture 24: Defn: Topological conjugacy: Given Z + d (resp, Zd ), actions T, S a topological conjugacy from T to S is a homeomorphism φ : M N s.t. φ T = S φ i.e., φ T n = S n φ for all n Z + d (resp, Zd
More informationMultiple Choice Questions for Review
Equivalence and Order Multiple Choice Questions for Review In each case there is one correct answer (given at the end of the problem set). Try to work the problem first without looking at the answer. Understand
More informationAbstract Algebra: Chapters 16 and 17
Study polynomials, their factorization, and the construction of fields. Chapter 16 Polynomial Rings Notation Let R be a commutative ring. The ring of polynomials over R in the indeterminate x is the set
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationInterpolation and Approximation
Interpolation and Approximation The Basic Problem: Approximate a continuous function f(x), by a polynomial p(x), over [a, b]. f(x) may only be known in tabular form. f(x) may be expensive to compute. Definition:
More informationICS141: Discrete Mathematics for Computer Science I
ICS141: Discrete Mathematics for Computer Science I Dept. Information & Computer Sci., Jan Stelovsky based on slides by Dr. Baek and Dr. Still Originals by Dr. M. P. Frank and Dr. J.L. Gross Provided by
More informationHYPERBOLIC SETS WITH NONEMPTY INTERIOR
HYPERBOLIC SETS WITH NONEMPTY INTERIOR TODD FISHER, UNIVERSITY OF MARYLAND Abstract. In this paper we study hyperbolic sets with nonempty interior. We prove the folklore theorem that every transitive hyperbolic
More information1. Find the domain of the function f(x) = (A) D = (-, 2) (2, ) (B) D = (-, -2) (-2, ) (C) D = (2, ) (D) D = (-2, ) (E) none of the above
1. Find the domain of the function f(x) = (A) D = (-, ) (, ) (B) D = (-, -) (-, ) (C) D = (, ) (D) D = (-, ) x. x + 1). A cell phone company offers a plan of $39 per month for 400 minutes. Each additional
More informationREAL VARIABLES: PROBLEM SET 1. = x limsup E k
REAL VARIABLES: PROBLEM SET 1 BEN ELDER 1. Problem 1.1a First let s prove that limsup E k consists of those points which belong to infinitely many E k. From equation 1.1: limsup E k = E k For limsup E
More informationMultivariate Topological Data Analysis
Cleveland State University November 20, 2008, Duke University joint work with Gunnar Carlsson (Stanford), Peter Kim and Zhiming Luo (Guelph), and Moo Chung (Wisconsin Madison) Framework Ideal Truth (Parameter)
More information25.1 Ergodicity and Metric Transitivity
Chapter 25 Ergodicity This lecture explains what it means for a process to be ergodic or metrically transitive, gives a few characterizes of these properties (especially for AMS processes), and deduces
More informationComputational Homology in Topological Dynamics
1 Computational Homology in Topological Dynamics ACAT School, Bologna, Italy MAY 26, 2012 Marian Mrozek Jagiellonian University, Kraków Outline 2 Dynamical systems Rigorous numerics of dynamical systems
More informationModeling with Exponential Functions
CHAPTER Modeling with Exponential Functions A nautilus is a sea creature that lives in a shell. The cross-section of a nautilus s shell, with its spiral of ever-smaller chambers, is a natural example of
More information7. Shortest Path Problems and Deterministic Finite State Systems
7. Shortest Path Problems and Deterministic Finite State Systems In the next two lectures we will look at shortest path problems, where the objective is to find the shortest path from a start node to an
More informationDynamical Systems in Biology
Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems
More informationPersistence and Stationary Distributions of Biochemical Reaction Networks
Persistence and Stationary Distributions of Biochemical Reaction Networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Discrete Models in Systems Biology SAMSI December
More informationConley Theory: A combinatorial approach to dynamics
Conley Theory: A combinatorial approach to dynamics Konstantin Mischaikow November 21, 2010 Contents 1 Why? 5 1.1 Physics............................... 5 1.2 Biology............................... 6 1.3
More informationNetworks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationWritten Exam 15 December Course name: Introduction to Systems Biology Course no
Technical University of Denmark Written Exam 15 December 2008 Course name: Introduction to Systems Biology Course no. 27041 Aids allowed: Open book exam Provide your answers and calculations on separate
More informationNBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011
NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets
More informationAn Introduction to Entropy and Subshifts of. Finite Type
An Introduction to Entropy and Subshifts of Finite Type Abby Pekoske Department of Mathematics Oregon State University pekoskea@math.oregonstate.edu August 4, 2015 Abstract This work gives an overview
More informationThe method of topological sections in the rigorous numerics of dynamical systems
1 Marian Mrozek Jagiellonian University, Kraków The method of topological sections in the rigorous numerics of dynamical systems Fourth International Workshop on Taylor Methods Boca Raton, December 17th,
More information