Source author record

Pieter Trapman

Pieter Trapman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2020arXiv

Key Questions for Modelling COVID-19 Exit Strategies

Combinations of intense non-pharmaceutical interventions ('lockdowns') were introduced in countries worldwide to reduce SARS-CoV-2 transmission. Many governments have begun to implement lockdown exit strategies that allow restrictions to be relaxed while attempting to control the risk of a surge in cases. Mathematical modelling has played a central role in guiding interventions, but the challenge of designing optimal exit strategies in the face of ongoing transmission is unprecedented. Here, we report discussions from the Isaac Newton Institute 'Models for an exit strategy' workshop (11-15 May 2020). A diverse community of modellers who are providing evidence to governments worldwide were asked to identify the main questions that, if answered, will allow for more accurate predictions of the effects of different exit strategies. Based on these questions, we propose a roadmap to facilitate the development of reliable models to guide exit strategies. The roadmap requires a global collaborative effort from the scientific community and policy-makers, and is made up of three parts: i) improve estimation of key epidemiological parameters; ii) understand sources of heterogeneity in populations; iii) focus on requirements for data collection, particularly in Low-to-Middle-Income countries. This will provide important information for planning exit strategies that balance socio-economic benefits with public health.

preprint2020arXiv

The disease-induced herd immunity level for Covid-19 is substantially lower than the classical herd immunity level

Most countries are suffering severely from the ongoing covid-19 pandemic despite various levels of preventive measures. A common question is if and when a country or region will reach herd immunity $h$. The classical herd immunity level $h_C$ is defined as $h_C=1-1/R_0$, where $R_0$ is the basic reproduction number, for covid-19 estimated to lie somewhere in the range 2.2-3.5 depending on country and region. It is shown here that the disease-induced herd immunity level $h_D$, after an outbreak has taken place in a country/region with a set of preventive measures put in place, is actually substantially smaller than $h_C$. As an illustration we show that if $R_0=2.5$ in an age-structured community with mixing rates fitted to social activity studies, and also categorizing individuals into three categories: low active, average active and high active, and where preventive measures affect all mixing rates proportionally, then the disease-induced herd immunity level is $h_D=43\%$ rather than $h_C=1-1/2.5=60\%$. Consequently, a lower fraction infected is required for herd immunity to appear. The underlying reason is that when immunity is induced by disease spreading, the proportion infected in groups with high contact rates is greater than that in groups with low contact rates. Consequently, disease-induced immunity is stronger than when immunity is uniformly distributed in the community as in the classical herd immunity level.

preprint2020arXiv

The duration of an $SIR$ epidemic on a configuration model

We consider the spread of a supercritical stochastic SIR (Susceptible, Infectious, Recovered) epidemic on a configuration model random graph. We mainly focus on the final stages of a large outbreak and provide limit results for the duration of the entire epidemic, while we allow for non-exponential distributions of the infectious period and for both finite and infinite variance of the asymptotic degree distribution in the graph. Our analysis relies on the analysis of some subcritical continuous time branching processes and on ideas from first-passage percolation. As an application we investigate the effect of vaccination with an all-or-nothing vaccine on the duration of the epidemic. We show that if vaccination fails to prevent the epidemic, it often -- but not always -- increases the duration of the epidemic.

preprint2016arXiv

A Dynamic Erdős-Rényi Graph Model

In this article we introduce a dynamic Erdős-Rényi graph model, in which, independently for each vertex pair, edges appear and disappear according to a Markov on-off process. In studying the dynamic graph we present two main results. The first being on how long it takes for the graph to reach stationarity. We give an explicit expression for this time, as well as proving that this is the fastest time to reach stationarity among all strong stationary times. The second result concerns the time it takes for the dynamic graph to reach a certain number of edges. We give an explicit expression for the expected value of such a time, as well as study its asymptotic behavior. This time is related to the first time the dynamic Erdős-Rényi graph contains a cluster exceeding a certain size.

preprint2016arXiv

Inferring $R_0$ in emerging epidemics - the effect of common population structure is small

