Linear Algebra of the Quasispecies Model
|
|
- Alicia McGee
- 5 years ago
- Views:
Transcription
1 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 Novozhilov (NDSU) August 6, / 1
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) August 6, / 1
3 Quasispecies theory: Manfred Eigen, born 1927 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: Artem Novozhilov (NDSU) August 6, / 1
4 Mathematical side of the story: Ref: M. Eigen, J. McCaskill, P. Schuster, J Phys Chem, 92(24), 1988: Artem Novozhilov (NDSU) August 6, / 1
5 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) August 6, / 1
6 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) August 6, / 1
7 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) August 6, / 1
8 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 = N m i p i. 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. i=0 Artem Novozhilov (NDSU) August 6, / 1
9 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) August 6, / 1
10 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) August 6, / 1
11 Known results: The error threshold Consider the single peaked fitness landscape M = diag(m 0,m 1,...,m 1 ), m 0 > m 1. Figure : Error threshold in the quasispecies model with the single peaked fitness landscape (m = (m 0,m 1,...,m 1), m 0 > m 1). The parameters are N = 30, m 0 = 20, m 1 = 1. (a) The mean population fitness m(µ) versus the mutation rate; (b) the stationary quasispecies distribution versus the mutation rate Artem Novozhilov (NDSU) August 6, / 1
12 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) August 6, / 1
13 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) August 6, / 1
14 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. Recall that N N N M = diag(m 0,...,m N ), Q N = 0 N 1 N N N Artem Novozhilov (NDSU) August 6, / 1
15 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) August 6, / 1
16 The key idea of the proof: Let Consider the linear operator N N N Q N = 0 N 1 N N N Q N : P(t) (1 t 2 )P (t) N(1 t)p(t), acting on the (N +1)-dimensional vector space R N [t] of the polynomials of degree less or equal N. Direct calculations show that matrix Q N is the matrix of Q N in the standard basis {1,t,...,t N } of R N [t]. Artem Novozhilov (NDSU) August 6, / 1
17 Behavior for µ : Let x = C 1 p, where C = (c ij ) is the matrix composed of the eigenvectors of Q, N ( )( ) k N k c ik = ( 1) j, j i j j=0 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) August 6, / 1
18 Parametric solution: Consider fitness landscape such that m j > 0 for some j, and all the rest m i = 0 for i j. Then a parametric solution to the eigenvalue problem can be found in the form Here x = C 1 p. (M +µq)p = mp x k (s) = 2 N c kj 1+ks, k = 0,...,N, p i (s) = 2 N N k=0 c ik c kj 1+ks, i = 0,...,N, m(s) = m j p j (s), µ = s 2 m(s). Artem Novozhilov (NDSU) August 6, / 1
19 Example 1: Let N = 2A and m = (0,...,0,N,0,...,0), where N is exactly on the 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 µ. Figure : Comparison of the limit distribution (black) with the numerical computations (gray, dashed). (a) The sequence length is N = 100, (b) N = 200 Artem Novozhilov (NDSU) August 6, / 1
20 Example 1: Let N = 2A and m = (0,...,0,N,0,...,0), where N is exactly on the A-th position. Consider a scaled fitness landscape Nr = m and r(µ) = m(µ)/n. Denoting r = µ 2 +1 µ we can prove that ( ) k lim p 1 r A±k(µ) = r, k = 0,...,A. N 1+r Figure : Comparison of the limit distribution (black) with the numerical computations (gray, dashed) for N = 100 (a) and N = 200 (b) Artem Novozhilov (NDSU) August 6, / 1
21 Example 2: Single peaked landscape Consider the single peaked landscape, defined as m = (N,0,...,0) = Nr = N(1,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 it can be proved 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 Ref: Saakian, D., and C.-K. Hu. Physical Review E 69.4 (2004) Ref: Bratus, A., Novozhilov, A., and Y. Semenov, Mathematical Biosciences, (2014), submitted Artem Novozhilov (NDSU) August 6, / 1
22 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 Figure : Comparison of numerical solutions of the quasispecies model (gray, dashed) with the single peaked fitness landscape for N = 100 with the geometric limit distribution (black). (a) The expressions for the leading eigenvalue; (b) the quasispecies distribution Artem Novozhilov (NDSU) August 6, / 1
23 Eigenvalue problem for probability generating functions (pgf): Recall that we solve the eigenvalue problem (M +µq)p = mp, and matrix Q is the representation of the linear differential operator Q N : P(t) (1 t 2 )P (t) N(1 t)p(t), acting on the space of polynomials, in the standard basis {1,t,...,t N }. Let P(t) be the pgf of p, i.e., P(t) = N k=0 p kt k. Then the original problem can be replaced with the eigenvalue problem Nr P(t)+Q N P(t) = NrP(t), where r = m/n and r P(t) = N k=0 r kp k t k. Artem Novozhilov (NDSU) August 6, / 1
24 Eigenvalue problem for pgf: Hence we have the problem to find a pgf P(t) satisfying where Nr P(t)+µ(1 t 2 )P (t) µn(1 t)p(t) = NrP(t), r = r(µ) = r P(1). This problem, due to Perron Frobenius theorem, has a unique solution P(t). After dividing by N and assuming that we obtain the limit problem µ N (1 t2 )P (t) 0, N, µ(1 t)p(t)+r P(t) = r P(t). Artem Novozhilov (NDSU) August 6, / 1
25 Example: The single peaked landscape: Start with Let Then Therefore, µ(1 t)p(t)+r P(t) = r P(t). r = (1,0,...,0,...). r P(t) = P(0). µ(1 t)p(t)+p(0) = r P(t). Plugging t = 0, we find r = 1 µ. Which implies that is the pgf of the geometric distribution. P(t) = 1 µ 1 µt, Artem Novozhilov (NDSU) August 6, / 1
26 Acknowledgements: Startup grant from Department of Mathematics, NDSU ND EPSCoR and NSF grant # EPS Artem Novozhilov (NDSU) August 6, / 1
27 Thank you for your attention! site: novozhil/ References: Bratus, A. S., Novozhilov, A. S., & Semenov, Y. S. (2014). Linear algebra of the permutation invariant Crow Kimura model of prebiotic evolution. Submitted to Math Biosciences. Semenov, Y. S., Bratus, A. S., & Novozhilov, A. S. (2014). On the behavior of the leading eigenvalue of the Eigen evolutionary matrices. arxiv: Semenov, Y. S., & Novozhilov, A. S. (2014) Exact solutions for the selection-mutation equilibrium in the Crow-Kimura evolutionary model. arxiv. Artem Novozhilov (NDSU) August 6, / 1
Linear Algebra of Eigen s Quasispecies Model
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,
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 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 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 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 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 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 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 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 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 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 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 information22m:033 Notes: 7.1 Diagonalization of Symmetric Matrices
m:33 Notes: 7. Diagonalization of Symmetric Matrices Dennis Roseman University of Iowa Iowa City, IA http://www.math.uiowa.edu/ roseman May 3, Symmetric matrices Definition. A symmetric matrix is a matrix
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 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 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 informationMa/CS 6b Class 20: Spectral Graph Theory
Ma/CS 6b Class 20: Spectral Graph Theory By Adam Sheffer Recall: Parity of a Permutation S n the set of permutations of 1,2,, n. A permutation σ S n is even if it can be written as a composition of an
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 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 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 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 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 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 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 informationMa/CS 6b Class 20: Spectral Graph Theory
Ma/CS 6b Class 20: Spectral Graph Theory By Adam Sheffer Eigenvalues and Eigenvectors A an n n matrix of real numbers. The eigenvalues of A are the numbers λ such that Ax = λx for some nonzero vector x
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 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 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 information18.600: Lecture 32 Markov Chains
18.600: Lecture 32 Markov Chains Scott Sheffield MIT Outline Markov chains Examples Ergodicity and stationarity Outline Markov chains Examples Ergodicity and stationarity Markov chains Consider a sequence
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 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 informationEigenvalues of the Adjacency Tensor on Products of Hypergraphs
Int. J. Contemp. Math. Sciences, Vol. 8, 201, no. 4, 151-158 HIKARI Ltd, www.m-hikari.com Eigenvalues of the Adjacency Tensor on Products of Hypergraphs Kelly J. Pearson Department of Mathematics and Statistics
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 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 informationDiagonalization of Matrix
of Matrix King Saud University August 29, 2018 of Matrix Table of contents 1 2 of Matrix Definition If A M n (R) and λ R. We say that λ is an eigenvalue of the matrix A if there is X R n \ {0} such that
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 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 informationLinear Algebra in Actuarial Science: Slides to the lecture
Linear Algebra in Actuarial Science: Slides to the lecture Fall Semester 2010/2011 Linear Algebra is a Tool-Box Linear Equation Systems Discretization of differential equations: solving linear equations
More informationMath113: Linear Algebra. Beifang Chen
Math3: Linear Algebra Beifang Chen Spring 26 Contents Systems of Linear Equations 3 Systems of Linear Equations 3 Linear Systems 3 2 Geometric Interpretation 3 3 Matrices of Linear Systems 4 4 Elementary
More information1 - Systems of Linear Equations
1 - Systems of Linear Equations 1.1 Introduction to Systems of Linear Equations Almost every problem in linear algebra will involve solving a system of equations. ü LINEAR EQUATIONS IN n VARIABLES We are
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 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 informationMath Assignment 11
Math 2280 - Assignment 11 Dylan Zwick Fall 2013 Section 8.1-2, 8, 13, 21, 25 Section 8.2-1, 7, 14, 17, 32 Section 8.3-1, 8, 15, 18, 24 1 Section 8.1 - Introduction and Review of Power Series 8.1.2 - Find
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 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 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 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 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 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 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 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 informationLinear Algebra Methods for Data Mining
Linear Algebra Methods for Data Mining Saara Hyvönen, Saara.Hyvonen@cs.helsinki.fi Spring 2007 1. Basic Linear Algebra Linear Algebra Methods for Data Mining, Spring 2007, University of Helsinki Example
More information8 Autonomous first order ODE
8 Autonomous first order ODE There are different ways to approach differential equations. Prior to this lecture we mostly dealt with analytical methods, i.e., with methods that require a formula as a final
More informationPREPRINT 2007:25. Multitype Galton-Watson processes escaping extinction SERIK SAGITOV MARIA CONCEIÇÃO SERRA
PREPRINT 2007:25 Multitype Galton-Watson processes escaping extinction SERIK SAGITOV MARIA CONCEIÇÃO SERRA Department of Mathematical Sciences Division of Mathematical Statistics CHALMERS UNIVERSITY OF
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 information(VI.C) Rational Canonical Form
(VI.C) Rational Canonical Form Let s agree to call a transformation T : F n F n semisimple if there is a basis B = { v,..., v n } such that T v = λ v, T v 2 = λ 2 v 2,..., T v n = λ n for some scalars
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 informationHigh-dimensional two-sample tests under strongly spiked eigenvalue models
1 High-dimensional two-sample tests under strongly spiked eigenvalue models Makoto Aoshima and Kazuyoshi Yata University of Tsukuba Abstract: We consider a new two-sample test for high-dimensional data
More informationSPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint
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 informationDetailed Proof of The PerronFrobenius Theorem
Detailed Proof of The PerronFrobenius Theorem Arseny M Shur Ural Federal University October 30, 2016 1 Introduction This famous theorem has numerous applications, but to apply it you should understand
More informationThe Perron Frobenius theorem and the Hilbert metric
The Perron Frobenius theorem and the Hilbert metric Vaughn Climenhaga April 7, 03 In the last post, we introduced basic properties of convex cones and the Hilbert metric. In this post, we loo at how these
More information18.175: Lecture 30 Markov chains
18.175: Lecture 30 Markov chains Scott Sheffield MIT Outline Review what you know about finite state Markov chains Finite state ergodicity and stationarity More general setup Outline Review what you know
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 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 informationSubadditivity of entropy for stochastic matrices
Subadditivity of entropy for stochastic matrices Wojciech S lomczyński 1 PREPRIT 2001/16 Abstract The entropy of irreducible stochastic matrix measures its mixing properties. We show that this quantity
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 informationOn rate of convergence in distribution of asymptotically normal statistics based on samples of random size
Annales Mathematicae et Informaticae 39 212 pp. 17 28 Proceedings of the Conference on Stochastic Models and their Applications Faculty of Informatics, University of Debrecen, Debrecen, Hungary, August
More informationThe number of simple modules of a cellular algebra
Science in China Ser. A Mathematics 2005 Vol.?? No. X: XXX XXX 1 The number of simple modules of a cellular algebra LI Weixia ( ) & XI Changchang ( ) School of Mathematical Sciences, Beijing Normal University,
More informationEdexcel GCE A Level Maths Further Maths 3 Matrices.
Edexcel GCE A Level Maths Further Maths 3 Matrices. Edited by: K V Kumaran kumarmathsweebly.com kumarmathsweebly.com 2 kumarmathsweebly.com 3 kumarmathsweebly.com 4 kumarmathsweebly.com 5 kumarmathsweebly.com
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 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 informationStochastic Processes
Stochastic Processes 8.445 MIT, fall 20 Mid Term Exam Solutions October 27, 20 Your Name: Alberto De Sole Exercise Max Grade Grade 5 5 2 5 5 3 5 5 4 5 5 5 5 5 6 5 5 Total 30 30 Problem :. True / False
More informationMean-field dual of cooperative reproduction
The mean-field dual of systems with cooperative reproduction joint with Tibor Mach (Prague) A. Sturm (Göttingen) Friday, July 6th, 2018 Poisson construction of Markov processes Let (X t ) t 0 be a continuous-time
More informationViral evolution model with several time scales
Viral evolution model with several time scales AA Archibasov Samara National Research University, 4 Moskovskoe Shosse, 4486, Samara, Russia Abstract In this paper a viral evolution model with specific
More informationCompetitive Exclusion in a Discrete-time, Size-structured Chemostat Model
Competitive Exclusion in a Discrete-time, Size-structured Chemostat Model Hal L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 1804, USA E-mail: halsmith@asu.edu Xiao-Qiang Zhao
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 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 informationChapter 1: Systems of Linear Equations
Chapter : Systems of Linear Equations February, 9 Systems of linear equations Linear systems Lecture A linear equation in variables x, x,, x n is an equation of the form a x + a x + + a n x n = b, where
More informationLecture 32: Infinite-dimensional/Functionvalued. Functions and Random Regressions. Bruce Walsh lecture notes Synbreed course version 11 July 2013
Lecture 32: Infinite-dimensional/Functionvalued Traits: Covariance Functions and Random Regressions Bruce Walsh lecture notes Synbreed course version 11 July 2013 1 Longitudinal traits Many classic quantitative
More informationSemi-Relativistic Reflection and Transmission Coefficients for Two Spinless Particles Separated by a Rectangular-Shaped Potential Barrier
Commun. Theor. Phys. 66 (2016) 389 395 Vol. 66, No. 4, October 1, 2016 Semi-Relativistic Reflection and Transmission Coefficients for Two Spinless Particles Separated by a Rectangular-Shaped Potential
More informationLINEAR MAPS PRESERVING NUMERICAL RADIUS OF TENSOR PRODUCTS OF MATRICES
LINEAR MAPS PRESERVING NUMERICAL RADIUS OF TENSOR PRODUCTS OF MATRICES AJDA FOŠNER, ZEJUN HUANG, CHI-KWONG LI, AND NUNG-SING SZE Abstract. Let m, n 2 be positive integers. Denote by M m the set of m m
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 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 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 informationRegularized Estimation of High Dimensional Covariance Matrices. Peter Bickel. January, 2008
Regularized Estimation of High Dimensional Covariance Matrices Peter Bickel Cambridge January, 2008 With Thanks to E. Levina (Joint collaboration, slides) I. M. Johnstone (Slides) Choongsoon Bae (Slides)
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 informationSection 1.7: Properties of the Leslie Matrix
Section 1.7: Properties of the Leslie Matrix Definition: A matrix A whose entries are nonnegative (positive) is called a nonnegative (positive) matrix, denoted as A 0 (A > 0). Definition: A square m m
More informationJUST THE MATHS SLIDES NUMBER 9.6. MATRICES 6 (Eigenvalues and eigenvectors) A.J.Hobson
JUST THE MATHS SLIDES NUMBER 96 MATRICES 6 (Eigenvalues and eigenvectors) by AJHobson 96 The statement of the problem 962 The solution of the problem UNIT 96 - MATRICES 6 EIGENVALUES AND EIGENVECTORS 96
More informationGeneralized hashing and applications to digital fingerprinting
Generalized hashing and applications to digital fingerprinting Noga Alon, Gérard Cohen, Michael Krivelevich and Simon Litsyn Abstract Let C be a code of length n over an alphabet of q letters. An n-word
More informationDecomposition of supercritical branching processes with countably many types
Pomorie 26 June 2012 Serik Sagitov and Altynay Shaimerdenova 20 slides Chalmers University and University of Gothenburg Al-Farabi Kazakh National University Decomposition of supercritical branching processes
More informationChapter 2 Heisenberg s Matrix Mechanics
Chapter 2 Heisenberg s Matrix Mechanics Abstract The quantum selection rule and its generalizations are capable of predicting energies of the stationary orbits; however they should be obtained in a more
More information