K - College of Natural Resources - University of California, Berkeley
K - College of Natural Resources - University of California, Berkeley
K - College of Natural Resources - University of California, Berkeley
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
fv(u). Demographic stochasticity in transmission is then represented by a Poisson<br />
process, as is standard in branching process treatments <strong>of</strong> epidemics (Diekmann and<br />
Heesterbeek 2000). This yields the following pgf for a Poisson distribution with mean<br />
ν distributed as fv(u):<br />
If ν is a constant, R0, then the pgf is:<br />
g<br />
( )<br />
∫ ∞<br />
− u 1−s<br />
( s)<br />
= e fν<br />
( u)<br />
du<br />
0<br />
g(<br />
s)<br />
= e<br />
−R0<br />
( 1−s<br />
)<br />
If ν is exponentially distributed with mean R0, the resulting <strong>of</strong>fspring<br />
distribution is geometric with mean R0 (Taylor and Karlin 1998) and pgf:<br />
( ( ) ) 1<br />
1<br />
−<br />
+ R −<br />
) ( g s = s<br />
If ν is gamma distributed, with mean R0 and shape parameter k, the resulting<br />
<strong>of</strong>fspring distribution is negative binomial, also with mean R0 and shape parameter k<br />
(Taylor and Karlin 1998), with pgf:<br />
0 1<br />
⎛ R0<br />
g(<br />
s)<br />
= ⎜1+<br />
1<br />
⎝ k<br />
( − s)<br />
This expression was applied in all <strong>of</strong> the general branching process results<br />
shown above to derive our results. The expression q=g(q) was solved numerically to<br />
generate Figures 3B and S2B, showing the dependence <strong>of</strong> the extinction probability on<br />
R0 and k. The negative binomial pgf itself is plotted in Figure S2A, showing how the<br />
probability <strong>of</strong> infecting zero others (p0) increases sharply with k for a given R0. The<br />
probability <strong>of</strong> extinction in the n th generation (Figure S2C) was calculated using<br />
gn(0)−gn−1(0). These numerical solutions match the averaged output <strong>of</strong> many<br />
138<br />
⎞<br />
⎟<br />
⎠<br />
−k