Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Comparative estimation of the reproduction number for pandemic influenza from daily case notification data</title>
<author><name sortKey="Chowell, Gerardo" sort="Chowell, Gerardo" uniqKey="Chowell G" first="Gerardo" last="Chowell">Gerardo Chowell</name>
<affiliation><nlm:aff id="aff1"><institution>Theoretical Division (MS B284), Los Alamos National Laboratory</institution>
<addr-line>Los Alamos, NM 87545, USA</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Nishiura, Hiroshi" sort="Nishiura, Hiroshi" uniqKey="Nishiura H" first="Hiroshi" last="Nishiura">Hiroshi Nishiura</name>
<affiliation><nlm:aff id="aff2"><institution>Department of Medical Biometry, University of Tübingen</institution>
<addr-line>Westbahnhofstrasse 55, Tübingen 72070, Germany</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><institution>Research Centre for Tropical Infectious Diseases, Nagasaki University Institute of Tropical Medicine</institution>
<addr-line>1-12-4 Sakamoto, Nagasaki 852-8523, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bettencourt, Luis M A" sort="Bettencourt, Luis M A" uniqKey="Bettencourt L" first="Luís M. A" last="Bettencourt">Luís M. A Bettencourt</name>
<affiliation><nlm:aff id="aff1"><institution>Theoretical Division (MS B284), Los Alamos National Laboratory</institution>
<addr-line>Los Alamos, NM 87545, USA</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">17254982</idno>
<idno type="pmc">2358966</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2358966</idno>
<idno type="RBID">PMC:2358966</idno>
<idno type="doi">10.1098/rsif.2006.0161</idno>
<date when="2006">2006</date>
<idno type="wicri:Area/Pmc/Corpus">000521</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000521</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Comparative estimation of the reproduction number for pandemic influenza from daily case notification data</title>
<author><name sortKey="Chowell, Gerardo" sort="Chowell, Gerardo" uniqKey="Chowell G" first="Gerardo" last="Chowell">Gerardo Chowell</name>
<affiliation><nlm:aff id="aff1"><institution>Theoretical Division (MS B284), Los Alamos National Laboratory</institution>
<addr-line>Los Alamos, NM 87545, USA</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Nishiura, Hiroshi" sort="Nishiura, Hiroshi" uniqKey="Nishiura H" first="Hiroshi" last="Nishiura">Hiroshi Nishiura</name>
<affiliation><nlm:aff id="aff2"><institution>Department of Medical Biometry, University of Tübingen</institution>
<addr-line>Westbahnhofstrasse 55, Tübingen 72070, Germany</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><institution>Research Centre for Tropical Infectious Diseases, Nagasaki University Institute of Tropical Medicine</institution>
<addr-line>1-12-4 Sakamoto, Nagasaki 852-8523, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bettencourt, Luis M A" sort="Bettencourt, Luis M A" uniqKey="Bettencourt L" first="Luís M. A" last="Bettencourt">Luís M. A Bettencourt</name>
<affiliation><nlm:aff id="aff1"><institution>Theoretical Division (MS B284), Los Alamos National Laboratory</institution>
<addr-line>Los Alamos, NM 87545, USA</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Journal of the Royal Society Interface</title>
<idno type="ISSN">1742-5689</idno>
<idno type="eISSN">1742-5662</idno>
<imprint><date when="2006">2006</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>The reproduction number, <italic>R</italic>
, defined as the average number of secondary cases generated by a primary case, is a crucial quantity for identifying the intensity of interventions required to control an epidemic. Current estimates of the reproduction number for seasonal influenza show wide variation and, in particular, uncertainty bounds for <italic>R</italic>
for the pandemic strain from 1918 to 1919 have been obtained only in a few recent studies and are yet to be fully clarified. Here, we estimate <italic>R</italic>
using daily case notifications during the autumn wave of the influenza pandemic (Spanish flu) in the city of San Francisco, California, from 1918 to 1919. In order to elucidate the effects from adopting different estimation approaches, four different methods are used: estimation of <italic>R</italic>
using the early exponential-growth rate (Method 1), a simple susceptible–exposed–infectious–recovered (SEIR) model (Method 2), a more complex SEIR-type model that accounts for asymptomatic and hospitalized cases (Method 3), and a stochastic susceptible–infectious–removed (SIR) with Bayesian estimation (Method 4) that determines the effective reproduction number <inline-formula><inline-graphic xlink:href="rsif20060161e18.jpg" alternate-form-of="M18"></inline-graphic>
<mml:math id="M18"><mml:mrow><mml:msub><mml:mi mathvariant="script">R</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
at a given time <italic>t</italic>
. The first three methods fit the initial exponential-growth phase of the epidemic, which was explicitly determined by the goodness-of-fit test. Moreover, Method 3 was also fitted to the whole epidemic curve. Whereas the values of <italic>R</italic>
obtained using the first three methods based on the initial growth phase were estimated to be 2.98 (95% confidence interval (CI): 2.73, 3.25), 2.38 (2.16, 2.60) and 2.20 (1.55, 2.84), the third method with the entire epidemic curve yielded a value of 3.53 (3.45, 3.62). This larger value could be an overestimate since the goodness-of-fit to the initial exponential phase worsened when we fitted the model to the entire epidemic curve, and because the model is established as an autonomous system without time-varying assumptions. These estimates were shown to be robust to parameter uncertainties, but the theoretical exponential-growth approximation (Method 1) shows wide uncertainty. Method 4 provided a maximum-likelihood effective reproduction number 2.10 (1.21, 2.95) using the first 17 epidemic days, which is consistent with estimates obtained from the other methods and an estimate of 2.36 (2.07, 2.65) for the entire autumn wave. We conclude that the reproduction number for pandemic influenza (Spanish flu) at the city level can be robustly assessed to lie in the range of 2.0–3.0, in broad agreement with previous estimates using distinct data.</p>
</div>
</front>
</TEI>
<pmc article-type="research-article" xml:lang="EN"><pmc-comment>The publisher of this article does not allow downloading of the full text in XML form.</pmc-comment>
<front><journal-meta><journal-id journal-id-type="nlm-ta">J R Soc Interface</journal-id>
<journal-id journal-id-type="publisher-id">RSIF</journal-id>
<journal-title>Journal of the Royal Society Interface</journal-title>
<issn pub-type="ppub">1742-5689</issn>
<issn pub-type="epub">1742-5662</issn>
<publisher><publisher-name>The Royal Society</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">17254982</article-id>
<article-id pub-id-type="pmc">2358966</article-id>
<article-id pub-id-type="publisher-id">rsif20060161</article-id>
<article-id pub-id-type="doi">10.1098/rsif.2006.0161</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Comparative estimation of the reproduction number for pandemic influenza from daily case notification data</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Chowell</surname>
<given-names>Gerardo</given-names>
</name>
<xref ref-type="aff" rid="aff1">1</xref>
<xref ref-type="corresp" rid="cor1">*</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Nishiura</surname>
<given-names>Hiroshi</given-names>
</name>
<xref ref-type="aff" rid="aff2">2</xref>
<xref ref-type="aff" rid="aff3">3</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Bettencourt</surname>
<given-names>Luís M.A</given-names>
</name>
<xref ref-type="aff" rid="aff1">1</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<institution>Theoretical Division (MS B284), Los Alamos National Laboratory</institution>
<addr-line>Los Alamos, NM 87545, USA</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<institution>Department of Medical Biometry, University of Tübingen</institution>
<addr-line>Westbahnhofstrasse 55, Tübingen 72070, Germany</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<institution>Research Centre for Tropical Infectious Diseases, Nagasaki University Institute of Tropical Medicine</institution>
<addr-line>1-12-4 Sakamoto, Nagasaki 852-8523, Japan</addr-line>
</aff>
<author-notes><corresp id="cor1"><label>*</label>
Author for correspondence (<email>chowell@lanl.gov</email>
)</corresp>
</author-notes>
<pub-date pub-type="epub"><day>12</day>
<month>10</month>
<year>2006</year>
</pub-date>
<pub-date pub-type="ppub"><day>22</day>
<month>2</month>
<year>2007</year>
</pub-date>
<volume>4</volume>
<issue>12</issue>
<fpage>155</fpage>
<lpage>166</lpage>
<history><date date-type="received"><day>11</day>
<month>8</month>
<year>2006</year>
</date>
<date date-type="accepted"><day>8</day>
<month>9</month>
<year>2006</year>
</date>
</history>
<permissions><copyright-statement>© 2006 The Royal Society</copyright-statement>
<copyright-year>2006</copyright-year>
</permissions>
<abstract xml:lang="EN"><p>The reproduction number, <italic>R</italic>
, defined as the average number of secondary cases generated by a primary case, is a crucial quantity for identifying the intensity of interventions required to control an epidemic. Current estimates of the reproduction number for seasonal influenza show wide variation and, in particular, uncertainty bounds for <italic>R</italic>
for the pandemic strain from 1918 to 1919 have been obtained only in a few recent studies and are yet to be fully clarified. Here, we estimate <italic>R</italic>
using daily case notifications during the autumn wave of the influenza pandemic (Spanish flu) in the city of San Francisco, California, from 1918 to 1919. In order to elucidate the effects from adopting different estimation approaches, four different methods are used: estimation of <italic>R</italic>
using the early exponential-growth rate (Method 1), a simple susceptible–exposed–infectious–recovered (SEIR) model (Method 2), a more complex SEIR-type model that accounts for asymptomatic and hospitalized cases (Method 3), and a stochastic susceptible–infectious–removed (SIR) with Bayesian estimation (Method 4) that determines the effective reproduction number <inline-formula><inline-graphic xlink:href="rsif20060161e18.jpg" alternate-form-of="M18"></inline-graphic>
<mml:math id="M18"><mml:mrow><mml:msub><mml:mi mathvariant="script">R</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
at a given time <italic>t</italic>
. The first three methods fit the initial exponential-growth phase of the epidemic, which was explicitly determined by the goodness-of-fit test. Moreover, Method 3 was also fitted to the whole epidemic curve. Whereas the values of <italic>R</italic>
obtained using the first three methods based on the initial growth phase were estimated to be 2.98 (95% confidence interval (CI): 2.73, 3.25), 2.38 (2.16, 2.60) and 2.20 (1.55, 2.84), the third method with the entire epidemic curve yielded a value of 3.53 (3.45, 3.62). This larger value could be an overestimate since the goodness-of-fit to the initial exponential phase worsened when we fitted the model to the entire epidemic curve, and because the model is established as an autonomous system without time-varying assumptions. These estimates were shown to be robust to parameter uncertainties, but the theoretical exponential-growth approximation (Method 1) shows wide uncertainty. Method 4 provided a maximum-likelihood effective reproduction number 2.10 (1.21, 2.95) using the first 17 epidemic days, which is consistent with estimates obtained from the other methods and an estimate of 2.36 (2.07, 2.65) for the entire autumn wave. We conclude that the reproduction number for pandemic influenza (Spanish flu) at the city level can be robustly assessed to lie in the range of 2.0–3.0, in broad agreement with previous estimates using distinct data.</p>
</abstract>
<kwd-group><kwd>Spanish flu</kwd>
<kwd>pandemic</kwd>
<kwd>influenza</kwd>
<kwd>reproduction number</kwd>
<kwd>San Francisco</kwd>
</kwd-group>
</article-meta>
</front>
<floats-wrap><fig id="fig1" position="float"><label>Figure 1</label>
<caption><p>Daily number of influenza notifications in San Francisco, California during the 1918–1919 influenza pandemic (<xref ref-type="bibr" rid="bib35">Department of Hygiene 1922</xref>
).</p>
</caption>
<graphic xlink:href="rsif20060161f01"></graphic>
</fig>
<fig id="fig2" position="float"><label>Figure 2</label>
<caption><p>The course of the outbreak can be visualized in an epidemic time-delay diagram of new cases Δ<italic>C</italic>
at consecutive times (black dots). For data that are not too stochastic, this provides a very simple method to estimate <italic>R</italic>
, via the tangent at the origin (dashed lines) of the initial growth trajectory (grey arrows). Jumps in case numbers (indicated for 22–23 Oct) lead to greater uncertainty in the estimation of the reproduction number.</p>
</caption>
<graphic xlink:href="rsif20060161f02"></graphic>
</fig>
<fig id="fig3" position="float"><label>Figure 3</label>
<caption><p>The <italic>Χ</italic>
<sup>2</sup>
density provided by three different models (Methods 1–3) to the initial epidemic growth phase of the cumulative number of influenza notifications as a function of the length of the initial epidemic phase. Using the goodness-of-fit statistic, the initial growth phase is predicted to be 5 days by Method 1 and 17 days by Methods 2 and 3.</p>
</caption>
<graphic xlink:href="rsif20060161f03"></graphic>
</fig>
<fig id="fig4" position="float"><label>Figure 4</label>
<caption><p>Model fits, residuals plots and the resulting distributions of the reproduction number obtained after fitting the simple SEIR epidemic model (Method 2) to the initial phase of the autumn influenza wave using 5 and 17 epidemic days of the Spanish flu pandemic in San Francisco, California. In the top panel, the epidemic data of the cumulative number of reported influenza cases are the circles, the solid line is the model best-fit and the solid grey lines are 1000 realizations of the model fit to the data obtained through parametric bootstrapping as explained in the text.</p>
</caption>
<graphic xlink:href="rsif20060161f04"></graphic>
</fig>
<fig id="fig5" position="float"><label>Figure 5</label>
<caption><p>Model fits, residuals plots and the resulting distributions of the reproduction number and the proportion of the clinical reporting obtained after fitting the complex SEIR epidemic model (Method 3) to the initial phase of the autumn influenza wave using 17 epidemic days of the Spanish flu pandemic in San Francisco, California. In the top panel, the epidemic data of the cumulative number of reported influenza cases are the circles, the solid line is the model best-fit and the solid grey lines are 1000 realizations of the model fit to the data obtained through parametric bootstrapping as explained in the text.</p>
</caption>
<graphic xlink:href="rsif20060161f05"></graphic>
</fig>
<fig id="fig6" position="float"><label>Figure 6</label>
<caption><p>Sequential Bayesian estimation of the full distribution of <inline-formula><inline-graphic xlink:href="rsif20060161e03.jpg" alternate-form-of="M3"></inline-graphic>
<mml:math id="M3"><mml:mrow><mml:msub><mml:mi mathvariant="script">R</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
leads to the estimation of its maximum-likelihood value (grey dots) and 95% CIs (solid black lines). Uncertainty, measured by the width of the CI, decreases with more case observations. The estimates eventually lead to smaller <inline-formula><inline-graphic xlink:href="rsif20060161e04.jpg" alternate-form-of="M4"></inline-graphic>
<mml:math id="M4"><mml:mrow><mml:msub><mml:mi mathvariant="script">R</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
owing to depletion of susceptibles. At late times, <inline-formula><inline-graphic xlink:href="rsif20060161e05.jpg" alternate-form-of="M5"></inline-graphic>
<mml:math id="M5"><mml:mrow><mml:msub><mml:mi mathvariant="script">R</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>
as a result of averaging periods in which the epidemic grows and declines.</p>
</caption>
<graphic xlink:href="rsif20060161f06"></graphic>
</fig>
<fig id="fig7" position="float"><label>Figure 7</label>
<caption><p>The complex SEIR model fit (solid line) to the entire epidemic curve (circles) and the simple SEIR model fit (dashed line) calibrated using the initial exponential phase (17 days) are shown for comparison (cumulative cases are shown in logarithmic scale). Solid grey lines are 1000 realizations of the complex SEIR model fit to the data obtained through parametric bootstrapping as explained in the text.</p>
</caption>
<graphic xlink:href="rsif20060161f07"></graphic>
</fig>
<table-wrap id="tbl1" position="float"><label>Table 1</label>
<caption><p>Estimates of the reproduction number for the autumn wave of the Spanish flu pandemic in San Francisco, California. n.a., not applicable. The number of data points is smaller than the number of parameters being estimated (seven parameters for the complex SEIR model). Note that the stochastic SIR method provides the effective reproduction number at time <italic>t</italic>
, while the other methods estimate the reproduction number by fitting the models to a specified number of epidemic days of data. The number of degrees of freedom (d.f.) is different by method. Initial growth rate, simple SEIR and complex SEIR estimate 1, 2 and 7 parameters, respectively. d.f. was determined by the difference between the observed number of epidemic days, <italic>n</italic>
, and the number of parameters to be estimated (e.g. for the complex SEIR, d.f. at 17 days was <italic>n</italic>
−7=10).</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th valign="bottom" rowspan="1" colspan="1"></th>
<th valign="bottom" colspan="2" align="left" rowspan="1">initial growth rate</th>
<th valign="bottom" colspan="2" align="left" rowspan="1">simple SEIR</th>
<th valign="bottom" colspan="2" align="left" rowspan="1">complex SEIR</th>
<th valign="bottom" colspan="2" align="left" rowspan="1">stochastic SIR</th>
</tr>
<tr><th valign="bottom" rowspan="1" colspan="1"></th>
<th valign="bottom" colspan="2" rowspan="1"><hr></hr>
</th>
<th valign="bottom" colspan="2" rowspan="1"><hr></hr>
</th>
<th valign="bottom" colspan="2" rowspan="1"><hr></hr>
</th>
<th valign="bottom" colspan="2" rowspan="1"><hr></hr>
</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">epidemic days</td>
<td rowspan="1" colspan="1"><italic>R</italic>
</td>
<td rowspan="1" colspan="1"><italic>R</italic>
95% CI</td>
<td rowspan="1" colspan="1"><italic>R</italic>
</td>
<td rowspan="1" colspan="1"><italic>R</italic>
95% CI</td>
<td rowspan="1" colspan="1"><italic>R</italic>
</td>
<td rowspan="1" colspan="1"><italic>R</italic>
95% CI</td>
<td rowspan="1" colspan="1"><inline-formula><inline-graphic xlink:href="rsif20060161e12.jpg" alternate-form-of="M12"></inline-graphic>
<mml:math id="M12"><mml:mrow><mml:msub><mml:mi mathvariant="script">R</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><inline-formula><inline-graphic xlink:href="rsif20060161e13.jpg" alternate-form-of="M13"></inline-graphic>
<mml:math id="M13"><mml:mrow><mml:msub><mml:mi mathvariant="script">R</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
95% CI</td>
</tr>
<tr><td rowspan="1" colspan="1">5</td>
<td align="char" char="." rowspan="1" colspan="1">5.78</td>
<td align="char" char="." rowspan="1" colspan="1">(3.80, 8.15)</td>
<td align="char" char="." rowspan="1" colspan="1">3.72</td>
<td align="char" char="." rowspan="1" colspan="1">(2.01, 5.44)</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="char" char="." rowspan="1" colspan="1">1.96</td>
<td align="char" char="." rowspan="1" colspan="1">(0.83, 3.09)</td>
</tr>
<tr><td rowspan="1" colspan="1">17</td>
<td align="char" char="." rowspan="1" colspan="1">2.98</td>
<td align="char" char="." rowspan="1" colspan="1">(2.73, 3.25)</td>
<td align="char" char="." rowspan="1" colspan="1">2.38</td>
<td align="char" char="." rowspan="1" colspan="1">(2.16, 2.60)</td>
<td align="char" char="." rowspan="1" colspan="1">2.20</td>
<td align="char" char="." rowspan="1" colspan="1">(1.55, 2.84)</td>
<td align="char" char="." rowspan="1" colspan="1">2.10</td>
<td align="char" char="." rowspan="1" colspan="1">(1.21, 2.95)</td>
</tr>
</tbody>
</table>
</table-wrap>
</floats-wrap>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/PandemieGrippaleV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000521 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000521 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= PandemieGrippaleV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.34. | ![]() |