Real-time Estimates in Early Detection of SARS
Identifieur interne : 000F80 ( Pmc/Corpus ); précédent : 000F79; suivant : 000F81Real-time Estimates in Early Detection of SARS
Auteurs : Simon Cauchemez ; Pierre-Yves Boëlle ; Christl A. Donnelly ; Neil M. Ferguson ; Guy Thomas ; Gabriel M. Leung ; Anthony J. Hedley ; Roy M. Anderson ; Alain-Jacques ValleronSource :
- Emerging Infectious Diseases [ 1080-6040 ] ; 2006.
Abstract
A statistical method can be used for early monitoring of the effect of disease control measures.
Url:
DOI: 10.3201/eid1201.050593
PubMed: 16494726
PubMed Central: 3293464
Links to Exploration step
PMC:3293464Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Real-time Estimates in Early Detection of SARS</title>
<author><name sortKey="Cauchemez, Simon" sort="Cauchemez, Simon" uniqKey="Cauchemez S" first="Simon" last="Cauchemez">Simon Cauchemez</name>
<affiliation><nlm:aff id="aff1">Institut National de la Santé et de la Recherche Médicale, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2">Université Pierre et Marie Curie, Paris, France;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Boelle, Pierre Yves" sort="Boelle, Pierre Yves" uniqKey="Boelle P" first="Pierre-Yves" last="Boëlle">Pierre-Yves Boëlle</name>
<affiliation><nlm:aff id="aff1">Institut National de la Santé et de la Recherche Médicale, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2">Université Pierre et Marie Curie, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Assistance Publique–Hôpitaux de Paris, Paris, France;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Donnelly, Christl A" sort="Donnelly, Christl A" uniqKey="Donnelly C" first="Christl A." last="Donnelly">Christl A. Donnelly</name>
<affiliation><nlm:aff id="aff4">Imperial College, London, United Kingdom;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ferguson, Neil M" sort="Ferguson, Neil M" uniqKey="Ferguson N" first="Neil M" last="Ferguson">Neil M. Ferguson</name>
<affiliation><nlm:aff id="aff4">Imperial College, London, United Kingdom;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Thomas, Guy" sort="Thomas, Guy" uniqKey="Thomas G" first="Guy" last="Thomas">Guy Thomas</name>
<affiliation><nlm:aff id="aff1">Institut National de la Santé et de la Recherche Médicale, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2">Université Pierre et Marie Curie, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Assistance Publique–Hôpitaux de Paris, Paris, France;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Leung, Gabriel M" sort="Leung, Gabriel M" uniqKey="Leung G" first="Gabriel M." last="Leung">Gabriel M. Leung</name>
<affiliation><nlm:aff id="aff5">University of Hong Kong, Hong Kong Special Administrative Region, People's Republic of China</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hedley, Anthony J" sort="Hedley, Anthony J" uniqKey="Hedley A" first="Anthony J" last="Hedley">Anthony J. Hedley</name>
<affiliation><nlm:aff id="aff5">University of Hong Kong, Hong Kong Special Administrative Region, People's Republic of China</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Anderson, Roy M" sort="Anderson, Roy M" uniqKey="Anderson R" first="Roy M." last="Anderson">Roy M. Anderson</name>
<affiliation><nlm:aff id="aff4">Imperial College, London, United Kingdom;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Valleron, Alain Jacques" sort="Valleron, Alain Jacques" uniqKey="Valleron A" first="Alain-Jacques" last="Valleron">Alain-Jacques Valleron</name>
<affiliation><nlm:aff id="aff1">Institut National de la Santé et de la Recherche Médicale, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2">Université Pierre et Marie Curie, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Assistance Publique–Hôpitaux de Paris, Paris, France;</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">16494726</idno>
<idno type="pmc">3293464</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3293464</idno>
<idno type="RBID">PMC:3293464</idno>
<idno type="doi">10.3201/eid1201.050593</idno>
<date when="2006">2006</date>
<idno type="wicri:Area/Pmc/Corpus">000F80</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000F80</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Real-time Estimates in Early Detection of SARS</title>
<author><name sortKey="Cauchemez, Simon" sort="Cauchemez, Simon" uniqKey="Cauchemez S" first="Simon" last="Cauchemez">Simon Cauchemez</name>
<affiliation><nlm:aff id="aff1">Institut National de la Santé et de la Recherche Médicale, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2">Université Pierre et Marie Curie, Paris, France;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Boelle, Pierre Yves" sort="Boelle, Pierre Yves" uniqKey="Boelle P" first="Pierre-Yves" last="Boëlle">Pierre-Yves Boëlle</name>
<affiliation><nlm:aff id="aff1">Institut National de la Santé et de la Recherche Médicale, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2">Université Pierre et Marie Curie, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Assistance Publique–Hôpitaux de Paris, Paris, France;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Donnelly, Christl A" sort="Donnelly, Christl A" uniqKey="Donnelly C" first="Christl A." last="Donnelly">Christl A. Donnelly</name>
<affiliation><nlm:aff id="aff4">Imperial College, London, United Kingdom;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ferguson, Neil M" sort="Ferguson, Neil M" uniqKey="Ferguson N" first="Neil M" last="Ferguson">Neil M. Ferguson</name>
<affiliation><nlm:aff id="aff4">Imperial College, London, United Kingdom;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Thomas, Guy" sort="Thomas, Guy" uniqKey="Thomas G" first="Guy" last="Thomas">Guy Thomas</name>
<affiliation><nlm:aff id="aff1">Institut National de la Santé et de la Recherche Médicale, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2">Université Pierre et Marie Curie, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Assistance Publique–Hôpitaux de Paris, Paris, France;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Leung, Gabriel M" sort="Leung, Gabriel M" uniqKey="Leung G" first="Gabriel M." last="Leung">Gabriel M. Leung</name>
<affiliation><nlm:aff id="aff5">University of Hong Kong, Hong Kong Special Administrative Region, People's Republic of China</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hedley, Anthony J" sort="Hedley, Anthony J" uniqKey="Hedley A" first="Anthony J" last="Hedley">Anthony J. Hedley</name>
<affiliation><nlm:aff id="aff5">University of Hong Kong, Hong Kong Special Administrative Region, People's Republic of China</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Anderson, Roy M" sort="Anderson, Roy M" uniqKey="Anderson R" first="Roy M." last="Anderson">Roy M. Anderson</name>
<affiliation><nlm:aff id="aff4">Imperial College, London, United Kingdom;</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Valleron, Alain Jacques" sort="Valleron, Alain Jacques" uniqKey="Valleron A" first="Alain-Jacques" last="Valleron">Alain-Jacques Valleron</name>
<affiliation><nlm:aff id="aff1">Institut National de la Santé et de la Recherche Médicale, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2">Université Pierre et Marie Curie, Paris, France;</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Assistance Publique–Hôpitaux de Paris, Paris, France;</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Emerging Infectious Diseases</title>
<idno type="ISSN">1080-6040</idno>
<idno type="eISSN">1080-6059</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>A statistical method can be used for early monitoring of the effect of disease control measures.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Emerg Infect Dis</journal-id>
<journal-id journal-id-type="publisher-id">EID</journal-id>
<journal-title-group><journal-title>Emerging Infectious Diseases</journal-title>
</journal-title-group>
<issn pub-type="ppub">1080-6040</issn>
<issn pub-type="epub">1080-6059</issn>
<publisher><publisher-name>Centers for Disease Control and Prevention</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">16494726</article-id>
<article-id pub-id-type="pmc">3293464</article-id>
<article-id pub-id-type="publisher-id">05-0593</article-id>
<article-id pub-id-type="doi">10.3201/eid1201.050593</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research</subject>
</subj-group>
<subj-group subj-group-type="article-type"><subject>Research</subject>
</subj-group>
</article-categories>
<title-group><article-title>Real-time Estimates in Early Detection of SARS</article-title>
<alt-title alt-title-type="running-head">Real-Time Estimates in Early Detection of SARS</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author" corresp="yes"><name><surname>Cauchemez</surname>
<given-names>Simon</given-names>
</name>
<xref ref-type="aff" rid="aff1">*</xref>
<xref ref-type="aff" rid="aff2">†</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Boëlle</surname>
<given-names>Pierre-Yves</given-names>
</name>
<xref ref-type="aff" rid="aff1">*</xref>
<xref ref-type="aff" rid="aff2">†</xref>
<xref ref-type="aff" rid="aff3">‡</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Donnelly</surname>
<given-names>Christl A.</given-names>
</name>
<xref ref-type="aff" rid="aff4">§</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ferguson</surname>
<given-names>Neil M</given-names>
</name>
<xref ref-type="aff" rid="aff4">§</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Thomas</surname>
<given-names>Guy</given-names>
</name>
<xref ref-type="aff" rid="aff1">*</xref>
<xref ref-type="aff" rid="aff2">†</xref>
<xref ref-type="aff" rid="aff3">‡</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Leung</surname>
<given-names>Gabriel M.</given-names>
</name>
<xref ref-type="aff" rid="aff5">¶</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Hedley</surname>
<given-names>Anthony J</given-names>
</name>
<xref ref-type="aff" rid="aff5">¶</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Anderson</surname>
<given-names>Roy M.</given-names>
</name>
<xref ref-type="aff" rid="aff4">§</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Valleron</surname>
<given-names>Alain-Jacques</given-names>
</name>
<xref ref-type="aff" rid="aff1">*</xref>
<xref ref-type="aff" rid="aff2">†</xref>
<xref ref-type="aff" rid="aff3">‡</xref>
</contrib>
<aff id="aff1"><label>*</label>
Institut National de la Santé et de la Recherche Médicale, Paris, France;</aff>
<aff id="aff2"><label>†</label>
Université Pierre et Marie Curie, Paris, France;</aff>
<aff id="aff3"><label>‡</label>
Assistance Publique–Hôpitaux de Paris, Paris, France;</aff>
<aff id="aff4"><label>§</label>
Imperial College, London, United Kingdom;</aff>
<aff id="aff5"><label>¶</label>
University of Hong Kong, Hong Kong Special Administrative Region, People's Republic of China</aff>
</contrib-group>
<author-notes><corresp id="cor1">Address for correspondence: Simon Cauchemez, Department of Infectious Disease Epidemiology, St Mary's Campus, Norfolk Place, London W2 1PG, United Kingdom; fax: 44-20-0759-43282; email: <email xlink:href="s.cauchemez@imperial.ac.uk">s.cauchemez@imperial.ac.uk</email>
</corresp>
</author-notes>
<pub-date pub-type="ppub"><month>1</month>
<year>2006</year>
</pub-date>
<volume>12</volume>
<issue>1</issue>
<fpage>110</fpage>
<lpage>113</lpage>
<abstract abstract-type="toc"><p>A statistical method can be used for early monitoring of the effect of disease control measures.</p>
</abstract>
<abstract><p>We propose a Bayesian statistical framework for estimating the reproduction number <italic>R</italic>
early in an epidemic. This method allows for the yet-unrecorded secondary cases if the estimate is obtained before the epidemic has ended. We applied our approach to the severe acute respiratory syndrome (SARS) epidemic that started in February 2003 in Hong Kong. Temporal patterns of <italic>R</italic>
estimated after 5, 10, and 20 days were similar. Ninety-five percent credible intervals narrowed when more data were available but stabilized after 10 days. Using simulation studies of SARS-like outbreaks, we have shown that the method may be used for early monitoring of the effect of control measures.</p>
</abstract>
<kwd-group kwd-group-type="author"><title>Keywords: </title>
<kwd>Severe acute respiratory syndrome</kwd>
<kwd>communicable diseases</kwd>
<kwd>emerging</kwd>
<kwd>disease outbreaks</kwd>
<kwd>epidemiologic methods</kwd>
<kwd>population surveillance</kwd>
<kwd>reproduction number</kwd>
<kwd>Hong Kong</kwd>
<kwd>research</kwd>
</kwd-group>
</article-meta>
</front>
<body><p>The reproduction number <italic>R</italic>
of an epidemic (the mean number of secondary cases infected by a single infectious case) is a key parameter for the analysis of infectious diseases because it summarizes the potential transmissibility of the disease and indicates whether an epidemic is under control (<italic>R</italic>
<1). Up to now, this parameter has only been estimated retrospectively for periods from which all secondary cases had been detected. In terms of policy development and evaluation during the epidemic, obtaining estimates of the temporal trends in the reproduction number relating to as recent a time as possible would be critical.</p>
<p>If all incident cases could be traced to their index cases, estimating the reproduction number would simply be a matter of counting secondary cases. However, if tracing information is incomplete or ambiguous, modeling or statistical approaches are required. For example, a mathematical model for disease transmission fitted to available data can provide estimates of <italic>R</italic>
(<xref ref-type="bibr" rid="R1"><italic>1</italic>
</xref>
). An approach requiring fewer assumptions has been proposed by Wallinga and Teunis (<xref ref-type="bibr" rid="R2"><italic>2</italic>
</xref>
), in which the distribution of the generation interval of the disease and the epidemic curve are directly analyzed and suffice to provide estimates. For an ongoing epidemic, this method could be used to estimate the number of secondary cases infected by a primary case-patient, but only for periods from which all secondary cases would have been detected. For severe acute respiratory syndrome (SARS), the required lag would be on the order of 15 days (95th percentile of the distribution of the generation interval described by Lipsitch et al.) (<xref ref-type="bibr" rid="R3"><italic>3</italic>
</xref>
).</p>
<p>In this report, we show how to estimate the reproduction number in an ongoing epidemic, which will account for yet unobserved secondary cases. The method is applied to data from the 2003 SARS outbreak in Hong Kong (<xref ref-type="bibr" rid="R4"><italic>4</italic>
</xref>
). Using simulated data, we demonstrate how the method may be used for early detection of the effect of control measures.</p>
<sec sec-type="materials|methods"><title>Materials and Methods</title>
<sec><title>Statistical Framework</title>
<p>We propose a Bayesian statistical framework for real-time inference on the temporal pattern of the reproduction number of an epidemic. Here, the reproduction number <italic>R<sub>t</sub>
</italic>
for day <italic>t</italic>
will be defined as the mean number of secondary cases infected by a case with symptom onset at day <italic>t</italic>
. Denoting <italic>n<sub>t</sub>
</italic>
as the number of cases with symptom onset at day <italic>t</italic>
and <italic>X<sub>t</sub>
</italic>
as the number of secondary cases they infected, the reproduction number <italic>R<sub>t</sub>
</italic>
is the ratio <italic>X<sub>t</sub>
</italic>
/<italic>n<sub>t</sub>
</italic>
, defined for <italic>n<sub>t</sub>
</italic>
>0.</p>
<p>Assume that we would like to compute the daily values <italic>R<sub>t</sub>
</italic>
from day 0 to present day <italic>T</italic>
, before the epidemic has ended. Although daily incident case counts can be known up to day <italic>T</italic>
, provided no delay in reporting occurs, the corresponding counts of secondary cases <italic>X<sub>t</sub>
</italic>
cannot. Secondary case-patients infected before day <italic>T</italic>
, whose illness had a long incubation time, may have clinical onset only after day <italic>T</italic>
. Furthermore, since the exact chain of transmission is seldom observed in practice, attributing secondary cases to previous cases is difficult. Focusing on these 2 issues, we show that the daily counts of symptom onset available until day <italic>T</italic>
are sufficient to estimate <italic>R<sub>t</sub>
</italic>
.</p>
<p>A 3-step construct is necessary. We first predict the eventual number of late secondary cases (as yet unobserved), for cases reported at day <italic>t</italic>
, assuming the number of early secondary cases (reported before day <italic>T</italic>
) is known. The method described by Wallinga and Teunis (<xref ref-type="bibr" rid="R2"><italic>2</italic>
</xref>
) is then used to estimate the number of early secondary cases from the daily counts of symptom onsets. These 2 steps are finally combined and yield an estimate of the predictive distribution of <italic>R<sub>t</sub>
</italic>
. Technical details are given in the Appendix. The estimation procedure depends on 3 assumptions: 1) ascertainment of patients whose symptoms appear before day <italic>T</italic>
is complete, 2) transmission events are independent, and 3) the generation interval, the time from symptom onset in a primary case to symptom onset in a secondary case, has a known frequency distribution.</p>
</sec>
<sec><title>Data from Hong Kong</title>
<p>The method was retrospectively used to analyze the SARS outbreak in Hong Kong. The data consisted of the dates of symptom onset of the 1,755 case-patients who were detected in Hong Kong in 2003 (<xref ref-type="bibr" rid="R4"><italic>4</italic>
</xref>
).</p>
</sec>
<sec><title>Simulated Data</title>
<p>Using simulations, we explored the ability of the method to quickly detect the effect of control measures. Five hundred epidemics were simulated with the following characteristics. During the first 20 days of the epidemics, the theoretical reproduction number was 3. Control measures were implemented at day 20. In a first scenario, control measures were completely effective (no transmission occurred after day 20). In a second scenario, the theoretical reproduction number after control measures were implemented was 0.7. Details on the simulations are available from the corresponding author.</p>
<p>In a simulation study, the bias and precision of the real-time estimator were investigated in situations in which the theoretical reproduction number remained constant with time. We also evaluated the effect of the length of the generation interval on the results. Detailed information can be obtained from the corresponding author.</p>
</sec>
</sec>
<sec><title>Results</title>
<sec><title>Application to Hong Kong SARS Data</title>
<p><xref ref-type="fig" rid="F1">Figure 1A</xref>
shows the dates of symptom onset of the 1,755 SARS patients detected in Hong Kong in 2003. <xref ref-type="fig" rid="F1">Figure 1B–F</xref>
shows the expectation and 95% credible intervals of the predictive distribution of <italic>R<sub>t</sub>
</italic>
based on data available at the end of the epidemic and after a lag of 2, 5, 10, and 20 days.</p>
<fig id="F1" fig-type="figure" position="float"><label>Figure 1</label>
<caption><p>Application of real-time estimation to the severe acute respiratory syndrome outbreak in Hong Kong. A) Data. B–F) Expectation (solid lines) and 95% credible intervals (dashed lines) of the real-time estimator of <italic>R<sub>t</sub>
</italic>
were calculated at the end of the epidemic (B) and after a lag of 2 (C), 5 (D), 10 (E), and 20 (F) days. The gray zones indicate that <italic>R</italic>
is <1.</p>
</caption>
<graphic xlink:href="05-0593-F1"></graphic>
</fig>
<p>After a lag of 2 days, the 95% credible intervals were wide and displayed an undesirable feature: they sharply decreased to 0 as soon as no cases had been observed for 2 consecutive days (<xref ref-type="fig" rid="F1">Figure 1C</xref>
; note especially days 1–4 and 13). After a 5-day lag, this undesirable feature had vanished (<xref ref-type="fig" rid="F1">Figure 1D</xref>
).</p>
<p>With lags >5 days, the trends of expected values were relatively similar, with a peak around day 20, a decreasing trend after this date, and the expectation of <italic>R<sub>t</sub>
</italic>
decreasing to <1 around day 40. These observations suggest that after a lag of only 5 days, the temporal trends in the expectation of <italic>R<sub>t</sub>
</italic>
are well captured. For a lag of 5 days, the credible interval of <italic>R<sub>t</sub>
</italic>
was wide when <20 cases were detected (periods 0<<italic>t</italic>
<20 and <italic>t</italic>
>63), but was relatively narrow when more cases were detected (period 21<<italic>t</italic>
<62). As expected, the width of the credible interval narrowed as the lag increased and more complete data were available. The expectations and credible intervals were very similar for lags of 10 and 20 days, 67.8th and 99.7th percentiles, respectively, of the distribution of the SARS generation interval described by Lipsitch et al. (<xref ref-type="bibr" rid="R3"><italic>3</italic>
</xref>
). No difference was detected between retrospective and 20-day estimates.</p>
</sec>
<sec><title>Detecting the Effect of Control Measures</title>
<p>In <xref ref-type="fig" rid="F2">Figure 2</xref>
, the method is used to estimate the impact of control measures implemented on day 20 in the simulated datasets with completely effective or limited control measures. The curves show the temporal pattern of <italic>R<sub>t</sub>
</italic>
based on an average over the 500 simulated datasets as a function of <italic>T</italic>
. Even when control measures are completely effective, based on data available up to day 21, the average expectation of <italic>R<sub>20</sub>
</italic>
is ≈3. Based on data available up to day 25, a downward trend is apparent, whereas based on data available up to day 29, the average expectation of <italic>R<sub>t</sub>
</italic>
is <1 from <italic>t</italic>
= 27 days. Based on data available up to day 40 (20 days after the implementation of the control measures), the estimates indicate that the threshold value 1 is crossed at day 22, which is 2 days after control measures were implemented. With limited control measures, the observed changes are qualitatively the same, although slightly more time is required for <italic>R<sub>t</sub>
</italic>
estimates to decrease to <1.</p>
<fig id="F2" fig-type="figure" position="float"><label>Figure 2</label>
<caption><p>Average expectation of the temporal pattern of <italic>R<sub>t</sub>
</italic>
after implementation of control measures according to the day <italic>T</italic>
of the last observation. A) Completely effective control measures. B) Limited control measures. Simulation values of <italic>R</italic>
are also given: before day 20, <italic>R</italic>
= 3; after day 20 <italic>R</italic>
= 0 (A) and <italic>R</italic>
= 0.7 (B). The gray zone indicates that <italic>R</italic>
is <1. Information that the average expectation of <italic>R</italic>
has passed <1 was obtained 9 (A) and 12 (B) days after control measures were implemented.</p>
</caption>
<graphic xlink:href="05-0593-F2"></graphic>
</fig>
</sec>
</sec>
<sec><title>Discussion</title>
<p>Our statistical framework provided real-time estimates of the reproduction number of an epidemic, and thus quickly showed the impact of control measures. In simulations of SARS-like diseases, the derived estimator detected the decrease of <italic>R<sub>t</sub>
</italic>
only 5 days after control measures were implemented. Furthermore, the average estimate had crossed the threshold value of 1 only 9 days after control measures were implemented.</p>
<p>In theory, the method could be applied to communicable diseases with the following characteristics: 1) no asymptomatic cases; 2) no underreporting; 3) knowledge of the generation interval. The list of communicable diseases that could be monitored is therefore relatively large, although it does not include diseases such as influenza, for which the proportion of asymptomatic or unreported cases may be large. In practice, the delay until estimates of the reproduction number become reliable will depend critically on the generation interval distribution. For SARS, when the reproduction number was constant over time, our real-time estimates were almost unbiased after only 1 day. With the original estimator of Wallinga and Teunis (<xref ref-type="bibr" rid="R2"><italic>2</italic>
</xref>
), which is not intended for real-time estimation, downward bias would be a concern for at least 2 weeks after observation. However, real-time estimates obtained for recent days displayed wide 95% credible intervals and zero-width intervals when no cases had been observed or reported for a few days. Here, owing to the relatively short generation interval of SARS (mean 8.4 days) (<xref ref-type="bibr" rid="R3"><italic>3</italic>
</xref>
), reliable estimates were obtained after only 5 days, albeit with wide credible intervals, and they were consolidated after 10 days. These lags corresponded to the 20th and 70th percentiles of the SARS generation interval (<xref ref-type="bibr" rid="R3"><italic>3</italic>
</xref>
). When the generation interval doubled, the time delay required to detect the effect of control measures implementation or to consolidate estimates roughly doubled.</p>
<p>We assumed that the distribution of the generation interval was known and remained unchanged during the course of the outbreak. In practice, however, this distribution is derived from a subset of traced cases. If the subset is small, e.g., the case at the beginning of an emerging disease outbreak, uncertainty will be large. Furthermore, the generation interval may decrease during the course of the outbreak because of quicker interventions, leading to possible bias in the estimates of <italic>R</italic>
(<xref ref-type="bibr" rid="R2"><italic>2</italic>
</xref>
). Further developments of the method could take these issues into consideration. For example, one could use information on traced cases as it accrues to sequentially estimate the generation interval. Depending on how cases are traced during the epidemic, changes in the generation interval could also be monitored.</p>
<p>The approach smoothed the temporal pattern of the reproduction number, leading to overestimation of <italic>R</italic>
in the week after control measures were implemented. We are trying to find a correction factor for this bias in ongoing research.</p>
<p>The method has a natural real-time implementation in which 1) a first estimate of the reproduction number is available after a lag that depends on the generation interval, and 2) while the epidemic goes on, the estimate is consolidated, and its credible interval narrows. Incorporation of such a statistical estimation framework into real-time surveillance of future infectious disease outbreaks would enhance the ability of epidemiologists to provide timely advice to public health policymakers.</p>
</sec>
<sec><title>Appendix</title>
<sec><title>Statistical Framework</title>
<p>Denoting <italic>n<sub>t</sub>
</italic>
the number of cases with onset at day <italic>t</italic>
and <italic>X<sub>t</sub>
</italic>
the number of cases they infected, the reproduction number <italic>R<sub>t</sub>
</italic>
is simply the ratio <italic>X<sub>t</sub>
/n<sub>t</sub>
</italic>
defined for <italic>n<sub>t</sub>
</italic>
>0. Here, we define a method to obtain the predictive distribution of <italic>R<sub>t</sub>
</italic>
given the available data at day <italic>T</italic>
, where data <italic>I(T) = {n<sub>t</sub>
}<sub>0< t </sub>
</italic>
are the daily counts of incident onsets, assuming that the density <italic>w(.)</italic>
of the generation interval is known. We will make use of the decomposition <italic>X<sub>t</sub>
</italic>
= <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
+ <italic>X<sub>t</sub>
<sup>+</sup>
(T)</italic>
, where the number of secondary cases <italic>X<sub>t</sub>
</italic>
from cases with onset at day <italic>t</italic>
has been split in those with onset before <italic>T</italic>
(<italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
) (early secondary cases), and those with onset after <italic>T</italic>
(<italic>X<sub>t</sub>
<sup>+</sup>
(T)</italic>
) (late secondary cases).</p>
<p>The construction of a global estimator is carried out in 3 stages. First, we consider the problem of right censoring, under the assumption that the exact chain of transmission has been observed until day <italic>T</italic>
. In this situation, <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
is observed while <italic>X<sub>t</sub>
<sup>+</sup>
(T)</italic>
is censored and must be predicted, conditional on <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
and <italic>n<sub>t</sub>
</italic>
, to allow computation of the predictive distribution of <italic>R<sub>t</sub>
</italic>
. Second, when the exact chain of transmission has not been observed, the number of early secondary cases, <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
, is not available. Following the recommendations of Wallinga and Teunis (<xref ref-type="bibr" rid="R1"><italic>1</italic>
</xref>
), we show that it is possible to compute the distribution of <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
given <italic>I(T)</italic>
. Finally, the conditional distributions of <italic>X<sub>t</sub>
<sup>+</sup>
(T)</italic>
given <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
, <italic>n<sub>t</sub>
</italic>
and <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
given <italic>I(T)</italic>
are combined to derive the distribution of the reproduction number <italic>R<sub>t</sub>
</italic>
conditional on <italic>I(T)</italic>
. All distributions presented are conditional to the number <italic>n<sub>t</sub>
</italic>
of symptom onsets at day <italic>t</italic>
, but notation is omitted for the sake of clarity.</p>
</sec>
<sec><title>Distribution of X<sub>t</sub>
<sup>+</sup>
(T) | X<sub>t</sub>
<sup>-</sup>
(T)</title>
<p>We assume that <italic>X<sub>t</sub>
</italic>
is Poisson distributed with mean <italic>n<sub>t</sub>
λ<sub>t</sub>
</italic>
and choose a vague gamma <italic>prior</italic>
distribution for λ<italic><sub>t</sub>
</italic>
with shape parameter α = 10<sup>-5</sup>
and rate β = 10<sup>-5</sup>
.</p>
<p>Conditional on <italic>X<sub>t</sub>
</italic>
, the number <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
of early secondary cases is binomial with parameters <italic>X<sub>t</sub>
</italic>
, <italic>W<sub>tT</sub>
</italic>
, where <italic>W<sub>tT</sub>
</italic>
is the probability that the generation interval is <<italic>T – t</italic>
. It follows that <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
| λ<italic><sub>t</sub>
</italic>
is Poisson distributed with mean <italic>n<sub>t</sub>
λ<sub>t</sub>
W<sub>tT</sub>
</italic>
. The same argument would show that <italic>X<sub>t</sub>
<sup>+</sup>
(T)</italic>
| λ<italic><sub>t</sub>
</italic>
is Poisson with mean <italic>n<sub>t</sub>
λ<sub>t</sub>
(1 – W<sub>tT</sub>
)</italic>
. Given λ<italic><sub>t</sub>
</italic>
, <italic>X<sub>t</sub>
<sup>+</sup>
(T)</italic>
and <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
are independent so that<disp-formula id="E"><inline-formula><inline-graphic xlink:href="05-0593-M1.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
With Bayes' theorem,<disp-formula id="E___1"><inline-formula><inline-graphic xlink:href="05-0593-M2.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
where<disp-formula id="E___2"><inline-formula><inline-graphic xlink:href="05-0593-M3.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
Eventually, we obtain that the distribution <italic>X<sub>t</sub>
<sup>+</sup>
(T)</italic>
| <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
= <italic>y</italic>
is negative binomial with parameters <italic>p = (n<sub>t</sub>
W<sub>tT</sub>
+ β)/(n<sub>t</sub>
+ β</italic>
), <italic>m</italic>
= <italic>y + α</italic>
and probability</p>
<disp-formula id="E___3"><inline-formula><inline-graphic xlink:href="05-0593-M4.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
</sec>
<sec><title>Distribution of X<sub>t</sub>
<sup>-</sup>
(T) | I(T)</title>
<p>In practice, the exact realization of <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
is unknown, and inference must be based on <italic>I(T)</italic>
alone. Wallinga and Teunis (<xref ref-type="bibr" rid="R5"><italic>5</italic>
</xref>
) have shown that the probability that a case detected at day <italic>k</italic>
<<italic>T</italic>
has been infected by a case detected at day <italic>t</italic>
<<italic>T</italic>
is<disp-formula id="E___4"><inline-formula><inline-graphic xlink:href="05-0593-M5.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
where <italic>1{.}</italic>
is the indicator function. The distribution of <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
given <italic>I(T)</italic>
is a sum of independent binomial distributions<disp-formula id="E___5"><italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
| <italic>I(T)</italic>
~ ∑<italic><sub>k </sub>
</italic>
Bin(<italic>n<sub>k</sub>
</italic>
,<italic>p<sub>tk</sub>
</italic>
)</disp-formula>
This probability distribution may be determined numerically.</p>
</sec>
<sec><title>Distribution of R<sub>t</sub>
| I(T)</title>
<p>Using the decomposition in early and late secondary cases, we obtain<disp-formula id="E___6"><inline-formula><inline-graphic xlink:href="05-0593-M6.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
After calculation, we find that the expectation and variance of <italic>X<sub>t</sub>
</italic>
| <italic>I(T)</italic>
are functions of the expectation and variance of <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
| <italic>I(T)</italic>
alone, derived with the method of Wallinga and Teunis:<disp-formula id="E___7"><inline-formula><inline-graphic xlink:href="05-0593-M7.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
<disp-formula id="E___8"><inline-formula><inline-graphic xlink:href="05-0593-M8.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
For the vague <italic>prior</italic>
we specified for λ<italic><sub>t</sub>
</italic>
, we obtain</p>
<disp-formula id="E___9"><inline-formula><inline-graphic xlink:href="05-0593-M9.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
<disp-formula id="E___10"><inline-formula><inline-graphic xlink:href="05-0593-M10.jpg"></inline-graphic>
</inline-formula>
</disp-formula>
<p>As expected, the average proportion of secondary cases detected before <italic>T</italic>
is <italic>W<sub>tT</sub>
</italic>
. The first term of the variance is related to our imperfect knowledge of the realization of <italic>X<sub>t</sub>
<sup>-</sup>
(T)</italic>
while the second term is related to the natural randomness of <italic>X<sub>t</sub>
<sup>+</sup>
(T)</italic>
. We stress that when the lag between day <italic>t</italic>
and day <italic>T</italic>
is large (i.e., <italic>W<sub>tT</sub>
</italic>
≈1), our estimates are similar to those of Wallinga and Teunis for complete epidemics.</p>
<p>Given <italic>X<sub>t</sub>
</italic>
| <italic>I(T)</italic>
, the derivation of the predictive distribution of the reproduction number <italic>R<sub>t</sub>
</italic>
is straightforward considering the deterministic relation <italic>R<sub>t</sub>
</italic>
= <italic>X<sub>t</sub>
</italic>
/<italic>n<sub>t</sub>
</italic>
.</p>
</sec>
</sec>
</body>
<back><fn-group><fn fn-type="citation"><p><italic>Suggested citation for this article</italic>
: Cauchemez S, Boëlle P-Y, Donnelly CA, Ferguson NM, Thomas G, Leung GM, et al. Real-time estimates in early detection of SARS. Emerg Infect Dis [serial on the Internet]. 2006 Jan [<italic>date cited</italic>
]. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.3201/eid1201.050593">http://dx.doi.org/10.3201/eid1201.050593</ext-link>
</p>
</fn>
</fn-group>
<ack><p>The work in Hong Kong was supported in part by a commissioned grant from the Research Fund for the Control of Infectious Diseases of the Health, Welfare and Food Bureau of the Hong Kong SAR Government.</p>
<p>Mr Cauchemez is a doctoral student at University Pierre et Marie Curie, Paris, France. He develops statistical methods to analyze transmission of infectious diseases using incomplete information.</p>
</ack>
<ref-list><title>References</title>
<ref id="R1"><label>1. </label>
<mixed-citation publication-type="journal"><string-name><surname>Riley</surname>
<given-names>S</given-names>
</string-name>
, <string-name><surname>Fraser</surname>
<given-names>C</given-names>
</string-name>
, <string-name><surname>Donnelly</surname>
<given-names>CA</given-names>
</string-name>
, <string-name><surname>Ghani</surname>
<given-names>AC</given-names>
</string-name>
, <string-name><surname>Abu-Raddad</surname>
<given-names>LJ</given-names>
</string-name>
, <string-name><surname>Hedley</surname>
<given-names>AJ</given-names>
</string-name>
, <etal></etal>
<article-title>Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions.</article-title>
<source>Science</source>
. <year>2003</year>
;<volume>300</volume>
:<fpage>1961</fpage>
–<lpage>6</lpage>
<pub-id pub-id-type="doi">10.1126/science.1086478</pub-id>
<pub-id pub-id-type="pmid">12766206</pub-id>
</mixed-citation>
</ref>
<ref id="R2"><label>2. </label>
<mixed-citation publication-type="journal"><string-name><surname>Wallinga</surname>
<given-names>J</given-names>
</string-name>
, <string-name><surname>Teunis</surname>
<given-names>P</given-names>
</string-name>
<article-title>Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures.</article-title>
<source>Am J Epidemiol</source>
. <year>2004</year>
;<volume>160</volume>
:<fpage>509</fpage>
–<lpage>16</lpage>
<pub-id pub-id-type="doi">10.1093/aje/kwh255</pub-id>
<pub-id pub-id-type="pmid">15353409</pub-id>
</mixed-citation>
</ref>
<ref id="R3"><label>3. </label>
<mixed-citation publication-type="journal"><string-name><surname>Lipsitch</surname>
<given-names>M</given-names>
</string-name>
, <string-name><surname>Cohen</surname>
<given-names>T</given-names>
</string-name>
, <string-name><surname>Cooper</surname>
<given-names>B</given-names>
</string-name>
, <string-name><surname>Robins</surname>
<given-names>JM</given-names>
</string-name>
, <string-name><surname>Ma</surname>
<given-names>S</given-names>
</string-name>
, <string-name><surname>James</surname>
<given-names>L</given-names>
</string-name>
, <etal></etal>
<article-title>Transmission dynamics and control of severe acute respiratory syndrome.</article-title>
<source>Science</source>
. <year>2003</year>
;<volume>300</volume>
:<fpage>1966</fpage>
–<lpage>70</lpage>
<pub-id pub-id-type="doi">10.1126/science.1086616</pub-id>
<pub-id pub-id-type="pmid">12766207</pub-id>
</mixed-citation>
</ref>
<ref id="R4"><label>4. </label>
<mixed-citation publication-type="journal"><string-name><surname>Leung</surname>
<given-names>GM</given-names>
</string-name>
, <string-name><surname>Hedley</surname>
<given-names>AJ</given-names>
</string-name>
, <string-name><surname>Ho</surname>
<given-names>LM</given-names>
</string-name>
, <string-name><surname>Chau</surname>
<given-names>P</given-names>
</string-name>
, <string-name><surname>Wong</surname>
<given-names>IO</given-names>
</string-name>
, <string-name><surname>Thach</surname>
<given-names>TQ</given-names>
</string-name>
, <etal></etal>
<article-title>The epidemiology of severe acute respiratory syndrome in the 2003 Hong Kong epidemic: an analysis of all 1,755 patients.</article-title>
<source>Ann Intern Med</source>
. <year>2004</year>
;<volume>141</volume>
:<fpage>662</fpage>
–<lpage>73</lpage>
<pub-id pub-id-type="pmid">15520422</pub-id>
</mixed-citation>
</ref>
<ref id="R5"><label>5. </label>
<mixed-citation publication-type="journal"><string-name><surname>Wallinga</surname>
<given-names>J</given-names>
</string-name>
, <string-name><surname>Teunis</surname>
<given-names>P</given-names>
</string-name>
<article-title>Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures.</article-title>
<source>Am J Epidemiol</source>
. <year>2004</year>
;<volume>160</volume>
:<fpage>509</fpage>
–<lpage>16</lpage>
<pub-id pub-id-type="doi">10.1093/aje/kwh255</pub-id>
<pub-id pub-id-type="pmid">15353409</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/SrasV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000F80 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000F80 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= SrasV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:3293464 |texte= Real-time Estimates in Early Detection of SARS }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:16494726" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a SrasV1
![]() | This area was generated with Dilib version V0.6.33. | ![]() |