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.

RA<br />

The variance-to-mean ratio <strong>of</strong> Z can be calculated from<br />

2<br />

( g ′ () + g′<br />

() 1 − ( g′<br />

() 1 ) ) g′<br />

( 1)<br />

RA 1 RA RA<br />

RA<br />

c<br />

′ (Harris 1989) and shown to equal 1+R0/k+cR0,<br />

which increases monotonically as c increases. For direct comparison with other<br />

<strong>of</strong>fspring distributions, this composite distribution under RA control can be<br />

RA<br />

RA RA<br />

approximated by a new negative binomial distribution, Z ~NegB( R , k ) where<br />

RA<br />

Rc<br />

RA<br />

is given above and k is estimated using the proportion <strong>of</strong> zeros method as the<br />

solution to<br />

p<br />

RA<br />

0<br />

0<br />

c<br />

( ) ( ) RA<br />

RA RA −k<br />

c<br />

1−<br />

p = 1+<br />

R k<br />

= p + c<br />

. This approximation yields better<br />

0<br />

c<br />

than 95% overlap with the exact distribution for k§1, and better than 85% overlap for<br />

almost all <strong>of</strong> parameter space (Figure S3A). (The proportion <strong>of</strong> overlap is calculated as<br />

1<br />

⎛<br />

⎜<br />

⎝<br />

− ∑ ∞<br />

i=<br />

0<br />

c<br />

⎞<br />

Z i,<br />

exact − Zi<br />

, approx ⎟ 2 , which scales from 0 to 1 as the two distributions go from<br />

⎠<br />

completely non-overlapping to identical.) The approximation approaches exactness for<br />

cØ0 and cØ1, and is least accurate for large values <strong>of</strong> k because it is unable to mimic<br />

RA<br />

the bimodal distribution <strong>of</strong> Z (Figure S3B). The approximated shape parameter<br />

c<br />

decreases monotonically as control effort c increases (Figure S3C).<br />

Relative efficacy <strong>of</strong> control policies<br />

For HP control, with all individuals’ transmission reduced by a factor c, the<br />

HP<br />

<strong>of</strong>fspring distribution is Z ~NegB((1−c)R0,k) and has pgf:<br />

c<br />

g<br />

HP<br />

⎛<br />

( s)<br />

= ⎜1+<br />

⎝<br />

R<br />

k<br />

0 ( 1−<br />

c)<br />

( 1−<br />

s)<br />

.<br />

141<br />

c<br />

⎞<br />

⎟<br />

⎠<br />

−k<br />

c<br />

c<br />

RA<br />

kc

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

Saved successfully!

Ooh no, something went wrong!