When controlling an emerging outbreak of an infectious disease it is essential to know the key epidemiological parameters, such as the basic reproduction number $R_0$ and the control effort required to prevent a large outbreak. These parameters are estimated from the observed incidence of new cases and information about the infectious contact structures of the population in which the disease spreads. However, the relevant infectious contact structures for new, emerging infections are often unknown or hard to obtain. Here we show that for many common true underlying heterogeneous contact structures, the simplification to neglect such structures and instead assume that all contacts are made homogeneously in the whole population, results in conservative estimates for $R_0$ and the required control effort. This means that robust control policies can be planned during the early stages of an outbreak, using such conservative estimates of the required control effort.

preprint2015arXiv

An epidemic in a dynamic population with importation of infectives

Consider a large uniformly mixing dynamic population, which has constant birth rate and exponentially distributed lifetimes, with mean population size $n$. A Markovian SIR (susceptible $\to$ infective $\to$ recovered) infectious disease, having importation of infectives, taking place in this population is analysed. The main situation treated is where $n\to\infty$, keeping the basic reproduction number $R_0$ as well as the importation rate of infectives fixed, but assuming that the quotient of the average infectious period and the average lifetime tends to 0 faster than $1/\log n$. It is shown that, as $ n \to \infty$, the behaviour of the 3-dimensional process describing the evolution of the fraction of the population that are susceptible, infective and recovered, is encapsulated in a 1-dimensional regenerative process $S=\{ S(t);t\ge 0\}$ describing the limiting fraction of the population that are susceptible. The process $S$ grows deterministically, except at one random time point per regenerative cycle, where it jumps down by a size that is completely determined by the waiting time since the previous jump. Properties of the process $S$, including the jump size and stationary distributions, are determined.

preprint2015arXiv

Reproduction numbers for epidemic models with households and other social structures II: comparisons and implications for vaccination

In this paper we consider epidemic models of directly transmissible SIR (susceptible $\to$ infective $\to$ recovered) and SEIR (with an additional latent class) infections in fully-susceptible populations with a social structure, consisting either of households or of households and workplaces. We review most reproduction numbers defined in the literature for these models, including the basic reproduction number $R_0$ introduced in the companion paper of this, for which we provide a simpler, more elegant derivation. Extending previous work, we provide a complete overview of the inequalities among these reproduction numbers and resolve some open questions. Special focus is put on the exponential-growth-associated reproduction number $R_r$, which is loosely defined as the estimate of $R_0$ based on the observed exponential growth of an emerging epidemic obtained when the social structure is ignored. We show that for the vast majority of the models considered in the literature $R_r \geq R_0$ when $R_0 \ge 1$ and $R_r \leq R_0$ when $R_0 \le 1$. We show that, in contrast to models without social structure, vaccination of a fraction $1-1/R_0$ of the population, chosen uniformly at random, with a perfect vaccine is usually insufficient to prevent large epidemics. In addition, we provide significantly sharper bounds than the existing ones for bracketing the critical vaccination coverage between two analytically tractable quantities, which we illustrate by means of extensive numerical examples.

preprint2015arXiv

Stochastic SIR epidemics in a population with households and schools

We study the spread of stochastic SIR (Susceptible $\to$ Infectious $\to$ Recovered) epidemics in two types of structured populations, both consisting of schools and households. In each of the types, every individual is part of one school and one household. In the independent partition model, the partitions of the population into schools and households are independent of each other. This model corresponds to the well-studied household-workplace model. In the hierarchical model which we introduce here, members of the same household are also members of the same school. We introduce computable branching process approximations for both types of populations and use these to compare the probabilities of a large outbreak. The branching process approximation in the hierarchical model is novel and of independent interest. We prove by a coupling argument that if all households and schools have the same size, an epidemic spreads easier (in the sense that the number of individuals infected is stochastically larger) in the independent partition model. We also show by example that this result does not necessarily hold if households and/or schools do not all have the same size.

preprint2014arXiv

Epidemics on random intersection graphs

