Researcher profile

Pieter Trapman

Pieter Trapman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 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.

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&#39;)$ 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&#39;$ 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

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.