The reproduction number Rt in structured and nonstructured populations

Tom Burr, Gerardo Chowell

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

Using daily counts of newly infected individuals, Wallinga and Teunis (WT) introduced a conceptually simple method to estimate the number of secondary cases per primary case (Rt) for a given day. The method requires an estimate of the generation interval probability density function (pdf), which specifies the probabilities for the times between symptom onset in a primary case and symptom onset in a corresponding secondary case. Other methods to estimate Rt are based on explicit models such as the SIR model; therefore, one might expect the WT method to be more robust to departures from SIR- type behavior. This paper uses simulated data to compare the quality of daily Rt estimates based on a SIR model to those using the WT method for both structured (classical SIR assumptions are violated) and nonstructured (classical SIR assumptions hold) populations. By using detailed simulations that record the infection day of each new infection and the donor-recipient identities, the true Rt and the generation interval pdf is known with negligible error. We find that the generation interval pdf is time dependent in all cases, which agrees with recent results reported elsewhere. We also find that the WT method performs essentially the same in the structured populations (except for a spatial network) as it does in the nonstructured population. And, the WT method does as well or better than a SIR-model based method in three of the four structured populations. Therefore, even if the contact patterns are heterogeneous as in the structured populations evaluated here, the WT method provides reasonable estimates of Rt, as does the SIR method.

Original languageEnglish (US)
Pages (from-to)239-259
Number of pages21
JournalMathematical Biosciences and Engineering
Volume6
Issue number2
DOIs
StatePublished - Apr 2009

Fingerprint

Reproduction number
Reproduction
Probability density function
Population
Interval Probability
SIR Model
Structured Populations
methodology
Estimate
Infection
Spatial Networks
infection
Count
Contact
Model-based

Keywords

  • Generation interval
  • Reproduction number
  • Structured population

ASJC Scopus subject areas

  • Applied Mathematics
  • Modeling and Simulation
  • Computational Mathematics
  • Agricultural and Biological Sciences(all)
  • Medicine(all)

Cite this

The reproduction number Rt in structured and nonstructured populations. / Burr, Tom; Chowell, Gerardo.

In: Mathematical Biosciences and Engineering, Vol. 6, No. 2, 04.2009, p. 239-259.

Research output: Contribution to journalArticle

Burr, Tom ; Chowell, Gerardo. / The reproduction number Rt in structured and nonstructured populations. In: Mathematical Biosciences and Engineering. 2009 ; Vol. 6, No. 2. pp. 239-259.
@article{1c7ea5b96a194b3a99dc9e0986220759,
title = "The reproduction number Rt in structured and nonstructured populations",
abstract = "Using daily counts of newly infected individuals, Wallinga and Teunis (WT) introduced a conceptually simple method to estimate the number of secondary cases per primary case (Rt) for a given day. The method requires an estimate of the generation interval probability density function (pdf), which specifies the probabilities for the times between symptom onset in a primary case and symptom onset in a corresponding secondary case. Other methods to estimate Rt are based on explicit models such as the SIR model; therefore, one might expect the WT method to be more robust to departures from SIR- type behavior. This paper uses simulated data to compare the quality of daily Rt estimates based on a SIR model to those using the WT method for both structured (classical SIR assumptions are violated) and nonstructured (classical SIR assumptions hold) populations. By using detailed simulations that record the infection day of each new infection and the donor-recipient identities, the true Rt and the generation interval pdf is known with negligible error. We find that the generation interval pdf is time dependent in all cases, which agrees with recent results reported elsewhere. We also find that the WT method performs essentially the same in the structured populations (except for a spatial network) as it does in the nonstructured population. And, the WT method does as well or better than a SIR-model based method in three of the four structured populations. Therefore, even if the contact patterns are heterogeneous as in the structured populations evaluated here, the WT method provides reasonable estimates of Rt, as does the SIR method.",
keywords = "Generation interval, Reproduction number, Structured population",
author = "Tom Burr and Gerardo Chowell",
year = "2009",
month = "4",
doi = "10.3934/mbe.2009.6.239",
language = "English (US)",
volume = "6",
pages = "239--259",
journal = "Mathematical Biosciences and Engineering",
issn = "1547-1063",
publisher = "Arizona State University",
number = "2",

}

TY - JOUR

T1 - The reproduction number Rt in structured and nonstructured populations

AU - Burr, Tom

AU - Chowell, Gerardo

PY - 2009/4

Y1 - 2009/4

N2 - Using daily counts of newly infected individuals, Wallinga and Teunis (WT) introduced a conceptually simple method to estimate the number of secondary cases per primary case (Rt) for a given day. The method requires an estimate of the generation interval probability density function (pdf), which specifies the probabilities for the times between symptom onset in a primary case and symptom onset in a corresponding secondary case. Other methods to estimate Rt are based on explicit models such as the SIR model; therefore, one might expect the WT method to be more robust to departures from SIR- type behavior. This paper uses simulated data to compare the quality of daily Rt estimates based on a SIR model to those using the WT method for both structured (classical SIR assumptions are violated) and nonstructured (classical SIR assumptions hold) populations. By using detailed simulations that record the infection day of each new infection and the donor-recipient identities, the true Rt and the generation interval pdf is known with negligible error. We find that the generation interval pdf is time dependent in all cases, which agrees with recent results reported elsewhere. We also find that the WT method performs essentially the same in the structured populations (except for a spatial network) as it does in the nonstructured population. And, the WT method does as well or better than a SIR-model based method in three of the four structured populations. Therefore, even if the contact patterns are heterogeneous as in the structured populations evaluated here, the WT method provides reasonable estimates of Rt, as does the SIR method.

AB - Using daily counts of newly infected individuals, Wallinga and Teunis (WT) introduced a conceptually simple method to estimate the number of secondary cases per primary case (Rt) for a given day. The method requires an estimate of the generation interval probability density function (pdf), which specifies the probabilities for the times between symptom onset in a primary case and symptom onset in a corresponding secondary case. Other methods to estimate Rt are based on explicit models such as the SIR model; therefore, one might expect the WT method to be more robust to departures from SIR- type behavior. This paper uses simulated data to compare the quality of daily Rt estimates based on a SIR model to those using the WT method for both structured (classical SIR assumptions are violated) and nonstructured (classical SIR assumptions hold) populations. By using detailed simulations that record the infection day of each new infection and the donor-recipient identities, the true Rt and the generation interval pdf is known with negligible error. We find that the generation interval pdf is time dependent in all cases, which agrees with recent results reported elsewhere. We also find that the WT method performs essentially the same in the structured populations (except for a spatial network) as it does in the nonstructured population. And, the WT method does as well or better than a SIR-model based method in three of the four structured populations. Therefore, even if the contact patterns are heterogeneous as in the structured populations evaluated here, the WT method provides reasonable estimates of Rt, as does the SIR method.

KW - Generation interval

KW - Reproduction number

KW - Structured population

UR - http://www.scopus.com/inward/record.url?scp=67650499116&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=67650499116&partnerID=8YFLogxK

U2 - 10.3934/mbe.2009.6.239

DO - 10.3934/mbe.2009.6.239

M3 - Article

VL - 6

SP - 239

EP - 259

JO - Mathematical Biosciences and Engineering

JF - Mathematical Biosciences and Engineering

SN - 1547-1063

IS - 2

ER -