Search

We study the effectiveness of regional lockdown strategies to mitigate the spread of a pathogen across regional units, in the following called cities, within a country or region for a single infection wave. The heterogeneity in the epidemically relevant connectivity is defined via a random network model with cities as nodes, where the city's sizes determine their connectivity via a gravity type kernel function. Isolation of a whole city is initiated when infection numbers surpass defined thresholds. We consider two basic strategies for the lockdowns. Strategy~$(P)$ isolates cities based on a proportional threshold of infections, while stra\-tegy~\((U)\) uses a uniform infection threshold for all cities. Given the heavy-tailed distribution of city sizes, strategy \((P)\) can potentially result in more secondary infections from larger cities than strategy \((U)\). As an efficiency measure we use the ratio of individuals under lockdown and the number of infected individuals. Additionally, we analytically derive formulas for the basic reproduction numbers and prevalences. Our model is fitted to mobility data from France, Japan, and Poland, and validated through simulations. The findings indicate that while the model aligns well with data from France and Poland, it deviates in Japan, highlighting the importance of geographical nuances in pathogen spread modeling. Furthermore, it suggests that for France (and Japan) both strategies perform equally well, while for Poland strategy \((U)\) outperforms strategy \((P)\)., Comment: main text: 48 pages, supplementary material: 15p, 109 figures (including the supplementary material)

We study the long-time behavior of an additive functional that takes into account the jumps of a symmetric Markov process. This process is assumed to be observed through a biased observation scheme that includes the survival to events of extinction and the Feynman-Kac weight by another similar additive functional. Under conditioning for the convergence to a quasi-stationary distribution and for two-sided estimates of the Feynmac-Kac semigroup to be obtained, we shall discuss general assumptions on the symmetric Markov process. For the law of additive functionals, we will prove a quasi-ergodic theorem, namely a conditional version of the ergodic theorem and a conditional functional weak law of large numbers. As an application, we also establish a large deviation principle for the mean ratio of additive functionals., Comment: 25 pages, no figure

Starting from a stochastic individual-based description of an SIS epidemic spreading on a random network, we study the dynamics when the size $n$ of the network tends to infinity. We recover in the limit an infinite-dimensional integro-differential equation studied by Delmas, Dronnier and Zitt (2022) for an SIS epidemic propagating on a graphon. Our work covers the case of dense and sparse graphs, provided that the number of edges grows faster than $n$, but not the case of very sparse graphs with $O(n)$ edges. In order to establish our limit theorem, we have to deal with both the convergence of the random graphs to the graphon and the convergence of the stochastic process spreading on top of these random structures: in particular, we propose a coupling between the process of interest and an epidemic that spreads on the complete graph but with a modified infection rate. Keywords: Random graph, mathematical models of epidemics, measure-valued process, large network limit, limit theorem, graphon., Comment: This work was financed by the Labex B\'ezout (ANR-10-LABX-58) and the COCOON grant (ANR-22-CE48-0011), and by the platform MODCOV19 of the CNRS. V.C.T. is partly financed by the Chaire "Mod\'elisation Math\'ematique et Biodiversit\'e''. The research leading to this article was largely performed while A. Velleret was a researcher at LAMA and at GIPSA-Lab

We prove the existence and uniqueness of a quasi-stationary distribution for three stochastic processes derived from the model of Muller's ratchet. This model was invented with the aim of evaluating the limitations of an asexual reproduction mode in preventing the accumulation of deleterious mutations through natural selection alone. The main considered model is non-classical, as it is a stochastic diffusion evolving on an irregular set of infinite dimension with hard killing on an hyperplane. We are nonetheless able to prove exponential convergence in total variation to the quasi-stationary distribution even in this case. The parameters in this last convergence result are directly related to the core parameters of Muller's ratchet. The speed of convergence to the quasi-stationary distribution is deduced both for the infinite dimensional model and for approximations with a large yet finite number of potential mutations. Likewise, we give uniform moment estimates of the empirical distribution of mutations in the population under quasi-stationarity., Comment: 74 pages, no figure

The purpose of this paper is to ensure the conditions of G\"artner-Ellis Theorem for evaluations of the empirical measure. We show that up-to-date conditions for ensuring the convergence to a quasi-stationary distribution can be applied efficiently. By this mean, we are able to prove Large Deviation results even with a conditioning that the process is not extinct at the end of the evaluation. The domain on which these Large Deviation results apply is implicitly given by the range of penalization for which one can prove the above-mentioned results of quasi-stationarity. We propose a way to relate the range of controlled deviations to the range of admissible penalization. Central Limit Theorems are deduced from these results of Large Deviations. As an application, we consider the empirical measure of position and jumps of a continuous-time process on aunbounded domain of $R^d$. This model is inspired by the adaptation of a population to a changing environment. Jumps makes it possible for the process to face a deterministic dynamics leading to high extinction areas., Comment: The argument in the proof of Proposition 1.3.4 is wrong and I cannot give simple conditions ensuring that the function is strictly increasing. This makes the statements much less clear if we cannot ensure it, so I hope I will be able to handle it. Besides, the limit in Lemma 3.3.2 is wrong. Nonetheless, the proof that the variance goes to 0 is true

