## ON EXPECTED DURATIONS OF BIRTH-DEATH PROCESSES, WITH APPLICATIONS TO BRANCHING PROCESSES AND SIS EPIDEMICS

FRANK BALL, ∗ University of Nottingham

TOM BRITTON, ∗∗ Stockholm University

PETER NEAL, ∗∗∗ Lancaster University

∗ Postal address: The Mathematical Sciences Building, University Park, Nottingham, NG7 2RD, UK ∗∗ Postal address: Department of Mathematics, Stockholm University, SE-106 91 Stockholm, Sweden. ∗∗∗ Postal address: Department of Mathematics and Statistics, Fylde College, Lancaster University,

LA1 4YF, United Kingdom

## Abstract

We study continuous-time birth-death type processes, where individuals have independent and identically distributed lifetimes, according to a random variable Q , with E[ Q ] = 1, and where the birth rate if the population is currently in state (has size) n is α n ( ). We focus on two important examples, namely α n ( ) = λn being a branching process, and α n ( ) = λn N ( -n /N ) which corresponds to an SIS (susceptible → infective → susceptible) epidemic model in a homogeneously mixing community of fixed size N . The processes are assumed to start with a single individual, i.e. in state 1. Let T , A n , C and S denote the (random) time to extinction, the total time spent in state n , the total number of individuals ever alive and the sum of the lifetimes of all individuals in the birth-death process, respectively. The main results of the paper give expressions for the expectation of all these quantities, and shows that these expectations are insensitive to the distribution of Q . We also derive an asymptotic expression for the expected time to extinction of the SIS epidemic, but now starting at the endemic state, which is not independent of the distribution of Q . The results are also applied to the household SIS epidemic, showing that its threshold parameter R ∗ is insensitive to the distribution of Q , contrary to the household SIR (susceptible → infective → recovered) epidemic, for which R ∗ does depend on Q .

Keywords: Birth-death process; branching processes; SIS epidemics; insensitivity results.

2010 Mathematics Subject Classification: Primary 60J80

Secondary 60G10;92D30

## 1. Introduction

Akey question for population processes of a birth-death type, for example, branching processes and epidemic processes (with infection and recovery corresponding to birth and death, respectively), is what effect does the lifetime distribution have on key quantities of scientific interest? For example, consider a single-type branching process, where individuals have independent and identically distributed ( iid ) lifetimes according to a random variable Q having an arbitrary, but specified, distribution and, whilst alive, give birth at the points of a homogeneous Poisson point process with rate λ . The basic reproduction number, R 0 = λ E[ Q ], the mean number of offspring produced by an

individual during its lifetime, only depends upon Q through its mean E[ Q ]. The mean total size of a subcritical branching process ( R &lt; 0 1) with one ancestor is 1 / (1 -R 0 ), which is again independent of the distribution of Q . However, other quantities of interest such as the probability of extinction and the Malthusian parameter of the branching process depend upon the distributional form of Q . Thus, in the language of stochastic networks, R 0 can be viewed as an insensitivity result in that it depends on Q only through its mean, see, for example, Zachary (2007). Without loss of generality throughout the paper we assume that E[ Q ] = 1.

Insensitivity results for stochastic networks are well known, see for example, Sevast'yanov (1957), Whittle (1985) and Zachary (2007). In particular, in Zachary (2007), Theorem 1, it is shown that for a wide class of queueing networks, where arrivals (births) into the system are Poissonian with rate depending upon the total number of individuals in the system and each arrival has an iid workload, the stationary distribution of the total number of individuals in the system is insensitive to the distribution of Q . It then follows automatically that, for example, the mean duration of a busy period of the network (at least one individual in the system) is insensitive to the distribution of Q .

Given the similarities between queueing networks and birth-death type models, arrivals equating to births and workload equating to lifetime, we seek in this paper to explore insensitivity results for birth-death type processes with particular emphasis upon branching processes and SIS (susceptible → infective → susceptible) epidemic models. In many cases, Zachary (2007), Theorem 1, cannot be applied directly to birth-death processes, as many birth-death processes do not exhibit stationary behaviour. For example, a branching process will either go extinct or grow exponentially. However, we can exploit Zachary (2007), Theorem 1, for birth-death type processes whose mean time to extinction is finite by introducing a regeneration step ( cf. Hern´ andez-Su´rez and Castillo-Chavez (1999)) whenever the population goes exa tinct. That is, whenever the population goes extinct, it spends an exponential length of time in state 0 (no individuals) before a new individual is introduced into the population (regeneration). The birth-death type process with regenerations then fits into the framework of Zachary (2007), provided that the birth rate is Poissonian and depends upon the population only through its size. Insensitivity results are then easy to

obtain for the regenerative process, and also for the original birth-death type process. In particular, we obtain that the mean duration of the birth-death type process is insensitive to the distribution of Q .

The paper is structured as follows. In Section 2, we formally introduce the generic birth-death type process with arbitrary birth rate α n ( ), where n denotes population size, and introduce regeneration. We identify key insensitivity results for birth-death type processes, namely, that the mean duration, the mean time with n individuals alive ( n = 1 2 , , . . . ) and the mean total number of individuals ever alive in the process are insensitive to the distribution of Q . In Section 3, we focus on three special cases of the birth-death type process, namely, branching processes with constant birth rate, and homogeneously mixing and household SIS epidemic models. In Section 3.1, we prove a conjecture of Neal (2014), that for a subcritical branching process, the mean time with n ( n = 1 2 , , . . . ) individuals alive is insensitive to Q and, using Lambert (2011), Lemma 3.1, give a corresponding insensitivity result for critical and supercritical branching processes. In Section 3.2, we apply the insensitivity results to homogeneously mixing SIS epidemics and obtain a simple approximation for the mean duration of the epidemic starting from a single infective. Moreover, we show that for a supercritical epidemic ( R 0 &gt; 1), the mean duration of the epidemic starting from the quasi-endemic equilibrium does depend upon the distribution of Q and we give a simple asymptotic expression for this quantity. Finally, in Section 3.3 we exploit the results obtained for the homogeneously mixing SIS epidemic to show that both the threshold parameter R ∗ and the quasi-endemic equilibrium of the household SIS epidemic are insensitive to the distribution of Q . These are interesting findings, as in the household SIR (susceptible → infective → recovered) epidemic both R ∗ and the fraction of the population ultimately recovered if the epidemic takes off do depend upon the distribution of Q .

## 2. Generic model

The generic birth-death type process is defined as follows. The process is initiated at time t = 0 with one individual. All individuals, including the initial individual, have iid lifetimes according to an arbitrary, but specified, positive random variable Q with finite