In this paper we consider a model for the spread of a stochastic SIR (Susceptible $\to$ Infectious $\to$ Recovered) epidemic on a network of individuals described by a random intersection graph. Individuals belong to a random number of cliques, each of random size, and infection can be transmitted between two individuals if and only if there is a clique they both belong to. Both the clique sizes and the number of cliques an individual belongs to follow mixed Poisson distributions. An infinite-type branching process approximation (with type being given by the length of an individual's infectious period) for the early stages of an epidemic is developed and made fully rigorous by proving an associated limit theorem as the population size tends to infinity. This leads to a threshold parameter $R_*$, so that in a large population an epidemic with few initial infectives can give rise to a large outbreak if and only if $R_*>1$. A functional equation for the survival probability of the approximating infinite-type branching process is determined; if $R_*\le1$, this equation has no nonzero solution, while if $R_*>1$, it is shown to have precisely one nonzero solution. A law of large numbers for the size of such a large outbreak is proved by exploiting a single-type branching process that approximates the size of the susceptibility set of a typical individual.

preprint2013arXiv

Stochastic epidemics in growing populations

Consider a uniformly mixing population which grows as a super-critical linear birth and death process. At some time an infectious disease (of SIR or SEIR type) is introduced by one individual being infected from outside. It is shown that three different scenarios may occur: 1) an epidemic never takes off, 2) an epidemic gets going and grows but at a slower rate than the community thus still being negligible in terms of population fractions, or 3) an epidemic takes off and grows quicker than the community eventually leading to an endemic equilibrium. Depending on the parameter values, either scenario 1 is the only possibility, both scenario 1 and 2 are possible, or scenario 1 and 3 are possible.

preprint2012arXiv

Inferring global network properties from egocentric data with applications to epidemics

Social networks are rarely observed in full detail. In many situations properties are known for only a sample of the individuals in the network and it is desirable to induce global properties of the full social network from this "egocentric" network data. In the current paper we study a few different types of egocentric data, and show what global network properties are consistent with those egocentric data. Two global network properties are considered: the size of the largest connected component in the network (the giant), and secondly, the possible size of an epidemic outbreak taking place on the network, in which transmission occurs only between network neighbours, and with probability $p$. The main conclusion is that in most cases, egocentric data allow for a large range of possible sizes of the giant and the outbreak. However, there is an upper bound for the latter. For the case that the network is selected uniformly among networks with prescribed egocentric data (satisfying some conditions), the asymptotic size of the giant and the outbreak is characterised.

preprint2011arXiv

Splitting trees stopped when the first clock rings and Vervaat's transformation

