No Arabic abstract
The dynamics of a population undergoing selection is a central topic in evolutionary biology. This question is particularly intriguing in the case where selective forces act in opposing directions at two population scales. For example, a fast-replicating virus strain outcompetes slower-replicating strains at the within-host scale. However, if the fast-replicating strain causes host morbidity and is less frequently transmitted, it can be outcompeted by slower-replicating strains at the between-host scale. Here we consider a stochastic ball-and-urn process which models this type of phenomenon. We prove the weak convergence of this process under two natural scalings. The first scaling leads to a deterministic nonlinear integro-partial differential equation on the interval $[0,1]$ with dependence on a single parameter, $lambda$. We show that the fixed points of this differential equation are Beta distributions and that their stability depends on $lambda$ and the behavior of the initial data around $1$. The second scaling leads to a measure-valued Fleming-Viot process, an infinite dimensional stochastic process that is frequently associated with a population genetics.
We formulate a mathematical model for daily activities of a cow (eating, lying down, and standing) in terms of a piecewise affine dynamical system. We analyze the properties of this bovine dynamical system representing the single animal and develop an exact integrative form as a discrete-time mapping. We then couple multiple cow oscillators together to study synchrony and cooperation in cattle herds. We comment on the relevant biology and discuss extensions of our model. With this abstract approach, we not only investigate equations with interesting dynamics but also develop interesting biological predictions. In particular, our model illustrates that it is possible for cows to synchronize emph{less} when the coupling is increased.
We discuss various limits of a simple random exchange model that can be used for the distribution of wealth. We start from a discrete state space - discrete time version of this model and, under suitable scaling, we show its functional convergence to a continuous space - discrete time model. Then, we show a thermodynamic limit of the empirical distribution to the solution of a kinetic equation of Boltzmann type. We solve this equation and we show that the solutions coincide with the appropriate limits of the invariant measure for the Markov chain. In this way we complete Boltzmanns program of deriving kinetic equations from random dynamics for this simple model. Three families of invariant measures for the mean field limit are discovered and we show that only two of those families can be obtained as limits of the discrete system and the third is extraneous. Finally, we cast our results in the framework of integer partitions and strengthen some results already available in the literature.
A number of discrete time, finite population size models in genetics describing the dynamics of allele frequencies are known to converge (subject to suitable scaling) to a diffusion process in the infinite population limit, termed the Wright-Fisher diffusion. In this article we show that the diffusion is ergodic uniformly in the selection and mutation parameters, and that the measures induced by the solution to the stochastic differential equation are uniformly locally asymptotically normal. Subsequently these two results are used to analyse the statistical properties of the Maximum Likelihood and Bayesian estimators for the selection parameter, when both selection and mutation are acting on the population. In particular, it is shown that these estimators are uniformly over compact sets consistent, display uniform in the selection parameter asymptotic normality and convergence of moments over compact sets, and are asymptotically efficient for a suitable class of loss functions.
We consider a one-dimensional traffic model with a slow-to-start rule. The initial position of the cars in $mathbb R$ is a Poisson process of parameter $lambda$. Cars have speed 0 or 1 and travel in the same direction. At time zero the speed of all cars is 0; each car waits an exponential time to switch speed from $0$ to $1$ and stops when it collides with a stopped car. When the car is no longer blocked, it waits a new exponential time to assume speed one, and so on. We study the emergence of condensation for the saturated regime $lambda>1$ and the critical regime $lambda=1$, showing that in both regimes all cars collide infinitely often and each car has asymptotic mean velocity $1/lambda$. In the saturated regime the moving cars form a point process whose intensity tends to 1. The remaining cars condensate in a set of points whose intensity tends to zero as $1/sqrt t$. We study the scaling limit of the traffic jam evolution in terms of a collection of coalescing Brownian motions.
In this paper, a mathematical model is proposed to analyze the dynamic behavior of COVID-19. Based on inter-city networked coupling effects, a fractional-order SEIHDR system with the real-data from 23 January to 18 March, 2020 of COVID-19 is discussed. Meanwhile, hospitalized individuals and the mortality rates of three types of individuals (exposed, infected and hospitalized) are firstly taken into account in the proposed model. And infectivity of individuals during incubation is also considered in this paper. By applying least squares method and predictor-correctors scheme, the numerical solutions of the proposed system in the absence of the inter-city network and with the inter-city network are stimulated by using the real-data from 23 January to $18-m$ March, 2020 where $m$ is equal to the number of prediction days. Compared with integer-order system ($alpha=0$), the fractional-order model without network is validated to have a better fitting of the data on Beijing, Shanghai, Wuhan, Huanggang and other cities. In contrast to the case without network, the results indicate that the inter-city network system may be not a significant case to virus spreading for China because of the lock down and quarantine measures, however, it may have an impact on cities that have not adopted city closure. Meanwhile, the proposed model better fits the data from 24 February to 31, March in Italy, and the peak number of confirmed people is also predicted by this fraction-order model. Furthermore, the existence and uniqueness of a bounded solution under the initial condition are considered in the proposed system. Afterwards, the basic reproduction number $R_0$ is analyzed and it is found to hold a threshold: the disease-free equilibrium point is locally asymptotically stable when $R_0le 1$, which provides a theoretical basis for whether COVID-19 will become a pandemic in the future.