26.10.2012 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!