We consider a branching population where individuals have i.i.d.\ life lengths (not necessarily exponential) and constant birth rate. We let $N_t$ denote the population size at time $t$. %(called homogeneous, binary Crump--Mode--Jagers process). We further assume that all individuals, at birth time, are equipped with independent exponential clocks with parameter $δ$. We are interested in the genealogical tree stopped at the first time $T$ when one of those clocks rings. This question has applications in epidemiology, in population genetics, in ecology and in queuing theory. We show that conditional on $\{T<\infty\}$, the joint law of $(N_T, T, X^{(T)})$, where $X^{(T)}$ is the jumping contour process of the tree truncated at time $T$, is equal to that of $(M, -I_M, Y_M')$ conditional on $\{M\not=0\}$, where : $M+1$ is the number of visits of 0, before some single independent exponential clock $\mathbf{e}$ with parameter $δ$ rings, by some specified L{é}vy process $Y$ without negative jumps reflected below its supremum; $I_M$ is the infimum of the path $Y_M$ defined as $Y$ killed at its last 0 before $\mathbf{e}$; $Y_M'$ is the Vervaat transform of $Y_M$. This identity yields an explanation for the geometric distribution of $N_T$ \cite{K,T} and has numerous other applications. In particular, conditional on $\{N_T=n\}$, and also on $\{N_T=n, T<a\}$, the ages and residual lifetimes of the $n$ alive individuals at time $T$ are i.i.d.\ and independent of $n$. We provide explicit formulae for this distribution and give a more general application to outbreaks of antibiotic-resistant bacteria in the hospital.

preprint2011arXiv

The growth of the infinite long-range percolation cluster

We consider long-range percolation on $\mathbb{Z}^d$, where the probability that two vertices at distance $r$ are connected by an edge is given by $p(r)=1-\exp[-λ(r)]\in(0,1)$ and the presence or absence of different edges are independent. Here, $λ(r)$ is a strictly positive, nonincreasing, regularly varying function. We investigate the asymptotic growth of the size of the $k$-ball around the origin, $|\mathcal{B}_k|$, that is, the number of vertices that are within graph-distance $k$ of the origin, for $k\to\infty$, for different $λ(r)$. We show that conditioned on the origin being in the (unique) infinite cluster, nonempty classes of nonincreasing regularly varying $λ(r)$ exist, for which, respectively: $\bullet$ $|\mathcal{B}_k|^{1/k}\to\infty$ almost surely; $\bullet$ there exist $1<a_1<a_2<\infty$ such that $\lim_{k\to \infty}\mathbb{P}(a_1<|\mathcal{B}_k|^{1/k}<a_2)=1$; $\bullet$ $|\mathcal{B}_k|^{1/k}\to1$ almost surely. This result can be applied to spatial SIR epidemics. In particular, regimes are identified for which the basic reproduction number, $R_0$, which is an important quantity for epidemics in unstructured populations, has a useful counterpart in spatial epidemics.

preprint2010arXiv

Bounding basic characteristics of spatial epidemics with a new percolation model

We introduce a new percolation model to describe and analyze the spread of an epidemic on a general directed and locally finite graph. We assign a two-dimensional random weight vector to each vertex of the graph in such a way that the weights of different vertices are i.i.d., but the two entries of the vector assigned to a vertex need not be independent. The probability for an edge to be open depends on the weights of its end vertices, but conditionally on the weights, the states of the edges are independent of each other. In an epidemiological setting, the vertices of a graph represent the individuals in a (social) network and the edges represent the connections in the network. The weights assigned to an individual denote its (random) infectivity and susceptibility, respectively. We show that one can bound the percolation probability and the expected size of the cluster of vertices that can be reached by an open path starting at a given vertex from above and below by the corresponding quantities for respectively independent bond and site percolation with certain densities; this generalizes a result of Kuulasmaa. Many models in the literature are special cases of our general model.

preprint2010arXiv

Long-range percolation on the hierarchical lattice

We study long-range percolation on the hierarchical lattice of order $N$, where any edge of length $k$ is present with probability $p_k=1-\exp(-β^{-k} α)$, independently of all other edges. For fixed $β$, we show that the critical value $α_c(β)$ is non-trivial if and only if $N < β< N^2$. Furthermore, we show uniqueness of the infinite component and continuity of the percolation probability and of $α_c(β)$ as a function of $β$. This means that the phase diagram of this model is well understood.

preprint2010arXiv

Maximizing the size of the giant

We consider two classes of random graphs: $(a)$ Poissonian random graphs in which the $n$ vertices in the graph have i.i.d.\ weights distributed as $X$, where $\mathbb{E}(X) = μ$. Edges are added according to a product measure and the probability that a vertex of weight $x$ shares and edge with a vertex of weight $y$ is given by $1-e^{-xy/(μn)}$. $(b)$ A thinned configuration model in which we create a ground-graph in which the $n$ vertices have i.i.d.\ ground-degrees, distributed as $D$, with $\mathbb{E}(D) = μ$. The graph of interest is obtained by deleting edges independently with probability $1-p$. In both models the fraction of vertices in the largest connected component converges in probability to a constant $1-q$, where $q$ depends on $X$ or $D$ and $p$. We investigate for which distributions $X$ and $D$ with given $μ$ and $p$, $1-q$ is maximized. We show that in the class of Poissonian random graphs, $X$ should have all its mass at 0 and one other real, which can be explicitly determined. For the thinned configuration model $D$ should have all its mass at 0 and two subsequent positive integers.

preprint2010arXiv

Threshold behaviour and final outcome of an epidemic on a random network with household structure

This paper considers a stochastic SIR (susceptible$\to$infective$\to$removed) epidemic model in which individuals may make infectious contacts in two ways, both within `households' (which for ease of exposition are assumed to have equal size) and along the edges of a random graph describing additional social contacts. Heuristically-motivated branching process approximations are described, which lead to a threshold parameter for the model and methods for calculating the probability of a major outbreak, given few initial infectives, and the expected proportion of the population who are ultimately infected by such a major outbreak. These approximate results are shown to be exact as the number of households tends to infinity by proving associated limit theorems. Moreover, simulation studies indicate that these asymptotic results provide good approximations for modestly sized finite populations. The extension to unequal sized households is discussed briefly.