This article is a presentation of specific recent results describing scaling limits of individual-based models. Thanks to them, we wish to relate the time-scales typical of demographic dynamics and natural selection to the parameters of the individual-based models. Although these results are by no means exhaustive, both on the mathematical and the biological level, they complement each other. Indeed, they provide a viewpoint for many classical time-scales. Namely, they encompass the timescale typical of the life-expectancy of a single individual, the longer one wherein a population can be characterized through its demographic dynamics, and at least four interconnected ones wherein selection occurs. The limiting behavior is generally deterministic. Yet, since there are selective effects on randomness in the history of lineages, probability theory is shown to be a key factor in understanding the results. Besides, randomness can be maintained in the limiting dynamics, for instance to model rare mutations fixing in the population., Comment: Review on several papers based on the presentations given during the session Stochastic Processes and Biology at the conference Journ\'{e}es MAS 2018 (Mod\'{e}lisation Al\'{e}atoire et Stochastiques). arXiv admin note: text overlap with arXiv:1310.6274, arXiv:1711.10732, arXiv:1505.02421, arXiv:1507.00397 by other authors. additional note of the author : also arXiv:1903.10161

In a view for a simple model where natural selection at the individual level is confronted to selection effects at the group level, we consider some individual-based models of some large population subdivided into a large number of groups. We then obtain the convergence to the law of a stochastic process with some Feynman-Kac penalization. To analyze the limiting behavior of this law, we exploit a recent approach, designed for the convergence to quasi-stationary distributions. We are able to deal with the fixation of the stochastic process and relate the convergence to equilibrium to the one where fixation implies extinction. We notably establish different regimes of convergence. Besides the case of an exponential rate (the rate being uniform over the initial condition), critical regimes with convergence in 1/t are also to notice. We finally address the relevance of such limiting behaviors to predict the long-time behavior of the individual-based model and describe more specifically the cases of weak selection. Consequences in term of evolutive dynamics are also derived, where such competition is assumed to occur repeatedly at each de novo mutation., Comment: 38 pages, 13 figures

We analyze the long-term stability of a stochastic model designed to illustrate the adaptation of a population to variation in its environment. A piecewise-deterministic process modeling adaptation is coupled to a Feller logistic diffusion modeling population size. As the individual features in the population become further away from the optimal ones, the growth rate declines, making population extinction more likely. Assuming that the environment changes deterministically and steadily in a constant direction, we obtain the existence and uniqueness of the quasi-stationary distribution, the associated survival capacity and the Q-process. Our approach also provides several exponential convergence results (in total variation for the measures). From this synthetic information, we can characterize the efficiency of internal adaptation (i.e. population turnover from mutant invasions). When the latter is lacking, there is still stability, but because of the high level of population extinction. Therefore, such a characterization must be based on specific features of this quasi-ergodic regime., Comment: 48 pages, 15 figure

This paper tackles the issue of establishing a lower-bound on the asymptotic ratio of survival probabilities between two different initial conditions, asymptotically in time for a given Markov process with extinction. Such a comparison is a crucial step in recent techniques for proving exponential convergence to a quasi-stationary distribution. We introduce a weak form of the Harnack inequality as the essential ingredient for such a comparison. This property is actually a consequence of the convergence property we intend to prove. Its complexity appears as the price to pay for the level of flexibility required by our applications. We show in our illustrations how simply and efficiently it can be used nonetheless. As illustrations, we consider two continuous-time processes on R^d that do not satisfy the classical Harnack inequalities, even in a local version. The first one is a piecewise deterministic process while the second is a pure jump process with restrictions on the directions of its jumps., Comment: Keywords : continuous-time and continuous-space Markov process; jumps; quasistationary distribution; survival capacity; Q-process; Harris recurrence

We establish sufficient conditions for exponential convergence to a unique quasi-stationary distribution in the total variation norm. These conditions also ensure the existence and exponential ergodicity of the Q-process, the process conditionned upon never being absorbed. The technique relies on a coupling procedure that is related to Harris recurrence (for Markov Chains). It applies to general continuous-time and continuous-space Markov processes. The main novelty is that we modulate each coupling step depending both on a final horizon of time (for survival) and on the initial distribution. By this way, we could notably include in the convergence a dependency on the initial condition. As an illustration, we consider a continuous-time birth-death process with catastrophes and a diffusion process describing a (localized) population adapting to its environment., Comment: 2 figures