Linear Algebra of Eigen s Quasispecies Model
|
|
- Rosemary Jasmine Wade
- 5 years ago
- Views:
Transcription
1 Linear Algebra of Eigen s Quasispecies Model Artem Novozhilov Department of Mathematics North Dakota State University Midwest Mathematical Biology Conference, University of Wisconsin La Crosse May 17, 2014 Artem Novozhilov (NDSU) May 17, / 29
2 Collaborators: Yuri S. Semenov Moscow State University of Railway Engineering, Moscow, Russia Alexander S. Bratus Lomonosov Moscow State University Moscow, Russia Artem Novozhilov (NDSU) May 17, / 29
3 Outline: Motivation and historic outline Mathematical problem Known results New results Artem Novozhilov (NDSU) May 17, / 29
4 M. Eigen, Naturwisenschaften, 58(10), 1971: Artem Novozhilov (NDSU) May 17, / 29
5 Quasispecies theory: M. Eigen, Naturwisenschaften, 58(10), 1971: M. Eigen, P. Schuster, The Hypercycle, Springer, 1979 M. Eigen, J. McCaskill, P. Schuster, J Phys Chem, 92(24), 1988: Manfred Eigen, born 1927 Artem Novozhilov (NDSU) May 17, / 29
6 Mathematical side of the story: Ref: M. Eigen, J. McCaskill, P. Schuster, J Phys Chem, 92(24), 1988: Artem Novozhilov (NDSU) May 17, / 29
7 Model statement: Consider a population of sequences of fixed length N composed of two-letter alphabet, say, {0,1}, therefore 2 N different sequences. The population is subject to two evolutionary forces. First evolutionary force is selection, which is included in the system through the Malthusian fitness, defined here for simplicity as m(particular sequence σ) = m(h σ ), where H σ is the Hamming norm of this sequence, i.e., number of 1s in sequence σ. In this way we do not distinguish between sequences with the same number of 1s and hence reduce the dimensionality of the problem from 2 N 2 N to (N +1) (N +1). Hence, we consider at this point only permutation invariant fitness landscapes M = diag(m 0,...,m N ) or m = (m 0,...,m N ). Artem Novozhilov (NDSU) May 17, / 29
8 Model statement: The second evolutionary force is mutation. In particular, assuming N +1 classes of sequences, we have that the mutations µ ij (i.e., the mutation rate from class j to class i) can be described by the matrix N N N N 1 N M = (µ ij ) = µq = 0 0 N 2 N , N N N where µ is the mutation rate per site per sequence per replication event. Artem Novozhilov (NDSU) May 17, / 29
9 Model statement: Let p(t) denote the vector of frequencies of different classes of sequences, then, assuming uncoupled reproduction and mutation events, we arrive at where is the mean population fitness. ṗ(t) = (M +µq)p(t) m(t)p(t), m(t) = m p(t) = N m i p i (t) This model is often called a paramuse of Crow Kimura quasispecies model with permutation invariant fitness landscape. Ref: Baake and Gabriel, Annual Reviews of Computational Physics VII, 1999: Ref: Crow and Kimura, An introduction to population genetics theory, 1970 i=0 Artem Novozhilov (NDSU) May 17, / 29
10 Elementary results: The asymptotic behavior of the quasispecies model is determined by the equilibrium p = lim t p(t), which solves the eigenvalue problem where (M +µq)p = mp, m = m p. By Perron Frobenius theorem it follows that there is a unique positive solution p > 0, which is the right eigenvector of M +µq corresponding to the simple real dominant eigenvalue λ = m. This vector p was called by Eigen the quasispecies. It is globally stable for the quasispecies system. We are mostly interested in properties of m and p depending on the fitness landscape M and mutation rate µ, therefore, we use the notation m = m(µ) and p = p(µ) for the mean fitness and equilibrium distribution. Artem Novozhilov (NDSU) May 17, / 29
11 Known results: Thompson and McBride, Math Biosciences, 21: , 1974 The quasispecies model is essentially linear. Rumschitzki, J of Math Biol, 24: , 1987 For zero epistasis fitness landscape the spectral properties of the Eigen evolutionary matrices can be inferred with the representation of M and M using tensor products. Swetina and Schuster, Bioph Chem, 16: , 1982 Numerical analysis of single peaked fitness landscape yields the error threshold. Artem Novozhilov (NDSU) May 17, / 29
12 Known results: The error threshold Ref: Swetina and Schuster, Bioph Chem, 16: , 1982 Consider the single peaked fitness landscape M = diag(m 0,0,...,0), m 0 > 0. Artem Novozhilov (NDSU) May 17, / 29
13 Known results: The error threshold Consider the single peaked fitness landscape M = diag(m 0,0,...,0), m 0 > p µ Artem Novozhilov (NDSU) May 17, / 29
14 Known results: Statistical Physics Leuthäusser, I. (1986). An exact correspondence between Eigen s evolution model and a two-dimensional Ising system. The Journal of Chemical Physics, 84(3), Tarazona, P. (1992). Error thresholds for molecular quasispecies as phase transitions: From simple landscapes to spin-glass models. Physical Review A, 45(8), Baake, E., Baake, M., & Wagner, H. (1997). Ising quantum chain is equivalent to a model of biological evolution. Physical Review Letters, 78(3), Galluccio, S. (1997). Exact solution of the quasispecies model in a sharply peaked fitness landscape. Physical Review E, 56(4), Baake, E., & Wagner, H. (2001). Mutation-selection models solved exactly with methods of statistical mechanics. Genetical Research, 78(01), Saakian, D. B., & Hu, C. K. (2006). Exact solution of the Eigen model with general fitness functions and degradation rates. Proceedings of the National Academy of Sciences of the USA, 103(13), Artem Novozhilov (NDSU) May 17, / 29
15 Known results: Maximum principle Hermisson, J., Redner, O., Wagner, H., & Baake, E. (2002). Mutation-selection balance: Ancestry, load, and maximum principle. Theoretical Population Biology, 62(1), Baake, E., & Georgii, H. O. (2007). Mutation, selection, and ancestry in branching models: a variational approach. Journal of Mathematical Biology, 54(2), Assume that m i = Nr i = Nr(x i ), x i = i N [0,1] and define g(x) = µ ( 1 2 x(1 x) ). Then the mean fitness m(µ) = Nr is given by ( ) r r = sup r(x) g(x). x [0,1] r(x) may have only finite number of discontinuities and be either left or right continuous at every point. Artem Novozhilov (NDSU) May 17, / 29
16 Main idea: We consider the eigenvalue problem (M +µq)p = mp, m = m p, where p = p(µ), m = m(µ) with a fixed fitness landscape m. We claim that this problem simplifies in the coordinates of the basis composed of the eigenvectors of the matrix Q = Q N. Artem Novozhilov (NDSU) May 17, / 29
17 Proposition: For the matrix Q = Q N : 1. The eigenvalues of Q N are simple (all have algebraic multiplicities one) and given by q k = 2k, k = 0,...,N. 2. Let v k = (c 0k,...,c Nk ) be the right eigenvector of Q N corresponding to q k and normalized such that c 0k = 1, C = C N = (c ik ) (N+1) (N+1) be the matrix composed of v k (v k is the k-th column of C N). Then the generating function for the elements of the k-th column has the form P k (t) = N c ik t i = (1 t) k (1+t) N k, k = 0,...,N. i=0 3. C 2 = 2 N I, where I is the identity matrix, or, equivalently, C 1 = 2 N C norm of C is C 1 = max 0 k N N c ik = 2 N. i=0 Artem Novozhilov (NDSU) May 17, / 29
18 Behavior for µ : Let x = C 1 p and ˆx = 2 N (1,0,...,0). Then x ˆx 1 1 µ M 1, and Therefore, x (µ) 1 p(µ) ˆp = 2 N (( N 0 2N 2µ (2 N+1 +1) M 1. ) ( )) N,...,. N Artem Novozhilov (NDSU) May 17, / 29
19 Parametric solution: Consider fitness landscape such that m j > 0 for some j, and all the rest m i = 0 for i j. Then x k (s) = 2 N c kj 1+ks, p i (s) = 2 N N k=0 m(s) = m j p j (s), µ = s 2 m(s). c ik c kj 1+ks, Artem Novozhilov (NDSU) May 17, / 29
20 Example 1: Let N = 2A and m = (0,...,0,N,0,...,0), where N is exactly on A-th position. Consider a scaled fitness landscape Nr = m and r(µ) = m(µ)/n. Then, using the properties of the mutation matrix Q and the parametric formulas from the previous slide, it can be proved that r r = lim N r(µ) = µ 2 +1 µ. 1 N N 200 r Μ r Μ 0 3 Μ 0 3 Μ Artem Novozhilov (NDSU) May 17, / 29
21 Example 1: Let N = 2A and m = (0,...,0,N,0,...,0), where N is exactly on A-th position. Consider a scaled fitness landscape Nr = m and r(µ) = m(µ)/n. Denoting r = µ 2 +1 µ we can prove that lim N p A±k(µ) = r ( ) k 1 r, k = 0,...,A. 1+r p p µ µ Artem Novozhilov (NDSU) May 17, / 29
22 Example 2: Single peaked landscape Let m = (N,0,...,0). Consider a scaled fitness landscape Nr = m and r(µ) = m(µ)/n. Then, using the properties of the mutation matrix Q and the parametric formulas we can prove that if µ < 1 and p j = 0 for all j if µ 1. r = lim p 0 = 1 µ, N lim p i = (1 µ)µ i, i 1. N Artem Novozhilov (NDSU) May 17, / 29
23 Example 2: Single peaked landscape Comparison with numerical computations for N = 100. r = lim N p 0 = 1 µ, lim p i = (1 µ)µ i, i 1 N r(µ) p µ µ Artem Novozhilov (NDSU) May 17, / 29
24 Epsilon stabilization: Definition: Dominant eigenvalue m(µ) admits epsilon stabilization if for a small enough ε > 0 there exists constant m ε and µ ε such that for all µ > µ ε we have m(µ) m ǫ < ε, m (µ) < ε. We know that epsilon stabilization is a sure event, due to the previous. The question is how to determine µ ε given ε > 0. Artem Novozhilov (NDSU) May 17, / 29
25 Epsilon stabilization: Let m = (m 0,m 1,...,m N ) such that m 0 > m 1 m 2... m N. Then we have µ = m 0 m 1, (Classical Result) N ( ) µ ε = (m 0 m 1 ) 1 1 2(m 0 1), (m 0 m 1 )N N ( ) (Our result 1) N if δ = 2 N (m k 1) < ε, k k=0 m(µ ) = m min +2µ. (Our result 2) Artem Novozhilov (NDSU) May 17, / 29
26 p Example: m Μ Μ Μ Figure : Error threshold in the quasispecies model with the single picked fitness landscape Estimates from the previous slide: µ 1 = 0.633, µ 2 = 0.644, µ 3 = Artem Novozhilov (NDSU) May 17, / 29
27 Example: p 0.15 m Μ Μ Μ Figure : Error threshold in the quasispecies model with m k = k α log(1 s) Estimates from the previous slide: µ 1 = 0.09, µ 2 = 0.19, µ 3 = 0.33 Artem Novozhilov (NDSU) May 17, / 29
28 Acknowledgements: Startup grant from Department of Mathematics, NDSU ND EPSCoR and NSF grant # EPS Artem Novozhilov (NDSU) May 17, / 29
29 Thank you for your attention! site: novozhil/ References: Bratus, A. S., Novozhilov, A. S., & Semenov, Y. S. (2013). Linear algebra of the permutation invariant Crow Kimura model of prebiotic evolution. arxiv: Semenov, Y. S., Bratus, A. S., & Novozhilov, A. S. (2014). On the behavior of the leading eigenvalue of the Eigen evolutionary matrices, in preparation. Artem Novozhilov (NDSU) May 17, / 29
Linear Algebra of the Quasispecies Model
Linear Algebra of the Quasispecies Model Artem Novozhilov Department of Mathematics North Dakota State University SIAM Conference on the Life Sciences, Charlotte, North Caroline, USA August 6, 2014 Artem
More informationarxiv: v1 [q-bio.pe] 11 Dec 2017
arxiv:1712.03855v1 [q-bio.pe] 11 Dec 2017 Rigorous mathematical analysis of the quasispecies model: From Manfred Eigen to the recent developments Alexander S. Bratus 1,2, Artem S. Novozhilov 3,, Yuri S.
More informationarxiv: v1 [q-bio.pe] 3 May 2017
arxiv:1705.01568v1 [q-bio.pe] 3 May 017 Adaptive Fitness Landscape for Replicator Systems: To Maximize or not to Maximize Alexander S. Bratus 1,,, Artem S. Novozhilov 3,, Yuri S. Semenov, 1 Faculty of
More informationRelaxation Property and Stability Analysis of the Quasispecies Models
Commun. Theor. Phys. (Beijing, China) 52 (2009) pp. 726 734 c Chinese Physical Society and IOP Publishing Ltd Vol. 52, No. 4, October 15, 2009 Relaxation Property and Stability Analysis of the Quasispecies
More informationQuasispecies distribution of Eigen model
Vol 16 No 9, September 2007 c 2007 Chin. Phys. Soc. 1009-1963/2007/16(09)/2600-08 Chinese Physics and IOP Publishing Ltd Quasispecies distribution of Eigen model Chen Jia( ), Li Sheng( ), and Ma Hong-Ru(
More informationarxiv:physics/ v1 [physics.bio-ph] 27 Jun 2001
Maternal effects in molecular evolution Claus O. Wilke Digital Life Laboratory, Mail Code 36-93, Caltech, Pasadena, CA 925 wilke@caltech.edu (Printed: May 3, 27) arxiv:physics/693v [physics.bio-ph] 27
More informationError thresholds on realistic fitness landscapes
Error thresholds on realistic fitness landscapes Peter Schuster Institut für Theoretische Chemie, Universität Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA Evolutionary Dynamics:
More informationMolecular models of biological evolution have attracted
Exact solution of the Eigen model with general fitness functions and degradation rates David B. Saakian* and Chin-Kun Hu* *Institute of Physics, Academia Sinica, ankang, Taipei 11529, Taiwan; and Yerevan
More informationA Method for Estimating Mean First-passage Time in Genetic Algorithms
A Method for Estimating Mean First-passage Time in Genetic Algorithms Hiroshi Furutani Department of Information Science, Kyoto University of Education, Fushimi-ku, Kyoto, 612-8522 Japan This paper presents
More informationError threshold in simple landscapes
arxiv:cond-mat/9610028v1 [cond-mat.dis-nn] 4 Oct 1996 Error threshold in simple landscapes Silvio Franz ICTP, Strada Costiera 10, I-34100 Trieste (Italy) and Luca Peliti Groupe de Physico-Chimie Théorique,
More informationConvergence of a Moran model to Eigen s quasispecies model
Convergence of a Moran odel to Eigen s quasispecies odel Joseba Dalau Université Paris Sud and ENS Paris June 6, 2018 arxiv:1404.2133v1 [q-bio.pe] 8 Apr 2014 Abstract We prove that a Moran odel converges
More informationPerron Frobenius Theory
Perron Frobenius Theory Oskar Perron Georg Frobenius (1880 1975) (1849 1917) Stefan Güttel Perron Frobenius Theory 1 / 10 Positive and Nonnegative Matrices Let A, B R m n. A B if a ij b ij i, j, A > B
More informationDEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular
form) Given: matrix C = (c i,j ) n,m i,j=1 ODE and num math: Linear algebra (N) [lectures] c phabala 2016 DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix
More informationEvolving on a Single-peak Fitness Landscape. Lionel Barnett. Centre for the Study of Evolution. University of Sussex.
The Eects of Recombination on a Haploid Quasispecies Evolving on a Single-peak Fitness Landscape Lionel Barnett Centre for the Study of Evolution Centre for Computational Neuroscience and Robotics School
More informationExploration of population fixed-points versus mutation rates for functions of unitation
Exploration of population fixed-points versus mutation rates for functions of unitation J Neal Richter 1, Alden Wright 2, John Paxton 1 1 Computer Science Department, Montana State University, 357 EPS,
More informationThe Tightness of the Kesten-Stigum Reconstruction Bound for a Symmetric Model With Multiple Mutations
The Tightness of the Kesten-Stigum Reconstruction Bound for a Symmetric Model With Multiple Mutations City University of New York Frontier Probability Days 2018 Joint work with Dr. Sreenivasa Rao Jammalamadaka
More informationLinear Algebra Practice Problems
Math 7, Professor Ramras Linear Algebra Practice Problems () Consider the following system of linear equations in the variables x, y, and z, in which the constants a and b are real numbers. x y + z = a
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationFIXED POINT ITERATIONS
FIXED POINT ITERATIONS MARKUS GRASMAIR 1. Fixed Point Iteration for Non-linear Equations Our goal is the solution of an equation (1) F (x) = 0, where F : R n R n is a continuous vector valued mapping in
More informationHigh-dimensional Markov Chain Models for Categorical Data Sequences with Applications Wai-Ki CHING AMACL, Department of Mathematics HKU 19 March 2013
High-dimensional Markov Chain Models for Categorical Data Sequences with Applications Wai-Ki CHING AMACL, Department of Mathematics HKU 19 March 2013 Abstract: Markov chains are popular models for a modelling
More informationThe Miracles of Tropical Spectral Theory
The Miracles of Tropical Spectral Theory Emmanuel Tsukerman University of California, Berkeley e.tsukerman@berkeley.edu July 31, 2015 1 / 110 Motivation One of the most important results in max-algebra
More informationAn estimation of the spectral radius of a product of block matrices
Linear Algebra and its Applications 379 (2004) 267 275 wwwelseviercom/locate/laa An estimation of the spectral radius of a product of block matrices Mei-Qin Chen a Xiezhang Li b a Department of Mathematics
More informationMath 4153 Exam 3 Review. The syllabus for Exam 3 is Chapter 6 (pages ), Chapter 7 through page 137, and Chapter 8 through page 182 in Axler.
Math 453 Exam 3 Review The syllabus for Exam 3 is Chapter 6 (pages -2), Chapter 7 through page 37, and Chapter 8 through page 82 in Axler.. You should be sure to know precise definition of the terms we
More informationChemistry on the Early Earth
Chemistry on the Early Earth Peter Schuster Institut für Theoretische Chemie, Universität Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA Germany-Japan Round Table Heidelberg, 01. 03.11.2011
More information1 Introduction This work follows a paper by P. Shields [1] concerned with a problem of a relation between the entropy rate of a nite-valued stationary
Prexes and the Entropy Rate for Long-Range Sources Ioannis Kontoyiannis Information Systems Laboratory, Electrical Engineering, Stanford University. Yurii M. Suhov Statistical Laboratory, Pure Math. &
More informationPermutation transformations of tensors with an application
DOI 10.1186/s40064-016-3720-1 RESEARCH Open Access Permutation transformations of tensors with an application Yao Tang Li *, Zheng Bo Li, Qi Long Liu and Qiong Liu *Correspondence: liyaotang@ynu.edu.cn
More informationIn particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with
Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient
More informationSeries Solutions Near a Regular Singular Point
Series Solutions Near a Regular Singular Point MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background We will find a power series solution to the equation:
More informationMarkov Chains and Stochastic Sampling
Part I Markov Chains and Stochastic Sampling 1 Markov Chains and Random Walks on Graphs 1.1 Structure of Finite Markov Chains We shall only consider Markov chains with a finite, but usually very large,
More informationEvolution on Realistic Landscapes
Evolution on Realistic Landscapes Peter Schuster Institut für Theoretische Chemie, Universität Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA Santa Fe Institute Seminar Santa Fe, 22.05.2012
More informationTHE POINT SPECTRUM OF FROBENIUS-PERRON AND KOOPMAN OPERATORS
PROCEEDINGS OF THE MERICN MTHEMTICL SOCIETY Volume 126, Number 5, May 1998, Pages 1355 1361 S 0002-9939(98)04188-4 THE POINT SPECTRUM OF FROBENIUS-PERRON ND KOOPMN OPERTORS J. DING (Communicated by Palle
More informationIs the Concept of Error Catastrophy Relevant for Viruses? Peter Schuster
Is the Concept of Error Catastrophy Relevant for Viruses? Quasispecies and error thresholds on realistic landscapes Peter Schuster Institut für Theoretische Chemie, Universität Wien, Austria and The Santa
More informationSupporting Information
Supporting Information I. INFERRING THE ENERGY FUNCTION We downloaded a multiple sequence alignment (MSA) for the HIV-1 clade B Protease protein from the Los Alamos National Laboratory HIV database (http://www.hiv.lanl.gov).
More informationRecent progress on the classical and quantum Kuramoto synchronization
Recent progress on the classical and quantum Kuramoto synchronization Seung Yeal Ha Department of Mathematical Sciences Seoul National University March 23rd, 2017 Outline PROLOGUE What is synchronization?
More informationNon white sample covariance matrices.
Non white sample covariance matrices. S. Péché, Université Grenoble 1, joint work with O. Ledoit, Uni. Zurich 17-21/05/2010, Université Marne la Vallée Workshop Probability and Geometry in High Dimensions
More informationMarkov Chains, Stochastic Processes, and Matrix Decompositions
Markov Chains, Stochastic Processes, and Matrix Decompositions 5 May 2014 Outline 1 Markov Chains Outline 1 Markov Chains 2 Introduction Perron-Frobenius Matrix Decompositions and Markov Chains Spectral
More informationVast Volatility Matrix Estimation for High Frequency Data
Vast Volatility Matrix Estimation for High Frequency Data Yazhen Wang National Science Foundation Yale Workshop, May 14-17, 2009 Disclaimer: My opinion, not the views of NSF Y. Wang (at NSF) 1 / 36 Outline
More information05. Multiparticle Systems
05. Multiparticle Systems I. -Particle Product Spaces Suppose: Particle and particle are represented by vector spaces V and W. Then: The composite -particle system is represented by a product vector space
More informationarxiv:q-bio/ v1 [q-bio.bm] 8 Jun 2005
Mutation model for nucleotide sequences based on crystal basis arxiv:q-bio/0506010v1 [q-bio.bm] 8 Jun 2005 C. Minichini, A. Sciarrino Dipartimento di Scienze Fisiche, Università di Napoli Federico II and
More informationInvertibility and stability. Irreducibly diagonally dominant. Invertibility and stability, stronger result. Reducible matrices
Geršgorin circles Lecture 8: Outline Chapter 6 + Appendix D: Location and perturbation of eigenvalues Some other results on perturbed eigenvalue problems Chapter 8: Nonnegative matrices Geršgorin s Thm:
More informationMaximizing the numerical radii of matrices by permuting their entries
Maximizing the numerical radii of matrices by permuting their entries Wai-Shun Cheung and Chi-Kwong Li Dedicated to Professor Pei Yuan Wu. Abstract Let A be an n n complex matrix such that every row and
More informationSpectral Properties of Matrix Polynomials in the Max Algebra
Spectral Properties of Matrix Polynomials in the Max Algebra Buket Benek Gursoy 1,1, Oliver Mason a,1, a Hamilton Institute, National University of Ireland, Maynooth Maynooth, Co Kildare, Ireland Abstract
More informationPopulation Genetics: a tutorial
: a tutorial Institute for Science and Technology Austria ThRaSh 2014 provides the basic mathematical foundation of evolutionary theory allows a better understanding of experiments allows the development
More informationLab 8: Measuring Graph Centrality - PageRank. Monday, November 5 CompSci 531, Fall 2018
Lab 8: Measuring Graph Centrality - PageRank Monday, November 5 CompSci 531, Fall 2018 Outline Measuring Graph Centrality: Motivation Random Walks, Markov Chains, and Stationarity Distributions Google
More informationTo appear in Monatsh. Math. WHEN IS THE UNION OF TWO UNIT INTERVALS A SELF-SIMILAR SET SATISFYING THE OPEN SET CONDITION? 1.
To appear in Monatsh. Math. WHEN IS THE UNION OF TWO UNIT INTERVALS A SELF-SIMILAR SET SATISFYING THE OPEN SET CONDITION? DE-JUN FENG, SU HUA, AND YUAN JI Abstract. Let U λ be the union of two unit intervals
More informationError Propagation in the Hypercycle
Error Propagation in the Hypercycle Paulo R. A. Campos J. F. Fontanari Peter F. Stadler SFI WORKING PAPER: 1999-9-63 SFI Working Papers contain accounts of scientific work of the author(s) and do not necessarily
More informationOBSERVABILITY AND OBSERVERS IN A FOOD WEB
OBSERVABILITY AND OBSERVERS IN A FOOD WEB I. LÓPEZ, M. GÁMEZ, AND S. MOLNÁR Abstract. The problem of the possibility to recover the time-dependent state of a whole population system out of the observation
More informationSection 3.9. Matrix Norm
3.9. Matrix Norm 1 Section 3.9. Matrix Norm Note. We define several matrix norms, some similar to vector norms and some reflecting how multiplication by a matrix affects the norm of a vector. We use matrix
More informationSWITCHED SYSTEMS THAT ARE PERIODICALLY STABLE MAY BE UNSTABLE
SWITCHED SYSTEMS THAT ARE PERIODICALLY STABLE MAY BE UNSTABLE VINCENT D. BLONDEL, JACQUES THEYS, AND ALEXANDER A. VLADIMIROV 1. Introduction In this contribution, we prove the existence of switched linear
More informationReflections in Hilbert Space III: Eigen-decomposition of Szegedy s operator
Reflections in Hilbert Space III: Eigen-decomposition of Szegedy s operator James Daniel Whitfield March 30, 01 By three methods we may learn wisdom: First, by reflection, which is the noblest; second,
More informationReplica Condensation and Tree Decay
Replica Condensation and Tree Decay Arthur Jaffe and David Moser Harvard University Cambridge, MA 02138, USA Arthur Jaffe@harvard.edu, David.Moser@gmx.net June 7, 2007 Abstract We give an intuitive method
More informationNote that in the example in Lecture 1, the state Home is recurrent (and even absorbing), but all other states are transient. f ii (n) f ii = n=1 < +
Random Walks: WEEK 2 Recurrence and transience Consider the event {X n = i for some n > 0} by which we mean {X = i}or{x 2 = i,x i}or{x 3 = i,x 2 i,x i},. Definition.. A state i S is recurrent if P(X n
More informationNotes on Linear Algebra and Matrix Theory
Massimo Franceschet featuring Enrico Bozzo Scalar product The scalar product (a.k.a. dot product or inner product) of two real vectors x = (x 1,..., x n ) and y = (y 1,..., y n ) is not a vector but a
More informationMutation, selection, and ancestry in branching models: a variational approach
J. Math. Biol. (2007) 54:257 303 DOI 10.1007/s00285-006-0039-5 Mathematical Biology Mutation, selection, and ancestry in branching models: a variational approach Ellen Baake Hans-Otto Georgii Received:
More informationThe Particle Filter. PD Dr. Rudolph Triebel Computer Vision Group. Machine Learning for Computer Vision
The Particle Filter Non-parametric implementation of Bayes filter Represents the belief (posterior) random state samples. by a set of This representation is approximate. Can represent distributions that
More informationTotal Least Squares Approach in Regression Methods
WDS'08 Proceedings of Contributed Papers, Part I, 88 93, 2008. ISBN 978-80-7378-065-4 MATFYZPRESS Total Least Squares Approach in Regression Methods M. Pešta Charles University, Faculty of Mathematics
More informationMATH3200, Lecture 31: Applications of Eigenvectors. Markov Chains and Chemical Reaction Systems
Lecture 31: Some Applications of Eigenvectors: Markov Chains and Chemical Reaction Systems Winfried Just Department of Mathematics, Ohio University April 9 11, 2018 Review: Eigenvectors and left eigenvectors
More informationResolvent Algebras. An alternative approach to canonical quantum systems. Detlev Buchholz
Resolvent Algebras An alternative approach to canonical quantum systems Detlev Buchholz Analytical Aspects of Mathematical Physics ETH Zürich May 29, 2013 1 / 19 2 / 19 Motivation Kinematics of quantum
More informationIntroduction to Machine Learning CMU-10701
Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos & Aarti Singh Contents Markov Chain Monte Carlo Methods Goal & Motivation Sampling Rejection Importance Markov
More informationLecture 02 Linear Algebra Basics
Introduction to Computational Data Analysis CX4240, 2019 Spring Lecture 02 Linear Algebra Basics Chao Zhang College of Computing Georgia Tech These slides are based on slides from Le Song and Andres Mendez-Vazquez.
More informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationComputer Vision Group Prof. Daniel Cremers. 11. Sampling Methods: Markov Chain Monte Carlo
Group Prof. Daniel Cremers 11. Sampling Methods: Markov Chain Monte Carlo Markov Chain Monte Carlo In high-dimensional spaces, rejection sampling and importance sampling are very inefficient An alternative
More informationLeslie Matrix II Formal Demography Stanford Spring Workshop in Formal Demography May 2008
Leslie Matrix II Formal Demography Stanford Spring Workshop in Formal Demography May 2008 James Holland Jones Department of Anthropology Stanford University May 2, 2008 1 Outline 1. Spectral Decomposition
More informationA Single-letter Upper Bound for the Sum Rate of Multiple Access Channels with Correlated Sources
A Single-letter Upper Bound for the Sum Rate of Multiple Access Channels with Correlated Sources Wei Kang Sennur Ulukus Department of Electrical and Computer Engineering University of Maryland, College
More information4 Insect outbreak model
4 Insect outbreak model In this lecture I will put to a good use all the mathematical machinery we discussed so far. Consider an insect population, which is subject to predation by birds. It is a very
More informationRNA From Mathematical Models to Real Molecules
RNA From Mathematical Models to Real Molecules 3. Optimization and Evolution of RNA Molecules Peter Schuster Institut für Theoretische hemie und Molekulare Strukturbiologie der Universität Wien IMPA enoma
More informationSARAH P. OTTO and TROY DAY
A Biologist's Guide to Mathematical Modeling in Ecology and Evolution SARAH P. OTTO and TROY DAY Univsr?.ltats- und Lender bibliolhek Darmstadt Bibliothek Biotogi Princeton University Press Princeton and
More informationPerron-Frobenius theorem for nonnegative multilinear forms and extensions
Perron-Frobenius theorem for nonnegative multilinear forms and extensions Shmuel Friedland Univ. Illinois at Chicago ensors 18 December, 2010, Hong-Kong Overview 1 Perron-Frobenius theorem for irreducible
More informationIntroduction. Genetic Algorithm Theory. Overview of tutorial. The Simple Genetic Algorithm. Jonathan E. Rowe
Introduction Genetic Algorithm Theory Jonathan E. Rowe University of Birmingham, UK GECCO 2012 The theory of genetic algorithms is beginning to come together into a coherent framework. However, there are
More informationCritical Mutation Rate Has an Exponential Dependence on Population Size
Critical Mutation Rate Has an Exponential Dependence on Population Size Alastair Channon 1, Elizabeth Aston 1, Charles Day 1, Roman V. Belavkin 2 and Christopher G. Knight 3 1 Research Institute for the
More informationDeterminants. Samy Tindel. Purdue University. Differential equations and linear algebra - MA 262
Determinants Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. Determinants Differential equations
More informationWeierstraß-Institut für Angewandte Analysis und Stochastik
Weierstraß-Institut fr Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISS 0946-8633 An asymptotic maximum principle for essentially linear evolution models Ellen Baake, Michael
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More information6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities
6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities 1 Outline Outline Dynamical systems. Linear and Non-linear. Convergence. Linear algebra and Lyapunov functions. Markov
More informationON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES
Volume 10 (2009), Issue 4, Article 91, 5 pp. ON THE HÖLDER CONTINUITY O MATRIX UNCTIONS OR NORMAL MATRICES THOMAS P. WIHLER MATHEMATICS INSTITUTE UNIVERSITY O BERN SIDLERSTRASSE 5, CH-3012 BERN SWITZERLAND.
More informationMath 113 Homework 5. Bowei Liu, Chao Li. Fall 2013
Math 113 Homework 5 Bowei Liu, Chao Li Fall 2013 This homework is due Thursday November 7th at the start of class. Remember to write clearly, and justify your solutions. Please make sure to put your name
More informationarxiv:quant-ph/ v1 29 Mar 2003
Finite-Dimensional PT -Symmetric Hamiltonians arxiv:quant-ph/0303174v1 29 Mar 2003 Carl M. Bender, Peter N. Meisinger, and Qinghai Wang Department of Physics, Washington University, St. Louis, MO 63130,
More informationPerron-Frobenius theorem
Perron-Frobenius theorem V.S. Sunder Institute of Mathematical Sciences sunder@imsc.res.in Unity in Mathematics Lecture Chennai Math Institute December 18, 2009 Overview The aim of the talk is to describe
More informationConvex Optimization CMU-10725
Convex Optimization CMU-10725 Simulated Annealing Barnabás Póczos & Ryan Tibshirani Andrey Markov Markov Chains 2 Markov Chains Markov chain: Homogen Markov chain: 3 Markov Chains Assume that the state
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationWhat is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A =
STUDENT S COMPANIONS IN BASIC MATH: THE ELEVENTH Matrix Reloaded by Block Buster Presumably you know the first part of matrix story, including its basic operations (addition and multiplication) and row
More informationDiscrete Math Class 19 ( )
Discrete Math 37110 - Class 19 (2016-12-01) Instructor: László Babai Notes taken by Jacob Burroughs Partially revised by instructor Review: Tuesday, December 6, 3:30-5:20 pm; Ry-276 Final: Thursday, December
More informationThermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model
Thermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model Nikolaj Veniaminov (in collaboration with Frédéric Klopp) CEREMADE, University of Paris IX Dauphine
More informationELEC633: Graphical Models
ELEC633: Graphical Models Tahira isa Saleem Scribe from 7 October 2008 References: Casella and George Exploring the Gibbs sampler (1992) Chib and Greenberg Understanding the Metropolis-Hastings algorithm
More informationSystem theory and system identification of compartmental systems Hof, Jacoba Marchiena van den
University of Groningen System theory and system identification of compartmental systems Hof, Jacoba Marchiena van den IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF)
More informationLecture 11: Finish Gaussian elimination and applications; intro to eigenvalues and eigenvectors (1)
Lecture 11: Finish Gaussian elimination and applications; intro to eigenvalues and eigenvectors (1) Travis Schedler Tue, Oct 18, 2011 (version: Tue, Oct 18, 6:00 PM) Goals (2) Solving systems of equations
More informationUnderstanding MCMC. Marcel Lüthi, University of Basel. Slides based on presentation by Sandro Schönborn
Understanding MCMC Marcel Lüthi, University of Basel Slides based on presentation by Sandro Schönborn 1 The big picture which satisfies detailed balance condition for p(x) an aperiodic and irreducable
More informationLinear-fractional branching processes with countably many types
Branching processes and and their applications Badajoz, April 11-13, 2012 Serik Sagitov Chalmers University and University of Gothenburg Linear-fractional branching processes with countably many types
More informationMatrices: 2.1 Operations with Matrices
Goals In this chapter and section we study matrix operations: Define matrix addition Define multiplication of matrix by a scalar, to be called scalar multiplication. Define multiplication of two matrices,
More informationEntropy and Ergodic Theory Notes 21: The entropy rate of a stationary process
Entropy and Ergodic Theory Notes 21: The entropy rate of a stationary process 1 Sources with memory In information theory, a stationary stochastic processes pξ n q npz taking values in some finite alphabet
More informationA Tropical Extremal Problem with Nonlinear Objective Function and Linear Inequality Constraints
A Tropical Extremal Problem with Nonlinear Objective Function and Linear Inequality Constraints NIKOLAI KRIVULIN Faculty of Mathematics and Mechanics St. Petersburg State University 28 Universitetsky Ave.,
More information1 Planar rotations. Math Abstract Linear Algebra Fall 2011, section E1 Orthogonal matrices and rotations
Math 46 - Abstract Linear Algebra Fall, section E Orthogonal matrices and rotations Planar rotations Definition: A planar rotation in R n is a linear map R: R n R n such that there is a plane P R n (through
More informationGeometric ρ-mixing property of the interarrival times of a stationary Markovian Arrival Process
Author manuscript, published in "Journal of Applied Probability 50, 2 (2013) 598-601" Geometric ρ-mixing property of the interarrival times of a stationary Markovian Arrival Process L. Hervé and J. Ledoux
More informationA PARAMETRIC MODEL FOR DISCRETE-VALUED TIME SERIES. 1. Introduction
tm Tatra Mt. Math. Publ. 00 (XXXX), 1 10 A PARAMETRIC MODEL FOR DISCRETE-VALUED TIME SERIES Martin Janžura and Lucie Fialová ABSTRACT. A parametric model for statistical analysis of Markov chains type
More informationField Theory of a Reaction- Diffusion Model of Quasispecies Dynamics
Field Theory of a Reaction- Diffusion Model of Quasispecies Dynamics Romualdo Pastor-Satorras Ricard V. Solé SFI WORKING PAPER: 2001-05-024 SFI Working Papers contain accounts of scientific work of the
More informationSynchronization between different motifs
Synchronization between different motifs Li Ying( ) a) and Liu Zeng-Rong( ) b) a) College of Information Technology, Shanghai Ocean University, Shanghai 201306, China b) Institute of Systems Biology, Shanghai
More informationPhylogenetic Algebraic Geometry
Phylogenetic Algebraic Geometry Seth Sullivant North Carolina State University January 4, 2012 Seth Sullivant (NCSU) Phylogenetic Algebraic Geometry January 4, 2012 1 / 28 Phylogenetics Problem Given a
More informationOn the Spread of Epidemics in a Closed Heterogeneous Population
On the Spread of Epidemics in a Closed Heterogeneous Population Artem Novozhilov Applied Mathematics 1 Moscow State University of Railway Engineering (MIIT) the 3d Workshop on Mathematical Models and Numerical
More informationMarkov Chains, Random Walks on Graphs, and the Laplacian
Markov Chains, Random Walks on Graphs, and the Laplacian CMPSCI 791BB: Advanced ML Sridhar Mahadevan Random Walks! There is significant interest in the problem of random walks! Markov chain analysis! Computer
More informationFerromagnetic Ordering of Energy Levels for XXZ Spin Chains
Ferromagnetic Ordering of Energy Levels for XXZ Spin Chains Bruno Nachtergaele and Wolfgang L. Spitzer Department of Mathematics University of California, Davis Davis, CA 95616-8633, USA bxn@math.ucdavis.edu
More information