A MARKOV PROCESS OF GENE FREQUENCY CHANGE IN A
|
|
- Ursula Holmes
- 6 years ago
- Views:
Transcription
1 A MARKOV PROCESS OF GENE FREQUENCY CHANGE IN A GEOGRAPHICALLY STRUCTURED POPULATION TAKE0 MARUYAMA National Instituie of Genetics, Mishima, Japan and Cenier for Demographic and Population Genetics, University of Texas at Houston, Texas Manuscript received February 9, 1972 Revised copy received September 26, 1973 ABSTRACT A Markov process (chain) of gene frequency change is derived for a geographically-structured model of a population. The population consists of colonies which are connected by migration. Selection operates in each colony independently. It is shown that there exists a stochastic clock that transforms the originally complicated process of gene frequency change to a random walk which is independent of the geographical structure of the population. The time parameter is a local random time that is dependent on the sample path. In fact, if the alleles are selectively neutral, the time parameter is exactly equal to the sum of the average local genetic variation appearing in the population, and otherwise they are approximately equal. The Kolmogorov forward and backward equations of the process are obtained. As a limit of large population size, a diffusion process is derived. The transition probabilities of the Markov chain and of the diffusion process are obtained explicitly. Certain quantities of biological interest are shown to be independent of the population structure. The quantities are the fixation probability of a mutant, the sum of the average local genetic variation and the variation summed over the generations in which the gene frequency in the whole population assumes a specified value. HE stochastic model of gene frequency change put forward by FISHER (1930) WRIGHT (1931 and later) has become increasingly important in population genetics. For reviews, see MORAN (1962); KIMURA (1964); EWENS (1969); WRIGHT (1969); CROW and KIMURA (1970); CAVALLI-SFORZA and BODMER (1971); KIMURA and OHTA (1971). However, most of the theory has dealt with random mating populations, while hardly any natural population meets this assumption. When we take into consideration the effect of geographical structure which prevents random mating, the mathematics becomes very difficult. WRIGHT ( 1931, 1943, 1951 ) pioneered in the analysis of structured populations. Later MALBCOT (1948, 1967) developed a new model and contributed greatly in this area. KIMURA and WEISS (1964), and BODMER and CAVALLI- SFORZA (1968) have developed still new models. Although different models and different methods are used, the essential point in these works is to investigate the genetic similarity between different localities. WRIGHT (1931,1937,1938) and KIMURA (1955) have dealt with the stochastic Genetm 76: February, 1974
2 368 T. MARUYAMA process of gece frequency change and of gene frequency distribution in a panmictic population of finite size. WRIGHT obtained the stationary distribution when mutation in both directions is assumed, and KIMURA obtained the time-dependent distribution of gene frequency change under the assumption of no mutation. EWENS (1963) obtained the mean of the sojourn time, that is, the duration of time during which the population assumes a particular gene frequency. Although the geographical structure of the population may not necessarily alter significantly the theory based on the assumption of random mating, it is certainly worthwhile to investigate the situation. The purposes of this paper is to develop a Markov process of gene frequency change in a special case and to investigate its properties. Throughout this paper, we consider one locus where two alleles, A and a, are segregating. The organism is haploid. The fitnesses of A and a are 1 + s and 1 respectively, and the fitnesses are constant in time and in the geographical structure. No mutation or migration from the outside occurs during the time under consideration. 2. DERIVATION OF THE MARKOV PROCESS The population consists of colonies. Migration occurs and connects colonies genetically. Each individual has a negative exponential life time distribution with unit mean. When an individual dies, it is immediately replaced by a new individual born in the same colony. This assumption implies that the total population size is constant. We denote the size of the total population by N. When an individual dies, each preexisting individual in the colony has a chance of being chosen to reproduce which is proportional to that individual's fitness. This is a geographically-structured version of the MORAN (1958) model. Let Nt,k be the size of colony IC at time t and Mt be the number of A genes in the entire population at time t. Then the random variable {Mt} is a stochastic process whose state space is [0, 1,2,..., NI. This process with time parameter t can be very complicated, but we will introduce a new time measure which is sample-path-dependent local time based on an additive functional and which makes the process into a Markov process, actually a random walk. First, note that, in any short time interval, At, two or more death-birth events occur with probability of order (At) z. We may therefore assume that the random variable Mt changes by at most one at any moment. Now consider a particular moment when a birth-death occurs, and ask what is the probability, given Mt = i immediately before this event, that the value of Mt changes from i to if+ 1 through this event? Let ik be the number of A genes in colony k at this moment, i.e., i = z ik. k The probability that the death-birth occurs in colony k whose size is Nt,k, is Nt,k/N, and given this, the probability, qi, +1, that one a in the colony dies and ' I it is replaced by one A, 'is The probability, qik,i,-;, that one A dies and replaced by one a, is
3 STRUCTURED POPULATION MARKOV PROCESS 369 The probability, qib,{&, that no change in the number of A occurs, is Therefore the probability, H (t), that the value of Mt changes through the deathbirth event is H(t) can be defined for all t. We assume that migration takes place instantaneously and it is independent of death-birth events. The value of H(t) changes discontinuously as migration or a death-birth event occurs, and otherwise H(t) remains constant. For example, assume that the population consists of two colonies, say 1 and 2, and that each of the two colonies has two A genes and two a genes and s = 0. Under this situation, H(t) = (2/8) (4 x 2 ~2 X 2/ X 2 X 2/16) = x. Now suppose that one A gene moves from colony 1 to colony 2, then H(t) = (2/8) (3 X 1 x 2/9 + 5 x 3 x 2/25) = 7/15. It is important that the effects of migration and of the population structure are incorporated in H(t). Note that the quantity in the square bracket is the probability that two randomly chosen gametes from a colony with replacement are different types. We call this quantity the local genetic variation. Therefore, H(t) is equal to the average local genetic variation if the alleles are selectively neutral, and otherwise they are approximately equal. We now ask, if Mt is changed, what is the probability that Mt is changed from i to i + l? This conditional probability is equal to Through simple algebra, this becomes Ifs -- 2+s. Similarly, the conditional probability that Mt is changed from i to i - 1 is 1 2+s We should note two properties here: (1) these conditional probabilities are independent of state i; (2) they are independent of the geographical structure of the population. Property (1) was shown by MORAN (1959, 1962) "for a panmictic case, and it was used to obtain the ultimate fixation probability of a mutant gene.
4 3 70 T. MARUYAMA Consider a short time interval At. Then the probability of occurrence of one death-birth event is NAt + o(4t). Let Rt be the probability that Mt is unchanged for time interval t. then The first term on the right side of the above equation is the probability that Mt remains unchanged until time t and no death-birth occurs in (t, t + At). The second term is the probability that Mt remains until t and one death-birth occurs in (t, t '+ At) but Mt is unchanged. From this equation we have and therefore drt --NH(t), Rtdt- Et t e -xj H(t)dt (2-2) 0 Consider the stochastic process as a collection of sample path {U}, and let where w indicates a particular sample path and H(w, t) is the same quantity as H([) of (2-1) associated with sample path.u. We assume that no part of the population is completely separated from the other. Consider this T(0,t) = T as a new time measure. The T is a local time based on the functional H(t) of the original process. Because the T depends on the history of the sample path, it is often called a stochastic clock. For a full account of this kind of time substitution, the reader is referred to ITO and MCKFAN (1965, chapters 2,5, and 6). From here on, we shall use the T as the time parameter and construct a Markov process. Let Mr be the number of A genes in the entire population at time T. The M7 is different from the Mt and it is not to be confused with the M,. Let RI be the probability that M7 is unchanged for time interval 7. Then, from (2-2) together with (2-3), we have R 7 = e-n7. (2-4) This process Mr is characterized as follows: (1) The state space is [0, 1,2,... A'] ; (2) The sojourn time at each state at each visit is exponentially distributed with mean 1/N, provided 0 < i < N. This follows from (2-4): (3) The conditional probability that Mr is changed from i to i + i is (1 + s)/(2 + s), and the conditional probability that M is changed from i to i - 1 is 1/(2 4- s), provided 0 < i < N. These conditional probabilities are independent of the state and of the geographical structure of the population. Let be the probability, given that Mr is i at time 0, that it is j at time T. From the properties (1)-(3) above, we have
5 STRUCTURED POPULATION MARKOV PROCESS 371 q..= z q q,. for O < ~ < T, 7,%,1 r&&k &h,~ and M7 is a Markov process. Then we have, for small AT, q7-47,irj = e-n*,q7,i,j + (1- e-n4t) +o(ar), O<i,j< N chqr,i-l,j + W,,i+l,j 1 where A = (1 4- s)/(2 4- s) and p = 1/(2 + s). From this equation, we derive dqr,i,j --- -N[I-q7,i,j+Xq7,i-l,j++q,,i+l,j I, o<i,i<n (2-5) dr the Kolmogorov backward equation. Using the additive functional (2-3), this equation can be also derived from DYNKIN S formula, (cf. DYNKIN 1965 vol. 1, p. 133). We can derive the Kolmogorov forward equation similarly: These differential equations characterize the Markov process, and all the information concerning the process can be obtained from these equations. It is often convenient to express the system of equation in matrix notation. We define the following matrices: Q, = [q,,i,jl (N-WN-~) 0 < i, i < N ; A= A-IpO A O h - l (N-l)S (N-1) where x = (1*+ s)/(2:+ s) and p = 1/(2 + s). Then the system of differential equations (2-5) becomes and (2-6) becomes where AS is the transpose of A. dq NAQ, dt (2-7) -- dq --NASQ (2-8) dt
6 3 752 T. MARUYAMA We shall next derive a diffusion process as a limit for large N. The backward equations (2-5) can be rewritten as dqr,i,j ~ - - ds 2 N[q7,i-i,j - 2q7,i,j + qr,i+i,j 1 Let T be a new time measure such that NT = T, and Ns = S where we assume S is a constant. Let qn (T, i, i -) i = q... Then (2-9) can be written as N N 7,573 Therefore as N + a, this difference equation converges to - --mi, %m(t,x,y) at 2 8x2 2 2X - 1 8% (T x Y ) + s a%j(t,%y). (2-10) The derivation of this equation is not mathematically rigorous, but rigorous derivation is possible (cf. KAC 1947). The Markov process governed by (2-10) is a Brownian motion. 3. SOME SOLUTIONS Based on the stochastic clock, whose time-rate depends on the total amount of the average local genetic variation which has appeared in the sample path, we have derived a Markov process of the gene frequency change. It is governed by differential equations (2-5),(2-6): (2-7),(2-8),or(2-10). Here we shall deal with (2-7) and (2-lo), and investigate biologically interesting questions. These equations together with definitions (2-1) and (2-3) tell us that, if we measure the time by 7 the stochastic process becomes a Markov process that is independent of the geographical structure of the population and furthermore, that it is a random walk (or approximately a Brownian motion) on the finite interval [ 1,2,..., N - I], (or (0,l)) with exit boundaries at both ends. A standard procedure in dealing with equation (2-7) is to obtain a complete spectral analysis of matrix A. The right eigenvectors of A are the column vectors - 37rk P N-1, (N-l)7rk * (A) sin ), k=1,2,..., N-I. E.l N P (3-1)
7 STRUCTURED POPULATION MARKOV PROCESS 373 where a* indicates the transpose of any row vector a. It can be easily shown that where Now let - 2 dk = (($) sin, h N-l/N (+) sin N Aek = hkek xk h,=-112~rpcos-. (3-2) N rrk (YP sin (N - 1)xk. 2xk, (5) ), k=l,2,..., N-I. These row vectors are the left eigenvectors of A, i.e., dha = Xkdk 3xk 3/2 sin,..., where the eigenvalues are the same as (3-2). Note that {ek}, or {d,}, themselves are not orthogonal systems, but {ek} and {d,} are biorthogonal, i.e., where Skk = 1 and Skl = 0 if k f 1. With these eigenvectors and eigenvalues, we can easily obtain the fundamental solution of the process of equation (2-7), that is the transition probability: where 0 indicates the direct product of two vectors. This can be written as The fundamental solutions of (2-10) are m qm(t, 2, y) = 2e-S(~9)/2 2 e-%((k29+wz))t/4 sin kxx sin kxy. k=1 The Markov process governed by qa (T, 5, y) is independent of the geographical structure and it is a Brownian motion. The ultimate fixation probability, ui, of allele A in the whole population is given by the solution of equation Au=O where U = (ul, u2,..., un-l) *, in which * indicates the transpose:
8 3 74 T. MARUYAMA where i is the initial number of A in the whole population. This is the same as MORAN (1959) in which the formula was derived for a panmictic population. Therefore the fixation probability is independent of the geographical structure proposed in this paper. If we replace si and sn of (3-3) by Sx and S respectively, the right side in (3-3) is the solution of the differential equation with boundary conditions f (0) = 0 and f (1) = 1 ( KIMURA 1962). The independence of the fixation probability here was derived under the assumption of constant size of the total population and local random mating, but importantly the structure of the population was not assumed to be fixed. Even if population structure depends on the genetic constitution, as would occur if genetically similar individuals tend to live closer to one another, or the structure depends on the gene frequency, the fixation probability is unaltered. This independence of the fixation probability was suggested by MARUYAMA (1970), based on an approximation. However it was not clear whether the independence is exact or an approximation. The first exit time of the process from [I, 2,..., N - 11 is the solution of equation Affl=O where 1 = (1,1,..., 1) * and 0 = (0, 0, -.., O)*. The i-th entry, fi, of f is the mean duration of time, measured by T, until allele A is fixed in the population or A is lost from the population. (For the first exit time in the case of a panmictic population, the reader is referred I to EWENS 1963.) If we let f (-) fi and let N become large while S = Ns stays N constant, f (5) becomes the solution of differential equation 1 d2f S df l=o: 2 dx2 2 dx The biological meaning of fi is that it is the sum of the average local genetic variation that appears in the population while the two alleles are segregating, provided that initially there are i A genes. This is exact if s = 0 and is approximately so if s# 0, (cf. (2-3)). The second moment, f i2), of the first exit time measured by T is the solution of Af@) + 2f = 0
9 STRUCTURED POPULATION MARKOV PROCESS 3 75 where f is the mean obtained above and f(z) = ( f?), fi2),..., fe!l). The higher moments are similarly obtained: AfW) + nf(n-1) = 0. In the above calculations, we included both kinds of sample paths in which A is established in the whole population and in which A is lost. We can make a distinction between the two kinds, and calculate the quantities for those paths in which A is fixed, (or in which A is lost). Let f i be the mean exit time in the path in which A is fixed. Then g = (ufl, uzfz,. *., un-lj -l) *, is the solution of equation Ag+u=O where U is the vector of the fixation probabilities, U+ Biologically, f+ is the expectation of the sum of the average local genetic variation given that A is established in the population. The solution is uj sin- 1 N-1 N-1 A7. irk. irk sm-, foro<i<n where ui is the fixation probability given in (3-3). The higher moments are the solution of Ag( ) f ng(n-1) = 0 where g(o) = U the fixation probability vector. As N becomes large, while S = Ns remains constant, fl ($) = fi and U (&) = ui satisfy the differential equation Finally we shall obtain the sojourn time at a given state, or in a given range of gene frequency. For 0 < i, j < N, let ai5 be the sojourn time, measured by 7, at state i while A and a are segregating in the population, provided the process starts from state i. Let Q1ij be the sojourn time among the paths in which the fixation of A occurs. Let = [air] (N-l) (N-l) and a1 = [alijui] (N-l) (N-l). Then these matrices are and where U = [64+i]. a = - A-1 Ql = - A- U
10 376 T. MARUYAMA As N becomes large, aij z- (1- e-zs(n-i) ) (1- elzsi) if j>i s( 1 - (3-5) 1.3 == 2i(N - j ) N if j>i a.. = ~- 2j(N - i) 23 if j<i N Formula (3-6) is a well-known result for a panmictic population, (cf. EWENS i i 1963). If we y) -) = aij, this satisfies the following differ- N'N entia1 equation wherewhere 6 (.) is the Dirac delta function. For large N, alii can be obtained as the solution of Approximation formulas (3-4), (3-5), and (3-6) were obtained in MARUYAMA (1972), where differential equation (2-10) and the method described in DYNKIN (1965, vol. 2, pp ) were used. Before concluding this paper, I would like to emphasize that the sojourn time, measured by the number of generations, the time until fixation of an allele, the distribution of the gene frequencies and the rate of decay of genetic variability, depends upon the population structure. I thank the referees for many criticism and suggestions, and for correcting the English. I also thank PROFESSOR T. Umo of the Mathematics Department, University of Tokyo, and PROFESSOR T. G. KURTZ of University of Wisconsin, for helpful conversation connected with the stochastic clock used in the paper. This work was supported in part by U.S. Public Health Service Research Grants GM and GM LITERATURE CITED BODMER, W. F. and L. L. CAVALLI-SFORW, 1968 A migration matrix model for the study of random genetic drift. Genetics 59:
11 ~, STRUCTURED POPULATION MARKOV PROCESS 377 CAVALLI-SFORZA, L. L. and W. F. BODMER, 1971 The genetics of human population. Freeman, San Francisco. CROW, J. F. and M. KIMURA, 1970 An introduction to population genetics theory. Harper and Row, New York. DYNKIN, E. B., 1965 Markov Processes, Vols. 1 and 2. Academic Press, New York. EWENS, W. J., 1963 The diffusion equation and a pseudo-distribution in genetics. J. Roy. Stat. Soc. B25: , 1969 Population Genetics. Methuen, London. FISHER, R. A., 1930 The genetical theory of natural selection. Clarendon Press, Oxford, Rev. Ed., 1958, New York, Dover Press. ITO, K. and H. P. MCKEAN, JR., 1965 Diffusion processes and their sample paths. Springer- Verlag, Berlin-Heidelberg-New York. Gc, M., 1947 Random walk and the theory of Brownian motion. Am. Math. Monthly 54: KIMURA, M., 1955 Solution of a process of random genetic drift with a continuous model. Proc. Natl. Acad. Sci. US. 41: I& On the probability of fixation of mutant genes in a population. Gsnetics 47: , 1964 Diffusion models in populahon genetics. Journal of Applied Probabiility 1 : KIMURA, M. and T. OHTA, 1971 Theoretical aspects of populaiion genetics. Princeton University Press, Princeton, New Jersey. KIMURA, M. and G. H. WEIS, 1964 The stepping stone model of population structure and the decrease of genetic correlation with distance. Genetics 49: MALBCOT, G., 1948 Les Mathkmatiques de l'mr6ditl. Masson et Cie, Paris., 1967 Identical loci and relationship. Proc. Fifth Berkeley Symp. Math. Stat. Prob. 4: MARUYAMA, T., 1970 On the fixation probability of mutant genes in a subdivided population. Genet. Res. 15: , 1972 Some invariant properties of a geographically structured finite population. Genet. Res. 20 : MORAN, P. A. P., 1958 Random processes in genetics. Proc. Camb. Phil. Soc. 54: , 1959 The survival of a mutant gene under selection. Aust. Math. Soc. 1: , 1962 The Statistical Processes of Evolutionary Theory. Clarendon Press, Oxford. WRIGHT, S., 1931 Evolution in mendelian populations. Genetics 16: , 1937 The distribution of gene frequencies in populations. Proc. Natl. Acad. Sci. 23: , 1938 The distribution of gene frequencies under irreversible mutation. Proc. Natl. Acad. Sci. US. 24: , 1943 Isolation by distance. Genetics 28: , 1951 The genetical structure of populations. Annals of Eugenics 15: , 1969 Evolution and the Genetics of Populations. Vol. 2. The Theory of Gene Frequencies. University of Chicago Press, Chicago. Corresponding Editor: T. PROUT
NATURAL SELECTION FOR WITHIN-GENERATION VARIANCE IN OFFSPRING NUMBER JOHN H. GILLESPIE. Manuscript received September 17, 1973 ABSTRACT
NATURAL SELECTION FOR WITHIN-GENERATION VARIANCE IN OFFSPRING NUMBER JOHN H. GILLESPIE Department of Biology, University of Penmyluania, Philadelphia, Pennsyluania 19174 Manuscript received September 17,
More informationThe Wright-Fisher Model and Genetic Drift
The Wright-Fisher Model and Genetic Drift January 22, 2015 1 1 Hardy-Weinberg Equilibrium Our goal is to understand the dynamics of allele and genotype frequencies in an infinite, randomlymating population
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research education use, including for instruction at the authors institution
More informationI of a gene sampled from a randomly mating popdation,
Copyright 0 1987 by the Genetics Society of America Average Number of Nucleotide Differences in a From a Single Subpopulation: A Test for Population Subdivision Curtis Strobeck Department of Zoology, University
More informationHow robust are the predictions of the W-F Model?
How robust are the predictions of the W-F Model? As simplistic as the Wright-Fisher model may be, it accurately describes the behavior of many other models incorporating additional complexity. Many population
More informationAARMS Homework Exercises
1 For the gamma distribution, AARMS Homework Exercises (a) Show that the mgf is M(t) = (1 βt) α for t < 1/β (b) Use the mgf to find the mean and variance of the gamma distribution 2 A well-known inequality
More informationURN MODELS: the Ewens Sampling Lemma
Department of Computer Science Brown University, Providence sorin@cs.brown.edu October 3, 2014 1 2 3 4 Mutation Mutation: typical values for parameters Equilibrium Probability of fixation 5 6 Ewens Sampling
More informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More informationDISTRIBUTION OF NUCLEOTIDE DIFFERENCES BETWEEN TWO RANDOMLY CHOSEN CISTRONS 1N A F'INITE POPULATION'
DISTRIBUTION OF NUCLEOTIDE DIFFERENCES BETWEEN TWO RANDOMLY CHOSEN CISTRONS 1N A F'INITE POPULATION' WEN-HSIUNG LI Center for Demographic and Population Genetics, University of Texas Health Science Center,
More informationLECTURE NOTES: Discrete time Markov Chains (2/24/03; SG)
LECTURE NOTES: Discrete time Markov Chains (2/24/03; SG) A discrete time Markov Chains is a stochastic dynamical system in which the probability of arriving in a particular state at a particular time moment
More informationPopulation Genetics and Evolution III
Population Genetics and Evolution III The Mechanisms of Evolution: Drift São Paulo / January 2019 SMRI (Italy) luca@peliti.org 1 Drift 2 The Population Genetics Triad Sewall Wright Ronald A. Fisher Motoo
More informationERGODICITY OF DIFFUSION AND TEMPORAL UNIFORMITY OF DIFFUSION APPROXIMATION. 1. Introduction
ERGODICITY OF DIFFUSION AND TEMPORAL UNIFORMITY OF DIFFUSION APPROXIMATION J. Appl. Prob. 14, 399-404 (1977) Printed in Israel 0 Applied Probability Trust 1977 M. FRANK NORMAN, University of Pennsylvania
More informationEvolution in a spatial continuum
Evolution in a spatial continuum Drift, draft and structure Alison Etheridge University of Oxford Joint work with Nick Barton (Edinburgh) and Tom Kurtz (Wisconsin) New York, Sept. 2007 p.1 Kingman s Coalescent
More informationEvolutionary dynamics on graphs
Evolutionary dynamics on graphs Laura Hindersin May 4th 2015 Max-Planck-Institut für Evolutionsbiologie, Plön Evolutionary dynamics Main ingredients: Fitness: The ability to survive and reproduce. Selection
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 informationFrequency Spectra and Inference in Population Genetics
Frequency Spectra and Inference in Population Genetics Although coalescent models have come to play a central role in population genetics, there are some situations where genealogies may not lead to efficient
More informationStochastic Demography, Coalescents, and Effective Population Size
Demography Stochastic Demography, Coalescents, and Effective Population Size Steve Krone University of Idaho Department of Mathematics & IBEST Demographic effects bottlenecks, expansion, fluctuating population
More informationPROBABILITY OF FIXATION OF A MUTANT GENE IN A FINITE POPULATION WHEN SELECTIVE ADVANTAGE DECREASES WITH TIME1
PROBABILITY OF FIXATION OF A MUTANT GENE IN A FINITE POPULATION WHEN SELECTIVE ADVANTAGE DECREASES WITH TIME1 MOT00 KIMURA AND TOMOKO OHTA National Institute of Genetics, Mishima, Japan Received December
More informationLIST OF MATHEMATICAL PAPERS
LIST OF MATHEMATICAL PAPERS 1961 1999 [1] K. Sato (1961) Integration of the generalized Kolmogorov-Feller backward equations. J. Fac. Sci. Univ. Tokyo, Sect. I, Vol. 9, 13 27. [2] K. Sato, H. Tanaka (1962)
More informationLectures on Probability and Statistical Models
Lectures on Probability and Statistical Models Phil Pollett Professor of Mathematics The University of Queensland c These materials can be used for any educational purpose provided they are are not altered
More informationGENETICS. On the Fixation Process of a Beneficial Mutation in a Variable Environment
GENETICS Supporting Information http://www.genetics.org/content/suppl/2/6/6/genetics..24297.dc On the Fixation Process of a Beneficial Mutation in a Variable Environment Hildegard Uecker and Joachim Hermisson
More informationSIMILAR MARKOV CHAINS
SIMILAR MARKOV CHAINS by Phil Pollett The University of Queensland MAIN REFERENCES Convergence of Markov transition probabilities and their spectral properties 1. Vere-Jones, D. Geometric ergodicity in
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 informationLINKAGE DISEQUILIBRIUM, SELECTION AND RECOMBINATION AT THREE LOCI
Copyright 0 1984 by the Genetics Society of America LINKAGE DISEQUILIBRIUM, SELECTION AND RECOMBINATION AT THREE LOCI ALAN HASTINGS Defartinent of Matheinntics, University of California, Davis, Calijornia
More informationThe fate of beneficial mutations in a range expansion
The fate of beneficial mutations in a range expansion R. Lehe (Paris) O. Hallatschek (Göttingen) L. Peliti (Naples) February 19, 2015 / Rutgers Outline Introduction The model Results Analysis Range expansion
More informationSo in terms of conditional probability densities, we have by differentiating this relationship:
Modeling with Finite State Markov Chains Tuesday, September 27, 2011 1:54 PM Homework 1 due Friday, September 30 at 2 PM. Office hours on 09/28: Only 11 AM-12 PM (not at 3 PM) Important side remark about
More informationThe Combinatorial Interpretation of Formulas in Coalescent Theory
The Combinatorial Interpretation of Formulas in Coalescent Theory John L. Spouge National Center for Biotechnology Information NLM, NIH, DHHS spouge@ncbi.nlm.nih.gov Bldg. A, Rm. N 0 NCBI, NLM, NIH Bethesda
More informationWright-Fisher Models, Approximations, and Minimum Increments of Evolution
Wright-Fisher Models, Approximations, and Minimum Increments of Evolution William H. Press The University of Texas at Austin January 10, 2011 1 Introduction Wright-Fisher models [1] are idealized models
More informationMarkov Chains Handout for Stat 110
Markov Chains Handout for Stat 0 Prof. Joe Blitzstein (Harvard Statistics Department) Introduction Markov chains were first introduced in 906 by Andrey Markov, with the goal of showing that the Law of
More informationComputational Systems Biology: Biology X
Bud Mishra Room 1002, 715 Broadway, Courant Institute, NYU, New York, USA Human Population Genomics Outline 1 2 Damn the Human Genomes. Small initial populations; genes too distant; pestered with transposons;
More informationSMSTC (2007/08) Probability.
SMSTC (27/8) Probability www.smstc.ac.uk Contents 12 Markov chains in continuous time 12 1 12.1 Markov property and the Kolmogorov equations.................... 12 2 12.1.1 Finite state space.................................
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
More informationA Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit
A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit David Vener Department of Mathematics, MIT May 5, 3 Introduction In 8.366, we discussed the relationship
More informationMathematical Population Genetics. Introduction to the Stochastic Theory. Lecture Notes. Guanajuato, March Warren J Ewens
1 Mathematical Population Genetics Introduction to the Stochastic Theory Lecture Notes Guanajuato, March 2009 Warren J Ewens 2 Preface These notes should, ideally, be read before the Guanajuato meeting
More informationGenetic Variation in Finite Populations
Genetic Variation in Finite Populations The amount of genetic variation found in a population is influenced by two opposing forces: mutation and genetic drift. 1 Mutation tends to increase variation. 2
More informationElementary Applications of Probability Theory
Elementary Applications of Probability Theory With an introduction to stochastic differential equations Second edition Henry C. Tuckwell Senior Research Fellow Stochastic Analysis Group of the Centre for
More informationMUTATIONS that change the normal genetic system of
NOTE Asexuals, Polyploids, Evolutionary Opportunists...: The Population Genetics of Positive but Deteriorating Mutations Bengt O. Bengtsson 1 Department of Biology, Evolutionary Genetics, Lund University,
More informationSurfing genes. On the fate of neutral mutations in a spreading population
Surfing genes On the fate of neutral mutations in a spreading population Oskar Hallatschek David Nelson Harvard University ohallats@physics.harvard.edu Genetic impact of range expansions Population expansions
More informationStochastic processes and
Stochastic processes and Markov chains (part II) Wessel van Wieringen w.n.van.wieringen@vu.nl wieringen@vu nl Department of Epidemiology and Biostatistics, VUmc & Department of Mathematics, VU University
More informationThe Moran Process as a Markov Chain on Leaf-labeled Trees
The Moran Process as a Markov Chain on Leaf-labeled Trees David J. Aldous University of California Department of Statistics 367 Evans Hall # 3860 Berkeley CA 94720-3860 aldous@stat.berkeley.edu http://www.stat.berkeley.edu/users/aldous
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 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 informationLecture 18 : Ewens sampling formula
Lecture 8 : Ewens sampling formula MATH85K - Spring 00 Lecturer: Sebastien Roch References: [Dur08, Chapter.3]. Previous class In the previous lecture, we introduced Kingman s coalescent as a limit of
More informationContrasts for a within-species comparative method
Contrasts for a within-species comparative method Joseph Felsenstein, Department of Genetics, University of Washington, Box 357360, Seattle, Washington 98195-7360, USA email address: joe@genetics.washington.edu
More informationSegregation versus mitotic recombination APPENDIX
APPENDIX Waiting time until the first successful mutation The first time lag, T 1, is the waiting time until the first successful mutant appears, creating an Aa individual within a population composed
More informationA Sampling Theory of Selectively Neutral Alleles in a Subdivided Population
Copyright 0 1988 by the Genetics Society of America A Sampling Theory of Selectively Neutral Alleles in a Subdivided Population Elisabeth R. Tillier and G. Brian Golding Department of Biology, York University,
More informationG(t) := i. G(t) = 1 + e λut (1) u=2
Note: a conjectured compactification of some finite reversible MCs There are two established theories which concern different aspects of the behavior of finite state Markov chains as the size of the state
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 1982 NOTES ON MATRIX METHODS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 198 NOTES ON MATRIX METHODS 1. Matrix Algebra Margenau and Murphy, The Mathematics of Physics and Chemistry, Chapter 10, give almost
More informationreciprocal altruism by kin or group selection can be analyzed by using the same approach (6).
Proc. Nati. Acad. Sci. USA Vol. 81, pp. 6073-6077, October 1984 Evolution Group selection for a polygenic behavioral trait: Estimating the degree of population subdivision (altruism/kin selection/population
More informationProblems for 3505 (2011)
Problems for 505 (2011) 1. In the simplex of genotype distributions x + y + z = 1, for two alleles, the Hardy- Weinberg distributions x = p 2, y = 2pq, z = q 2 (p + q = 1) are characterized by y 2 = 4xz.
More informationThe Kemeny constant of a Markov chain
The Kemeny constant of a Markov chain Peter Doyle Version 1.0 dated 14 September 009 GNU FDL Abstract Given an ergodic finite-state Markov chain, let M iw denote the mean time from i to equilibrium, meaning
More informationStochastic Modelling Unit 1: Markov chain models
Stochastic Modelling Unit 1: Markov chain models Russell Gerrard and Douglas Wright Cass Business School, City University, London June 2004 Contents of Unit 1 1 Stochastic Processes 2 Markov Chains 3 Poisson
More informationStochastic processes and Markov chains (part II)
Stochastic processes and Markov chains (part II) Wessel van Wieringen w.n.van.wieringen@vu.nl Department of Epidemiology and Biostatistics, VUmc & Department of Mathematics, VU University Amsterdam, The
More informationThe mathematical challenge. Evolution in a spatial continuum. The mathematical challenge. Other recruits... The mathematical challenge
The mathematical challenge What is the relative importance of mutation, selection, random drift and population subdivision for standing genetic variation? Evolution in a spatial continuum Al lison Etheridge
More informationTHE genetic consequences of population structure simple migration models with no selection (cited above),
Copyright 2003 by the Genetics Society of America A Diffusion Approximation for Selection and Drift in a Subdivided Population Joshua L. Cherry 1 and John Wakeley Department of Organismic and Evolutionary
More information4452 Mathematical Modeling Lecture 16: Markov Processes
Math Modeling Lecture 16: Markov Processes Page 1 4452 Mathematical Modeling Lecture 16: Markov Processes Introduction A stochastic model is one in which random effects are incorporated into the model.
More informationVARIANCE AND COVARIANCE OF HOMOZYGOSITY IN A STRUCTURED POPULATION
Copyright 0 1983 by the Genetics Society of America VARIANCE AND COVARIANCE OF HOMOZYGOSITY IN A STRUCTURED POPULATION G. B. GOLDING' AND C. STROBECK Deportment of Genetics, University of Alberta, Edmonton,
More informationGenetic hitch-hiking in a subdivided population
Genet. Res., Camb. (1998), 71, pp. 155 160. With 3 figures. Printed in the United Kingdom 1998 Cambridge University Press 155 Genetic hitch-hiking in a subdivided population MONTGOMERY SLATKIN* AND THOMAS
More information1.3 Forward Kolmogorov equation
1.3 Forward Kolmogorov equation Let us again start with the Master equation, for a system where the states can be ordered along a line, such as the previous examples with population size n = 0, 1, 2,.
More informationMathematical Population Genetics II
Mathematical Population Genetics II Lecture Notes Joachim Hermisson June 9, 2018 University of Vienna Mathematics Department Oskar-Morgenstern-Platz 1 1090 Vienna, Austria Copyright (c) 2013/14/15/18 Joachim
More informationPopulation Genetics I. Bio
Population Genetics I. Bio5488-2018 Don Conrad dconrad@genetics.wustl.edu Why study population genetics? Functional Inference Demographic inference: History of mankind is written in our DNA. We can learn
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationRare Alleles and Selection
Theoretical Population Biology 59, 8796 (001) doi:10.1006tpbi.001.153, available online at http:www.idealibrary.com on Rare Alleles and Selection Carsten Wiuf Department of Statistics, University of Oxford,
More informationA Method for Evaluating the Distribution of the Total Cost of a Random Process over its Lifetime
A Method for Evaluating the Distribution of the Total Cost of a Random Process over its Lifetime P.K. Pollett a and V.T. Stefanov b a Department of Mathematics, The University of Queensland Queensland
More informationWXML Final Report: Chinese Restaurant Process
WXML Final Report: Chinese Restaurant Process Dr. Noah Forman, Gerandy Brito, Alex Forney, Yiruey Chou, Chengning Li Spring 2017 1 Introduction The Chinese Restaurant Process (CRP) generates random partitions
More informationGENETIC VARIABILITY AND RATE OF GENE SUBSTITUTION IN A FINITE POPULATION UNDER MUTATION AND FLUCTUATING SELECTION* NAOYUKI TAKAHATA ABSTRACT
GENETC VARABLTY AND RATE OF GENE SUBSTTUTON N A FNTE POPULATON UNDER MUTATON AND FLUCTUATNG SELECTON* NAOYUK TAKAHATA National nstitute of Genetics, Mishima, 411 Japan ABSTRACT By using a numerical method
More informationMathematical Population Genetics II
Mathematical Population Genetics II Lecture Notes Joachim Hermisson March 20, 2015 University of Vienna Mathematics Department Oskar-Morgenstern-Platz 1 1090 Vienna, Austria Copyright (c) 2013/14/15 Joachim
More informationProcesses of Evolution
15 Processes of Evolution Forces of Evolution Concept 15.4 Selection Can Be Stabilizing, Directional, or Disruptive Natural selection can act on quantitative traits in three ways: Stabilizing selection
More informationEvolutionary Theory. Sinauer Associates, Inc. Publishers Sunderland, Massachusetts U.S.A.
Evolutionary Theory Mathematical and Conceptual Foundations Sean H. Rice Sinauer Associates, Inc. Publishers Sunderland, Massachusetts U.S.A. Contents Preface ix Introduction 1 CHAPTER 1 Selection on One
More informationPITMAN S 2M X THEOREM FOR SKIP-FREE RANDOM WALKS WITH MARKOVIAN INCREMENTS
Elect. Comm. in Probab. 6 (2001) 73 77 ELECTRONIC COMMUNICATIONS in PROBABILITY PITMAN S 2M X THEOREM FOR SKIP-FREE RANDOM WALKS WITH MARKOVIAN INCREMENTS B.M. HAMBLY Mathematical Institute, University
More informationAn Introduction to Stochastic Modeling
F An Introduction to Stochastic Modeling Fourth Edition Mark A. Pinsky Department of Mathematics Northwestern University Evanston, Illinois Samuel Karlin Department of Mathematics Stanford University Stanford,
More informationMemoirs of My Research on Stochastic Analysis
Memoirs of My Research on Stochastic Analysis Kiyosi Itô Professor Emeritus, Kyoto University, Kyoto, 606-8501 Japan It is with great honor that I learned of the 2005 Oslo Symposium on Stochastic Analysis
More informationGene Pool Recombination in Genetic Algorithms
Gene Pool Recombination in Genetic Algorithms Heinz Mühlenbein GMD 53754 St. Augustin Germany muehlenbein@gmd.de Hans-Michael Voigt T.U. Berlin 13355 Berlin Germany voigt@fb10.tu-berlin.de Abstract: A
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More information2 Discrete-Time Markov Chains
2 Discrete-Time Markov Chains Angela Peace Biomathematics II MATH 5355 Spring 2017 Lecture notes follow: Allen, Linda JS. An introduction to stochastic processes with applications to biology. CRC Press,
More informationThe Transition Probability Function P ij (t)
The Transition Probability Function P ij (t) Consider a continuous time Markov chain {X(t), t 0}. We are interested in the probability that in t time units the process will be in state j, given that it
More informationTutorial on Theoretical Population Genetics
Tutorial on Theoretical Population Genetics Joe Felsenstein Department of Genome Sciences and Department of Biology University of Washington, Seattle Tutorial on Theoretical Population Genetics p.1/40
More informationSTAT 536: Genetic Statistics
STAT 536: Genetic Statistics Frequency Estimation Karin S. Dorman Department of Statistics Iowa State University August 28, 2006 Fundamental rules of genetics Law of Segregation a diploid parent is equally
More informationIEOR 6711: Stochastic Models I Professor Whitt, Thursday, November 29, Weirdness in CTMC s
IEOR 6711: Stochastic Models I Professor Whitt, Thursday, November 29, 2012 Weirdness in CTMC s Where s your will to be weird? Jim Morrison, The Doors We are all a little weird. And life is a little weird.
More informationOutline of lectures 3-6
GENOME 453 J. Felsenstein Evolutionary Genetics Autumn, 007 Population genetics Outline of lectures 3-6 1. We want to know what theory says about the reproduction of genotypes in a population. This results
More informationIntroduction to population genetics & evolution
Introduction to population genetics & evolution Course Organization Exam dates: Feb 19 March 1st Has everybody registered? Did you get the email with the exam schedule Summer seminar: Hot topics in Bioinformatics
More informationHadamard matrices of order 32
Hadamard matrices of order 32 H. Kharaghani a,1 B. Tayfeh-Rezaie b a Department of Mathematics and Computer Science, University of Lethbridge, Lethbridge, Alberta, T1K3M4, Canada b School of Mathematics,
More informationIntroductory seminar on mathematical population genetics
Exercises Sheets Introductory seminar on mathematical population genetics WS 20/202 Kristan Schneider, Ada Akerman Ex Assume a single locus with alleles A and A 2 Denote the frequencies of the three (unordered
More informationExact results for fixation probability of bithermal evolutionary graphs.
Exact results for fixation probability of bithermal evolutionary graphs. Bahram Houchmandzadeh, Marcel Vallade To cite this version: Bahram Houchmandzadeh, Marcel Vallade. Exact results for fixation probability
More informationProblems on Evolutionary dynamics
Problems on Evolutionary dynamics Doctoral Programme in Physics José A. Cuesta Lausanne, June 10 13, 2014 Replication 1. Consider the Galton-Watson process defined by the offspring distribution p 0 =
More informationBustamante et al., Supplementary Nature Manuscript # 1 out of 9 Information #
Bustamante et al., Supplementary Nature Manuscript # 1 out of 9 Details of PRF Methodology In the Poisson Random Field PRF) model, it is assumed that non-synonymous mutations at a given gene are either
More informationOutline of lectures 3-6
GENOME 453 J. Felsenstein Evolutionary Genetics Autumn, 009 Population genetics Outline of lectures 3-6 1. We want to know what theory says about the reproduction of genotypes in a population. This results
More informationON THE DIFFUSION OPERATOR IN POPULATION GENETICS
J. Appl. Math. & Informatics Vol. 30(2012), No. 3-4, pp. 677-683 Website: http://www.kcam.biz ON THE DIFFUSION OPERATOR IN POPULATION GENETICS WON CHOI Abstract. W.Choi([1]) obtains a complete description
More informationMathematical Biology. Sabin Lessard
J. Math. Biol. DOI 10.1007/s00285-008-0248-1 Mathematical Biology Diffusion approximations for one-locus multi-allele kin selection, mutation and random drift in group-structured populations: a unifying
More informationLecture 4a: Continuous-Time Markov Chain Models
Lecture 4a: Continuous-Time Markov Chain Models Continuous-time Markov chains are stochastic processes whose time is continuous, t [0, ), but the random variables are discrete. Prominent examples of continuous-time
More informationLecture 1 Hardy-Weinberg equilibrium and key forces affecting gene frequency
Lecture 1 Hardy-Weinberg equilibrium and key forces affecting gene frequency Bruce Walsh lecture notes Introduction to Quantitative Genetics SISG, Seattle 16 18 July 2018 1 Outline Genetics of complex
More informationDepartment of Applied Mathematics Faculty of EEMCS. University of Twente. Memorandum No Birth-death processes with killing
Department of Applied Mathematics Faculty of EEMCS t University of Twente The Netherlands P.O. Box 27 75 AE Enschede The Netherlands Phone: +3-53-48934 Fax: +3-53-48934 Email: memo@math.utwente.nl www.math.utwente.nl/publications
More informationA simple genetic model with non-equilibrium dynamics
J. Math. Biol. (1998) 36: 550 556 A simple genetic model with non-equilibrium dynamics Michael Doebeli, Gerdien de Jong Zoology Institute, University of Basel, Rheinsprung 9, CH-4051 Basel, Switzerland
More informationFollow links Class Use and other Permissions. For more information, send to:
COPYRIGHT NOTICE: Sarah P. Otto & Troy Day: A Biologist's Guide to Mathematical Modeling in Ecology and Evolution is published by Princeton University Press and copyrighted, 2007, by Princeton University
More informationUNIVERSITY OF CALIFORNIA, MERCED. Genetic Correlation in the One-Dimensional Stepping Stone Model of Population Structure. by Elizabeth Owens
UNIVERSITY OF CALIFORNIA, MERCED MASTER S THESIS Genetic Correlation in the One-Dimensional Stepping Stone Model of Population Structure by Elizabeth Owens A technical report submitted in partial fulfillment
More informationP(I -ni < an for all n > in) = 1 - Pm# 1
ITERATED LOGARITHM INEQUALITIES* By D. A. DARLING AND HERBERT ROBBINS UNIVERSITY OF CALIFORNIA, BERKELEY Communicated by J. Neyman, March 10, 1967 1. Introduction.-Let x,x1,x2,... be a sequence of independent,
More informationModeling Evolution DAVID EPSTEIN CELEBRATION. John Milnor. Warwick University, July 14, Stony Brook University
Modeling Evolution John Milnor Stony Brook University DAVID EPSTEIN CELEBRATION Warwick University, July 14, 2007 A First Model for Evolution: The Phylogenetic Tree, U time V W X A B C D E F G 1500 1000
More informationMATHEMATICAL MODELS OF LIFE SUPPORT SYSTEMS Vol. II - Population Models - V.P. Passekov
POPULATION MODELS V.P. Passekov Division of Modeling Problems, Computing Center of Russian Academy of Sciences, Russia Keywords: Population Models, Discrete population Models, Continuous Population Models,
More informationINVARIANT SUBSETS OF THE SEARCH SPACE AND THE UNIVERSALITY OF A GENERALIZED GENETIC ALGORITHM
INVARIANT SUBSETS OF THE SEARCH SPACE AND THE UNIVERSALITY OF A GENERALIZED GENETIC ALGORITHM BORIS MITAVSKIY Abstract In this paper we shall give a mathematical description of a general evolutionary heuristic
More informationGENETIC DRIFT IN A CLINE MONTGOMERY SLATKIN
GENETIC DRIFT IN A CLINE MONTGOMERY SLATKIN Department of Biophysics ana! Theoretical Biology, The University of Chicago, 920 East 58th Street, Chicago: Illinois 60637 and TAKE0 MARUYAMA National Institute
More information