mean. At the end of its lifetime an individual dies and is removed from the population. New individuals are born and enter the population at the points of an independent inhomogeneous Poisson point process with rate α n ( ) ≥ 0, where n denotes the total number of individuals in the population. Without loss of generality, we assume that E[ Q ] = 1, since otherwise we can simply rescale time by dividing Q and multiplying α n ( ) by E[ Q ]. The special cases of a branching process with individuals giving birth at the points of a homogeneous Poisson point process with rate λ and the homogeneously mixing SIS epidemic (see, for example, Kryscio and Lef`vre (1989)) in a population e of size N with infection rate λ , correspond to α n ( ) = nλ and α n ( ) = nλ N ( -n /N ) , respectively.

The birth-death type process is similar to the single-class networks studied in Zachary (2007), Section 2. In Zachary (2007), it is assumed that new individuals enter the system (births) at the points of a Poisson process with state-dependent rate α n ( ), where n is the total number of individuals currently in the system. Individuals have independent and identically distributed workloads, according to a random variable Q with E[ Q ] = 1. While there are n individuals in the system, the total workload is reduced at rate β n ( ) ≥ 0, with β n ( ) &gt; 0 if and only if n &gt; 0. In a biological setting, where the workload Q associated with an individual is its lifetime, it only makes sense to take β n ( ) = n , so each individual's remaining lifetime decreases at a constant rate 1. In queueing terminology this corresponds to an infinite server queue.

In Zachary (2007), Theorem 1, it is shown that if the proper distribution π = ( π (0) , π (1) , . . . ) satisfies the detailed balance equations

$$\pi ( n + 1 ) \beta ( n + 1 ) = \pi ( n ) \alpha ( n ) \quad n = 0, 1, \dots,, \quad \quad ( 2. 1 )$$

and

$$\sum _ { n = 0 } ^ { \infty } \pi ( n ) \alpha ( n ) & < \infty, \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon$$

then π is the stationary distribution of the size of the system, irrespective of the distribution of Q . That is, if Z t denotes the number of individuals in the system at time t , then for n = 0 1 , , . . .

$$P ( Z _ { t } = n ) \to \pi ( n ) \quad \text{as } t \to \infty.$$

Note that the total number of individuals in the system is not a Markov process, unless Q has an exponential distribution.

For many biological systems, Zachary (2007) does not apply since often α (0) = 0. That is, the population can go extinct and then remains extinct forever. This is the case for branching processes and the homogeneously mixing SIS epidemic model. Moreover, for a branching process if λ ≤ 1 (subcritical/critical), the branching process goes extinct with probability 1, whereas for a supercritical branching process λ &gt; 1, the branching process either goes extinct or grows unboundedly. In either case, a stationary distribution for the total number of individuals alive does not exist. The solution to make Zachary (2007), Theorem 1, relevant to birth-death type processes is to follow Hern´ndez-Su´rez and Castillo-Chavez (1999) and consider a birth-death a a type process with regeneration. We introduce regeneration by setting α (0) = 1, leaving all other transition rates unchanged. This corresponds to the process, if it goes extinct, spending an exponentially distributed time, having mean 1, with no individual before a new individual enters the population leading to the process restarting (regeneration).

The key questions are, how to analyse the regenerative process and what does it tell us about the original birth-death type process? Firstly, if the regenerative process satisfies the detailed balance equation (2.1) and is non-explosive (satisfies condition (2.2)), then it has a stationary distribution π , which is insensitive to the distribution of Q . Secondly, the behaviour of the process between regenerations are independent and identically distributed copies of the original birth-death process. Therefore this gives us a way to explore characteristics of the original process, and of identifying quantities which are insensitive to the distribution of Q and also, later in Section 3.2, some which do depend upon the distribution of Q . Recall that in the processes we study β n ( ) = n . Then (2.1) implies that

$$\pi ( n ) \alpha ( n ) \ = \ ( n + 1 ) \pi ( n + 1 ) \quad \text{for } n = 0, 1, \dots. \quad \text{(2.4)}$$

Therefore, it follows that

$$\pi ( n ) \ = \ \pi ( 0 ) \prod _ { i = 1 } ^ { n - 1 } \frac { \alpha ( i ) } { i + 1 }, \quad n = 1, 2, \dots,$$

where the product is 1 when n = 1. It then follows from (2.4) and (2.5) that

$$\pi ( 0 ) \ = \ \left \{ 1 + \sum _ { n = 1 } ^ { \infty } \prod _ { i = 0 } ^ { n - 1 } \frac { \alpha ( i ) } { i + 1 } \right \} ^ { - 1 }. \quad \quad \ ( 2. 6 ) \\. \quad \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

Note that π being a proper distribution implicitly implies that π (0) &gt; 0 or, equivalently, that the sum in (2.6) is finite, and hence that the process is positive recurrent. Thus the regenerative process is not suitable for critical or supercritical branching processes. We discuss this in more detail in Section 3.1 below.

We complete this section by identifying a number of key quantities whose means are insensitive to the distribution of Q and are summarised in Theorem 2.1.

Theorem 2.1. Let Y t denote the total number of individuals in the birth-death process at time t . Let T = ∫ ∞ 0 1 { Y &gt; t 0 } dt denote the duration of the birth-death process and, for n = 1 2 , , . . . , let A n = ∫ ∞ 0 1 { Y t = n } dt denote the total time the birth-death process spends with n individuals alive. Then

$$\mathbf E [ T ] = \sum _ { n = 1 } ^ { \infty } \prod _ { i = 0 } ^ { n - 1 } \frac { \alpha ( i ) } { i + 1 },$$

and for n = 1 2 , , . . . ,

$$\mathbf E [ A _ { n } ] \ = \ \prod _ { i = 0 } ^ { n - 1 } \frac { \alpha ( i ) } { i + 1 }. \quad \quad \ \ ( 2. 8 ) \\ \cdot \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

Finally, let C and S denote the total number of individuals ever alive in the population and the sum of the lifetimes of all individuals ever alive in the birth-death process, respectively. Then E[ C ] = E[ S ] with

$$E [ S ] \ = \ \sum _ { k = 1 } ^ { \infty } k \prod _ { i = 1 } ^ { k - 1 } \frac { \alpha ( i ) } { i + 1 }. \quad \quad \ \ ( 2. 9 ) \\.. \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Proof. An immediate consequence of the above construction is that the mean time between regenerations is 1 /π (0), irrespective of the distribution of Q . On average one unit of time is spent with no individual in the population, so (see Ball and Milne (2004) for a formal justification)

$$E [ T ] = \frac { 1 } { \pi ( 0 ) } - 1 = \left \{ \sum _ { n = 1 } ^ { \infty } \prod _ { i = 0 } ^ { n - 1 } \frac { \alpha ( i ) } { i + 1 } \right \},$$

as required. Moreover, for n = 1 2 , , . . . ,

$$E [ A _ { n } ] \ = \ \pi ( n ) \frac { 1 } { \pi ( 0 ) } = \prod _ { i = 0 } ^ { n - 1 } \frac { \alpha ( i ) } { i + 1 },$$

using (2.5).

Using Fubini's theorem,

$$\text{bini} ^ { \text{s theorem} } \\ E [ S ] \ & = \ \text{E} \left [ \int _ { 0 } ^ { \infty } Y _ { t } \, d t \right ] \\ & = \ \text{E} \left [ \int _ { 0 } ^ { \infty } \sum _ { k = 1 } ^ { \infty } k 1 _ { \{ Y _ { t } = k \} } \, d t \right ] \\ & = \ \sum _ { k = 1 } ^ { \infty } k E \left [ \int _ { 0 } ^ { \infty } 1 _ { \{ Y _ { t } = k \} } \, d t \right ] \\ & = \ \sum _ { k = 1 } ^ { \infty } k E [ A _ { k } ] = \sum _ { k = 1 } ^ { \infty } k \prod _ { i = 1 } ^ { k - 1 } \frac { \alpha ( i ) } { i + 1 }. \quad & ( 2. 1 2 ) \\. \text{fter using} \ ( 2. 1 ) \, \text{with} \, \beta ( n ) = n, \, ( 2. 2 ) \, \text{ensures that} \, E [ S ] \, \text{is finite}.$$

Note that, after using (2.1) with β n ( ) = n , (2.2) ensures that E[ S ] is finite.

Given A k , the mean number of births whilst the process is in state k is α k A ( ) k . Therefore, including the initial ancestor and noting from (2.11) that E[ A 1 ] = 1, we have that

$$E [ C ] \ & = \ \ E \left [ 1 + \sum _ { k = 1 } ^ { \infty } \alpha ( k ) A _ { k } \right ] \\ & = \ \ 1 + \sum _ { k = 1 } ^ { \infty } \alpha ( k ) \frac { \pi ( k ) } { \pi ( 0 ) } \\ & = \ \ E [ A _ { 1 } ] + \sum _ { k = 1 } ^ { \infty } ( k + 1 ) \frac { \pi ( k + 1 ) } { \pi ( 0 ) } \\ & = \ \ E [ A _ { 1 } ] + \sum _ { k = 2 } ^ { \infty } k E [ A _ { k } ] \\ & = \ \ E [ S ].$$

## 3. Special cases

In this section we apply the results obtained in Section 2 to three special cases, namely branching processes, homogeneously mixing SIS epidemics and household SIS epidemics, yielding fresh insight into these models.

## 3.1. Branching process

As mentioned above, we consider branching processes where individuals have iid lifetimes according to Q (with E Q [ ] = 1) and whilst alive give birth at the points of a homogeneous Poisson point process with rate λ . Therefore we have that α n ( ) = nλ . The key result is Lemma 3.1 which is a generalization of Neal (2014), Conjecture 2.1.

Lemma 3.1. For n = 1 2 , , . . . ,

$$E [ A _ { n } ] & = \frac { \lambda ^ { n - 1 } } { n ( \max \{ 1, \lambda \} ) ^ { n } }. & ( 3. 1 ) \\ a \colon = ( \tilde { n } \circ \widetilde { n } \circ \widetilde { n } \circ \tilde { n } \dots.$$

The Lemma is proved in (3.6) and (3.9) below.

For λ &lt; 1, the branching process is subcritical and the results of Section 2 hold. The form of α n ( ) allows explicit expressions to be obtained. Firstly, it follows from (2.5) that

$$\pi ( n ) \ = \ \pi ( 0 ) \prod _ { i = 1 } ^ { n - 1 } \frac { i \lambda } { i + 1 } = \pi ( 0 ) \frac { \lambda ^ { n - 1 } } { n }, \quad n = 1, 2, \dots. \quad \quad ( 3. 2 ) \intertext { \text{hen follows immediately that}$$

It then follows immediately that

$$\pi ( 0 ) \ & = \ \left \{ 1 + \sum _ { n = 1 } ^ { \infty } \frac { \lambda ^ { n - 1 } } { n } \right \} ^ { - 1 } \\ & = \ \left \{ 1 - \log ( 1 - \lambda ) / \lambda \right \} ^ { - 1 }.$$

$$= \ \left \{ 1 - \log ( 1 - \lambda ) / \lambda \right \} ^ { - 1 }.$$

(Note that the sum in (3.3) diverges if λ ≥ 1.) Therefore, it follows from (3.4) and (2.7) that the mean duration of the branching process is

$$E [ T ] = - \log ( 1 - \lambda ) / \lambda,$$

a result which is well known for the linear birth-death process with Q ∼ Exp(1). Also Neal (2014), Conjecture 2.1, is proved in that, for n = 1 2 , , . . . ,

$$E [ A _ { n } ] \ = \ \frac { \lambda ^ { n - 1 } } { n }.$$

Finally, we obtain the classical result that the mean total number of individuals ever alive in the branching process is

$$E [ C ] \ = \ \sum _ { k = 1 } ^ { \infty } k E [ A _ { k } ] = \sum _ { k = 1 } ^ { \infty } \lambda ^ { k - 1 } = \frac { 1 } { 1 - \lambda }.$$

The approach taken in Section 2 is valid only in the above scenario. A key component of the proofs is that the birth rate depends only upon the total number of individuals in the population. Therefore the above insensitivity results do not hold for more general reproductive life histories. However, progress can be made in extending Neal (2014), Conjecture 2.1, and (3.6) to the critical and supercritical cases ( λ ≥ 1) by using Lambert (2011), Lemma 3.1. (Note that the mean duration and mean total number of individuals ever alive in the branching process is infinite in these cases.) Specifically, Lambert (2011), Lemma 3.1, shows that, for n = 1 2 , , . . . and t ≥ 0,

$$\Phi ( Y _ { t } = n ) = \left ( 1 - \frac { 1 } { W ( t ) } \right ) ^ { n - 1 } \frac { W ^ { \prime } ( t ) } { \lambda W ( t ) ^ { 2 } }, \quad \quad \quad ( 3. 8 ) \\ \therefore \quad. \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots$$

where W t ( ) solves Lambert (2011), equation (6). For most choices of Q , it is not possible to get an explicit expression for W t ( ) for all t ≥ 0. However, for any Q (E[ Q ] = 1) W (0) = 1 and for λ ≥ 1, W ( ∞ ) = ∞ . Therefore

$$1 ) \, W ( 0 ) = 1 \, \text{ and for } \lambda \geq 1, \, W ( \infty ) = \infty. \, \text{Therefore} \\ E [ A _ { n } ] \ & = \ \int _ { 0 } ^ { \infty } \, P ( Y _ { t } = n ) \, d t \\ & = \ \frac { 1 } { \lambda } \int _ { 0 } ^ { \infty } \left ( 1 - \frac { 1 } { W ( t ) } \right ) ^ { n - 1 } \frac { W ^ { \prime } ( t ) } { W ( t ) ^ { 2 } } \, d t \\ & = \ \frac { 1 } { \lambda } \left [ \frac { 1 } { n } \left ( 1 - \frac { 1 } { W ( t ) } \right ) ^ { n } \right ] _ { 0 } ^ { \infty } = \frac { 1 } { \lambda n }. \\ \text{I be noted that in the subcritical case } ( \lambda < 1 ) \, W ( \infty ) = 1 / ( 1 - \lambda ) \text{ and } ( 3. 9 )$$

It should be noted that in the subcritical case ( λ &lt; 1) W ( ∞ ) = 1 / (1 -λ ) and (3.9) can be used to obtain (3.6) directly.

## 3.2. SIS epidemic

3.2.1. Mean duration with one initial infective For a homogeneously mixing SIS epidemic in a population of size N , the rate at which new infections (births) occur, given that there are currently n infectives is α n ( ) = λn N ( -n /N ) , where λ is the rate at which a typical infective makes infectious contacts. Infectious contacts are assumed to be with individuals chosen independently and uniformly at random from the entire population. Thus, if there are y infectives, and hence N -y susceptibles, the probability that an infectious contact is with a suscpetible, and hence results in a new infective, is ( N -y /N ) . The infectious periods of infectives are independently and identically distributed according to Q and an infective becomes susceptible again as soon as its infectious period ends. Using (2.5), it is straightforward to show that,

for n = 1 2 , , . . . , N ,

$$\pi ^ { ( N ) } ( n ) = \frac { \pi ^ { ( N ) } ( 0 ) } { n } \frac { ( N - 1 )! } { ( N - n )! } \left ( \frac { \lambda } { N } \right ) ^ { n - 1 }, \quad \quad \quad ( 3. 1 0 ) \\. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdets \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \csrf \quad \,.$$

which has previously been obtained for the Markov case ( Q ∼ Exp(1)) by Hern´ndez-Su´rez and Castillo-Chavez a a (1999). Consequently, the mean duration of the epidemic starting from a single infective is

$$\mathrm E \left [ T ^ { ( N ) } \right ] \ = \ \sum _ { n = 1 } ^ { N } \frac { ( N - 1 )! } { n ( N - n )! } \left ( \frac { \lambda } { N } \right ) ^ { n - 1 }. \quad \quad \ \ ( 3. 1 1 ) \\ \dots \, \hat { \ }. \, \dots \, \hat { \ }. \, \hat { \ }. \, \hat { \ }. \, \dots \, \hat { \ }. \, \dots \, \hat { \ }. \, \dots \, \hat { \ }. \, \dots$$

Note that this follows directly from (2.7), on recalling that there α (0) = 1. We index statistics of interest by the total population size N to highlight the role played by N in the analysis below.

It is interesting to investigate the behaviour of E [ T ( N ) ] for large N . We summarise the results in Lemma 3.2.

## Lemma 3.2. For the subcritical case, λ &lt; 1 ,

$$E \left [ T ^ { ( N ) } \right ] \to - \frac { \log ( 1 - \lambda ) } { \lambda } \quad a s \ N \to \infty. \quad \quad ( 3. 1 2 ) \\ \dots \. \ \. \ \.$$

For the critical case, λ = 1 ,

$$E \left [ T ^ { ( N ) } \right ] & \sim \frac { 1 } { 2 } \log N, \\. \quad. \quad. \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \
 \,. \quad. \quad. \quad. \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,  \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \   \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad\, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad\ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \$$

where a ( N ) ∼ b ( N ) denotes that lim N →∞ a ( N ) /b ( N ) = 1 .

For the supercritical case, λ &gt; 1 ,

$$\dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

We give a heuristic proof of Lemma 3.2 here, deferring a formal proof to Appendix A.

For λ &lt; 1, we have that π ( N ) ( n ) ≈ π ( N ) (1) λ n -1 /n , a branching process approximation ( c.f. Whittle (1955), Ball and Donnelly (1995)) and E [ T ( N ) ] ≈ -log(1 -λ /λ ) ( c.f. Section 3.1).

For λ = 1, it is fruitful to recall the birthday problem. Balls are drawn uniformly at random with replacement from an urn containing N balls numbered 1 2 , , . . . , N . Let

M ( N ) be the number of draws required until a ball is drawn that had previously been drawn. Then (3.11) yields, with products being one if vacuous,

$$\begin{array} { l l l } \left ( 3. 1 1 \right ) \text{ yields, with products being one if vacuous,} \\ \\ E \left [ T ^ { ( N ) } \right ] & = & \sum _ { n = 1 } ^ { N } \frac { 1 } { n } \prod _ { i = 1 } ^ { n - 1 } \left ( 1 - \frac { i } { N } \right ) \\ & = & \sum _ { n = 1 } ^ { N } \frac { 1 } { n } P \left ( M ^ { ( N ) } > n \right ) \\ & = & E \left [ \sum _ { n = 1 } ^ { M ^ { ( N ) } } \frac { 1 } { n } \right ]. \\ \\ \text{ple, Aldous } (1985), page 96, \end{array}$$

Now, for example, Aldous (1985), page 96,

$$N ^ { - \frac { 1 } { 2 } } M ^ { ( N ) } \stackrel { D } { \longrightarrow } M \quad \text{as} \quad N \to \infty,$$

where D -→ denotes convergence in distribution and M is a random variable having probability density function f M ( x ) = x exp ( -x 2 2 ) ( x &gt; 0). Thus, when N is large, M ( N ) is concentrated on n -values of the form xN 1 2 , for which

$$\sum _ { i = 1 } ^ { n } \frac { 1 } { i } \approx \log \left ( x N ^ { \frac { 1 } { 2 } } \right ) = \log x + \frac { 1 } { 2 } \log N$$

and (3.13) follows.

For λ &gt; 1, rearranging (3.11) we get

$$E \left [ T ^ { ( N ) } \right ] = ( N - 1 )! ( \lambda / N ) ^ { N - 1 } \sum _ { j = 0 } ^ { N - 1 } \frac { ( N / \lambda ) ^ { j } } { ( N - j ) j! }.$$

Let us first consider the sum in (3.17). If we multiply by e -N/λ we are adding Poisson( N/λ ) probabilities multiplied by 1 / N ( -j ). When N is large the probability mass is concentrated on j -values near the mean, j ≈ N/λ , by the law of large numbers. Owing to the assumption λ &gt; 1, these terms are smaller than N and hence contained in the finite sum. From this it follows that the sum is asymptotically equivalent to

$$\sum _ { j = 0 } ^ { N - 1 } \frac { ( N / \lambda ) ^ { j } } { ( N - j ) j! } \ & = \ e ^ { N / \lambda } \sum _ { j = 0 } ^ { N - 1 } \frac { ( N / \lambda ) ^ { j } e ^ { - N / \lambda } } { ( N - j ) j! } \\ & \approx \ \frac { e ^ { N / \lambda } } { N ( 1 - 1 / \lambda ) } \sum _ { j = 0 } ^ { N - 1 } \frac { e ^ { - N / \lambda } ( N / \lambda ) ^ { j } } { j! } \approx \frac { e ^ { N / \lambda } } { N ( 1 - 1 / \lambda ) }. \quad ( 3. 1 8 )$$

The second approximation we use is Stirling's formula implying that ( N -1)! ∼ √ 2 π N ( -1)( N -1) ( N -1) e -( N -1) . Combining these two approximations yields

$$\L e _ { \pi ( \lambda \nu - 1 ) } \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1 ) \L ( \lambda \nu - 1. \ \end{cases} \L ( \Lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \L ( \lambda \nu \nu - 1 ) \\ \vdots \text{required.}$$

as required.

Let A ( N ) n denote the total amount of time that the SIS epidemic, initiated with a single infective, spends with n infectious individuals. Then, from (2.11) and (2.10),

$$E \left [ A _ { n } ^ { ( N ) } \right ] = \frac { \pi ^ { ( N ) } ( n ) } { \pi ^ { ( N ) } ( 0 ) } = \frac { \pi ^ { ( N ) } ( n ) } { 1 - \pi ^ { ( N ) } ( 0 ) } E \left [ T ^ { ( N ) } \right ]. \quad \text{(3.20)} \\ \intertext { 1 \times \text{the fact condition in } } \text{$(2.90)$ and $(2.10)$ to mean total number of infotimes}$$

Using (2.13), the first equation in (3.20) and (3.10), the mean total number of infectives during the course of a supercritical epidemic is

$$\text{uring} \, \text{the course or a supercritica} \, \text{epiaemic is} \\ \text{E} \left [ C ^ { ( N ) } \right ] = \sum _ { n = 1 } ^ { N } n E \left [ A _ { n } ^ { ( N ) } \right ] \ & = \ \sum _ { n = 1 } ^ { N } n \frac { \pi ^ { ( N ) } ( n ) } { \pi ^ { ( N ) } ( 0 ) } \\ & = \ \sum _ { n = 1 } ^ { N } \frac { ( N - 1 )! } { ( N - n )! } \left ( \frac { \lambda } { N } \right ) ^ { n - 1 } \quad & ( 3. 2 1 ) \\ & \sim \ \frac { \sqrt { 2 \pi } } { \lambda } \sqrt { N } \exp \left ( \{ \log \lambda - 1 + 1 / \lambda \} N \right ). \quad & ( 3. 2 2 ) \\ \text{The derivation of} \, ( 3. 2 2 ) \, \text{is similar to but simnler than that of} \, E \, \lceil T ^ { ( N ) } \rceil.$$

The derivation of (3.22) is similar to but simpler than that of E [ T ( N ) ] . Note that the second equation in (3.20) gives

$$\mathbf E \left [ C ^ { ( N ) } \right ] & = \left ( \sum _ { n = 1 } ^ { N } \frac { n \pi ^ { ( N ) } ( n ) } { 1 - \pi ^ { ( N ) } ( 0 ) } \right ) \mathbf E \left [ T ^ { ( N ) } \right ]. \\ \cdot \, \tau _ { N } \cup \quad \dots \tau _ { N } \cup \cdot \tau _ { N } \cup \cdot \quad \cdot \tau _ { N } \cup \cdot \quad \cdot \tau _ { N } \cup \cdot \quad \cdot$$

The distribution ˜ π ( N ) = (˜ π ( N ) 1 , π ˜ ( N ) 2 , . . . , π ˜ ( N ) N ), where ˜ π ( N ) n = π ( N ) ( n / ) (1 -π ( N ) (0)), gives a 'quasi-equilibrium' distribution for the SIS epidemic. Thus, the mean total number of infectives in the epidemic is given by the mean number of infectives in quasiequilibrium multiplied by the mean duration of the epidemic. When the epidemic is supercritical ( λ &gt; 1), the distribution of ˜ π ( N ) is concentrated on values close to N (1 -λ -1 ), as indicated previously, which explains the simple multiplicative relationship between the approximations (3.19) and (3.22).

3.2.2. Mean extinction time from quasi-endemic equilibrium The above calculations of E [ T ( N ) ] are insensitive to the distribution of Q . However, for supercritical SIS epidemics there is interest in the time to extinction of the epidemic starting from the quasi-endemic equilibrium of N (1 -1 /λ ) infectives. We outline how the mean time to extinction from the quasi-endemic equilibrium, E [ T ( N ) Q ] , does depend upon the distribution of Q . The epidemic initiated from a single infective either goes extinct very quickly or takes-off and reaches an endemic equilibrium of a proportion ( λ -1) /λ of the population infected, see Kryscio and Lef`vre (1989). e The epidemic then spends a long time fluctuating about the endemic equilibrium before eventually going extinct. This can be seen from π ( N ) , with most of the probability mass centred about ( λ -1) N/λ infectives. There has been considerable interest in investigating the distribution of the time to extinction from the endemic equilibrium, see for example Kryscio and Lef` evre (1989), Andersson and Djehiche (1998), N˚ asell (1999) and Britton and Neal (2010). This is a difficult problem on which to make analytical progress. In Andersson and Djehiche (1998), it was shown that, in the limit as N → ∞ , for Q ∼ Exp(1) the time to extinction divided by E [ T ( N ) Q ] converges in distribution to an exponential random variable with mean 1. Moreover, E [ T ( N ) Q ] ∼ √ 2 π/Nλ exp( N { log λ +1 /λ - } 1 ) / λ ( -1) 2 = µ ( N ) / (1 -1 /λ ), where µ ( N ) = E [ T ( N ) ] . It is conjectured that an exponential distribution for the time to extinction holds more generally than for Q ∼ Exp(1), but even computing E [ T ( N ) Q ] up to leading terms in N has proved difficult. By studying Gaussian approximations for the endemic equilibrium qualitative results on the time to extinction have been obtained, see N˚sell a (1999) and Britton and Neal (2010). Whilst, such approaches have given a qualitative understanding of extinction of SIS epidemics, the estimates obtained for the mean time to extinction are incorrect by orders of magnitude. Moreover, it is noted in Neal (2014) that simulation results suggest that the distribution of Q does affect the mean time to extinction from the quasi-endemic equilibrium, which is not predicted by using the qualitative Gaussian approximation.

Let p Q denote the extinction probability of a branching process, in which individuals have iid infectious periods according to Q and whilst alive, give birth at the points of a homogeneous Poisson point process with rate λ . We show in Lemma 3.3 that E [ T ( N ) Q ] depends upon the distribution Q through p Q .

## Lemma 3.3. For λ &gt; 1 and var Q ( ) &lt; ∞ ,

$$E \left [ T _ { Q } ^ { ( N ) } \right ] & \sim \frac { 1 } { 1 - p _ { Q } } \mu ^ { ( N ) }. \\ \intertext { f o f f o m m a } 3 \, 3 \, \text{ in } \Delta \tt n o n d i v \ R$$

We give a formal proof of Lemma 3.3 in Appendix B, where we also show that the mean extinction time of the SIS epidemic ∼ 1 1 -p Q µ ( N ) whenever it starts with a strictly positive fraction of the population infected. Here we present a heuristic proof of Lemma 3.3. The requirement that var Q ( ) &lt; ∞ is almost certainly not necessary but is assumed in the proof below. As noted above, the supercritical SIS epidemic will either quickly go extinct or will take-off and reach the endemic equilibrium. Let 1 -P ( N ) Q denote the probability that the total number of infectives in the epidemic equals /floorleft /epsilon1N /floorright at some point in time, for some 0 &lt; /epsilon1 &lt; ( λ -1) /λ . (Throughout the paper, /floorleft x /floorright denotes the greatest integer ≤ x and /ceilingleft x /ceilingright denotes the smallest integer ≥ x .) Then it is straightforward to show, using a branching process approximation (see, for example, Whittle (1955) and Ball and Donnelly (1995)), that P ( N ) Q → p Q as N →∞ . Then the mean duration of an epidemic, initiated from a single infective satisfies

$$\mu ^ { ( N ) } & = \frac { 1 } { \pi ^ { ( N ) } ( 1 ) } \ = \ P _ { Q } ^ { ( N ) } A _ { Q } ^ { ( N ) } + ( 1 - P _ { Q } ^ { ( N ) } ) \left \{ B _ { Q } ^ { ( N ) } + \mathrm E \left [ T _ { Q } ^ { ( N ) } \right ] \right \}, \quad ( 3. 2 3 ) \\, \quad \cdot ( N ) \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

where A ( N ) Q is the mean duration of an epidemic which never reaches /floorleft /epsilon1N /floorright infected (epidemic dies off quickly) and B ( N ) Q is the mean time take for the epidemic to reach the endemic equilibrium given it reaches /floorleft /epsilon1N /floorright infected. The definition of B ( N ) Q is imprecise and correspondingly we take B ( N ) Q to be the mean time to reach [( λ -1) N/λ ] infectives given that the epidemic takes off. For the case Q ∼ Exp(1), it is shown in Andersson and Djehiche (1998) that A ( N ) Q = O (1) and B ( N ) Q = O (log N ). Therefore assuming that for general Q A , ( N ) Q , B ( N ) Q = ( o µ ( N ) ), we have that

$$\mathbf E \left [ T _ { Q } ^ { ( N ) } \right ] \ = \ \frac { 1 } { 1 - P _ { Q } ^ { ( N ) } } \left \{ \frac { 1 } { \pi ^ { ( N ) } ( 1 ) } - P _ { Q } ^ { ( N ) } A _ { Q } ^ { ( N ) } - ( 1 - P _ { Q } ^ { ( N ) } ) B _ { Q } ^ { ( N ) } \right \} \\ \approx \ \frac { 1 } { 1 - p _ { Q } } \times \frac { 1 } { \pi ^ { ( N ) } ( 1 ) } = \frac { 1 } { 1 - p _ { Q } } \mu ^ { ( N ) }, \\ \mathbf T \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{cases}$$

Therefore we can immediately see the role of Q in E [ T ( N ) Q ] . Specifically, the greater the extinction probability p Q , the longer the epidemic will on average persist, given that it takes off and becomes established. Note that, subject to E[ Q ] = 1, the extinction probability is least when Q is constant ( i.e. P( Q = 1) = 1), so the model with a constant

infectious period has the shortest mean time to extinction starting from quasi-endemic equilibrium.

## 3.3. Household SIS epidemic

The final special case we consider is the household SIS epidemic model. Consider a fixed community consisting of m households which, for simplicity of exposition, all have the same size h , so the population size is N h = mh . Our results extend straightforwardly to the case where the household sizes are unequal. We are particularly interested in the case where m , and hence N h , is large. Infectious individuals have iid infectious periods according to Q , after which they become susceptible again. While infectious an individual makes two types of contacts: the individual makes global infectious contacts at rate λ G , each time the contacted person is selected independently and uniformly at random from the whole community, including individuals in the same household, and the individual makes local infectious contacts at rate λ L with any given individual in their household. Therefore, an individual makes infectious contacts at a total rate of λ G +( h -1) λ L . By examining the within-household dynamics of the SIS epidemic in the initial stages of the epidemic and at the quasi-endemic equilibrium, we obtain interesting, and perhaps unexpected, insensitivity results for the household SIS epidemic model.

For large m , the initial stages of the household SIS epidemic can be approximated by a branching process; see Ball (1999), where the approximation is made fully rigorous using a coupling argument. The branching process approximation is similar to that used for the household SIR epidemic, Ball et al. (1997), with individuals in the approximating branching process corresponding to within-household epidemic outbreaks in the epidemic. For large m , in the initial stages of the household SIS epidemic the probability that a global infectious contact is with an infectious household (a household containing at least one infective) is very small. Therefore, we assume that all global infectious contacts are with totally susceptible households and we consider the epidemic within a household, ignoring for the moment global infectious contacts, initiated by a single infective and without any additional global infections from outside.

Let S denote the total severity of such a within-household epidemic, where the severity is the sum of the infectious periods of all infectives during the course of the epi-

demic from the initial infective until the epidemic within the household ceases. Then, conditional upon S , the total number of global infectious contacts emanating from the household has a Poisson distribution with mean λ G S , so the basic reproduction number of the approximating branching process is R ∗ = λ G E[ S ]. The household SIS epidemic is said to be subcritical, critical or supercritical if R &lt; ∗ 1, R ∗ = 1 or R &gt; ∗ 1, respectively. The above expression for R ∗ holds also for the household SIR epidemic, where it is known that E[ S ], the mean severity of the within-household epidemic, depends upon the distribution of Q , as does both the size of a major outbreak and the distribution of the ultimate number of susceptibles in a typical household in the event of a major outbreak. The following lemma shows that all of the corresponding quantities for the households SIS epidemic are insensitive to the distribution of Q .

Lemma 3.4. For any Q , satisfying E[ Q ] = 1 ,

$$R _ { * } \ = \ \lambda _ { G } \sum _ { n = 1 } ^ { h } \frac { ( h - 1 )! } { ( h - n )! } \lambda _ { L } ^ { n - 1 }, \quad \quad ( 3. 2 5 ) \\.. \quad \widehat { \lambda }.$$

(cf. Ball (1999), equation (8)).

For R ∗ &gt; 1 , in the limit as m → ∞ , there exists an endemic equilibrium with a proportion z of the population infected, where z is the unique non-zero solution of

$$s = \sum _ { i = 0 } ^ { h } i \phi _ { i } ( s ),$$

with

$$\phi _ { 0 } ( s ) = \left \{ 1 + \sum _ { k = 1 } ^ { h } \prod _ { i = 1 } ^ { k } \frac { ( h + 1 - i ) ( \lambda _ { G } s + ( i - 1 ) \lambda _ { L } ) } { i } \right \} ^ { - 1 } \quad \quad ( 3. 2 7 ) \\. \. \. \. \. \.$$

and, for 1 ≤ k ≤ h ,

$$\phi _ { k } ( s ) = \frac { ( h + 1 - k ) ( \lambda _ { G } s + ( k - 1 ) \lambda _ { L } ) } { k } \phi _ { k - 1 } ( s ).$$

Further, in the limit as m → ∞ , at the endemic equilibrium, the distribution of the number of infectives in a typical household is given by ( φ 0 ( z , φ ) 1 ( z , . . . , φ ) h ( z )) .

Proof. The within-household epidemic without additional global infections is simply a homogeneously mixing SIS epidemic with N = h and λ/N = λ L , so the insensitivity

results for the homogeneously mixing SIS epidemic are applicable. Of primary interest, this means that

$$\mathbf E [ S ] \ = \ \sum _ { n = 1 } ^ { h } \frac { ( h - 1 )! } { ( h - n )! } \lambda _ { L } ^ { n - 1 },$$

(recall (2.13) and (3.21)), whence (3.25) follows, regardless of the distribution of Q . It should be noted that the distribution of S does depend upon the distribution of Q , so the probability that the epidemic takes off, corresponding to the approximating branching process not going extinct, does depend upon the distribution of Q .

We turn our attention to the quasi-endemic equilibrium in the case R ∗ &gt; 1. For the Markov case, it is shown in Neal (2006), Section 4, that, in the limit as m → ∞ , there exists a stable endemic equilibrium satisfying (3.26) to (3.28). We again consider a within-household epidemic but now with a constant global force of infection from outside the household. Given that a proportion z of the population is infected, each individual in a household receives global infectious contacts at the points of a homogeneous Poisson point process with rate zλ G . Therefore letting α n ( ) = ( h -n ) { zλ G + nλ L } , new infections take place within the household at rate α n ( ) if there are currently n infectives and h -n susceptibles. Again individuals within the household have iid infectious periods distributed according to Q . The within-household epidemic satisfies the generic framework of Section 2 and it follows from (3.28) that φ satisfies the detailed balance equation (2.1). Thus φ is the stationary distribution of the withinhousehold epidemic, regardless of the distribution of Q , so at the endemic equilibrium of the household SIS epidemic, both the proportion of the population infected and the distribution of the number infected in a typical household are insensitive to the distribution of Q .

## Appendix A. Proof of Lemma 3.2

First note from (3.11) that

$$\mathbf E \left [ T ^ { ( N ) } \right ] = \sum _ { n = 1 } ^ { N } \prod _ { i = 1 } ^ { n - 1 } \left ( 1 - \frac { i } { N } \right ) \frac { \lambda ^ { n - 1 } } { n }.$$

Suppose that λ &lt; 1. Then, for any k ∈ N ,

$$\liminf _ { N \to \infty } E \left [ T ^ { ( N ) } \right ] & \geq \liminf _ { N \to \infty } \sum _ { n = 1 } ^ { k } \prod _ { i = 1 } ^ { n - 1 } \left ( 1 - \frac { i } { N } \right ) \frac { \lambda ^ { n - 1 } } { n } = \sum _ { n = 1 } ^ { k } \frac { \lambda ^ { n - 1 } } { n }, \\ \intertext { - a \i \i \dots \dots - t \i } \,.$$

and letting k →∞ yields

$$\liminf _ { N \to \infty } E \left [ T ^ { ( N ) } \right ] \geq - \frac { \log ( 1 - \lambda ) } { \lambda }. \quad \quad \quad ( A. 2 ) \\. \text{nr}$$

Also, for any k ∈ N ,

$$\lim _ { N \to \infty } & \text{E} \left [ T ^ { ( N ) } \right ] \leq \lim _ { N \to \infty } \sum _ { n = 1 } ^ { k } \prod _ { i = 1 } ^ { n - 1 } \left ( 1 - \frac { i } { N } \right ) \frac { \lambda ^ { n - 1 } } { n } + \sum _ { n = k + 1 } ^ { \infty } \lambda ^ { n - 1 } = \sum _ { n = 1 } ^ { k } \frac { \lambda ^ { n - 1 } } { n } + \frac { \lambda ^ { k } } { 1 - \lambda }, \\ & \dots \dots \cdot \cdot$$

and letting k →∞ yields

$$\lim _ { N \to \infty } E \left [ T ^ { ( N ) } \right ] \leq - \frac { \log ( 1 - \lambda ) } { \lambda }. \quad \quad \quad ( A. 3 ) \\ \dots \, 1 / 2 / 2 \, \partial \, \dots \, 1 / 2 \, \partial \, \dots$$

Combining (A.2) and (A.3) yields (3.12).

Suppose that λ = 1. Then, setting λ = 1 in (A.1) and noting that

$$\prod _ { i = 1 } ^ { n - 1 } \left ( 1 - \frac { i } { N } \right ) & \leq \prod _ { i = 1 } ^ { n - 1 } \exp \left ( - \frac { i } { N } \right ) = \exp \left ( - \frac { n ( n - 1 ) } { 2 N } \right ), \\. \, \cdots \, \ e \quad \ r \cdots \, \ e$$

yields that, for any L &gt; 0,

$$& \ell \left [ T ^ { ( N ) } \right ] \leq \sum _ { n = 1 } ^ { \prime } \frac { \ell } { n } + \sum _ { \left [ L \sqrt { N } \right ] + 1 } ^ { N } \frac { 1 } { n } \exp \left ( - \frac { n ( n - 1 ) } { 2 N } \right ) \leq \sum _ { n = 1 } ^ { \left [ L \sqrt { N } \right ] } \frac { 1 } { n } + \exp \left ( - L ^ { 2 } / 2 \right ) \sum _ { \left [ L \sqrt { N } \right ] + 1 } ^ { N } \frac { 1 } { n }. \\ & \mu _ { \dots \dots }$$

Hence,

$$\lim _ { N \to \infty } \frac { E \left [ T ^ { ( N ) } \right ] } { \frac { 1 } { 2 } \log N } \leq 1 + \exp ( - L ^ { 2 } / 2 ). \quad \quad ( A. 4 ) \\ \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

Setting λ = 1 in (A.1) yields that, for any K &gt; 0,

$$& \lambda \left [ T ^ { ( N ) } \right ] \geq \sum _ { n = 1 } ^ { \left \langle \, \left \langle \, \left \langle \, \right \rangle } \frac { \cdot \, \cdot \, \cdot } { n } \prod _ { i = 1 } ^ { \cdot \, \cdot \, \cdot } \left ( 1 - \frac { i } { N } \right ) \geq \left ( \sum _ { n = 1 } ^ { \left \langle \, \left \langle \, \left \langle \, \right \rangle } \sum _ { n } ^ { \left \langle \, \left \langle \, \right \rangle } \frac { \cdot \, \cdot } { n } \right ) \prod _ { i = 1 } ^ { \left \langle \, \left \langle \, \right \rangle } \left ( 1 - \frac { i } { N } \right ). \\ & \text{Now (3.16) implies that}$$

Now (3.16) implies that

$$\lim _ { N \to \infty } \prod _ { i = 1 } ^ { \left \lceil K \sqrt { N } \right \rceil } \left ( 1 - \frac { i } { N } \right ) = \exp ( - K ^ { 2 } / 2 ),$$

whence

$$\liminf _ { N \to \infty } \frac { \text{E} \left [ T ^ { ( N ) } \right ] } { \frac { 1 } { 2 } \log N } \geq \exp ( - K ^ { 2 } / 2 ). \quad \quad \quad ( \text{A.5} )$$

Letting L →∞ in (A.4) and K ↓ 0 in (A.5) yields (3.13).

To prove (3.14) we show that

$$\lim _ { N \to \infty } N \sum _ { j = 0 } ^ { N - 1 } \frac { ( N / \lambda ) ^ { j } \text{e} ^ { - N / \lambda } } { ( N - j ) j! } = \frac { \lambda } { \lambda - 1 }. \quad \quad \quad ( A. 6 ) \\. \quad \dots \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

Note that this makes fully rigorous the approximation at (3.18) and (3.14) then follows using (3.17) and Stirling's formula, as at (3.19).

Fix /epsilon1 ∈ (0 , λ -1 ) and let A /epsilon1 1 = { j ∈ Z : 0 ≤ j &lt; N λ ( -1 -/epsilon1 ) } , A /epsilon1 2 = { j ∈ Z : N λ ( -1 -/epsilon1 ) ≤ j ≤ N λ ( -1 + ) /epsilon1 } and A /epsilon1 3 = { j ∈ Z : N λ ( -1 + ) /epsilon1 &lt; j ≤ N -1 } . Further, let X ( N ) denote a Poisson random variable with mean N/λ . Then, using Chebyshev's inequality, P ( X ( N ) ∈ A /epsilon1 1 ) → 0 and P ( X ( N ) ∈ A /epsilon1 2 ) → 1 as N → ∞ . Also, by large deviation theory, there exists a &gt; 0, independent of N , such that P ( X ( N ) &gt; N λ ( -1 + ) /epsilon1 ) ≤ e -aN . Now

$$N \sum _ { j \in A _ { 1 } ^ { \circ } } \frac { ( N / \lambda ) ^ { j } \text{e} ^ { - N / \lambda } } { ( N - j ) j! } \leq \frac { 1 } { \lambda ^ { - 1 } - \epsilon } P \left ( X ^ { ( N ) } \in A _ { 1 } ^ { \epsilon } \right ) \to 0 \ \text{ as } \ N \to \infty \quad ( A. 7 )$$

and

$$N \sum _ { j \in A _ { 3 } ^ { 5 } } \frac { ( N / \lambda ) ^ { j } e ^ { - N / \lambda } } { ( N - j ) j! } \leq N \text{P} \left ( X ^ { ( N ) } > N ( \lambda ^ { - 1 } + \epsilon ) \right ) \rightarrow 0 \quad \text{as} \quad N \rightarrow \infty. \quad ( \text{A.8} )$$

Also,

$$\frac { 1 } { 1 - \lambda ^ { - 1 } + \epsilon } \mathsf P \left ( X ^ { ( N ) } \in A _ { 2 } ^ { \epsilon } \right ) \leq N \sum _ { j \in A _ { 2 } ^ { \epsilon } } \frac { ( N / \lambda ) ^ { j } \mathsf e ^ { - N / \lambda } } { ( N - j ) j! } \leq \frac { 1 } { 1 - \lambda ^ { - 1 } - \epsilon } \mathsf P \left ( X ^ { ( N ) } \in A _ { 2 } ^ { \epsilon } \right ),$$

whence, using (A.7), (A.8) and lim N →∞ P ( X ( N ) ∈ A /epsilon1 2 ) = 1,

$$\liminf _ { N \to \infty } N \sum _ { j = 0 } ^ { N - 1 } \frac { ( N / \lambda ) ^ { j } \mathrm e ^ { - N / \lambda } } { ( N - j ) j! } \geq \frac { 1 } { 1 - \lambda ^ { - 1 } + \epsilon }$$

and

$$\lim \sup _ { N \to \infty } N \sum _ { j = 0 } ^ { N - 1 } \frac { ( N / \lambda ) ^ { j } e ^ { - N / \lambda } } { ( N - j ) j! } \leq \frac { 1 } { 1 - \lambda ^ { - 1 } - \epsilon }. \quad \quad ( A. 1 0 ) \\ \intertext { + \infty \, \infty } \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \end{array}$$

Letting /epsilon1 ↓ 0 in (A.9) and (A.10) yields (A.6), as required.

where and

## Appendix B. Proof of Lemma 3.3

For i = 0 1 , , . . . , N , let T ( N ) ( ) i = inf { t : Y ( N ) t = i } be the first time that the number of infectives Y ( N ) t in the SIS epidemic equals i , with the convention that T ( N ) ( ) i = ∞ if Y ( N ) t never reaches i . Further, let F N ( ) i = { T ( N ) ( ) i &lt; ∞} and F C N ( ) i = { T ( N ) ( ) i = ∞} . Also, let y ( N ) = /floorleft ( λ -1) N/λ /floorright and, for k ∈ (0 , ∞ ), let y ( N ) k = ⌊ ( λ -1) N/λ -k √ N ⌋ . Then

$$\mu ^ { ( N ) } = \mathbb { E } \left [ T ^ { ( N ) } \right ] = \mathbb { E } \left [ T ^ { ( N ) } 1 _ { F ^ { C } _ { N } ( y _ { k } ^ { ( N ) } ) } \right ] + \mathbb { E } \left [ T ^ { ( N ) } 1 _ { F _ { N } ( y _ { k } ^ { ( N ) } ) } \right ]. \\.,$$

Further,

$$& \mathrm E \left [ T ^ { ( N ) } 1 _ { F _ { N } ( y _ { k } ^ { ( N ) } ) } \right ] = \mathrm E \left [ T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) 1 _ { F _ { N } ( y _ { k } ^ { ( N ) } ) } \right ] + \mathrm E \left [ \left ( T ^ { ( N ) } - T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) \right ) 1 _ { F _ { N } ( y _ { k } ^ { ( N ) } ) } \right ]. \\ & \mathrm r n$$

Thus,

$$\mu ^ { ( N ) } = a _ { k } ^ { ( N ) } + b _ { k } ^ { ( N ) },$$

$$a _ { k } ^ { ( N ) } = \mathbf E \left [ T ^ { ( N ) } 1 _ { F _ { N } ^ { C } ( y _ { k } ^ { ( N ) } ) } \right ] + \mathbf E \left [ T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) 1 _ { F _ { N } ( y _ { k } ^ { ( N ) } ) } \right ]$$

$$b _ { k } ^ { ( N ) } & = E \left [ \left ( T ^ { ( N ) } - T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) \right ) 1 _ { F _ { N } ( y _ { k } ^ { ( N ) } ) } \right ]. & & \text{(B.3)} \\ \text{that } A _ { - } ^ { ( N ) } & \text{is the total time that the endemic sends with $n$ individuals}$$

$$0 \leq a _ { k } ^ { ( N ) } \leq \sum _ { n = 1 } ^ { y _ { k } ^ { ( N ) } } \mathrm E \left [ A _ { n } ^ { ( N ) } \right ] = \sum _ { n = 1 } ^ { y _ { k } ^ { ( N ) } } \frac { ( N - 1 )! } { n ( N - n )! } \left ( \frac { \lambda } { N } \right ) ^ { n - 1 }, \quad \quad \text{(B.4)}$$

Recall that A ( N ) n is the total time that the epidemic spends with n individuals infective. Note that if T ( N ) ( y ( N ) k ) = ∞ then T ( N ) ≤ ∑ y ( N ) k n =1 A ( N ) n , and if T ( N ) ( y ( N ) k ) &lt; ∞ then T ( N ) ( y ( N ) k ) ≤ ∑ y ( N ) k n =1 A ( N ) n , so using (2.11) and (3.10). Denote the right hand sum in (B.4) by c ( N ) k . Then rearranging as at (3.17) yields

$$c _ { k } ^ { ( N ) } = ( N - 1 )! ( \lambda / N ) ^ { N - 1 } e ^ { N / \lambda } \sum _ { j = z _ { k } ^ { ( N ) } } ^ { N - 1 } \frac { ( N / \lambda ) ^ { j } e ^ { - N / \lambda } } { ( N - j ) j! }, \quad \text{(B.5)}$$

where z ( N ) k = ⌈ Nλ -1 + k √ N ⌉ . Omitting the details, fixing /epsilon1 ∈ (0 , λ -1 ), splitting the sum in (B.5) into z ( N ) k ≤ j ≤ N λ ( -1 + ) and /epsilon1 N λ ( -1 + ) /epsilon1 &lt; j ≤ N -1, invoking the

central limit theorem and letting /epsilon1 ↓ 0 yields

$$\lim _ { N \to \infty } \frac { c _ { k } ^ { ( N ) } } { \mu ^ { ( N ) } } = \Phi ( - k \sqrt { \lambda } ), \quad \ \quad \ \ ( B. 6 ) \\. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

where Φ denotes the standard normal cumulative distribution function. Thus, recalling (A.6), for any k &gt; 0,

$$0 \leq \liminf _ { N \to \infty } \frac { a _ { k } ^ { ( N ) } } { \mu ^ { ( N ) } } \leq \limsup _ { N \to \infty } \frac { a _ { k } ^ { ( N ) } } { \mu ^ { ( N ) } } \leq \Phi ( - k \sqrt { \lambda } ). \quad \quad ( B. 7 ) \\ \intertext { - \cdots \cdots \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots - \cdots } }$$

(Note that a similar argument using (3.9) shows that the mean time a supercritical branching process takes to reach size n , given that it does not go extinct, is O (log n ).)

We next determine lim N →∞ P ( F N ( y ( N ) k ) ) ( k ∈ (0 , ∞ )). For λ &gt; 0, let B λ denote a branching process, with one ancestor, in which individuals have iid lifetimes according to Q and, whilst alive, give birth at the points of a homogeneous Poisson point process with rate λ . Let p Q ( λ ) denote the probability that B λ does go extinct and, for t ≥ 0, let Y t ( λ ) denote the number of individuals alive in B λ at time t . Observe that by Lemma 4.1 of Britton et al. (2007), p Q ( λ ′ ) → p Q ( λ ) as λ ′ → λ . Note that p Q = p Q ( λ ).

For i = 0 1 , , . . . , let T λ ( ) i = inf { t : Y t ( λ ) = i } , where T λ ( ) i = ∞ if Y t ( λ ) never reaches i . Fix δ ∈ [0 , 1). Then, whilst Y ( N ) t ≤ Nδ , the SIS epidemic is bounded below by the branching process B (1 -δ λ ) ( c.f. Whittle (1955)), so

$$1 - p _ { Q } ( ( 1 - \delta ) \lambda ) & \leq P \left ( T _ { ( 1 - \delta ) \lambda } ( \lfloor \delta N \rfloor ) < \infty \right ) \leq P \left ( T ^ { ( N ) } ( \lfloor \delta N \rfloor ) < \infty \right ). \quad \text{(B.8)} \\ \text{Further, willst } Y _ { t } ^ { ( N ) } & \leq y _ { k } ^ { ( N ) }, \text{ the SIS epidemic is bounded below by the branching}$$

Further, whilst Y ( N ) t ≤ y ( N ) k , the SIS epidemic is bounded below by the branching process B 1+ λkN -1 2 . Using for example equation (5.63) of Haccou et al. (2005),

$$& p _ { Q } \left ( 1 + \lambda k N ^ { - \frac { 1 } { 2 } } \right ) \leq 1 - \frac { 2 \lambda k } { N ^ { \frac { 1 } { 2 } } E [ Q ^ { 2 } ] ( 1 + \lambda k N ^ { - \frac { 1 } { 2 } } ) ^ { 2 } } \leq \exp \left ( - \frac { 2 \lambda k } { N ^ { \frac { 1 } { 2 } } E [ Q ^ { 2 } ] ( 1 + \lambda k N ^ { - \frac { 1 } { 2 } } ) ^ { 2 } } \right ). \\ & \intertext { T h e n e } \text{ for anv } \delta \in \mathcal { N } \ \mathcal { N } - 1 \mathcal { N } / \mathcal { N } \ \mathcal { N } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \mathcal { n } \ \cdots 0 }$$

Thus, for any δ ∈ (0 , ( λ -1) /λ ) and any k &gt; 0,

$$\i i n s, \, \i i n a \, \rightarrow \, \upsilon \in \{ \upsilon, \, \langle \, \Lambda ^ { - } - \i i \rangle \, \Lambda ^ { \prime } \, \Lambda \, \rightarrow \, \upsilon \,, \\ \, P \left ( T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) = \infty \, | \, T ^ { ( N ) } ( \lfloor \delta N \rfloor ) < \infty \right ) \, & \leq \, \left ( 1 - \frac { 2 \lambda k } { N ^ { \frac { 1 } { 2 } } E [ Q ^ { 2 } ] ( 1 + \lambda k N ^ { - \frac { 1 } { 2 } } ) ^ { 2 } } \right ) ^ { \lfloor \delta N \rfloor } ( B. 9 ) \\ & \leq \, \exp \left ( - \frac { 2 \lfloor \delta N \rfloor \, \lambda k } { N ^ { \frac { 1 } { 2 } } E [ Q ^ { 2 } ] ( 1 + \lambda k N ^ { - \frac { 1 } { 2 } } ) ^ { 2 } } \right ) \\ & \to \, \ 0 \quad \text{as } N \to \infty, \\ \text{which on combinum with } / R \, \lambda \, \text{and} \, \text{lottime} \, \delta \, | \, \i i n \, \text{valid} \,.$$

which on combining with (B.8) and letting δ ↓ 0 yields

$$\liminf _ { N \to \infty } P \left ( F _ { N } ( y _ { k } ^ { ( N ) } ) \right ) \geq 1 - p _ { Q } ( \lambda ).$$

(Note that the bound (B.9) may not hold if Q does not have an exponential distribution, since the excess lives of the /floorleft δN /floorright individuals alive at time T ( N ) ( /floorleft δN /floorright ) will not be distributed as Q . However, it is clear that there exist /epsilon1 1 , /epsilon1 2 &gt; 0 such that P ( Y ( N ) ( /epsilon1 1 ) ≥ /epsilon1 2 /floorleft δN /floorright ) → 1 as N → ∞ , where Y ( N ) ( /epsilon1 1 ) denotes the number of individuals at time T ( N ) ( /floorleft δN /floorright ) whose excess life is at least /epsilon1 1 . Conditioning on the number of offspring of those Y ( N ) ( /epsilon1 1 ) individuals and arguing as above shows that P ( T ( N ) ( y ( N ) k ) = ∞| T ( N ) ( /floorleft δN /floorright ) &lt; ∞ → ) 0 as N → ∞ . A similar comment applies elsewhere in the proof.)

The SIS epidemic is bounded above by B λ , so

$$\lim _ { N \to \infty } P \left ( F _ { N } ( y _ { k } ^ { ( N ) } ) \right ) & \leq \lim _ { N \to \infty } P \left ( T _ { \lambda } ( y _ { k } ^ { ( N ) } ) < \infty \right ) = 1 - p _ { Q } ( \lambda ), \\ \text{whenever, for any $k>0$}.$$

whence, for any k &gt; 0,

$$\lim _ { N \to \infty } P \left ( F _ { N } ( y _ { k } ^ { ( N ) } ) \right ) = 1 - p _ { Q } ( \lambda ). \quad \quad \quad ( B. 1 0 ) \\ \text{that $k>l>0$.\ Ar}\sin\sigma$ as above shows that}$$

Suppose now that k &gt; l &gt; 0. Arguing as above shows that whence

$$P \left ( T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) < \infty \, | \, T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) < \infty \right ) \to 1 \ \text{ as } \ N \to \infty, \quad \text{(B.11)} \\ \text{whence}$$

$$\frac { b _ { k } ^ { ( N ) } } { \mu ^ { ( N ) } } & = \frac { 1 } { \mu ^ { ( N ) } } \mathrm E \left [ \left ( T ^ { ( N ) } - T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) \right ) \, 1 _ { F _ { N } ( y _ { l } ^ { ( N ) } ) } \right ] \\ & + \frac { 1 } { \mu ^ { ( N ) } } \mathrm E \left [ \left ( T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) - T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) \right ) \, 1 _ { F _ { N } ( y _ { l } ^ { ( N ) } ) } \right ] + o ( 1 ) \quad \text{as } N \to \infty. \quad ( B. 1 2 ) \\ \text{Let } \, \theta ^ { ( N ) } & = 1 + \lambda I N ^ { - \frac { 1 } { 2 } }. \ \text{Whilst} \, y _ { k } ^ { ( N ) } \leq Y _ { t } ^ { ( N ) } \leq y _ { l } ^ { ( N ) }, \text{ the epidemic is bounded below}$$

Let θ ( N ) = 1 + λlN -1 2 . Whilst y ( N ) k ≤ Y ( N ) t ≤ y ( N ) l , the epidemic is bounded below by B θ ( N ) but now starting with y ( N ) k individuals. Hence, using (B.11),

$$\begin{array} { c } \quad \. \ \\ \frac { 1 } { \mu ^ { ( N ) } } \text{E} \left [ \left ( T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) - T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) \right ) 1 _ { F _ { N } ( y _ { l } ^ { ( N ) } ) } \right ] \\ \leq \frac { 1 } { \mu ^ { ( N ) } } \text{E} \left [ \left ( T _ { \theta ^ { ( N ) } } ( y _ { l } ^ { ( N ) } ) - T _ { \theta ^ { ( N ) } } ( y _ { k } ^ { ( N ) } ) \right ) 1 _ { F _ { N } ( k, l ) } \right ] + o ( 1 ) \quad \text{as $N\to\infty$,} \quad ( B. 1 3 ) \\ \text{where $F_{N}(k,l)=\{T_{\theta^{\left(N})}(y_{i}^{(N)})<\infty\}$, assuming that $\mathcal{B}_{\theta^{\left(N})$ reaches $y_{i}^{(N)}$.} \end{array}$$

where F N ( k, l ) = { T θ ( N ) ( y ( N ) l ) &lt; ∞} , assuming that B θ ( N ) reaches y ( N ) k .

Consider the branching process B λ . For n = 1 2 , , . . . , let A λ ( n ) be the total time that B λ spends with n individuals alive. (Thus, in the notation of Section 3.1, A λ ( n ) = A n .) Suppose that λ &gt; 1. Then, since

$$T _ { \lambda } ( n ) 1 _ { \{ T _ { \lambda } ( n ) < \infty \} } \leq \sum _ { i = 1 } ^ { n } A _ { \lambda } ( i ) 1 _ { \{ T _ { \lambda } ( n ) < \infty \} } \leq \sum _ { i = 1 } ^ { n } A _ { \lambda } ( i ),$$

it follows from (3.9) that

$$\mathbf E \left [ T _ { \lambda } ( n ) \, | \, T _ { \lambda } ( n ) < \infty \right ] \leq \frac { 1 } { \lambda ( 1 - p _ { Q } ( \lambda ) ) } \sum _ { i = 1 } ^ { n } \frac { 1 } { i }.$$

Return to the branching process B θ ( N ) . The expected time for it to reach y ( N ) l , starting from y ( N ) k individuals, given that it does so, is less than the expected time for B θ ( N ) to reach y ( N ) l , starting from one individual, again conditional upon it doing so. Thus, since θ ( N ) &gt; 1,

$$\mathbf E \left [ \left ( T _ { \theta ^ { ( N ) } } ( y _ { l } ^ { ( N ) } ) - T _ { \theta ^ { ( N ) } } ( y _ { k } ^ { ( N ) } ) \right ) 1 _ { F _ { N } ( k, l ) } \right ] \leq \left ( 1 - p _ { Q } ( \theta ^ { ( N ) } ) ^ { y _ { k } ^ { ( N ) } } \right ) \frac { 1 } { 1 - p _ { Q } ( \theta ^ { ( N ) } ) } \sum _ { i = 1 } ^ { y _ { l } ^ { ( N ) } } \frac { 1 } { i }. \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ 0 \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \csize \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \vdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots  \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ \\ 0 \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cd> \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \clots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots
 \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cmathopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopoppopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopcoropopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopopop
 \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cd: \cdots \cdots \cd$$

Theorem 5.5 of Haccou et al. (2005) yields 1 -p Q ( θ ( N ) ) ∼ 2 λlN -1 2 / E[ Q 2 ], whence, using (B.14), E [( T θ ( N ) ( y ( N ) l ) -T θ ( N ) ( y ( N ) k ) ) 1 F N ( k,l ) ] = O N ( 1 2 log N ) as N → ∞ . It then follows from (B.10) to (B.13) that

$$\frac { b _ { k } ^ { ( N ) } } { \mu ^ { ( N ) } } = \frac { 1 } { \mu ^ { ( N ) } } ( 1 { - } p _ { Q } ( \lambda ) ) E \left [ T ^ { ( N ) } - T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) \, | \, T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) < \infty \right ] + o ( 1 ) \quad \text{as } N \to \infty.$$

Dividing (B.1) by µ ( N ) and letting N →∞ yields, after using (B.7),

$$1 - \Phi ( - k \sqrt { \lambda } ) & \leq \liminf _ { N \to \infty } ( 1 - p _ { Q } ( \lambda ) ) \frac { E \left [ T ^ { ( N ) } - T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) \, | \, T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) < \infty \right ] } { \mu ^ { ( N ) } } \\ & \leq \lim _ { N \to \infty } ( 1 - p _ { Q } ( \lambda ) ) \frac { E \left [ T ^ { ( N ) } - T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) \, | \, T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) < \infty \right ] } { \mu ^ { ( N ) } } \leq 1, \\ \text{whenever, letting $\sigma$ $k \to \infty$}.$$

whence, letting k →∞ ,

$$\mathbf E \left [ T ^ { ( N ) } - T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) \, | \, T ^ { ( N ) } ( y _ { l } ^ { ( N ) } ) < \infty \right ] \sim \frac { 1 } { 1 - p _ { Q } ( \lambda ) } \mu ^ { ( N ) }.$$

The above shows that, for any k &gt; 0, the mean time to extinction from y ( N ) k = ⌊ ( λ -1) N/λ -k √ N ⌋ , µ ( N ) k say, satisfies µ ( N ) k ∼ 1 1 -p Q µ ( N ) , since p Q ( λ ) = p Q . Suppose that we start at the endemic level with y ( N ) = ( /floorleft λ -1) N/λ /floorright infectives. Then the mean time to extinction from the endemic level satisfies

$$\mu _ { Q } ^ { ( N ) } = D _ { Q, k } ^ { ( N ) } + \mu _ { k } ^ { ( N ) },$$

where D ( N ) Q,k denotes the mean time that the epidemic takes to reach y ( N ) k for the first time. We prove that D ( N ) Q,k = O N ( ) for 1 &lt; λ ≤ 2 and that D ( N ) Q,k = O N ( log N ) for λ &gt; 2.

Turning first to the case when 1 &lt; λ ≤ 2, consider a subcritical branching process with immigration, B ( N ) , where individuals immigrate into the population at the points of a homogeneous Poisson process with rate λN (1 -1 /λ ) 2 . The lifetimes of individuals in B ( N ) are independent and identically distributed according to Q , and whilst alive, individuals give birth at the points of a homogeneous Poisson process with rate 2 -λ . Let B ( N ) t denote the total number of individuals alive in B ( N ) at time t . Then for any t ≥ 0, Neal (2014), Corollary 2.1, gives that

$$\sqrt { N } ( B _ { t } ^ { ( N ) } / N - ( 1 - 1 / \lambda ) ) \stackrel { D } { \longrightarrow } N ( 0, 1 / \lambda ) \text{ \ as } N \to \infty. \text{ \quad \ } ( B. 1 7 )$$

It is shown in Neal (2014), Section 3, that when 1 &lt; λ ≤ 2, the SIS epidemic can be coupled to B ( N ) , such that for all t ≥ 0, Y ( N ) t ≤ B ( N ) t . Therefore D ( N ) Q,k ≤ ˜ D ( N ) Q,k , the mean time that B ( N ) takes to reach y ( N ) k for the first time starting with /floorleft ( λ -1) N /floorright individuals. Therefore we focus on computing an upper bound for ˜ D ( N ) Q,k .

Let T ( N ) 1 denote the total length of time it takes for the family trees of all individuals alive at time 0 to go extinct which is stochastically smaller than the sum of the extinction times of the family trees from each of the /floorleft ( λ -1) N /floorright individuals. Thus, using (3.5),

$$\mathbf E \left [ T _ { 1 } ^ { ( N ) } \right ] & \leq \frac { \lambda - 1 } { \lambda } N \left \{ \frac { - \log ( 1 - ( 2 - \lambda ) ) } { 2 - \lambda } \right \} = O ( N ) \quad \text{as } N \to \infty. \quad \text{(B.18)} \\ \mathfrak { r } \cdot \hat { n } ^ { ( N ) } \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon$$

Let ˆ B ( N ) 1 ,t denote the total progeny still alive at time t of the family trees of individuals who immigrate into the population in the interval [0 , T ( N ) 1 ]. Then ˆ B ( N ) 1 ,t is independent of the family trees of individuals who immigrate into the population prior to time 0 and B ( N ) T ( N ) 1 = ˆ B ( N ) 1 ,T ( N ) 1 , where regardless of the value of T ( N ) 1 , ˆ B ( N ) 1 ,T ( N ) 1 is stochastically smaller than B ( N ) ∗ , the total number of individuals alive at time 0 of an independent copy of B ( N ) started at time -∞ . Using (B.17) it is straightforward to show that

$$P ( \hat { B } ^ { ( N ) } _ { 1, T ^ { ( N ) } _ { 1 } } < y ^ { ( N ) } _ { k } ) & \ \geq \ \ P ( B ^ { ( N ) } _ { * } < y ^ { ( N ) } _ { k } ) \\ & \quad \to \ \Phi ( - k \sqrt { \lambda } ) \quad \text{as } N \to \infty. \quad \quad ( B. 1 9 ) \\. \quad \ \text{$\quad \mathbb{ }$}$$

If ˆ B ( N ) 1 ,T ( N ) 1 &lt; y ( N ) k , we know that the first time B ( N ) t reaches y ( N ) k is less than or equal

to T ( N ) 1 . Otherwise, we consider B ( N ) t at a sequence of times S ( N ) k = ∑ k i =1 T ( N ) i , where for i = 2 3 , , . . . , T ( N ) i denotes the total length of time (from S ( N ) i -1 ) it takes for the family trees of all individuals alive at time S ( N ) i -1 to go extinct. This is stochastically smaller than the sum of the extinction times of the family trees from each of the ˆ B ( N ) i -1 ,S ( N ) i -1 individuals alive at time S ( N ) i -1 , where ˆ B ( N ) i -1 ,t denotes the total progeny still alive at time t of the family trees of individuals who immigrate into the population in the interval [ S ( N ) i -2 , S ( N ) i -1 ]. ( S ( N ) 0 = T ( N ) 0 = 0 ) . Note that the ˆ B ( N ) i,S ( N ) i ( i = 1 2 , , . . . , T ( N ) i -1 ) are conditionally independent given { T ( N ) i } and, moreover, T ( N ) i depends upon { ( ˆ B ( N ) 1 ,S ( N ) 1 , T ( N ) 1 ) , . . . , ( ˆ B ( N ) i -1 ,S ( N ) i -1 , T ( N ) i -1 ) } only through ˆ B ( N ) i -1 ,S ( N ) i -1 . It is then straightforward using similar arguments to (B.18) and (B.19) to show that

$$\intertext { a c u a g u a t o w a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a u a } } E \left [ T _ { i } ^ { ( N ) } | \hat { B } _ { i - 1, S _ { i - 1 } ^ { ( N ) } } > y _ { k } ^ { ( N ) } \right ] & \quad = \quad E \left [ \hat { E } \left [ \hat { B } _ { i } ^ { ( N ) } | \hat { B } _ { i - 1, S _ { i - 1 } ^ { ( N ) } } ^ { ( N ) } > y _ { k } ^ { ( N ) } \right ] \frac { - \log ( 1 - ( 2 - \lambda ) ) } { 2 - \lambda } \\ & = \quad O ( N ) \quad \text{as $N\to\infty$,} \quad \text{(B.20)} \\ R ^ { ( N ) } \ - \ \hat { R } ^ { ( N ) } \quad \text{and}$$

B ( N ) S ( N ) i -1 = ˆ B ( N ) i -1 ,S ( N ) i -1 and

$$\mathbf P \left ( \hat { B } _ { i - 1, S _ { i$$

Let L ( N ) = min { l ; ˆ B ( N ) l,S ( N ) l &lt; y ( N ) k } . Note that by (B.18) and (B.20), there exists M &lt; ∞ such that for all i = 1 2 , , . . . and for all sufficiently large N , E[ T ( N ) i | ˆ B ( N ) i -1 ,S ( N ) i -1 &gt; y ( N ) k ] ≤ MN . Therefore

$$\leq M N. \text{ Therefore} \\ \mathrm E [ S ^ { ( N$$

since for any l ≥ 1 and lim sup N →∞ P( L ( N ) &gt; l ) ≤ (1 -(Φ( -k √ λ )) . l Hence ˜ D N Q,k ≤ E[ S ( N ) L ( N ) ] = O N ( ).

Before considering the case λ &gt; 2, it is fruitful to derive an upper bound for the expected time the SIS epidemic takes to reach its endemic level y ( N ) given that it is

currently above that level. Thus, suppose that the SIS epidemic starts with k &gt; y ( N ) infectives and define T ( N ) ( y ( N ) ) as before, but note that now P ( T ( N ) ( y ( N ) ) &lt; ∞ ) = 1. A simple calculation shows that, whilst Y ( N ) t &gt; y ( N ) , Y ( N ) t is bounded below by the subcritical branching process B η ( N ) , where η ( N ) = 1 -λ/N , starting from k individuals. It follows that T ( N ) ( y ( N ) ) is stochastically smaller than the extinction time of this branching process, which in turn is stochastically smaller than the sum of the extinction times of the family trees from each of the k initial individuals in B η ( N ) . Thus, using (3.5), for any k &gt; y ( N ) ,

$$E \left [ T ^ { ( N ) } ( y ^ { ( N ) } )$$

Fix λ &gt; 2 and k &gt; 0. For all N &gt; ( k/ ( 1 2 -1 λ )) 2 and all n &gt; y ( N ) k , the infection rate satisfies

$$\begin{array} {$$

where β N,k,λ ( ) = ( 1 -1 λ ) N + kλ √ N . For such N it follows that the SIS process { Y ( N ) t : t ≥ 0 } can be coupled to a birth-death type process { ˜ Y ( N ) t : t ≥ 0 } having birth rate given by α n ( ) = β N,k,λ ( ) if y ( N ) k &lt; n ≤ N and α n ( ) = 0 otherwise, such that if Y ( N ) 0 = ˜ Y ( N ) = y ( N ) and these y ( N ) initial individuals have the same excess lifetimes in the two processes, then Y ( N ) t ≤ ˜ Y ( N ) t for all 0 ≤ t ≤ T ( N ) ( y ( N ) k ).

Let T ( N ) 1 denote the total length of time it takes for the y ( N ) initial individuals in { Y ( N ) t : t ≥ 0 } (or { ˜ Y ( N ) t : t ≥ 0 } ) to all die. Then, conditional upon T ( N ) 1 , ˜ Y ( N ) T ( N ) 1 has a Poisson distribution with mean bounded above by

$$\beta ( N, k, \lambda ) \int _ { 0 } ^ { T _ { 1$$

since E[ Q ] = 1. Let X ( N ) denote a Poisson random variable with mean β N,k,λ ( ). Then, the above coupling implies that

$$\mathbf P \left ( T ^ { ( N ) } ( y _ { k } ^ { ( N ) } ) \leq T _ { 1 } ^ { ( N ) } \right ) \geq \mathbf P \left ( \tilde { Y } _ { T _ { 1 } ^ { ( N ) } } ^ { ( N ) } \leq y _ { k } ^ { ( N ) } \right ) \geq \mathbf P \left ( X ^ { ( N ) } \leq y _ { k } ^ { ( N ) } \right ).$$

Straightforward application of the central limit theorem yields that

$$\mathsf P \left ( X ^ { ( N ) } & \leq y _ { k } ^ { ( N ) } \right ) \to \Phi \left ( - k ( \lambda + 1 ) \sqrt { \frac { \lambda } { \lambda - 1 } } \right ) \quad \text{as} \quad N \to \infty.$$

Note also that T ( N ) 1 is less than the sum of the excess lifetimes of the y ( N ) individuals alive at time t = 0, so E [ T ( N ) 1 ] = O N ( ).

Finally, for i = 1 2 , , . . . , N , let T ( N ) i = min { t &gt; 0 : Y ( N ) t = 0 } be the extinction time of the SIS epidemic given that initially there are i infectives. Then, together with Lemma 3.3, (B.23) implies that for any δ &gt; 0,

If Y ( N ) T ( N ) 1 &gt; y ( N ) k , let U ( N ) 1 = min { u ≥ 0 : Y ( N ) T ( N ) 1 + u ≤ y ( N ) } . Note that U ( N ) 1 = 0 if Y ( N ) T ( N ) 1 ≤ y ( N ) and (B.23) implies that E [ U ( N ) 1 | Y ( N ) T ( N ) 1 &gt; y ( N ) ] = O N ( log N ). Thus, E [ U ( N ) 1 ] = O N ( log N ). Now let T ( N ) 2 be the total length of time it takes for the Y ( N ) T ( N ) 1 + U ( N ) 1 individuals alive in { Y ( N ) t : t ≥ 0 } at time T ( N ) 1 + U ( N ) 1 to all die. The argument now proceeds in a similar fashion to the case when 1 &lt; λ ≤ 2 and, omitting the details, it is easily seen that if λ &gt; 2 then ˜ D N Q,k = O N ( log N ) and that this holds for all λ &gt; 1. It then follows using (B.16) that µ ( N ) Q ∼ 1 1 -p Q µ ( N ) , as required.

$$\mathrm E \left [ T _ { i } ^ { ( N ) } \right ] \sim \frac { 1 } { 1 - p _ { Q } } \mu ^ { ( N ) } \quad \text{for any } i \geq N \delta.$$

## Acknowledgements

We would like to thank Amaury Lambert for bringing Lambert (2011) and the proof of (3.9) to our attention. Tom Britton is grateful to the Swedish Research Council for financial support.

## References

Aldous, D. J. (1985) Exchangeability and related topics. In Ecole d' ´ t´ de Probabilit´s E e e de Saint-Flour XIII, 1993 (Lecture Notes in Mathematics 1117 ), ed. Hennequin, P. L. Springer, Berlin.

Andersson, H. and Djehiche, B. (1998) A threshold limit theorem for the stochastic logistic epidemic. J. Appl. Prob. 35 , 662-670.

Asmussen, S. (1987) Applied Probability and Queues. John Wiley, Chichester.

Ball, F. (1999). Stochastic and deterministic models for SIS epidemics among a population partitioned into households. Math. Biosci. , 156 , 41-67.

- Ball, F.G. and Donnelly, P. (1995) Strong approximations for epidemic models. Stoch. Proc. Appl. , 55 , 1-21.
- Ball, F.G. and Milne, R.K. (2004) Applications of simple point process methods to superpositions of aggregated stationary processes. Aust. N. Z. J.Stat. , 46 , 181-196.
- Ball, F., Mollison, D. and Scalia-Tomba, G. (1997). Epidemics with two levels of mixing. Ann. Appl Probab. , 7 , 46-89.
- Britton, T., Janson, S. and Martin-L¨f, o A. (2007) Graphs with specified degree distributions, simple epidemics, and local vaccination strategies. Adv. Appl. Prob. , 39 , 922-948.
- Britton, T. and Neal, P. (2010) The time to extinction for an SIS-household-epidemic model. J. Math. Biol. , 61 , 763-779.
- Ethier, S.N. and Kurtz, T. (1986) Markov Processes: Characterization and Convergence. John Wiley, New York.
- Ghoshal, G. Sander, L. and Sokolov, I. (2004). SIS epidemics with household structure: the self-consistent field method. Math. Biosci. , 190 , 71-85.
- Haccou, P.,Jagers, P. and Vatutin, V. A. (2005) Branching Processes: Variation, Growth and Extinction of Populations. Cambridge University Press.
- Hern´ andez-Su´rez, C.M. and Castillo-Chavez, C. (1999) A basic result on the integral a for birth-death Markov processes. Math Biosci. , 161 , 95-104.
- Kryscio, R. and Lef` evre, C. (1989). On the extinction of the SIS stochastic logistic epidemic. J. Appl. Prob. , 26 , 685-694.
- Lambert, A. (2011) Species abundance distributions in neutral models with immigration or mutation and general lifetimes. J. Math. Biol. , 63 , 57-72.
- N˚ asell, I. (1999) On the time to extinction in recurrent epidemics. J. R. Stat. Soc. Ser. B Stat. Methodol. , 61 , 309-330.

Neal, P. (2006) Stochastic and deterministic analysis of SIS household epidemics. Adv. Appl. Prob. 38 , 943-968. (See also correction in Adv. Appl. Prob. 44 , 309310 (2012).)

- Neal, P. (2014) Endemic behaviour of SIS epidemics with general infectious period distributions. Adv. Appl. Prob. , 46 , 241-255.
- Sevast'yanov, B.A. (1957) An ergodic theory for Markov processes and its application to a telephone system with refusals. Theory Prob. Appl. , 2 , 104-112.

Whittle, P. (1955) The outcome of a stochastic epidemic - a note on Bailey's paper. Biometrika , 42 , 116-122.

Whittle, P. (1985) Partial balance and insensitivity. J. Appl. Probab. , 22 , 168-176.

Zachary, S. (2007). A note on insensitivity in stochastic networks. J. Appl. Probab. , 44 , 238-248.
Frank Ball Tom Britton Peter Neal
On expected durations of birth-death processes, with applications to branching processes and SIS epidemics
We study continuous-time birth-death type processes, where individuals have
independent and identically distributed lifetimes, according to a random
variable Q, with E[Q]=1, and where the birth rate if the population is
currently in state (has size) n is \alpha(n). We focus on two important
examples, namely \alpha(n)=\lambda n being a branching process, and
\alpha(n)=\lambda n(N-n)/N which corresponds to an SIS epidemic model in a
homogeneously mixing community of fixed size N. The processes are assumed to
start with a single individual, i.e. in state 1. Let T, A_n, C and S denote the
(random) time to extinction, the total time spent in state $n$, the total
number of individuals ever alive and the sum of the lifetimes of all
individuals in the birth-death process, respectively. The main results of the
paper give expressions for the expectation of all these quantities, and shows
that these expectations are insensitive to the distribution of Q. We also
derive an asymptotic expression for the expected time to extinction of the SIS
epidemic, but now starting at the endemic state, which is not independent of
the distribution of Q. The results are also applied to the household SIS
epidemic, showing that its threshold parameter R_* is insensitive to the
distribution of Q, contrary to the household SIR epidemic, for which R_* does
depend on Q.