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.

original process is eliminated with probability c (regardless <strong>of</strong> the ν value <strong>of</strong> the source<br />

HP<br />

case) changes this to a marked Poisson process so that Z ~Poisson((1-c)ν) (Taylor<br />

and Karlin 1998). If uncontrolled individual reproductive numbers are gamma-<br />

distributed, ν~gamma(R0,k), then only the scale parameter <strong>of</strong> the resulting negative<br />

HP<br />

binomial distribution is affected by HP control and Z ~NegB((1−c)R0,k). The<br />

HP<br />

variance-to-mean ratio <strong>of</strong> Z is 1+(1−c)R/k, and decreases monotonically as control<br />

effort increases (Figure S3C).<br />

c<br />

Under random absolute (RA) control, each infected individual is controlled<br />

perfectly (such that they cause zero secondary infections) with probability c.<br />

Imposition <strong>of</strong> RA control influences transmission only for the fraction 1−p0 <strong>of</strong><br />

individuals whose natural Z value is greater than zero—<strong>of</strong> these a fraction c have<br />

RA<br />

Zc<br />

RA<br />

=0, while the remaining fraction 1−c are unaffected and have Z =Z. Under an<br />

RA<br />

RA control policy, therefore, the proportion <strong>of</strong> cases causing zero infections is p =<br />

p0+c(1−p0) and the population mean<br />

R<br />

RA<br />

c<br />

N<br />

N<br />

1<br />

1<br />

= ∑ Z i Pr( case i not controlled)<br />

= ∑<br />

N<br />

N<br />

i=<br />

1<br />

c<br />

( 1−<br />

c)<br />

Z = ( 1−<br />

c)<br />

R0<br />

i=<br />

1<br />

i<br />

c<br />

c<br />

. The exact<br />

RA<br />

RA<br />

RA<br />

RA<br />

distribution <strong>of</strong> Z is defined by Pr( Z =0)= p and Pr( Z =j)=(1−c)Pr(Z=j) for all<br />

c<br />

RA<br />

j>0, i.e. the distribution <strong>of</strong> Z has an expanded zero class relative to Z, while for non-<br />

zero values its density is simply reduced by a factor (1−c) from Z~NegB(R0,k). Hence,<br />

the <strong>of</strong>fspring distribution under RA control has pgf:<br />

g<br />

c<br />

RA<br />

c<br />

⎛ R<br />

( s)<br />

= c +<br />

1<br />

⎝ k<br />

0<br />

0<br />

( 1−<br />

c)<br />

⎜1+<br />

( − s)<br />

140<br />

⎞<br />

⎟<br />

⎠<br />

c<br />

−k<br />

0

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

Saved successfully!

Ooh no, something went wrong!