The Impact of Prior Information on Estimates of Disease Transmissibility Using Bayesian Tools
Identifieur interne : 001343 ( Pmc/Corpus ); précédent : 001342; suivant : 001344The Impact of Prior Information on Estimates of Disease Transmissibility Using Bayesian Tools
Auteurs : Carlee B. Moser ; Mayetri Gupta ; Brett N. Archer ; Laura F. WhiteSource :
- PLoS ONE [ 1932-6203 ] ; 2015.
Abstract
The basic reproductive number (R₀) and the distribution of the serial interval (SI) are often used to quantify transmission during an infectious disease outbreak. In this paper, we present estimates of R₀ and SI from the 2003 SARS outbreak in Hong Kong and Singapore, and the 2009 pandemic influenza A(H1N1) outbreak in South Africa using methods that expand upon an existing Bayesian framework. This expanded framework allows for the incorporation of additional information, such as contact tracing or household data, through prior distributions. The results for the R₀ and the SI from the influenza outbreak in South Africa were similar regardless of the prior information (
Url:
DOI: 10.1371/journal.pone.0118762
PubMed: 25793993
PubMed Central: 4368801
Links to Exploration step
PMC:4368801Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">The Impact of Prior Information on Estimates of Disease Transmissibility Using Bayesian Tools</title>
<author><name sortKey="Moser, Carlee B" sort="Moser, Carlee B" uniqKey="Moser C" first="Carlee B." last="Moser">Carlee B. Moser</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Biostatistics, Boston University School of Public Health, Boston University, Boston, Massachusetts, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gupta, Mayetri" sort="Gupta, Mayetri" uniqKey="Gupta M" first="Mayetri" last="Gupta">Mayetri Gupta</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Biostatistics, Boston University School of Public Health, Boston University, Boston, Massachusetts, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Archer, Brett N" sort="Archer, Brett N" uniqKey="Archer B" first="Brett N." last="Archer">Brett N. Archer</name>
<affiliation><nlm:aff id="aff002"><addr-line>National Institute for Communicable Diseases, National Health Laboratory Service, Johannesburg, South Africa</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="White, Laura F" sort="White, Laura F" uniqKey="White L" first="Laura F." last="White">Laura F. White</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Biostatistics, Boston University School of Public Health, Boston University, Boston, Massachusetts, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">25793993</idno>
<idno type="pmc">4368801</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4368801</idno>
<idno type="RBID">PMC:4368801</idno>
<idno type="doi">10.1371/journal.pone.0118762</idno>
<date when="2015">2015</date>
<idno type="wicri:Area/Pmc/Corpus">001343</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">001343</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">The Impact of Prior Information on Estimates of Disease Transmissibility Using Bayesian Tools</title>
<author><name sortKey="Moser, Carlee B" sort="Moser, Carlee B" uniqKey="Moser C" first="Carlee B." last="Moser">Carlee B. Moser</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Biostatistics, Boston University School of Public Health, Boston University, Boston, Massachusetts, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gupta, Mayetri" sort="Gupta, Mayetri" uniqKey="Gupta M" first="Mayetri" last="Gupta">Mayetri Gupta</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Biostatistics, Boston University School of Public Health, Boston University, Boston, Massachusetts, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Archer, Brett N" sort="Archer, Brett N" uniqKey="Archer B" first="Brett N." last="Archer">Brett N. Archer</name>
<affiliation><nlm:aff id="aff002"><addr-line>National Institute for Communicable Diseases, National Health Laboratory Service, Johannesburg, South Africa</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="White, Laura F" sort="White, Laura F" uniqKey="White L" first="Laura F." last="White">Laura F. White</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Biostatistics, Boston University School of Public Health, Boston University, Boston, Massachusetts, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint><date when="2015">2015</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>The basic reproductive number (R₀) and the distribution of the serial interval (SI) are often used to quantify transmission during an infectious disease outbreak. In this paper, we present estimates of R₀ and SI from the 2003 SARS outbreak in Hong Kong and Singapore, and the 2009 pandemic influenza A(H1N1) outbreak in South Africa using methods that expand upon an existing Bayesian framework. This expanded framework allows for the incorporation of additional information, such as contact tracing or household data, through prior distributions. The results for the R₀ and the SI from the influenza outbreak in South Africa were similar regardless of the prior information (<inline-formula id="pone.0118762.e001"><alternatives><graphic xlink:href="pone.0118762.e001.jpg" id="pone.0118762.e001g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M1"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 1.36–1.46, <inline-formula id="pone.0118762.e002"><alternatives><graphic xlink:href="pone.0118762.e002.jpg" id="pone.0118762.e002g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M2"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= 2.0–2.7, <inline-formula id="pone.0118762.e003"><alternatives><graphic xlink:href="pone.0118762.e003.jpg" id="pone.0118762.e003g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M3"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= mean of the SI). The estimates of R₀ and μ for the SARS outbreak ranged from 2.0–4.4 and 7.4–11.3, respectively, and were shown to vary depending on the use of contact tracing data. The impact of the contact tracing data was likely due to the small number of SARS cases relative to the size of the contact tracing sample.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Anderson, Rm" uniqKey="Anderson R">RM Anderson</name>
</author>
<author><name sortKey="May, M" uniqKey="May M">M May</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Svensson, A" uniqKey="Svensson A">A Svensson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name>
</author>
<author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name>
</author>
<author><name sortKey="Anderson, R" uniqKey="Anderson R">R Anderson</name>
</author>
<author><name sortKey="Ferguson, N" uniqKey="Ferguson N">N Ferguson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name>
</author>
<author><name sortKey="Fang, Vj" uniqKey="Fang V">VJ Fang</name>
</author>
<author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name>
</author>
<author><name sortKey="Malik, Peiris Js" uniqKey="Malik P">Peiris JS Malik</name>
</author>
<author><name sortKey="Leung, Gm" uniqKey="Leung G">GM Leung</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Donnelly, Ca" uniqKey="Donnelly C">CA Donnelly</name>
</author>
<author><name sortKey="Finelli, L" uniqKey="Finelli L">L Finelli</name>
</author>
<author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name>
</author>
<author><name sortKey="Olson, Sj" uniqKey="Olson S">SJ Olson</name>
</author>
<author><name sortKey="Doshi, S" uniqKey="Doshi S">S Doshi</name>
</author>
<author><name sortKey="Jackson, Ml" uniqKey="Jackson M">ML Jackson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="White, Lf" uniqKey="White L">LF White</name>
</author>
<author><name sortKey="Pagano, M" uniqKey="Pagano M">M Pagano</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Black, Aj" uniqKey="Black A">AJ Black</name>
</author>
<author><name sortKey="Ross, Jv" uniqKey="Ross J">JV Ross</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name>
</author>
<author><name sortKey="Carrat, F" uniqKey="Carrat F">F Carrat</name>
</author>
<author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name>
</author>
<author><name sortKey="Valleron, Aj" uniqKey="Valleron A">AJ Valleron</name>
</author>
<author><name sortKey="Boelle, Py" uniqKey="Boelle P">PY Boëlle</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chowell, G" uniqKey="Chowell G">G Chowell</name>
</author>
<author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L Simonsen</name>
</author>
<author><name sortKey="Towers, S" uniqKey="Towers S">S Towers</name>
</author>
<author><name sortKey="Miller, Ma" uniqKey="Miller M">MA Miller</name>
</author>
<author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcdonald, Sa" uniqKey="Mcdonald S">SA McDonald</name>
</author>
<author><name sortKey="Presanis, Am" uniqKey="Presanis A">AM Presanis</name>
</author>
<author><name sortKey="Angelis, D" uniqKey="Angelis D">D Angelis</name>
</author>
<author><name sortKey="Hoek, W" uniqKey="Hoek W">W Hoek</name>
</author>
<author><name sortKey="Hooiveld, M" uniqKey="Hooiveld M">M Hooiveld</name>
</author>
<author><name sortKey="Donker, G" uniqKey="Donker G">G Donker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shubin, M" uniqKey="Shubin M">M Shubin</name>
</author>
<author><name sortKey="Virtanen, M" uniqKey="Virtanen M">M Virtanen</name>
</author>
<author><name sortKey="Toikkanen, S" uniqKey="Toikkanen S">S Toikkanen</name>
</author>
<author><name sortKey="Lyytik Inen, O" uniqKey="Lyytik Inen O">O Lyytikäinen</name>
</author>
<author><name sortKey="Auranen, K" uniqKey="Auranen K">K Auranen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yang, F" uniqKey="Yang F">F Yang</name>
</author>
<author><name sortKey="Yuan, L" uniqKey="Yuan L">L Yuan</name>
</author>
<author><name sortKey="Tan, X" uniqKey="Tan X">X Tan</name>
</author>
<author><name sortKey="Huang, C" uniqKey="Huang C">C Huang</name>
</author>
<author><name sortKey="Feng, J" uniqKey="Feng J">J Feng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Becker, Ng" uniqKey="Becker N">NG Becker</name>
</author>
<author><name sortKey="Wang, D" uniqKey="Wang D">D Wang</name>
</author>
<author><name sortKey="Clements, M" uniqKey="Clements M">M Clements</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lunn, D" uniqKey="Lunn D">D Lunn</name>
</author>
<author><name sortKey="Spiegelhalter, D" uniqKey="Spiegelhalter D">D Spiegelhalter</name>
</author>
<author><name sortKey="Thomas, A" uniqKey="Thomas A">A Thomas</name>
</author>
<author><name sortKey="Best, N" uniqKey="Best N">N Best</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Thomas" uniqKey="Thomas">Thomas</name>
</author>
<author><name sortKey="O Hara, B" uniqKey="O Hara B">B O'Hara</name>
</author>
<author><name sortKey="Ligges, U" uniqKey="Ligges U">U Ligges</name>
</author>
<author><name sortKey="Sturtz, S" uniqKey="Sturtz S">S Sturtz</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Bernardo, Jm" uniqKey="Bernardo J">JM Bernardo</name>
</author>
<author><name sortKey="Ramon, Jm" uniqKey="Ramon J">JM Ramon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gelman, A" uniqKey="Gelman A">A Gelman</name>
</author>
<author><name sortKey="Carlin, Jb" uniqKey="Carlin J">JB Carlin</name>
</author>
<author><name sortKey="Stern, Hs" uniqKey="Stern H">HS Stern</name>
</author>
<author><name sortKey="Rubin, Db" uniqKey="Rubin D">DB Rubin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gelman, A" uniqKey="Gelman A">A. Gelman</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Kelly, D" uniqKey="Kelly D">D Kelly</name>
</author>
<author><name sortKey="Atwood, C" uniqKey="Atwood C">C Atwood</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhu, M" uniqKey="Zhu M">M Zhu</name>
</author>
<author><name sortKey="Lu, A" uniqKey="Lu A">A Lu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Spiegelhalter, Dj" uniqKey="Spiegelhalter D">DJ Spiegelhalter</name>
</author>
<author><name sortKey="Best, Ng" uniqKey="Best N">NG Best</name>
</author>
<author><name sortKey="Carlin, Bp" uniqKey="Carlin B">BP Carlin</name>
</author>
<author><name sortKey="Van Der Linde, A" uniqKey="Van Der Linde A">A Van Der Linde</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="White, Lf" uniqKey="White L">LF White</name>
</author>
<author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
<author><name sortKey="Finelli, L" uniqKey="Finelli L">L Finelli</name>
</author>
<author><name sortKey="Reed, C" uniqKey="Reed C">C Reed</name>
</author>
<author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name>
</author>
<author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Archer, Bn" uniqKey="Archer B">BN Archer</name>
</author>
<author><name sortKey="Cohen, C" uniqKey="Cohen C">C Cohen</name>
</author>
<author><name sortKey="Naidoo, D" uniqKey="Naidoo D">D Naidoo</name>
</author>
<author><name sortKey="Thomas, J" uniqKey="Thomas J">J Thomas</name>
</author>
<author><name sortKey="Makunga, C" uniqKey="Makunga C">C Makunga</name>
</author>
<author><name sortKey="Blumberg, L" uniqKey="Blumberg L">L Blumberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Archer, Bn" uniqKey="Archer B">BN Archer</name>
</author>
<author><name sortKey="Timothy, Ga" uniqKey="Timothy G">GA Timothy</name>
</author>
<author><name sortKey="Cohen, C" uniqKey="Cohen C">C Cohen</name>
</author>
<author><name sortKey="Tempia, S" uniqKey="Tempia S">S Tempia</name>
</author>
<author><name sortKey="Huma, M" uniqKey="Huma M">M Huma</name>
</author>
<author><name sortKey="Blumberg, L" uniqKey="Blumberg L">L Blumberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Archer, Bn" uniqKey="Archer B">BN Archer</name>
</author>
<author><name sortKey="Tempia, S" uniqKey="Tempia S">S Tempia</name>
</author>
<author><name sortKey="White, Lf" uniqKey="White L">LF White</name>
</author>
<author><name sortKey="Pagano, M" uniqKey="Pagano M">M Pagano</name>
</author>
<author><name sortKey="Cohen, C" uniqKey="Cohen C">C Cohen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcbryde, Es" uniqKey="Mcbryde E">ES McBryde</name>
</author>
<author><name sortKey="Bergeri, I" uniqKey="Bergeri I">I Bergeri</name>
</author>
<author><name sortKey="Van Gemert, C" uniqKey="Van Gemert C">C van Gemert</name>
</author>
<author><name sortKey="Rotty, J" uniqKey="Rotty J">J Rotty</name>
</author>
<author><name sortKey="Headley, Ej" uniqKey="Headley E">EJ Headley</name>
</author>
<author><name sortKey="Simpson, K" uniqKey="Simpson K">K Simpson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
<author><name sortKey="Teunis, P" uniqKey="Teunis P">P Teunis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name>
</author>
<author><name sortKey="Cohen, T" uniqKey="Cohen T">T Cohen</name>
</author>
<author><name sortKey="Cooper, B" uniqKey="Cooper B">B Cooper</name>
</author>
<author><name sortKey="Robins, J" uniqKey="Robins J">J Robins</name>
</author>
<author><name sortKey="Ma, S" uniqKey="Ma S">S Ma</name>
</author>
<author><name sortKey="James, L" uniqKey="James L">L James</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Opatowski, L" uniqKey="Opatowski L">L Opatowski</name>
</author>
<author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name>
</author>
<author><name sortKey="Griffin, J" uniqKey="Griffin J">J Griffin</name>
</author>
<author><name sortKey="De Silva, E" uniqKey="De Silva E">E de Silva</name>
</author>
<author><name sortKey="Van Kerkhove, Md" uniqKey="Van Kerkhove M">MD Van Kerkhove</name>
</author>
<author><name sortKey="Lyons, Ej" uniqKey="Lyons E">EJ Lyons</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kenah, E" uniqKey="Kenah E">E Kenah</name>
</author>
<author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name>
</author>
<author><name sortKey="Robins, Jm" uniqKey="Robins J">JM Robins</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Griffin, J" uniqKey="Griffin J">J Griffin</name>
</author>
<author><name sortKey="Garske, T" uniqKey="Garske T">T Garske</name>
</author>
<author><name sortKey="Ghani, A" uniqKey="Ghani A">A Ghani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hens, N" uniqKey="Hens N">N Hens</name>
</author>
<author><name sortKey="Calatayud, L" uniqKey="Calatayud L">L Calatayud</name>
</author>
<author><name sortKey="Kurkela, S" uniqKey="Kurkela S">S Kurkela</name>
</author>
<author><name sortKey="Tamme, T" uniqKey="Tamme T">T Tamme</name>
</author>
<author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Te Beest, De" uniqKey="Te Beest D">DE te Beest</name>
</author>
<author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
<author><name sortKey="Donker, T" uniqKey="Donker T">T Donker</name>
</author>
<author><name sortKey="Van Boven, M" uniqKey="Van Boven M">M van Boven</name>
</author>
</analytic>
</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">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group><journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, CA USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">25793993</article-id>
<article-id pub-id-type="pmc">4368801</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0118762</article-id>
<article-id pub-id-type="publisher-id">PONE-D-14-13446</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>The Impact of Prior Information on Estimates of Disease Transmissibility Using Bayesian Tools</article-title>
<alt-title alt-title-type="running-head">Disease Transmissibility Estimates with Bayesian Tools</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Moser</surname>
<given-names>Carlee B.</given-names>
</name>
<xref ref-type="aff" rid="aff001"><sup>1</sup>
</xref>
<xref rid="cor001" ref-type="corresp">*</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Gupta</surname>
<given-names>Mayetri</given-names>
</name>
<xref ref-type="aff" rid="aff001"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Archer</surname>
<given-names>Brett N.</given-names>
</name>
<xref ref-type="aff" rid="aff002"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>White</surname>
<given-names>Laura F.</given-names>
</name>
<xref ref-type="aff" rid="aff001"><sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff001"><label>1</label>
<addr-line>Department of Biostatistics, Boston University School of Public Health, Boston University, Boston, Massachusetts, United States of America</addr-line>
</aff>
<aff id="aff002"><label>2</label>
<addr-line>National Institute for Communicable Diseases, National Health Laboratory Service, Johannesburg, South Africa</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>van Boven</surname>
<given-names>Michiel</given-names>
</name>
<role>Academic Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>National Institute for Public Health and the Environment, NETHERLANDS</addr-line>
</aff>
<author-notes><fn fn-type="COI-statement" id="coi001"><p><bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con" id="contrib001"><p>Conceived and designed the experiments: CBM MG BNA LFW. Analyzed the data: CBM. Contributed reagents/materials/analysis tools: BNA LFW. Wrote the paper: CBM LFW.</p>
</fn>
<corresp id="cor001">* E-mail: <email>cmoser@sdac.harvard.edu</email>
</corresp>
</author-notes>
<pub-date pub-type="epub"><day>20</day>
<month>3</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="collection"><year>2015</year>
</pub-date>
<volume>10</volume>
<issue>3</issue>
<elocation-id>e0118762</elocation-id>
<history><date date-type="received"><day>27</day>
<month>3</month>
<year>2014</year>
</date>
<date date-type="accepted"><day>22</day>
<month>1</month>
<year>2015</year>
</date>
</history>
<permissions><copyright-statement>© 2015 Moser et al</copyright-statement>
<copyright-year>2015</copyright-year>
<copyright-holder>Moser et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:type="simple" xlink:href="pone.0118762.pdf"></self-uri>
<abstract><p>The basic reproductive number (R₀) and the distribution of the serial interval (SI) are often used to quantify transmission during an infectious disease outbreak. In this paper, we present estimates of R₀ and SI from the 2003 SARS outbreak in Hong Kong and Singapore, and the 2009 pandemic influenza A(H1N1) outbreak in South Africa using methods that expand upon an existing Bayesian framework. This expanded framework allows for the incorporation of additional information, such as contact tracing or household data, through prior distributions. The results for the R₀ and the SI from the influenza outbreak in South Africa were similar regardless of the prior information (<inline-formula id="pone.0118762.e001"><alternatives><graphic xlink:href="pone.0118762.e001.jpg" id="pone.0118762.e001g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M1"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 1.36–1.46, <inline-formula id="pone.0118762.e002"><alternatives><graphic xlink:href="pone.0118762.e002.jpg" id="pone.0118762.e002g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M2"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= 2.0–2.7, <inline-formula id="pone.0118762.e003"><alternatives><graphic xlink:href="pone.0118762.e003.jpg" id="pone.0118762.e003g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M3"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= mean of the SI). The estimates of R₀ and μ for the SARS outbreak ranged from 2.0–4.4 and 7.4–11.3, respectively, and were shown to vary depending on the use of contact tracing data. The impact of the contact tracing data was likely due to the small number of SARS cases relative to the size of the contact tracing sample.</p>
</abstract>
<funding-group><funding-statement>LFW was supported by the National Institute of General Medical Sciences (Grant U54GM088558). CBM, MG, and BNA received no specific funding for this work. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts><fig-count count="3"></fig-count>
<table-count count="5"></table-count>
<page-count count="16"></page-count>
</counts>
<custom-meta-group><custom-meta id="data-availability"><meta-name>Data Availability</meta-name>
<meta-value>Data used in these analyses are available upon request because they were obtained from a third party. Brett Archer (email: <email>brettarcher@hotmail.com</email>
) supplied the South Africa H1N1 data, Jacco Wallinga (email: <email>jacco.wallinga@rivm.nl</email>
), supplied the SARS data, and Laura F. White (email: <email>lfwhite@bu.edu</email>
) supplied the estimates of the SI distribution used for the simulations.</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
<notes><title>Data Availability</title>
<p>Data used in these analyses are available upon request because they were obtained from a third party. Brett Archer (email: <email>brettarcher@hotmail.com</email>
) supplied the South Africa H1N1 data, Jacco Wallinga (email: <email>jacco.wallinga@rivm.nl</email>
), supplied the SARS data, and Laura F. White (email: <email>lfwhite@bu.edu</email>
) supplied the estimates of the SI distribution used for the simulations.</p>
</notes>
</front>
<body><sec sec-type="intro" id="sec001"><title>Introduction</title>
<p>When an infectious disease outbreak occurs, public health officials need to understand the dynamics of disease transmission in order to launch an effective response. Two quantities that are often used to describe transmission are the basic reproductive number and the distribution of the serial interval (SI). The basic reproductive number (R<sub>0</sub>
) is the average number of secondary cases a primary case will infect, assuming a completely susceptible population [<xref rid="pone.0118762.ref001" ref-type="bibr">1</xref>
]. The reproductive number is always nonnegative; values less than one are indicative of an outbreak that will not continue to grow in the absence of imported cases. When R<sub>0</sub>
is larger than one, the magnitude of the value guides the types of control measures that are necessary to restrict transmission and control the outbreak. It is also essential to understand the timing between primary and secondary cases. For a given R<sub>0</sub>
, if secondary cases occur shortly after the primary cases, a rapidly growing outbreak will result, which can be more difficult to control than an outbreak with a longer time interval between cases. The timing of the secondary cases is most easily measured by the SI distribution, an observable quantity. The SI is defined as the time between symptom onset in successive cases in a chain of backward transmission. The SI is used as a surrogate measure for the generation interval, which is unobservable, and is defined as the time between consecutive infections in the chain of transmission [<xref rid="pone.0118762.ref002" ref-type="bibr">2</xref>
]. R<sub>0</sub>
and the SI distribution provide important information that is used to initiate an appropriate public health response to an infectious disease outbreak.</p>
<p>Many methods exist to quantify the R<sub>0</sub>
and the SI [<xref rid="pone.0118762.ref003" ref-type="bibr">3</xref>
]. Typically, the SI distribution is estimated using contact tracing or household data (see for example [<xref rid="pone.0118762.ref004" ref-type="bibr">4</xref>
–<xref rid="pone.0118762.ref005" ref-type="bibr">5</xref>
]); however, these studies are often small and subject to potential bias and errors in recall by participants. White and Pagano [<xref rid="pone.0118762.ref006" ref-type="bibr">6</xref>
] introduced a novel approach to simultaneously estimate the R<sub>0</sub>
and the SI using only data from the epidemic curve. In recent years Bayesian methods have been developed to estimate transmission parameters, and can be particularly useful in instances with sparse data or when prior data about an outbreak exists; however these types of models have often been limited to Bayesian evidence synthesis or compartmental models [<xref rid="pone.0118762.ref007" ref-type="bibr">7</xref>
–<xref rid="pone.0118762.ref012" ref-type="bibr">12</xref>
]. Becker et al. [<xref rid="pone.0118762.ref013" ref-type="bibr">13</xref>
] introduced a Bayesian framework to estimate R<sub>0</sub>
and the SI distribution by augmenting the likelihood function, introduced by White and Pagano, with independent observations of the SI from contact tracing data, and obtained posterior estimates through MCMC methods. They also made recommendations about the number of observations from the epidemic curve and contact tracing sample needed to obtain reliable estimates for R<sub>0</sub>
and the SI distribution.</p>
<p>In this paper we describe an extension of the Bayesian methods introduced by Becker et al. Our approach, like Becker et al., also allows for the inclusion of additional data, but does so through a different mechanism, as prior information via prior distributions. In what follows, we present the statistical model introduced by White and Pagano, and outline how to include additional information, such as contact tracing data, via the prior distributions. Details of a simulation study that examines our method are also discussed. Finally, we analyze data from the 2003 SARS outbreak in Hong Kong and Singapore, and the 2009 pandemic influenza A(H1N1) outbreak in South Africa with our method.</p>
</sec>
<sec sec-type="materials|methods" id="sec002"><title>Methods</title>
<sec id="sec003"><title>Statistical model</title>
<p>The method proposed in White and Pagano [<xref rid="pone.0118762.ref006" ref-type="bibr">6</xref>
] can simultaneously estimate the R<sub>0</sub>
and the SI by maximizing the likelihood shown in <xref rid="pone.0118762.e004" ref-type="disp-formula">equation 1</xref>
.</p>
<disp-formula id="pone.0118762.e004"><alternatives><graphic xlink:href="pone.0118762.e004.jpg" id="pone.0118762.e004g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M4"><mml:mi>L</mml:mi>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow><mml:mstyle displaystyle="false"><mml:msubsup><mml:mo>∏</mml:mo>
<mml:mrow><mml:mi>t</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mi>T</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mstyle>
<mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msup>
<mml:msubsup><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mspace width="0.25em"></mml:mspace>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mspace width="0.5em"></mml:mspace>
<mml:msub><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mstyle><mml:msubsup><mml:mi mathvariant="normal">∑</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mtext>min</mml:mtext>
<mml:mspace width="0.05em"></mml:mspace>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msubsup>
</mml:mstyle>
<mml:mrow><mml:msub><mml:mrow><mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:msub><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
<label>(1)</label>
</disp-formula>
<p>The number of new secondary cases at a given time <italic>t</italic>
is defined as N<sub>t</sub>
. For simplicity, we assume <italic>t</italic>
indexes days. Here R<sub>0</sub>
is the basic reproductive number and p<sub>i</sub>
describes the probability of a serial interval that is <italic>i</italic>
days long. The serial intervals are constrained to be no longer than <italic>k</italic>
day and, to follow a multinomial distribution, which is assumed to be stationary. We perform estimation using a Markov chain Monte Carlo (MCMC) method using OpenBUGS software via the BRugs package in R version 2.11.1 [<xref rid="pone.0118762.ref014" ref-type="bibr">14</xref>
–<xref rid="pone.0118762.ref016" ref-type="bibr">16</xref>
]. See <xref rid="pone.0118762.s004" ref-type="supplementary-material">S4 Appendix</xref>
for details.</p>
</sec>
<sec id="sec004"><title>Prior Distributions</title>
<p>To model an outbreak using a Bayesian method, prior distributions are necessary for each parameter of interest: in this instance R<sub>0</sub>
and <bold>p</bold>
= {p<sub>i</sub>
,…, p<sub>k</sub>
}. We choose the prior distribution for R<sub>0</sub>
to be log-normal, such that log(R<sub>0</sub>
) is distributed Gaussian with mean 0 and variance 1,000.</p>
<p>Because we model a discrete form of the serial interval distribution, like White and Pagano and Becker et al., a natural choice for a prior distribution for <bold>p</bold>
is the Dirichlet distribution, which is the conjugate prior for the multinomial distribution. The Dirichlet distribution,
<disp-formula id="pone.0118762.e005"><alternatives><graphic xlink:href="pone.0118762.e005.jpg" id="pone.0118762.e005g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M5"><mml:mi>f</mml:mi>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mi>p</mml:mi>
<mml:mo>|</mml:mo>
<mml:mi>α</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:mi mathvariant="normal">Γ</mml:mi>
<mml:mo>(</mml:mo>
<mml:mrow><mml:msub><mml:mo>∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:msub><mml:mrow><mml:mi>α</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mo>∏</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mi mathvariant="normal">Γ</mml:mi>
<mml:mo>(</mml:mo>
<mml:msub><mml:mrow><mml:mi>α</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:msubsup><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mi>α</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:msubsup><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mi>α</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>…</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mi>α</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>,</mml:mo>
</mml:math>
</alternatives>
<label>(2)</label>
</disp-formula>
is parameterized by a set <bold>α</bold>
, known as hyperparameters, with each α<sub>i</sub>
corresponding to a specific p<sub>i</sub>
. In the absence of prior information, we want to select hyperparameters for the Dirichlet prior such that the distribution will be noninformative [<xref rid="pone.0118762.ref017" ref-type="bibr">17</xref>
–<xref rid="pone.0118762.ref022" ref-type="bibr">22</xref>
]. We examined different priors such as Jeffrey’s prior, which assigns 0.5 for all values of α<sub>i</sub>
, and the Bayes-Laplace prior, which assigns 1 to all α<sub>i</sub>
values. We found that using hyperparameters that are less than 1 can result in a lack of convergence or poor mixing. Because of this, when no additional information about the SI distribution is available, we select the Bayes-Laplace prior. This choice is not noninformative, as it imposes a uniform distribution on the prior, but it is reasonable for outbreaks where there are no contact tracing samples. This prior will be referred to as the uniform prior in the following sections.</p>
<p>When additional information is available, such as contact tracing from initial cases, we can inform the prior distribution with these observations. The hyperparameters of the Dirichlet distribution can be updated with the number of observed serial intervals corresponding to each day of the SI. Previously, Becker et al. [<xref rid="pone.0118762.ref013" ref-type="bibr">13</xref>
] incorporated this information by augmenting it with the likelihood (in equation <xref rid="pone.0118762.e004" ref-type="disp-formula">1</xref>
) and used a Bayes-Laplace prior distribution for the SI (see <xref rid="pone.0118762.e006" ref-type="disp-formula">equation 3</xref>
). Our approach directly informs the SI prior distribution with the observed SIs, which enables us to guide the estimation of the SI by weighting each day of the SI according to the observed serial intervals. The Bayesian model considers the prior information in the estimation, but still puts emphasis on the epidemic curve (daily incidence data).</p>
<p>To further explain the differences between the approaches the following simplified model formulation can be used. The posterior distribution presented in Becker et al. can be generalized to: <italic>P</italic>
(<italic>R</italic>
<sub>0</sub>
,<italic>p</italic>
|<italic>N</italic>
,<italic>r</italic>
) ∝ <italic>P</italic>
(<italic>R</italic>
<sub>0</sub>
) * <italic>P</italic>
(<italic>p</italic>
) * <italic>L</italic>
<sub>1</sub>
(<italic>N</italic>
|<italic>R</italic>
<sub>0</sub>
,<italic>p</italic>
)* <italic>L</italic>
<sub>2</sub>
(<italic>r</italic>
|<italic>p</italic>
), where the posterior is proportional to the prior distributions and likelihood function. P(R<sub>0</sub>
) and P(p) represent the prior distributions for R<sub>0</sub>
and p, and L<sub>1</sub>
and L<sub>2</sub>
are the components of the likelihood, where L<sub>1</sub>
is from White and Pagano and L<sub>2</sub>
is the contact tracing sample that is used to augment the White and Pagano method. In this formulations, P(p) ∼ Dirichlet(<bold>α = 1</bold>
).</p>
<p>The posterior distribution that we propose has a slightly different formulation, in that the contact tracing data is incorporated in the SI prior: <italic>P</italic>
(<italic>R</italic>
<sub>0</sub>
,<italic>p</italic>
|<italic>N</italic>
,<italic>r</italic>
) ∝ <italic>P</italic>
(<italic>R</italic>
<sub>0</sub>
) * <italic>P</italic>
(<italic>p</italic>
|<italic>r</italic>
) * <italic>L</italic>
<sub>1</sub>
(<italic>N</italic>
|<italic>R</italic>
<sub>0</sub>
,<italic>p</italic>
). In our approach, the likelihood function only considers data from the epidemic curve, N, and the contact tracing sample is used to inform the hyperparameters of the SI prior, P(p|r) ∼ Dirichlet(<bold>α = r</bold>
). With this formulation we can consider instances when contact tracing is not available, P(p|r) ∼ Dirichlet(<bold>α = 1</bold>
), or we can weight the contact tracing sample as appropriate, P(p|r) ∼ Dirichlet(<bold>α = r</bold>
<sub><bold>w</bold>
</sub>
).</p>
<p>Our method is flexible because it does not require independent observations of the SI from the same population. For example, when a contact tracing sample is not available, a sample from another outbreak of the same or similar pathogen could be used to inform the prior. If the contact tracing is obtained from a population that is different from the current outbreak population or the contact tracing sample is very large relative to the current outbreak, the contact tracing could be down weighted in the prior, thus putting more emphasis on the current outbreak data as opposed to the outside source. It may not be appropriate to combine data sources from different populations via the likelihood function, as in Becker et al., and instances with large contact tracing samples relative to the outbreak size could bias the results.</p>
<p>If contact tracing samples are not available for an analysis, then the uniform prior (Dirichlet(<bold>α = 1</bold>
)) can be used instead; however, deciding the maximum length of the SI, k, may not be obvious. The deviance information criteria (DIC), a model assessment tool, can be used in these situations to select the best model. A range of plausible values for k should be selected and a model based on each k should be analyzed. The fit of each model can then be compared with the DIC; smaller DIC indicates better model fit. The DIC is useful to compare different models because it uses the posterior densities in its calculation, accounts for the complexity of the model and can easily be implemented in an MCMC simulation [<xref rid="pone.0118762.ref023" ref-type="bibr">23</xref>
].</p>
</sec>
<sec id="sec005"><title>Simulation Study</title>
<p>The following simulation study was designed to assess the validity and utility of our proposed approach and provide guidance in selecting prior distributions. A description of the simulations, including how the data are generated and summarized is included next.</p>
<sec id="sec006"><title>Description of Simulations</title>
<p>We generate outbreaks in the same manner as suggested by White and Pagano [<xref rid="pone.0118762.ref006" ref-type="bibr">6</xref>
], and model our simulated epidemics to be similar to influenza by setting the SI to have a maximum length 5 and mean of 2.74 (see <xref rid="pone.0118762.g001" ref-type="fig">Fig. 1</xref>
), with parameters selected based on estimates from CDC influenza data [<xref rid="pone.0118762.ref024" ref-type="bibr">24</xref>
]. R<sub>0</sub>
assumes the values of 1.25, 3, and 6, and we examine epidemic sizes of 50, 200 and 500 cases. For each scenario discussed, 300 epidemics are generated, with all cases on the first day assumed to be index cases. A description of the epidemic curves for each scenario is shown in <xref rid="pone.0118762.t001" ref-type="table">Table 1</xref>
.</p>
<fig id="pone.0118762.g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0118762.g001</object-id>
<label>Fig 1</label>
<caption><title>Serial Interval Distribution for Simulated Outbreaks.</title>
</caption>
<graphic xlink:href="pone.0118762.g001"></graphic>
</fig>
<table-wrap id="pone.0118762.t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0118762.t001</object-id>
<label>Table 1</label>
<caption><title>Means and Ranges of the Epidemic Lengths for each Simulation Scenario.</title>
</caption>
<alternatives><graphic id="pone.0118762.t001g" xlink:href="pone.0118762.t001"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="1" colspan="1"></th>
<th colspan="2" align="center" rowspan="1">R<sub>0</sub>
= 1.25</th>
<th colspan="2" align="center" rowspan="1">R<sub>0</sub>
= 3</th>
<th colspan="2" align="center" rowspan="1">R<sub>0</sub>
= 6</th>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">N = 500</td>
<td align="left" rowspan="1" colspan="1">50</td>
<td align="left" rowspan="1" colspan="1">25–123</td>
<td align="left" rowspan="1" colspan="1">13</td>
<td align="left" rowspan="1" colspan="1">11–19</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">7–11</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">N = 200</td>
<td align="left" rowspan="1" colspan="1">40</td>
<td align="left" rowspan="1" colspan="1">21–96</td>
<td align="left" rowspan="1" colspan="1">11</td>
<td align="left" rowspan="1" colspan="1">8–21</td>
<td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">6–9</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">N = 50</td>
<td align="left" rowspan="1" colspan="1">24</td>
<td align="left" rowspan="1" colspan="1">9–59</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">6–15</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">4–7</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="t001fn001"><p>There are 300 Epidemics Generated for each Scenario.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>The data from these nine scenarios are analyzed using the White and Pagano model (<xref rid="pone.0118762.e004" ref-type="disp-formula">equation 1</xref>
) with k = 5 and the prior distributions outlined in <xref rid="pone.0118762.t002" ref-type="table">Table 2</xref>
. We first include the uniform prior previously discussed (prior 1), which represents the situation when contact tracing data is not available and we only have daily incidence data (epidemic curve). Two informative prior distributions are also included, where the hyperparameters of the prior (observed SIs) match the SI distribution used to generate the data (<xref rid="pone.0118762.g001" ref-type="fig">Fig. 1</xref>
). The first informative prior (prior 2) has 20 observed SIs, which corresponds to 40%, 10%, and 4% of the epidemic sizes listed above; the second (prior 3) includes only 10 observed SIs, which corresponds to 20%, 5%, 2% of the epidemic sizes. These priors (priors 2 and 3) are considered the “gold standard” in our simulation because the distribution of SI observations matches the SI distribution used to generate the data. We also consider the situation when the contact tracing data used to inform the prior is biased, which is possible in a real life situation. The first misspecified prior (prior 4) is informed using an SI distribution that differs from the true distribution; however, the overall shape and mean of the SI distribution (μ = 2.8) are very similar to the distribution used to generate the data. This scenario could occur when contact tracing data is included from another study of the same or similar pathogen. The second misspecified prior (prior 5) has a SI distribution that differs dramatically from the distribution used to generate the data, with a heavy tail and a prior mean of four (μ = 4). Both priors 4 and 5 are composed of 20 observed SIs, which corresponds to 40%, 10% and 4% of the epidemic sizes.</p>
<table-wrap id="pone.0118762.t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0118762.t002</object-id>
<label>Table 2</label>
<caption><title>Dirichlet Prior Distributions for the Serial Interval for the Simulation Study.</title>
</caption>
<alternatives><graphic id="pone.0118762.t002g" xlink:href="pone.0118762.t002"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="1" colspan="1">Prior</th>
<th align="left" rowspan="1" colspan="1">Type</th>
<th align="left" rowspan="1" colspan="1">Hyperparameters</th>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">Uniform</td>
<td align="left" rowspan="1" colspan="1"><bold>α</bold>
= (1, 1, 1, 1, 1)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1"><italic>Ideally</italic>
Specified</td>
<td align="left" rowspan="1" colspan="1"><bold>α</bold>
= (3, 4, 10, 1, 2)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1"><italic>Ideally</italic>
Specified</td>
<td align="left" rowspan="1" colspan="1"><bold>α</bold>
= (1, 2, 5, 1, 1)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">Misspecified</td>
<td align="left" rowspan="1" colspan="1"><bold>α</bold>
= (1, 10, 4, 2, 3)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">Misspecified</td>
<td align="left" rowspan="1" colspan="1"><bold>α</bold>
= (1, 2, 3, 4, 10)</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>A second set of analyses examines in more detail the scenario when contact tracing data is not available and the investigator must choose the maximum length of serial interval. These analyses model the simulated data assuming various lengths of k, k = {5, 7, 10, 15, 20} and focus on outbreaks based on R<sub>0</sub>
= 1.25 and N = 500 and 200. We limit these analyses to only these outbreaks because of outbreak length limitations, that is, the maximum k considered (k = 20) is larger than some simulated outbreak lengths, which limits the number of observable generations of transmission (see <xref rid="pone.0118762.t001" ref-type="table">Table 1</xref>
). The uniform prior with corresponding k is utilized for each model.</p>
<p>Last, we analyze the simulated data using the Becker et al. method [<xref rid="pone.0118762.ref013" ref-type="bibr">13</xref>
]. The Becker likelihood function, shown in equation <xref rid="pone.0118762.e006" ref-type="disp-formula">3</xref>
, is based on the White and Pagano model (<xref rid="pone.0118762.e004" ref-type="disp-formula">equation 1</xref>
) and is augmented by an additional component composed of observed SIs from contact tracing data. The r<sub>j</sub>
r<sub>j</sub>
represents the number of observed serial intervals for each increment of the serial interval distribution.</p>
<disp-formula id="pone.0118762.e006"><alternatives><graphic xlink:href="pone.0118762.e006.jpg" id="pone.0118762.e006g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M6"><mml:mi>L</mml:mi>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow><mml:mstyle><mml:msubsup><mml:mo>∏</mml:mo>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mi>k</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mstyle>
<mml:mrow><mml:msubsup><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
<mml:mrow><mml:mstyle><mml:msubsup><mml:mo>∏</mml:mo>
<mml:mrow><mml:mi>t</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mi>T</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mstyle>
<mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msup>
<mml:msubsup><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>!</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mrow>
<mml:mspace width="0.25em"></mml:mspace>
<mml:mo>,</mml:mo>
<mml:mspace width="0.5em"></mml:mspace>
<mml:msub><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mstyle><mml:msubsup><mml:mi mathvariant="normal">∑</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mtext>min</mml:mtext>
<mml:mspace width="0.05em"></mml:mspace>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msubsup>
</mml:mstyle>
<mml:mrow><mml:msub><mml:mrow><mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:math>
</alternatives>
<label>(3)</label>
</disp-formula>
<p>For the analyses with the Becker method we also select prior distributions for the parameters, and as with our approach the prior distribution for the reproductive number is log-normal. The serial interval prior for the Becker method differs from our approach because it remains the same across all scenarios, and is set to be a Dirichlet distribution with <bold>α</bold>
= (1, 1, 1, 1, 1), as shown in their paper. The contact tracing data is incorporated in the likelihood function, and thus <xref rid="pone.0118762.t002" ref-type="table">Table 2</xref>
contains the information used to augment the Becker likelihood function (<xref rid="pone.0118762.e006" ref-type="disp-formula">equation 3</xref>
). We assume that these data are independent observations of the SI for the simulation analysis.</p>
<p>We analyze all the simulations via MCMC methods and provide posterior summaries. For each simulated outbreak, single chains were run for 20,000 iterations with a 10,000 iteration burn-in. Trace plots, density plots, and autocorrelation plots were examined to explore sampler convergence.</p>
</sec>
</sec>
<sec id="sec007"><title>Description of Outbreak Data</title>
<p>We apply our method to data from the 2009 pandemic influenza A(H1N1) outbreak in South Africa and to the 2003 SARS outbreak in Hong Kong and Singapore.</p>
<sec id="sec008"><title>South Africa Influenza A(H1N1)2009pdm</title>
<p>The National Institute for Communicable Diseases (NICD) of the National Health Laboratory (NHLS) in South Africa maintained a database of laboratory confirmed cases of H1N1 during the 2009 pandemic. Data were collected throughout the country beginning in April and lasting until October, and include basic demographic information, and spatial and temporal data for each case; date of symptom onset was imputed using a multiple imputation approach for some cases [<xref rid="pone.0118762.ref025" ref-type="bibr">25</xref>
–<xref rid="pone.0118762.ref027" ref-type="bibr">27</xref>
]. We randomly select one imputed outbreak, and because we are interested in estimating R<sub>0</sub>
, only the initial epidemic growth phase of the outbreak is examined. This includes 2,423 cases from days 15 to 35 of the epidemic. Large gaps in cases occurred during the first two weeks of the epidemic due to no sustained transmission from these presumably imported cases. Therefore, only cases occurring after sustained transmission was established are included. Three possible contact tracing samples are considered for informing the SI prior distribution, as shown in <xref rid="pone.0118762.g002" ref-type="fig">Fig. 2</xref>
. The first two samples are from contact tracing data collected in South Africa during the pandemic and consist of confirmed and probable influenza-like illness (ILI) secondary cases [<xref rid="pone.0118762.ref026" ref-type="bibr">26</xref>
]. The third set of data is empirically observed SIs from the influenza A(H1N1)2009pdm outbreak in Victoria, Australia [<xref rid="pone.0118762.ref028" ref-type="bibr">28</xref>
].</p>
<fig id="pone.0118762.g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0118762.g002</object-id>
<label>Fig 2</label>
<caption><title>Contact tracing data from influenza A(H1N1)2009pdm.</title>
<p>Black: Confirmed ILL; Light Gray: Probable ILL; Dark Gray: Australia A(H1N1) Data.</p>
</caption>
<graphic xlink:href="pone.0118762.g002"></graphic>
</fig>
</sec>
<sec id="sec009"><title>Severe Acute Respiratory Syndrome (SARS)</title>
<p>Severe acute respiratory syndrome (SARS) first appeared in Hong Kong and Singapore in 2003. Reported cases of SARS and the date of onset were collected and reported by the World Health Organization (WHO) for Hong Kong and Singapore [<xref rid="pone.0118762.ref029" ref-type="bibr">29</xref>
]. We focus our analysis on the initial growth phase in order to estimate R<sub>0</sub>
. Hong Kong experienced 659 cases in the first 39 days of the outbreak (February 15<sup>th</sup>
—March 25<sup>th</sup>
), while Singapore had 57 cases in the first 20 days of the outbreak (February 25<sup>th</sup>
—March 16<sup>th</sup>
) [<xref rid="pone.0118762.ref029" ref-type="bibr">29</xref>
]. The SI prior distribution is based on 179 (314% of the Singapore incidence data and 27% of the Hong Kong data) observed SIs in Singapore with a mean SI of 8.41 days [<xref rid="pone.0118762.ref030" ref-type="bibr">30</xref>
].</p>
</sec>
</sec>
</sec>
<sec sec-type="results" id="sec010"><title>Results</title>
<sec id="sec011"><title>Simulation Results</title>
<p>Detailed results from the simulations are provided in the supplement and are displayed in <xref rid="pone.0118762.s005" ref-type="supplementary-material">S1</xref>
and <xref rid="pone.0118762.s006" ref-type="supplementary-material">S2</xref>
Fig. and S1-S12 Tables in <xref rid="pone.0118762.s002" ref-type="supplementary-material">S2 Appendix</xref>
, and in <xref rid="pone.0118762.s001" ref-type="supplementary-material">S1 Appendix</xref>
.</p>
<p>In general, our simulations show that our approach provides good estimation of <inline-formula id="pone.0118762.e007"><alternatives><graphic xlink:href="pone.0118762.e007.jpg" id="pone.0118762.e007g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M7"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
and <inline-formula id="pone.0118762.e008"><alternatives><graphic xlink:href="pone.0118762.e008.jpg" id="pone.0118762.e008g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M8"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
when we include contact tracing data that is in agreement with the outbreak data, regardless of the outbreak size. However, the estimates for the reproductive number were affected by the number of observed generations of transmission. When outbreak sizes were small and reproductive numbers were large we observed short epidemic lengths, such as when N = 50 and R<sub>0</sub>
= 6 we only observe epidemic lengths of four to five days for an SI with maximum length five. We see for some scenarios with R<sub>0</sub>
= 6, that <inline-formula id="pone.0118762.e009"><alternatives><graphic xlink:href="pone.0118762.e009.jpg" id="pone.0118762.e009g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M9"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
is slightly underestimated, which is likely due to the limited number of observed generations of transmission and limited total outbreak lengths (see <xref rid="pone.0118762.t002" ref-type="table">Table 2</xref>
).</p>
<p>There were inconsistent results when using a misspecified prior based on biased contact tracing data. When informing the prior with a contact tracing sample that had a total number of observed SIs equaling more than 4% of the total outbreak size, the estimates were biased (i.e. simulations with N = 50 and N = 200). Because we will never know if the contact tracing data is in agreement with the epidemic curve, it is important to not weight the SI prior to be more than 4% of the total outbreak size. In general, if enough generations of transmission are observed, the prior weight has less impact; however, if the total epidemic length is relatively short and few generations of transmission have been observed, then a smaller weight of 1–2% may be more appropriate.</p>
<p>The comparison of the models with varying values of k, using uniform SI priors, shows that the DIC selects the model with k = 5 most often. To summarize the simulations with the DIC, the minimum DIC was calculated across models and then the frequency that each model was selected is depicted in <xref rid="pone.0118762.g003" ref-type="fig">Fig. 3</xref>
. Models with larger ks were sometimes selected; however, the differences between the DICs were small in these instances. The ranges of epidemic curve lengths for the two scenarios that were considered are 25–123 days (N = 500) and 21–96 days (N = 200). Some outbreak lengths are small relative to the largest maximum SI length (k = 20) that was considered, which could affect the estimates. However, we did observe that the DIC is a useful model selection tool.</p>
<fig id="pone.0118762.g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0118762.g003</object-id>
<label>Fig 3</label>
<caption><title>Summary of DIC-based Prior Selection for Simulated Outbreaks.</title>
<p>Gray scale from dark to light: k = 5, k = 7, k = 10, k = 15, k = 20.</p>
</caption>
<graphic xlink:href="pone.0118762.g003"></graphic>
</fig>
<sec id="sec012"><title>Simulation Conclusions</title>
<p>When contact tracing samples are available this data should be used to inform the SI prior distribution. Careful thought should be given to the size of the prior relative to the total outbreak size and length. For smaller outbreaks it is especially important to make sure that the contact tracing sample size is not too large, which could potentially bias the estimates. If enough generations of transmission have been observed to properly estimate a SI of a given maximum length k, then the SI prior should be weighted to be no more 4% of the total outbreak size; however, a range of weights could be considered to examine the robustness of the results. When total outbreak lengths are limited, weighting the SI prior to be 1–2% of the total outbreak size may be more appropriate. If contact tracing data is not available then a range of possible models with reasonable choices of k should be considered. The choice of k can be determined by prior knowledge of similar outbreaks or general knowledge about the pathogen. Each model can be analyzed using a uniform prior distribution for the SI and the models can be compared using the DIC (or other model assessment tool). Smallest DICs indicate best fit, which can be used to narrow down the best model and value for k. These estimates can also be updated and confirmed as the outbreak continues and more cases are observed.</p>
</sec>
</sec>
<sec id="sec013"><title>Outbreak Data</title>
<sec id="sec014"><title>South Africa Influenza A(H1N1)2009pdm</title>
<p>The results of the South African influenza A(H1N1)2009pdm outbreak are shown in <xref rid="pone.0118762.t003" ref-type="table">Table 3</xref>
. Five separate MCMC analyses were conducted, each with a different SI prior distribution, as well as two additional analyses using the White and Pagano method. The first three MCMC analyses were performed using the informative priors shown in <xref rid="pone.0118762.g002" ref-type="fig">Fig. 2</xref>
, and are composed of the actual number of observed serial intervals from the contact tracing samples, which corresponds to 0.8–1.5% of the total outbreak size. The confirmed influenza-like illness cases have a maximum SI length of five days, while the probable cases have a maximum length of six days. Two uniform prior distributions were included, with five and six day maximum SI lengths, chosen to match the lengths of contact tracing samples. To assess if these models (choices of k) are appropriate, additional models with varying values of k, using uniform priors, were also included. Values of k that were considered ranged from k = 5 to 20, with model comparison done via the DIC (<xref rid="pone.0118762.t004" ref-type="table">Table 4</xref>
).</p>
<table-wrap id="pone.0118762.t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0118762.t003</object-id>
<label>Table 3</label>
<caption><title>Means and 95% CIs for South Africa Influenza A(H1N1)2009pdm.</title>
</caption>
<alternatives><graphic id="pone.0118762.t003g" xlink:href="pone.0118762.t003"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="1" colspan="1"></th>
<th colspan="3" align="center" rowspan="1">Informative Prior</th>
<th colspan="2" align="center" rowspan="1">Uniform Prior</th>
<th colspan="2" align="center" rowspan="1">White and Pagano</th>
</tr>
<tr><th align="left" rowspan="1" colspan="1"></th>
<th align="left" rowspan="1" colspan="1">C ILI</th>
<th align="left" rowspan="1" colspan="1">P ILI</th>
<th align="left" rowspan="1" colspan="1">Australia</th>
<th align="left" rowspan="1" colspan="1">K = 5</th>
<th align="left" rowspan="1" colspan="1">K = 6</th>
<th align="left" rowspan="1" colspan="1">K = 5</th>
<th align="left" rowspan="1" colspan="1">K = 6</th>
</tr>
</thead>
<tbody><tr><td rowspan="2" align="left" colspan="1">R<sub>0</sub>
</td>
<td align="left" rowspan="1" colspan="1">1.36</td>
<td align="left" rowspan="1" colspan="1">1.41</td>
<td align="left" rowspan="1" colspan="1">1.43</td>
<td align="left" rowspan="1" colspan="1">1.38</td>
<td align="left" rowspan="1" colspan="1">1.46</td>
<td align="left" rowspan="1" colspan="1">1.36</td>
<td align="left" rowspan="1" colspan="1">1.47</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.27, 1.45</td>
<td align="left" rowspan="1" colspan="1">1.32, 1.51</td>
<td align="left" rowspan="1" colspan="1">1.35, 1.53</td>
<td align="left" rowspan="1" colspan="1">1.28, 1.49</td>
<td align="left" rowspan="1" colspan="1">1.34, 1.60</td>
<td align="left" rowspan="1" colspan="1">1.06, 2.65</td>
<td align="left" rowspan="1" colspan="1">1.08, 2.81</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">μ</td>
<td align="left" rowspan="1" colspan="1">2.00</td>
<td align="left" rowspan="1" colspan="1">2.33</td>
<td align="left" rowspan="1" colspan="1">2.42</td>
<td align="left" rowspan="1" colspan="1">2.19</td>
<td align="left" rowspan="1" colspan="1">2.66</td>
<td align="left" rowspan="1" colspan="1">2.07</td>
<td align="left" rowspan="1" colspan="1">2.70</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.62, 2.42</td>
<td align="left" rowspan="1" colspan="1">1.92, 2.81</td>
<td align="left" rowspan="1" colspan="1">2.06, 2.85</td>
<td align="left" rowspan="1" colspan="1">1.66, 2.67</td>
<td align="left" rowspan="1" colspan="1">2.03, 3.29</td>
<td align="left" rowspan="1" colspan="1">1.08, 3.70</td>
<td align="left" rowspan="1" colspan="1">1.35, 4.37</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">p<sub>1</sub>
</td>
<td align="left" rowspan="1" colspan="1">0.55</td>
<td align="left" rowspan="1" colspan="1">0.45</td>
<td align="left" rowspan="1" colspan="1">0.42</td>
<td align="left" rowspan="1" colspan="1">0.56</td>
<td align="left" rowspan="1" colspan="1">0.49</td>
<td align="left" rowspan="1" colspan="1">0.57</td>
<td align="left" rowspan="1" colspan="1">0.50</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">0.43, 0.67</td>
<td align="left" rowspan="1" colspan="1">0.38, 0.61</td>
<td align="left" rowspan="1" colspan="1">0.31, 0.52</td>
<td align="left" rowspan="1" colspan="1">0.43, 0.70</td>
<td align="left" rowspan="1" colspan="1">0.35, 0.63</td>
<td align="left" rowspan="1" colspan="1">0.08, 0.96</td>
<td align="left" rowspan="1" colspan="1">0.10, 0.85</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">p<sub>2</sub>
</td>
<td align="left" rowspan="1" colspan="1">0.19</td>
<td align="left" rowspan="1" colspan="1">0.15</td>
<td align="left" rowspan="1" colspan="1">0.21</td>
<td align="left" rowspan="1" colspan="1">0.15</td>
<td align="left" rowspan="1" colspan="1">0.15</td>
<td align="left" rowspan="1" colspan="1">0.22</td>
<td align="left" rowspan="1" colspan="1">0.20</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">0.09, 0.32</td>
<td align="left" rowspan="1" colspan="1">0.07, 0.25</td>
<td align="left" rowspan="1" colspan="1">0.12, 0.31</td>
<td align="left" rowspan="1" colspan="1">0.04, 0.29</td>
<td align="left" rowspan="1" colspan="1">0.04, 0.28</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.54</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.49</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">p<sub>3</sub>
</td>
<td align="left" rowspan="1" colspan="1">0.04</td>
<td align="left" rowspan="1" colspan="1">0.08</td>
<td align="left" rowspan="1" colspan="1">0.11</td>
<td align="left" rowspan="1" colspan="1">0.03</td>
<td align="left" rowspan="1" colspan="1">0.03</td>
<td align="left" rowspan="1" colspan="1">0.00</td>
<td align="left" rowspan="1" colspan="1">0.00</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">0.01, 0.11</td>
<td align="left" rowspan="1" colspan="1">0.03, 0.15</td>
<td align="left" rowspan="1" colspan="1">0.06, 0.19</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.10</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.09</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.42</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.38</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">p<sub>4</sub>
</td>
<td align="left" rowspan="1" colspan="1">0.12</td>
<td align="left" rowspan="1" colspan="1">0.14</td>
<td align="left" rowspan="1" colspan="1">0.11</td>
<td align="left" rowspan="1" colspan="1">0.05</td>
<td align="left" rowspan="1" colspan="1">0.04</td>
<td align="left" rowspan="1" colspan="1">0.00</td>
<td align="left" rowspan="1" colspan="1">0.00</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">0.04, 0.23</td>
<td align="left" rowspan="1" colspan="1">0.06, 0.24</td>
<td align="left" rowspan="1" colspan="1">0.05, 0.20</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.18</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.13</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.41</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.32</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">p<sub>5</sub>
</td>
<td align="left" rowspan="1" colspan="1">0.09</td>
<td align="left" rowspan="1" colspan="1">0.07</td>
<td align="left" rowspan="1" colspan="1">0.11</td>
<td align="left" rowspan="1" colspan="1">0.20</td>
<td align="left" rowspan="1" colspan="1">0.10</td>
<td align="left" rowspan="1" colspan="1">0.21</td>
<td align="left" rowspan="1" colspan="1">0.00</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">0.00, 0.22</td>
<td align="left" rowspan="1" colspan="1">0.01, 0.18</td>
<td align="left" rowspan="1" colspan="1">0.04, 0.20</td>
<td align="left" rowspan="1" colspan="1">0.04, 0.35</td>
<td align="left" rowspan="1" colspan="1">0.01, 0.26</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.57</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.39</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">p<sub>6</sub>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0.06</td>
<td align="left" rowspan="1" colspan="1">0.04</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0.18</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0.30</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0.00, 0.18</td>
<td align="left" rowspan="1" colspan="1">0.00, 0.14</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0.02, 0.35</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0.00, 0.61</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="t003fn001"><p>Australia: Australian Contact Trace data; C ILI: Confirmed Influenza-Like Illness; P ILI: Probable Influenza-Like Illness. Informative and Uniform priors are 95% credible intervals, White and Pagano Results are 95% Confidence Intervals.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="pone.0118762.t004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0118762.t004</object-id>
<label>Table 4</label>
<caption><title>DIC, Means and 95% CIs for South Africa Influenza A(H1N1)2009pdm for various uniform priors.</title>
</caption>
<alternatives><graphic id="pone.0118762.t004g" xlink:href="pone.0118762.t004"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="1" colspan="1"></th>
<th align="left" rowspan="1" colspan="1">DIC</th>
<th align="left" rowspan="1" colspan="1">R<sub>0</sub>
</th>
<th align="left" rowspan="1" colspan="1">μ</th>
</tr>
</thead>
<tbody><tr><td rowspan="2" align="left" colspan="1">K = 5</td>
<td rowspan="2" align="left" colspan="1">321.4</td>
<td align="left" rowspan="1" colspan="1">1.38</td>
<td align="left" rowspan="1" colspan="1">2.19</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.28, 1.49</td>
<td align="left" rowspan="1" colspan="1">1.66, 2.67</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1"><bold>K = 6</bold>
</td>
<td rowspan="2" align="left" colspan="1"><bold>315.9</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>1.46</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>2.66</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold>1.34, 1.60</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>2.03, 3.29</bold>
</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1"><bold>K = 7</bold>
</td>
<td rowspan="2" align="left" colspan="1"><bold>316</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>1.51</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>2.92</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold>1.37, 1.66</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>2.23, 3.61</bold>
</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 8</td>
<td rowspan="2" align="left" colspan="1">317.3</td>
<td align="left" rowspan="1" colspan="1">1.53</td>
<td align="left" rowspan="1" colspan="1">3.04</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.38, 1.69</td>
<td align="left" rowspan="1" colspan="1">2.33, 3.75</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 9</td>
<td rowspan="2" align="left" colspan="1">318.7</td>
<td align="left" rowspan="1" colspan="1">1.56</td>
<td align="left" rowspan="1" colspan="1">3.21</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.41, 1.72</td>
<td align="left" rowspan="1" colspan="1">2.50, 3.92</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 10</td>
<td rowspan="2" align="left" colspan="1">319.7</td>
<td align="left" rowspan="1" colspan="1">1.57</td>
<td align="left" rowspan="1" colspan="1">3.33</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.43, 1.75</td>
<td align="left" rowspan="1" colspan="1">2.60, 4.09</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 11</td>
<td rowspan="2" align="left" colspan="1">320.6</td>
<td align="left" rowspan="1" colspan="1">1.60</td>
<td align="left" rowspan="1" colspan="1">3.53</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.44, 1.78</td>
<td align="left" rowspan="1" colspan="1">2.75, 4.34</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 12</td>
<td rowspan="2" align="left" colspan="1">321.5</td>
<td align="left" rowspan="1" colspan="1">1.63</td>
<td align="left" rowspan="1" colspan="1">3.75</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.46, 1.83</td>
<td align="left" rowspan="1" colspan="1">2.89, 4.70</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 13</td>
<td rowspan="2" align="left" colspan="1">321.3</td>
<td align="left" rowspan="1" colspan="1">1.69</td>
<td align="left" rowspan="1" colspan="1">4.19</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.50, 1.92</td>
<td align="left" rowspan="1" colspan="1">3.19, 5.37</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 14</td>
<td rowspan="2" align="left" colspan="1">320.7</td>
<td align="left" rowspan="1" colspan="1">1.77</td>
<td align="left" rowspan="1" colspan="1">4.76</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.55, 2.04</td>
<td align="left" rowspan="1" colspan="1">3.54, 6.17</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 15</td>
<td rowspan="2" align="left" colspan="1">321</td>
<td align="left" rowspan="1" colspan="1">1.82</td>
<td align="left" rowspan="1" colspan="1">5.11</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.59, 2.10</td>
<td align="left" rowspan="1" colspan="1">3.82, 6.54</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 16</td>
<td rowspan="2" align="left" colspan="1">321.6</td>
<td align="left" rowspan="1" colspan="1">1.87</td>
<td align="left" rowspan="1" colspan="1">5.44</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.61, 2.17</td>
<td align="left" rowspan="1" colspan="1">4.06, 6.95</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 17</td>
<td rowspan="2" align="left" colspan="1">321.2</td>
<td align="left" rowspan="1" colspan="1">1.94</td>
<td align="left" rowspan="1" colspan="1">5.98</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.67, 2.27</td>
<td align="left" rowspan="1" colspan="1">4.49, 7.62</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 18</td>
<td rowspan="2" align="left" colspan="1">321.4</td>
<td align="left" rowspan="1" colspan="1">2.01</td>
<td align="left" rowspan="1" colspan="1">6.45</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.71, 2.36</td>
<td align="left" rowspan="1" colspan="1">4.82, 8.20</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 19</td>
<td rowspan="2" align="left" colspan="1">320.9</td>
<td align="left" rowspan="1" colspan="1">2.10</td>
<td align="left" rowspan="1" colspan="1">7.06</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.77, 2.50</td>
<td align="left" rowspan="1" colspan="1">5.27, 9.00</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">K = 20</td>
<td rowspan="2" align="left" colspan="1">320.3</td>
<td align="left" rowspan="1" colspan="1">2.19</td>
<td align="left" rowspan="1" colspan="1">7.71</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">1.82, 2.67</td>
<td align="left" rowspan="1" colspan="1">5.71, 9.92</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>R<sub>0</sub>
was estimated between 1.36–1.46 for when using the priors based on observed SIs from contact tracing, which is consistent with other estimates from South Africa [<xref rid="pone.0118762.ref027" ref-type="bibr">27</xref>
, <xref rid="pone.0118762.ref031" ref-type="bibr">31</xref>
]. The estimates for the average serial interval range from 2.00–2.66 days. Slightly larger posterior estimates for μ and R<sub>0</sub>
are observed for the analyses with a six-day maximum SI length, but overall the results are similar between approaches. The results from the White and Pagano maximum likelihood estimation (MLE) method are consistent with those from the Bayesian models. The MLEs for the five-day SI and six-day SI are <inline-formula id="pone.0118762.e010"><alternatives><graphic xlink:href="pone.0118762.e010.jpg" id="pone.0118762.e010g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M10"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 1.36 and 1.47, and <inline-formula id="pone.0118762.e011"><alternatives><graphic xlink:href="pone.0118762.e011.jpg" id="pone.0118762.e011g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M11"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= 2.07 and 2.70, respectively. The confidence intervals for White and Pagano are wider than the credible intervals obtained from the Bayesian approach, but otherwise they are similar.</p>
<p>The DIC assessment compared 16 possible models with differing values of k, each with a uniform prior. We observe the smallest DIC for models with k = 6 or k = 7, which supports our findings above and our choice of k. Additional exploratory analyses are included in <xref rid="pone.0118762.s003" ref-type="supplementary-material">S3 Appendix</xref>
, which examine the results when using a weighted contact tracing sample.</p>
</sec>
<sec id="sec015"><title>Severe Acute Respiratory Syndrome (SARS)</title>
<p>Singapore and Hong Kong SARS data are analyzed separately for each location and the results are shown in <xref rid="pone.0118762.t005" ref-type="table">Table 5</xref>
. For each location, two analyses are conducted: the first using a contact tracing informed prior distribution for the SI, and the second using a uniform prior distribution, both with a maximum length of 20 days. Additionally, both outbreaks are analyzed with the White and Pagano method. The R<sub>0</sub>
and μ posterior estimates vary depending upon if contact tracing is used as part of the prior. In both locations, there are larger estimates for R<sub>0</sub>
and μ from the analyses using the uniform prior compared to the informed prior (Singapore: <inline-formula id="pone.0118762.e012"><alternatives><graphic xlink:href="pone.0118762.e012.jpg" id="pone.0118762.e012g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M12"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 4.43 v. 3.86, <inline-formula id="pone.0118762.e013"><alternatives><graphic xlink:href="pone.0118762.e013.jpg" id="pone.0118762.e013g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M13"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= 10.68 v. 8.45; Hong Kong: <inline-formula id="pone.0118762.e014"><alternatives><graphic xlink:href="pone.0118762.e014.jpg" id="pone.0118762.e014g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M14"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 2.78 v. 2.02, <inline-formula id="pone.0118762.e015"><alternatives><graphic xlink:href="pone.0118762.e015.jpg" id="pone.0118762.e015g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M15"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= 11.33 v. 7.42). The results for the informative Singapore analysis are similar to those observed in Lipsitch et al. [<xref rid="pone.0118762.ref030" ref-type="bibr">30</xref>
], who observed a mean SI of 8.4 days. The epidemic curve data from Singapore is primarily composed of cases prior to March 12<sup>th</sup>
, the date in which WHO released their first global SARS alert. Lipsitch et al. noted that the mean SI was slightly larger before this date, with μ = 10 days, which is consistent with our results for the uniform prior. Our estimate of R<sub>0</sub>
, 3.86 (95% CI: 2.86–5.01), is similar to the estimate calculated by Wallinga and Teunis before the global health alert, <inline-formula id="pone.0118762.e016"><alternatives><graphic xlink:href="pone.0118762.e016.jpg" id="pone.0118762.e016g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M16"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 3.1 (95% CI: 2.3, 4.0), and to the R<sub>0</sub>
calculated by Lipsitch et al., <inline-formula id="pone.0118762.e017"><alternatives><graphic xlink:href="pone.0118762.e017.jpg" id="pone.0118762.e017g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M17"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 3.5 (95% CI: 1.5, 7.7). The Hong Kong results are slightly different than those from Singapore, but a similar trend is seen when including contact tracing data. Our estimates for R<sub>0</sub>
(2.02 for the informative prior; 2.78 for the uniform prior) are smaller than what Wallinga and Teunis calculated (<inline-formula id="pone.0118762.e018"><alternatives><graphic xlink:href="pone.0118762.e018.jpg" id="pone.0118762.e018g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M18"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 3.6, 95% CI: 3.1, 4.2). The estimates from White and Pagano are substantially larger, especially for Singapore, which is likely to be due to the small sample size and large number of parameters. The 95% confidence intervals estimates are also wide.</p>
<table-wrap id="pone.0118762.t005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0118762.t005</object-id>
<label>Table 5</label>
<caption><title>Means and 95% CIs for SARS outbreak.</title>
</caption>
<alternatives><graphic id="pone.0118762.t005g" xlink:href="pone.0118762.t005"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="1" colspan="1"></th>
<th colspan="2" align="center" rowspan="1">Singapore</th>
<th colspan="2" align="center" rowspan="1">Hong Kong</th>
</tr>
<tr><th align="left" rowspan="1" colspan="1"></th>
<th align="left" rowspan="1" colspan="1">R<sub>0</sub>
</th>
<th align="left" rowspan="1" colspan="1">μ</th>
<th align="left" rowspan="1" colspan="1">R<sub>0</sub>
</th>
<th align="left" rowspan="1" colspan="1">μ</th>
</tr>
</thead>
<tbody><tr><td rowspan="2" align="left" colspan="1">Informative Prior</td>
<td align="left" rowspan="1" colspan="1">3.86</td>
<td align="left" rowspan="1" colspan="1">8.45</td>
<td align="left" rowspan="1" colspan="1">2.02</td>
<td align="left" rowspan="1" colspan="1">7.42</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">2.56, 5.01</td>
<td align="left" rowspan="1" colspan="1">7.91, 9.01</td>
<td align="left" rowspan="1" colspan="1">1.83, 2.22</td>
<td align="left" rowspan="1" colspan="1">6.77, 8.11</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">Uniform</td>
<td align="left" rowspan="1" colspan="1">4.43</td>
<td align="left" rowspan="1" colspan="1">10.68</td>
<td align="left" rowspan="1" colspan="1">2.78</td>
<td align="left" rowspan="1" colspan="1">11.33</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">2.63, 7.19</td>
<td align="left" rowspan="1" colspan="1">8.23, 13.24</td>
<td align="left" rowspan="1" colspan="1">2.19, 3.55</td>
<td align="left" rowspan="1" colspan="1">9.27, 13.39</td>
</tr>
<tr><td rowspan="2" align="left" colspan="1">White and Pagano</td>
<td align="left" rowspan="1" colspan="1">23.95</td>
<td align="left" rowspan="1" colspan="1">16.86</td>
<td align="left" rowspan="1" colspan="1">4.93</td>
<td align="left" rowspan="1" colspan="1">15.73</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">16.67, 49.27</td>
<td align="left" rowspan="1" colspan="1">15.10, 18.39</td>
<td align="left" rowspan="1" colspan="1">3.54, 26.71</td>
<td align="left" rowspan="1" colspan="1">13.78, 18.55</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="t005fn001"><p>Informative and Uniform priors are 95% credible intervals, White and Pagano Results are 95% Confidence Intervals.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>The results from the analysis of the Singapore data should be interpreted with caution given the limited outbreak size and short total epidemic length. The Hong Kong outbreak was larger than Singapore, but was still limited in terms of total length and was not substantially larger than the contact tracing sample size; the contact tracing sample size was 27% of the total outbreak size. We would recommend down weighting the contact tracing sample in the prior in this instance, but due to the large k this may not be feasible. The simplest prior we consider is the uniform prior and with k = 20 this is already equivalent to 3% of the total outbreak size, implying that any informative prior would be weighted to be at least 3% of the total outbreak size. Because of this, the uniform prior may be most appropriate for this analysis.</p>
<p>The SARS outbreak in Hong Kong was analyzed using three models with varying lengths of k (k = 10, 15, and 20) and a uniform prior for the SI; Singapore was not included due to the short total outbreak length and larger values for k were not considered due to a limited epidemic curve length of 39 days in Hong Kong. The models were compared via the DIC and showed that the model with k = 20 was a better fit than k = 10 or k = 15. For k = 10, the DIC = 323.9, with <inline-formula id="pone.0118762.e019"><alternatives><graphic xlink:href="pone.0118762.e019.jpg" id="pone.0118762.e019g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M19"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 1.45 and <inline-formula id="pone.0118762.e020"><alternatives><graphic xlink:href="pone.0118762.e020.jpg" id="pone.0118762.e020g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M20"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= 3.40; for k = 15, the DIC = 289.6, with <inline-formula id="pone.0118762.e021"><alternatives><graphic xlink:href="pone.0118762.e021.jpg" id="pone.0118762.e021g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M21"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 2.02 and <inline-formula id="pone.0118762.e022"><alternatives><graphic xlink:href="pone.0118762.e022.jpg" id="pone.0118762.e022g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M22"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= 8.12; for k = 20, the DIC = 282.2, with <inline-formula id="pone.0118762.e023"><alternatives><graphic xlink:href="pone.0118762.e023.jpg" id="pone.0118762.e023g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M23"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
= 2.78 and <inline-formula id="pone.0118762.e024"><alternatives><graphic xlink:href="pone.0118762.e024.jpg" id="pone.0118762.e024g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M24"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
= 11.33. Modeling the SI with maximum length 20, we see that estimates for R<sub>0</sub>
= 2–3 and the mean SI = 8–11 days.</p>
</sec>
</sec>
</sec>
<sec sec-type="conclusions" id="sec016"><title>Discussion</title>
<p>In this paper, we modified the Bayesian framework introduced by Becker et al. to obtain robust estimates of the reproductive number and the SI distribution using the White and Pagano [<xref rid="pone.0118762.ref006" ref-type="bibr">6</xref>
] likelihood function. This framework allows for the inclusion of additional data sources beyond the epidemic curve through the prior distributions. We explored including a contact tracing sample or household study data with simulations and in the analysis of real outbreak data.</p>
<p>The South African H1N1 outbreak is an example of a large outbreak in which the epidemic curve contains many cases spread over time. We explored including the contact tracing sample with and without weighting the observed SIs in the prior. Using the contact tracing sample as is, which was equivalent to 0.8%-1.5% of the total outbreak size, led to similar results for the posterior estimates. When we increased the weight of the contact tracing sample to be more influential (27% of the total outbreak size), the estimates were impacted. The increase in prior information led to a slight increase in the estimates of the reproductive number and the mean of the serial interval (results shown in Table S13 in <xref rid="pone.0118762.s003" ref-type="supplementary-material">S3 Appendix</xref>
). The posterior mean of μ was almost identical to the prior mean, which indicates that a heavily weighted prior can bias the posterior results. For the SARS outbreak in Singapore and Hong Kong the contact tracing sample size was much larger relative to the number of cases in epidemic curve, which also impacted our results. Here, when we informed the SI prior distribution with contact tracing data, the posterior estimates of R<sub>0</sub>
and μ decreased, and the posterior estimate of μ was similar to the prior mean (mean of the contact tracing sample). This decrease in the estimates is likely related to the timing of the contact tracing sample collection and of the implementation of control measures, and is attributable to including a large contact tracing sample from a different phase of the outbreak. The results from the outbreak analyses support the work of Kenah et al. [<xref rid="pone.0118762.ref032" ref-type="bibr">32</xref>
], in which they show that when the prevalence of an infection is lower, the serial interval should be larger.</p>
<p>Contact tracing samples with a large number of observed SIs relative to the total outbreak size should be used with caution, especially if obtained from a different part of the outbreak. If the additional data source is substantially different from the true underlying process, the posterior estimates may be biased towards the contact tracing sample as seen with the SARS outbreak analysis, the South Africa sensitivity analysis, and prior 5 from our simulation study. We recommend analyzing the data with and without the contact tracing sample to see if the final estimates differ substantially and not using a SI prior that is weighted to be more than 4% of the total outbreak size. In instances when there are substantial differences, we recommend critical examination of the techniques used to obtain the samples and consideration of the total outbreak size and length in order to evaluate potential explanations for the inconsistencies. Ideally, one would be able to determine if the data are in conflict, or if the differences are due to an improvement in estimation because more information is being incorporated.</p>
<p>The SARS outbreaks highlighted the improvements of a Bayesian approach as compared to a frequentist approach. For these outbreaks, we observe more reasonable estimates for the posterior means compared to the MLEs, as well as narrower credible intervals. The White and Pagano method [<xref rid="pone.0118762.ref006" ref-type="bibr">6</xref>
] has been shown to overestimate R<sub>0</sub>
as R<sub>0</sub>
increases [<xref rid="pone.0118762.ref033" ref-type="bibr">33</xref>
], which is likely due to a flattening of the likelihood and the difficulty of finding a maximum using a numerical optimizer. If additional data (e.g. contact tracing) is not available for a given outbreak, in some instances, especially with a large number of parameters, a Bayesian approach with a uniform (Bayes-Laplace) prior for the SI shows improvement over the traditional White and Pagano approach.</p>
<p>When contact tracing data is not available or chosen to not be included in the analysis, then the Bayes-Laplace prior for the SI is recommended, which assigns ones to the hyperparameters of the Dirichlet distribution. One challenge, when there is little or no contact tracing data to inform the prior distribution of the serial interval, is to determine an appropriate value for k, the maximum SI length. We demonstrated through simulation and in the analysis of H1N1 in South Africa and SARS in Hong Kong, that the DIC can be utilized to select the best model and value of k. We recommend that a set of plausible values of k be selected, and then model comparison be done with the DIC; smaller values of DIC indicate better fit.</p>
<p>Our simulations also highlighted how posterior estimates can be affected by the number of observed generations of transmission. When outbreak sizes were small and reproductive numbers were large we observed short epidemic lengths, such as when N = 50 and R<sub>0</sub>
= 6 we only observe epidemic lengths of four to five days for an SI with maximum length five and saw that <inline-formula id="pone.0118762.e025"><alternatives><graphic xlink:href="pone.0118762.e025.jpg" id="pone.0118762.e025g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M25"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
is underestimated. For simulations based on larger values of the reproductive number the results were varied; the best results were observed for large epidemic sizes. Outbreaks based on smaller reproductive numbers resulted in better posterior estimates regardless of outbreak size. Griffin et al. [<xref rid="pone.0118762.ref033" ref-type="bibr">33</xref>
] notes that estimation is better when cases are spread over many generations, as opposed to having more individuals in fewer generations, as is likely to occur with a large reproductive number. We also observed that not only is the total number of observed cases important, but also when these cases appear [<xref rid="pone.0118762.ref013" ref-type="bibr">13</xref>
]. The issues with estimating the parameters with the larger R<sub>0</sub>
could be resolved by allowing the epidemic curve to cover more generations of data (larger N). In fact, one is unlikely to begin such an analysis for outbreak that has only been occurring over so few days. Additional information about the reproductive number that could be used to inform its prior distribution could potentially improve this issue as well. However, in practice, these scenarios where the SI is relatively short and the reproductive number is large are uncommon.</p>
<p>Realistically, the methods we have proposed are much more applicable to a setting where the reproductive number is relatively small. We are able to show when we have sufficient data from the epidemic curve and an accurate contact tracing sample that our approach is preferable to using only epidemic curve data, as demonstrated through the simulations. Our findings support those seen in Becker et al. [<xref rid="pone.0118762.ref013" ref-type="bibr">13</xref>
], in which they conclude that having only 10 additional observations on the SI can substantially improve estimation. We also showed through simulation that our method is comparable to the Becker approach. When no contact tracing data is available the methods are mathematically equivalent, and when we incorporate observations on the SI the two approaches are very similar. The primary difference between the two methods is in how the contact tracing data is incorporated in the estimation process, and the method we propose provides a more general framework in allowing outside data sources and different prior weights to be utilized.</p>
<p>Although our approach has been shown to be advantageous compared to previous approaches, there are some limitations and potential improvements that could be made. First, we assume a stationary distribution for the SI, which may not be appropriate in some settings. As previously discussed, changes in the prevalence of disease can affect transmission, which then affects the SI and R<sub>0</sub>
estimates. One potential solution would be to modify the model to obtain time-updated estimates of the SI. This could be done in phases as prevalence of the disease changes. In addition, one could consider using a hazard-based estimator, as suggested in Kenah et al., or other likelihood based approach [<xref rid="pone.0118762.ref032" ref-type="bibr">32</xref>
, <xref rid="pone.0118762.ref034" ref-type="bibr">34</xref>
]. Because our focus is in the initial exponential growth phase of the outbreak, clustering of cases and SI contraction may not have a substantial impact on estimates.</p>
<p>Another important limitation of our method is the use of a discrete SI distribution. In assuming a multinomial distribution for the SI, the time of infection is considered the date of report, although the actual symptom onset could have occurred any time between the prior day and the current date of report. te Beest et al. [<xref rid="pone.0118762.ref035" ref-type="bibr">35</xref>
], introduced an interval-censored approach that considers the time of symptom onset on the interval between consecutive days. We do not account for interval censoring in our analysis, which could lead to potential biases, particularly when the serial interval is short. An interesting extension of this work would be to use a continuous SI distribution, such as the gamma distribution shown in White and Pagano [<xref rid="pone.0118762.ref006" ref-type="bibr">6</xref>
] or other smoothness assumption and estimate the serial interval using the approach in te Beest et al. [<xref rid="pone.0118762.ref035" ref-type="bibr">35</xref>
]. The prior distributions for the parameters of the gamma distributed serial interval could then be informed from the SI estimate corrected for interval censoring. One could potentially further modify the likelihood function or overall model formulation to allow for interval-censoring to be directly incorporated in our estimation.</p>
<p>The Bayesian methods presented here offer a simple solution to improving estimation of R<sub>0</sub>
and the serial interval. By including contact tracing data via the prior distribution for the SI, we obtain better estimates for these measures when incidence data are sparse and comparable estimates with larger epidemic sizes as compared to the traditional frequentist approach. Our approach also allows for more flexibility when contact tracing samples contain a large number of observed SIs.</p>
</sec>
<sec sec-type="supplementary-material" id="sec017"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0118762.s001"><label>S1 Appendix</label>
<caption><title>Simulations Description.</title>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0118762.s001.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0118762.s002"><label>S2 Appendix</label>
<caption><title>Simulations Results Tables.</title>
<p>(DOC)</p>
</caption>
<media xlink:href="pone.0118762.s002.doc"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0118762.s003"><label>S3 Appendix</label>
<caption><title>South Africa Sensitivity Analysis Description and Results.</title>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0118762.s003.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0118762.s004"><label>S4 Appendix</label>
<caption><title>Bugs and R Code.</title>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0118762.s004.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0118762.s005"><label>S1 Fig</label>
<caption><title>Simulation results for <inline-formula id="pone.0118762.e026"><alternatives><graphic xlink:href="pone.0118762.e026.jpg" id="pone.0118762.e026g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M26"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
.</title>
<p>(DOC)</p>
</caption>
<media xlink:href="pone.0118762.s005.doc"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0118762.s006"><label>S2 Fig</label>
<caption><title>Simulation results for <inline-formula id="pone.0118762.e027"><alternatives><graphic xlink:href="pone.0118762.e027.jpg" id="pone.0118762.e027g" position="anchor" mimetype="image" orientation="portrait"></graphic>
<mml:math id="M27"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
.</title>
<p>(DOC)</p>
</caption>
<media xlink:href="pone.0118762.s006.doc"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><p>The authors thank Dr. Jacco Wallinga of for providing SARS outbreak data.</p>
</ack>
<ref-list><title>References</title>
<ref id="pone.0118762.ref001"><label>1</label>
<mixed-citation publication-type="book"><name><surname>Anderson</surname>
<given-names>RM</given-names>
</name>
, <name><surname>May</surname>
<given-names>M</given-names>
</name>
. <source>Infectious Disease of Humans: Dynamics and Control</source>
. <publisher-loc>Oxford</publisher-loc>
: <publisher-name>Oxford University Press</publisher-name>
; <year>1991</year>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref002"><label>2</label>
<mixed-citation publication-type="journal"><name><surname>Svensson</surname>
<given-names>A</given-names>
</name>
. <article-title>A note on generation times in epidemic models</article-title>
. <source>Math Biosci</source>
. <year>2007</year>
;<volume>208</volume>
(<issue>1</issue>
): <fpage>300</fpage>
–<lpage>311</lpage>
.
<pub-id pub-id-type="pmid">17174352</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref003"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Fraser</surname>
<given-names>C</given-names>
</name>
, <name><surname>Riley</surname>
<given-names>S</given-names>
</name>
, <name><surname>Anderson</surname>
<given-names>R</given-names>
</name>
, <name><surname>Ferguson</surname>
<given-names>N</given-names>
</name>
. <article-title>Factors that make an infectious disease outbreak controllable</article-title>
. <source>Proc Natl Acad Sci U S A</source>
. <year>2004</year>
;<volume>101</volume>
(<issue>16</issue>
): <fpage>6146</fpage>
–<lpage>6151</lpage>
.
<pub-id pub-id-type="pmid">15071187</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref004"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Cowling</surname>
<given-names>BJ</given-names>
</name>
, <name><surname>Fang</surname>
<given-names>VJ</given-names>
</name>
, <name><surname>Riley</surname>
<given-names>S</given-names>
</name>
, <name><surname>Malik</surname>
<given-names>Peiris JS</given-names>
</name>
, <name><surname>Leung</surname>
<given-names>GM</given-names>
</name>
. <article-title>Estimation of the serial interval of influenza</article-title>
. <source>Epidemiology</source>
<year>2009</year>
;<volume>20</volume>
(<issue>3</issue>
): <fpage>344</fpage>
–<lpage>347</lpage>
. <pub-id pub-id-type="doi">10.1097/EDE.0b013e31819d1092</pub-id>
<pub-id pub-id-type="pmid">19279492</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref005"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Donnelly</surname>
<given-names>CA</given-names>
</name>
, <name><surname>Finelli</surname>
<given-names>L</given-names>
</name>
, <name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <name><surname>Olson</surname>
<given-names>SJ</given-names>
</name>
, <name><surname>Doshi</surname>
<given-names>S</given-names>
</name>
, <name><surname>Jackson</surname>
<given-names>ML</given-names>
</name>
, <etal>et al</etal>
<article-title>Serial Intervals and the Temporal Distribution of Secondary Infectious within Households of 2009 Pandemic Influenza A (H1N1): Implications for Influenza Control Recommendations</article-title>
. <source>Clin Infec Dis</source>
. <year>2011</year>
;<volume>52</volume>
(<issue>Suppl 1</issue>
): <fpage>S123</fpage>
–<lpage>S130</lpage>
. <pub-id pub-id-type="doi">10.1093/cid/ciq028</pub-id>
<pub-id pub-id-type="pmid">21342883</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref006"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>White</surname>
<given-names>LF</given-names>
</name>
, <name><surname>Pagano</surname>
<given-names>M</given-names>
</name>
. <article-title>A likelihood-based method for real-time estimation of the serial interval and reproductive number of an epidemic</article-title>
. <source>Stat Med</source>
. <year>2008</year>
;<volume>27</volume>
(<issue>16</issue>
): <fpage>2999</fpage>
–<lpage>3016</lpage>
.
<pub-id pub-id-type="pmid">18058829</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref007"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Black</surname>
<given-names>AJ</given-names>
</name>
, <name><surname>Ross</surname>
<given-names>JV</given-names>
</name>
. <article-title>Estimating a Markovian epidemic model using household serial interval data from the early phase of an epidemic</article-title>
. <source>PLOS ONE</source>
. <year>2013</year>
;<volume>8</volume>
(<issue>8</issue>
): <fpage>e73420</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0073420</pub-id>
<pub-id pub-id-type="pmid">24023679</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref008"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <name><surname>Carrat</surname>
<given-names>F</given-names>
</name>
, <name><surname>Viboud</surname>
<given-names>C</given-names>
</name>
, <name><surname>Valleron</surname>
<given-names>AJ</given-names>
</name>
, <name><surname>Boëlle</surname>
<given-names>PY</given-names>
</name>
. <article-title>A Bayesian MCMC approach to study transmission of influenza: application to household longitudinal data</article-title>
. <source>Stat Med</source>
. <year>2004</year>
;<volume>23</volume>
(<issue>22</issue>
): <fpage>3469</fpage>
–<lpage>3487</lpage>
. <pub-id pub-id-type="doi">10.1002/sim.1912</pub-id>
<pub-id pub-id-type="pmid">15505892</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref009"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Chowell</surname>
<given-names>G</given-names>
</name>
, <name><surname>Simonsen</surname>
<given-names>L</given-names>
</name>
, <name><surname>Towers</surname>
<given-names>S</given-names>
</name>
, <name><surname>Miller</surname>
<given-names>MA</given-names>
</name>
, <name><surname>Viboud</surname>
<given-names>C</given-names>
</name>
. <article-title>Transmission potential of influenza A/H7N9, February to May 2013, China</article-title>
. <source>BMC medicine</source>
. <year>2013</year>
;<volume>11</volume>
(<issue>1</issue>
): <fpage>214</fpage>
.<pub-id pub-id-type="pmid">24083506</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref010"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>McDonald</surname>
<given-names>SA</given-names>
</name>
, <name><surname>Presanis</surname>
<given-names>AM</given-names>
</name>
, <name><surname>Angelis</surname>
<given-names>D</given-names>
</name>
, <name><surname>Hoek</surname>
<given-names>W</given-names>
</name>
, <name><surname>Hooiveld</surname>
<given-names>M</given-names>
</name>
, <name><surname>Donker</surname>
<given-names>G</given-names>
</name>
, <etal>et al</etal>
<article-title>An evidence synthesis approach to estimating the incidence of seasonal influenza in the Netherlands</article-title>
. <source>Influenza and Other Respi Viruses</source>
. <year>2013</year>
;<volume>8</volume>
(<issue>1</issue>
): <fpage>33</fpage>
–<lpage>41</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref011"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Shubin</surname>
<given-names>M</given-names>
</name>
, <name><surname>Virtanen</surname>
<given-names>M</given-names>
</name>
, <name><surname>Toikkanen</surname>
<given-names>S</given-names>
</name>
, <name><surname>Lyytikäinen</surname>
<given-names>O</given-names>
</name>
, <name><surname>Auranen</surname>
<given-names>K</given-names>
</name>
. <article-title>Estimating the burden of A (H1N1) pdm09 influenza in Finland during two seasons</article-title>
. <source>Epidemiol Infect</source>
. <year>2014</year>
;<volume>142</volume>
(<issue>5</issue>
): <fpage>964</fpage>
–<lpage>974</lpage>
. <pub-id pub-id-type="doi">10.1017/S0950268813002537</pub-id>
<pub-id pub-id-type="pmid">24139316</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref012"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Yang</surname>
<given-names>F</given-names>
</name>
, <name><surname>Yuan</surname>
<given-names>L</given-names>
</name>
, <name><surname>Tan</surname>
<given-names>X</given-names>
</name>
, <name><surname>Huang</surname>
<given-names>C</given-names>
</name>
, <name><surname>Feng</surname>
<given-names>J</given-names>
</name>
. <article-title>Bayesian estimation of the effective reproduction number for pandemic influenza A H1N1 in Guangdong Province, China</article-title>
. <source>Ann Epidemiol</source>
. <year>2013</year>
;<volume>23</volume>
(<issue>6</issue>
): <fpage>301</fpage>
–<lpage>306</lpage>
. <pub-id pub-id-type="doi">10.1016/j.annepidem.2013.04.005</pub-id>
<pub-id pub-id-type="pmid">23683708</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref013"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Becker</surname>
<given-names>NG</given-names>
</name>
, <name><surname>Wang</surname>
<given-names>D</given-names>
</name>
, <name><surname>Clements</surname>
<given-names>M</given-names>
</name>
. <article-title>Type and quantity of data needed for an early estimate of transmissibility when an infectious disease emerges</article-title>
. <source>Euro Surveill</source>
. <year>2010</year>
;<volume>14</volume>
(<issue>26</issue>
): pii = 19603.</mixed-citation>
</ref>
<ref id="pone.0118762.ref014"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Lunn</surname>
<given-names>D</given-names>
</name>
, <name><surname>Spiegelhalter</surname>
<given-names>D</given-names>
</name>
, <name><surname>Thomas</surname>
<given-names>A</given-names>
</name>
, <name><surname>Best</surname>
<given-names>N</given-names>
</name>
. <article-title>The BUGS project: Evolution, critique, and future directions</article-title>
. <source>Statist Med</source>
. <year>2009</year>
;<volume>28</volume>
: <fpage>3049</fpage>
–<lpage>3067</lpage>
. <pub-id pub-id-type="doi">10.1002/sim.3680</pub-id>
<pub-id pub-id-type="pmid">19630097</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref015"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Thomas</surname>
</name>
, <name><surname>O'Hara</surname>
<given-names>B</given-names>
</name>
, <name><surname>Ligges</surname>
<given-names>U</given-names>
</name>
, <name><surname>Sturtz</surname>
<given-names>S</given-names>
</name>
. <article-title>Making BUGS Open</article-title>
. <source>R News</source>
. <year>2006</year>
;<volume>6</volume>
(<issue>1</issue>
): <fpage>12</fpage>
–<lpage>17</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref016"><label>16</label>
<mixed-citation publication-type="book"><collab>R Development Core Team</collab>
. <chapter-title>R: A language and environment for statistical computing</chapter-title>
<publisher-name>R Foundation for Statistical Computing</publisher-name>
, <publisher-loc>Vienna, Austria</publisher-loc>
<year>2010</year>
;ISBN 3–900051–07–0, URL <ext-link ext-link-type="uri" xlink:href="http://www.R-project.org">http://www.R-project.org</ext-link>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref017"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Bernardo</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Ramon</surname>
<given-names>JM</given-names>
</name>
. <article-title>An Introduction to Bayesian Reference Analysis: Inference on the Ratio of Multinomial Parameters</article-title>
. <source>J R Stat Soc Series D</source>
. <year>1998</year>
;<volume>47</volume>
(<issue>1</issue>
): <fpage>101</fpage>
–<lpage>135</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref018"><label>18</label>
<mixed-citation publication-type="book"><name><surname>Gelman</surname>
<given-names>A</given-names>
</name>
, <name><surname>Carlin</surname>
<given-names>JB</given-names>
</name>
, <name><surname>Stern</surname>
<given-names>HS</given-names>
</name>
, <name><surname>Rubin</surname>
<given-names>DB</given-names>
</name>
. <source>Bayesian Data Analysis</source>
, <edition>Second Edition</edition>
<publisher-loc>Boca Raton, FL</publisher-loc>
: <publisher-name>Chapman and Hall/CRC</publisher-name>
; <year>2003</year>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref019"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Gelman</surname>
<given-names>A.</given-names>
</name>
<article-title>Prior distributions for variance parameters in hierarchical models</article-title>
. <source>Bayesian Anal</source>
. <year>2006</year>
;<volume>1</volume>
(<issue>3</issue>
): <fpage>515</fpage>
–<lpage>533</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref020"><label>20</label>
<mixed-citation publication-type="other">de Campos CP, Benavoli A. Inference with Multinomial Data: Why to Weaken the Prior Strength. IJCAI 2011, Proceedings-International Joint Conference on Artificial Intelligenc. 2011;22(3): 2107–2112.</mixed-citation>
</ref>
<ref id="pone.0118762.ref021"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>Kelly</surname>
<given-names>D</given-names>
</name>
, <name><surname>Atwood</surname>
<given-names>C</given-names>
</name>
. <article-title>Finding a minimally informative Dirichlet prior distribution using least squares</article-title>
. <source>Reliab Eng Sys Safe</source>
. <year>2011</year>
;<volume>96</volume>
: <fpage>398</fpage>
–<lpage>402</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref022"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Zhu</surname>
<given-names>M</given-names>
</name>
, <name><surname>Lu</surname>
<given-names>A</given-names>
</name>
. <article-title>The counter-intuitive non-informative prior for the Bernoulli family</article-title>
. <source>J. Stat. Educ.</source>
<year>2004</year>
;<volume>12</volume>
: <fpage>1</fpage>
–<lpage>10</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref023"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>Spiegelhalter</surname>
<given-names>DJ</given-names>
</name>
, <name><surname>Best</surname>
<given-names>NG</given-names>
</name>
, <name><surname>Carlin</surname>
<given-names>BP</given-names>
</name>
, <name><surname>Van Der Linde</surname>
<given-names>A</given-names>
</name>
. <article-title>Bayesian measures of model complexity and fit</article-title>
. <source>J R Stat Soc Series B</source>
. <year>2002</year>
;<volume>64</volume>
(<issue>4</issue>
): <fpage>583</fpage>
–<lpage>639</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0118762.ref024"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>White</surname>
<given-names>LF</given-names>
</name>
, <name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
, <name><surname>Finelli</surname>
<given-names>L</given-names>
</name>
, <name><surname>Reed</surname>
<given-names>C</given-names>
</name>
, <name><surname>Riley</surname>
<given-names>S</given-names>
</name>
, <name><surname>Lipsitch</surname>
<given-names>M</given-names>
</name>
, <etal>et al</etal>
<article-title>Estimation of the reproductive number and the serial interval in early phase of the 2009 influenza A/H1N1 pandemic in the USA</article-title>
. <source>Influenza and Other Respi Viruses</source>
. <year>2009</year>
;<volume>3</volume>
(<issue>6</issue>
): <fpage>267</fpage>
–<lpage>276</lpage>
<pub-id pub-id-type="doi">10.1111/j.1750-2659.2009.00106.x</pub-id>
<pub-id pub-id-type="pmid">19903209</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref025"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Archer</surname>
<given-names>BN</given-names>
</name>
, <name><surname>Cohen</surname>
<given-names>C</given-names>
</name>
, <name><surname>Naidoo</surname>
<given-names>D</given-names>
</name>
, <name><surname>Thomas</surname>
<given-names>J</given-names>
</name>
, <name><surname>Makunga</surname>
<given-names>C</given-names>
</name>
, <name><surname>Blumberg</surname>
<given-names>L</given-names>
</name>
, <etal>et al</etal>
<article-title>Interim report on pandemic H1N1 influenza virus infections in South Africa, April to October 2009: Epidemiology and factors associated with fatal cases</article-title>
. <source>Euro Surveill</source>
. <year>2009</year>
;<volume>14</volume>
(<issue>42</issue>
): pii = 19369.
<pub-id pub-id-type="pmid">19883551</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref026"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Archer</surname>
<given-names>BN</given-names>
</name>
, <name><surname>Timothy</surname>
<given-names>GA</given-names>
</name>
, <name><surname>Cohen</surname>
<given-names>C</given-names>
</name>
, <name><surname>Tempia</surname>
<given-names>S</given-names>
</name>
, <name><surname>Huma</surname>
<given-names>M</given-names>
</name>
, <name><surname>Blumberg</surname>
<given-names>L</given-names>
</name>
, <etal>et al</etal>
<article-title>Introduction of 2009 pandemic influenza A(H1N1) into South Africa: Clinical presentation, epidemiology and transmissibility of the first 100 cases</article-title>
. <source>J Infect Dis</source>
. <year>2012</year>
;<volume>206</volume>
(<issue>sup 1</issue>
): <fpage>S148</fpage>
–<lpage>S153</lpage>
. <pub-id pub-id-type="doi">10.1093/infdis/jis583</pub-id>
<pub-id pub-id-type="pmid">23169962</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref027"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Archer</surname>
<given-names>BN</given-names>
</name>
, <name><surname>Tempia</surname>
<given-names>S</given-names>
</name>
, <name><surname>White</surname>
<given-names>LF</given-names>
</name>
, <name><surname>Pagano</surname>
<given-names>M</given-names>
</name>
, <name><surname>Cohen</surname>
<given-names>C</given-names>
</name>
. <article-title>Reproductive Number and Serial Interval of the First Wave of Influenza A(H1N1)pdm09 Virus in South Africa</article-title>
. <source>PLOS ONE</source>
. <year>2012</year>
;<volume>7</volume>
(<issue>11</issue>
): <fpage>e49482</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0049482</pub-id>
<pub-id pub-id-type="pmid">23166682</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref028"><label>28</label>
<mixed-citation publication-type="journal"><name><surname>McBryde</surname>
<given-names>ES</given-names>
</name>
, <name><surname>Bergeri</surname>
<given-names>I</given-names>
</name>
, <name><surname>van Gemert</surname>
<given-names>C</given-names>
</name>
, <name><surname>Rotty</surname>
<given-names>J</given-names>
</name>
, <name><surname>Headley</surname>
<given-names>EJ</given-names>
</name>
, <name><surname>Simpson</surname>
<given-names>K</given-names>
</name>
, <etal>et al</etal>
<article-title>Early transmission characteristics of influenza A(H1N1)v in Australia: Victorian state, 16 May—3 June 2009</article-title>
. <source>Euro Surveill</source>
. <year>2009</year>
;<volume>14</volume>
(<issue>42</issue>
): pii = 19363.
<pub-id pub-id-type="pmid">19883551</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref029"><label>29</label>
<mixed-citation publication-type="journal"><name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
, <name><surname>Teunis</surname>
<given-names>P</given-names>
</name>
. <article-title>Different epidemic curves for Severe Acute Respiratory Syndrome</article-title>
. <source>Am J Epidemiol</source>
. <year>2004</year>
;<volume>160</volume>
(<issue>6</issue>
): <fpage>509</fpage>
–<lpage>516</lpage>
.
<pub-id pub-id-type="pmid">15353409</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref030"><label>30</label>
<mixed-citation publication-type="journal"><name><surname>Lipsitch</surname>
<given-names>M</given-names>
</name>
, <name><surname>Cohen</surname>
<given-names>T</given-names>
</name>
, <name><surname>Cooper</surname>
<given-names>B</given-names>
</name>
, <name><surname>Robins</surname>
<given-names>J</given-names>
</name>
, <name><surname>Ma</surname>
<given-names>S</given-names>
</name>
, <name><surname>James</surname>
<given-names>L</given-names>
</name>
, <etal>et al</etal>
<article-title>Transmission dynamics and control of Sever Acute Respiratory Syndrome</article-title>
. <source>Science</source>
. <year>2003</year>
;<volume>300</volume>
: <fpage>1966</fpage>
–<lpage>1970</lpage>
.
<pub-id pub-id-type="pmid">12766207</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref031"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>Opatowski</surname>
<given-names>L</given-names>
</name>
, <name><surname>Fraser</surname>
<given-names>C</given-names>
</name>
, <name><surname>Griffin</surname>
<given-names>J</given-names>
</name>
, <name><surname>de Silva</surname>
<given-names>E</given-names>
</name>
, <name><surname>Van Kerkhove</surname>
<given-names>MD</given-names>
</name>
, <name><surname>Lyons</surname>
<given-names>EJ</given-names>
</name>
, <etal>et al</etal>
<article-title>Transmission Characteristics of the 2009 H1N1 Influenza Pandemic: Comparison of 8 Southern Hemisphere Countries</article-title>
. <source>PLoS Pathog</source>
. <year>2011</year>
;<volume>7</volume>
(<issue>9</issue>
): <fpage>e1002225</fpage>
<pub-id pub-id-type="doi">10.1371/journal.ppat.1002225</pub-id>
<pub-id pub-id-type="pmid">21909272</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref032"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Kenah</surname>
<given-names>E</given-names>
</name>
, <name><surname>Lipsitch</surname>
<given-names>M</given-names>
</name>
, <name><surname>Robins</surname>
<given-names>JM</given-names>
</name>
. <article-title>Generation interval contraction and epidemic data analysis</article-title>
. <source>Math Biosci</source>
. <year>2008</year>
;<volume>213</volume>
(<issue>1</issue>
): <fpage>71</fpage>
–<lpage>79</lpage>
. <pub-id pub-id-type="doi">10.1016/j.mbs.2008.02.007</pub-id>
<pub-id pub-id-type="pmid">18394654</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref033"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>Griffin</surname>
<given-names>J</given-names>
</name>
, <name><surname>Garske</surname>
<given-names>T</given-names>
</name>
, <name><surname>Ghani</surname>
<given-names>A</given-names>
</name>
. <article-title>Joint estimation of the basic reproduction number and generation time parameters for infectious disease outbreaks</article-title>
. <source>Biostatistics</source>
. <year>2011</year>
;<volume>12</volume>
(<issue>2</issue>
): <fpage>303</fpage>
–<lpage>312</lpage>
. <pub-id pub-id-type="doi">10.1093/biostatistics/kxq058</pub-id>
<pub-id pub-id-type="pmid">20858771</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref034"><label>34</label>
<mixed-citation publication-type="journal"><name><surname>Hens</surname>
<given-names>N</given-names>
</name>
, <name><surname>Calatayud</surname>
<given-names>L</given-names>
</name>
, <name><surname>Kurkela</surname>
<given-names>S</given-names>
</name>
, <name><surname>Tamme</surname>
<given-names>T</given-names>
</name>
, <name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
. <article-title>Robust reconstruction and analysis of outbreak data: influenza A (H1N1)v transmission in a school-based population</article-title>
. <source>Am J Epidemiol</source>
. <year>2012</year>
;<volume>176</volume>
(<issue>3</issue>
): <fpage>196</fpage>
–<lpage>203</lpage>
. <pub-id pub-id-type="doi">10.1093/aje/kws006</pub-id>
<pub-id pub-id-type="pmid">22791742</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0118762.ref035"><label>35</label>
<mixed-citation publication-type="journal"><name><surname>te Beest</surname>
<given-names>DE</given-names>
</name>
, <name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
, <name><surname>Donker</surname>
<given-names>T</given-names>
</name>
, <name><surname>van Boven</surname>
<given-names>M</given-names>
</name>
. <article-title>Estimating the generation interval of influenza A (H1N1) in a range of social settings. Epidemiology</article-title>
. <year>2013</year>
;<volume>24</volume>
(<issue>2</issue>
): <fpage>244</fpage>
–<lpage>250</lpage>
.</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 001343 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 001343 | 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:4368801 |texte= The Impact of Prior Information on Estimates of Disease Transmissibility Using Bayesian Tools }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:25793993" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a SrasV1
![]() | This area was generated with Dilib version V0.6.33. | ![]() |