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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

integrated with theoretical modeling to demonstrate the universality and practical<br />

relevance <strong>of</strong> individual variation in infectiousness.<br />

2. Evidence <strong>of</strong> individual-level variation<br />

Let Z represent the number <strong>of</strong> secondary cases caused by a given infectious<br />

individual in an outbreak. Z is a random variable with an “<strong>of</strong>fspring distribution”<br />

Pr(Z=k) that incorporates influences <strong>of</strong> stochasticity and variation in the population.<br />

Stochastic effects in transmission are modeled using a Poisson process, as is<br />

conventional (Diekmann and Heesterbeek 2000), with intensity given by the individual<br />

reproductive number, ν i.e. Z~Poisson(ν). In conventional models neglecting individual<br />

variation, all individuals are characterized by the population mean, yielding<br />

Z~Poisson(R0). Another common approach, motivated by models with constant<br />

recovery or death rates and homogeneous transmission rates, is to assume that ν is<br />

exponentially distributed, yielding Z~geometric(R0). (For clarity we express all<br />

distributions using the mean as the scale parameter; see appendix A for correspondence<br />

with standard notation.) We introduce a more general formulation, in which ν follows a<br />

gamma distribution with mean R0 and shape parameter k, yielding Z~negative<br />

binomial(R0,k) (henceforth abbreviated Z~NegB(R0,k)) (Taylor and Karlin 1998). The<br />

negative binomial model includes the Poisson (k→∞) and geometric (k=1) models as<br />

special cases.<br />

Realized <strong>of</strong>fspring distributions can be determined from detailed contact tracing<br />

<strong>of</strong> particular outbreaks, or from surveillance data covering multiple introductions <strong>of</strong> a<br />

disease. The above candidate models can then be challenged with these data using<br />

99

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

Saved successfully!

Ooh no, something went wrong!