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