Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Early real-time estimation of the basic reproduction number of emerging or reemerging infectious diseases in a community with heterogeneous contact pattern: Using data from Hong Kong 2009 H1N1 Pandemic Influenza as an illustrative example</title><author><name sortKey="Kwok, Kin On" sort="Kwok, Kin On" uniqKey="Kwok K" first="Kin On" last="Kwok">Kin On Kwok</name><affiliation><nlm:aff id="aff001"><addr-line>School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, Hong Kong, People’s Republic of China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Davoudi, Bahman" sort="Davoudi, Bahman" uniqKey="Davoudi B" first="Bahman" last="Davoudi">Bahman Davoudi</name><affiliation><nlm:aff id="aff002"><addr-line>Division of Mathematical Modeling, British Columbia Centre for Disease Control, Vancouver, British Columbia, Canada</addr-line></nlm:aff></affiliation></author><author><name sortKey="Riley, Steven" sort="Riley, Steven" uniqKey="Riley S" first="Steven" last="Riley">Steven Riley</name><affiliation><nlm:aff id="aff001"><addr-line>School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, Hong Kong, People’s Republic of China</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff003"><addr-line>MRC Centre for Outbreak Analysis and Modelling, Department for Infectious Disease Epidemiology, Imperial College London, United Kingdom</addr-line></nlm:aff></affiliation></author><author><name sortKey="Pourbohloul, Babak" sort="Pourbohloul, Babak" uniqKey="Pourbohloul B" first="Babak" last="Pourbohloul">Babak Pourbohloul</name><affiliation><nlm:aff id="aff002"><addr-line>Division of Mathematical Modeling, British Columbia Centre for Disease Control, Vancouver, British Columbia, Canada</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff004"><addr-line>School of Population & Public Health, Faculty of Medicine, University of British Columbia, Vancouver, British Columbia, Canada</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">26372219</idno><idno type="pmc">4570805</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4570805</idno><idno type="RBID">PMC:4570805</idno><idno type="doi">10.1371/journal.pone.0137959</idno><date when="2015">2015</date><idno type="wicri:Area/Pmc/Corpus">001346</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">001346</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Early real-time estimation of the basic reproduction number of emerging or reemerging infectious diseases in a community with heterogeneous contact pattern: Using data from Hong Kong 2009 H1N1 Pandemic Influenza as an illustrative example</title><author><name sortKey="Kwok, Kin On" sort="Kwok, Kin On" uniqKey="Kwok K" first="Kin On" last="Kwok">Kin On Kwok</name><affiliation><nlm:aff id="aff001"><addr-line>School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, Hong Kong, People’s Republic of China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Davoudi, Bahman" sort="Davoudi, Bahman" uniqKey="Davoudi B" first="Bahman" last="Davoudi">Bahman Davoudi</name><affiliation><nlm:aff id="aff002"><addr-line>Division of Mathematical Modeling, British Columbia Centre for Disease Control, Vancouver, British Columbia, Canada</addr-line></nlm:aff></affiliation></author><author><name sortKey="Riley, Steven" sort="Riley, Steven" uniqKey="Riley S" first="Steven" last="Riley">Steven Riley</name><affiliation><nlm:aff id="aff001"><addr-line>School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, Hong Kong, People’s Republic of China</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff003"><addr-line>MRC Centre for Outbreak Analysis and Modelling, Department for Infectious Disease Epidemiology, Imperial College London, United Kingdom</addr-line></nlm:aff></affiliation></author><author><name sortKey="Pourbohloul, Babak" sort="Pourbohloul, Babak" uniqKey="Pourbohloul B" first="Babak" last="Pourbohloul">Babak Pourbohloul</name><affiliation><nlm:aff id="aff002"><addr-line>Division of Mathematical Modeling, British Columbia Centre for Disease Control, Vancouver, British Columbia, Canada</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff004"><addr-line>School of Population & Public Health, Faculty of Medicine, University of British Columbia, Vancouver, British Columbia, Canada</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>Emerging and re-emerging infections such as SARS (2003) and pandemic H1N1 (2009) have caused concern for public health researchers and policy makers due to the increased burden of these diseases on health care systems. This concern has prompted the use of mathematical models to evaluate strategies to control disease spread, making these models invaluable tools to identify optimal intervention strategies. A particularly important quantity in infectious disease epidemiology is the basic reproduction number, R<sub>0.</sub> Estimation of this quantity is crucial for effective control responses in the early phase of an epidemic. In our previous study, an approach for estimating the basic reproduction number in real time was developed. This approach uses case notification data and the structure of potential transmission contacts to accurately estimate R<sub>0</sub> from the limited amount of information available at the early stage of an outbreak. Based on this approach, we extend the existing methodology; the most recent method features intra- and inter-age groups contact heterogeneity. Given the number of newly reported cases at the early stage of the outbreak, with parsimony assumptions on removal distribution and infectivity profile of the diseases, experiments to estimate real time R<sub>0</sub> under different levels of intra- and inter-group contact heterogeneity using two age groups are presented. We show that the new method converges more quickly to the actual value of R<sub>0</sub> than the previous one, in particular when there is high-level intra-group and inter-group contact heterogeneity. With the age specific contact patterns, number of newly reported cases, removal distribution, and information about the natural history of the 2009 pandemic influenza in Hong Kong, we also use the extended model to estimate R<sub>0</sub> and age-specific R<sub>0</sub>.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name></author><author><name sortKey="Donnelly, Ca" uniqKey="Donnelly C">CA Donnelly</name></author><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name></author><author><name sortKey="Hanage, Wp" uniqKey="Hanage W">WP Hanage</name></author><author><name sortKey="Van Kerkhove, Md" uniqKey="Van Kerkhove M">MD Van Kerkhove</name></author><author><name sortKey="Hollingsworth, Td" uniqKey="Hollingsworth T">TD Hollingsworth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name></author><author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name></author><author><name sortKey="Donnelly, Ca" uniqKey="Donnelly C">CA Donnelly</name></author><author><name sortKey="Ghani, Ac" uniqKey="Ghani A">AC Ghani</name></author><author><name sortKey="Abu Raddad, Lj" uniqKey="Abu Raddad L">LJ Abu-Raddad</name></author><author><name sortKey="Hedley, Aj" uniqKey="Hedley A">AJ Hedley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, L D" uniqKey="Wang L">L-D Wang</name></author><author><name sortKey="Chen, H G" uniqKey="Chen H">H-G Chen</name></author><author><name sortKey="Guo, J G" uniqKey="Guo J">J-G Guo</name></author><author><name sortKey="Zeng, X J" uniqKey="Zeng X">X-J Zeng</name></author><author><name sortKey="Hong, X L" uniqKey="Hong X">X-L Hong</name></author><author><name sortKey="Xiong, J J" uniqKey="Xiong J">J-J Xiong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name></author><author><name sortKey="Van Boven, M" uniqKey="Van Boven M">M Van Boven</name></author><author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pourbohloul, B" uniqKey="Pourbohloul B">B Pourbohloul</name></author><author><name sortKey="Ahued, A" uniqKey="Ahued A">A Ahued</name></author><author><name sortKey="Davoudi, B" uniqKey="Davoudi B">B Davoudi</name></author><author><name sortKey="Meza, R" uniqKey="Meza R">R Meza</name></author><author><name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name></author><author><name sortKey="Skowronski, Dm" uniqKey="Skowronski D">DM Skowronski</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Anderson, Rm" uniqKey="Anderson R">RM Anderson</name></author><author><name sortKey="May, Rm" uniqKey="May R">RM May</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name></author><author><name sortKey="Pourbohloul, B" uniqKey="Pourbohloul B">B Pourbohloul</name></author><author><name sortKey="Newman, Me" uniqKey="Newman M">ME Newman</name></author><author><name sortKey="Skowronski, Dm" uniqKey="Skowronski D">DM Skowronski</name></author><author><name sortKey="Brunham, Rc" uniqKey="Brunham R">RC Brunham</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Donnelly, Ca" uniqKey="Donnelly C">CA Donnelly</name></author><author><name sortKey="Ghani, Ac" uniqKey="Ghani A">AC Ghani</name></author><author><name sortKey="Leung, Gm" uniqKey="Leung G">GM Leung</name></author><author><name sortKey="Hedley, Aj" uniqKey="Hedley A">AJ Hedley</name></author><author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name></author><author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wu, Jt" uniqKey="Wu J">JT Wu</name></author><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name></author><author><name sortKey="Lau, Ehy" uniqKey="Lau E">EHY Lau</name></author><author><name sortKey="M Ip, Dk" uniqKey="M Ip D">DK M Ip</name></author><author><name sortKey="Ho, Lm" uniqKey="Ho L">LM Ho</name></author><author><name sortKey="Tsang, T" uniqKey="Tsang T">T Tsang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Noel, Pa" uniqKey="Noel P">PA Noël</name></author><author><name sortKey="Davoudi, B" uniqKey="Davoudi B">B Davoudi</name></author><author><name sortKey="Brunham, Rc" uniqKey="Brunham R">RC Brunham</name></author><author><name sortKey="Dube, Lj" uniqKey="Dube L">LJ Dube</name></author><author><name sortKey="Pourbohloul, B" uniqKey="Pourbohloul B">B Pourbohloul</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Davoudi, B" uniqKey="Davoudi B">B Davoudi</name></author><author><name sortKey="Miller, Jc" uniqKey="Miller J">JC Miller</name></author><author><name sortKey="Meza, R" uniqKey="Meza R">R Meza</name></author><author><name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name></author><author><name sortKey="Earn, Djd" uniqKey="Earn D">DJD Earn</name></author><author><name sortKey="Pourbohloul, B" uniqKey="Pourbohloul B">B Pourbohloul</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Press, Wh" uniqKey="Press W">WH Press</name></author><author><name sortKey="Teukolsky, Sa" uniqKey="Teukolsky S">SA Teukolsky</name></author><author><name sortKey="Vetterling, Wt" uniqKey="Vetterling W">WT Vetterling</name></author><author><name sortKey="Flannery, Bp" uniqKey="Flannery B">BP Flannery</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tuite, Ar" uniqKey="Tuite A">AR Tuite</name></author><author><name sortKey="Greer, Al" uniqKey="Greer A">AL Greer</name></author><author><name sortKey="Whelan, M" uniqKey="Whelan M">M Whelan</name></author><author><name sortKey="Winter, Al" uniqKey="Winter A">AL Winter</name></author><author><name sortKey="Lee, B" uniqKey="Lee B">B Lee</name></author><author><name sortKey="Yan, P" uniqKey="Yan P">P Yan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name></author><author><name sortKey="Cummings, Dat" uniqKey="Cummings D">DAT Cummings</name></author><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name></author><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="Meeyai, A" uniqKey="Meeyai A">A Meeyai</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name></author><author><name sortKey="Nizam, A" uniqKey="Nizam A">A Nizam</name></author><author><name sortKey="Xu, S" uniqKey="Xu S">S Xu</name></author><author><name sortKey="Ungchusak, K" uniqKey="Ungchusak K">K Ungchusak</name></author><author><name sortKey="Hanshaoworakul, W" uniqKey="Hanshaoworakul W">W Hanshaoworakul</name></author><author><name sortKey="Cummings, Da" uniqKey="Cummings D">DA Cummings</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name></author><author><name sortKey="Kwok, Ko" uniqKey="Kwok K">KO Kwok</name></author><author><name sortKey="Wu, Km" uniqKey="Wu K">KM Wu</name></author><author><name sortKey="Ning, Dy" uniqKey="Ning D">DY Ning</name></author><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name></author><author><name sortKey="Wu, Jt" uniqKey="Wu J">JT Wu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kwok, Ko" uniqKey="Kwok K">KO Kwok</name></author><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name></author><author><name sortKey="Peiris, Js" uniqKey="Peiris J">JS Peiris</name></author><author><name sortKey="Wei, V" uniqKey="Wei V">V Wei</name></author><author><name sortKey="Wu, Km" uniqKey="Wu K">KM Wu</name></author><author><name sortKey="Ning, Dy" uniqKey="Ning D">DY Ning</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name></author><author><name sortKey="Lau, Ms" uniqKey="Lau M">MS Lau</name></author><author><name sortKey="Ho, Lm" uniqKey="Ho L">LM Ho</name></author><author><name sortKey="Chuang, Sk" uniqKey="Chuang S">SK Chuang</name></author><author><name sortKey="Tsang, T" uniqKey="Tsang T">T Tsang</name></author><author><name sortKey="Liu, Sh" uniqKey="Liu S">SH Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Glass, K" uniqKey="Glass K">K Glass</name></author><author><name sortKey="Kelly, H" uniqKey="Kelly H">H Kelly</name></author><author><name sortKey="Mercer, Gn" uniqKey="Mercer G">GN Mercer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kwok, Ko" uniqKey="Kwok K">KO Kwok</name></author><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name></author><author><name sortKey="Wei, V" uniqKey="Wei V">V Wei</name></author><author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name></author><author><name sortKey="Read, J" uniqKey="Read J">J Read</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">26372219</article-id><article-id pub-id-type="pmc">4570805</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0137959</article-id><article-id pub-id-type="publisher-id">PONE-D-13-29570</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Early real-time estimation of the basic reproduction number of emerging or reemerging infectious diseases in a community with heterogeneous contact pattern: Using data from Hong Kong 2009 H1N1 Pandemic Influenza as an illustrative example</article-title><alt-title alt-title-type="running-head">Early real time estimation of R0</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Kwok</surname><given-names>Kin On</given-names></name><xref ref-type="aff" rid="aff001"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Davoudi</surname><given-names>Bahman</given-names></name><xref ref-type="aff" rid="aff002"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Riley</surname><given-names>Steven</given-names></name><xref ref-type="aff" rid="aff001"><sup>1</sup></xref><xref ref-type="aff" rid="aff003"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Pourbohloul</surname><given-names>Babak</given-names></name><xref ref-type="aff" rid="aff002"><sup>2</sup></xref><xref ref-type="aff" rid="aff004"><sup>4</sup></xref><xref rid="cor001" ref-type="corresp">*</xref></contrib></contrib-group><aff id="aff001"><label>1</label><addr-line>School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, Hong Kong, People’s Republic of China</addr-line></aff><aff id="aff002"><label>2</label><addr-line>Division of Mathematical Modeling, British Columbia Centre for Disease Control, Vancouver, British Columbia, Canada</addr-line></aff><aff id="aff003"><label>3</label><addr-line>MRC Centre for Outbreak Analysis and Modelling, Department for Infectious Disease Epidemiology, Imperial College London, United Kingdom</addr-line></aff><aff id="aff004"><label>4</label><addr-line>School of Population & Public Health, Faculty of Medicine, University of British Columbia, Vancouver, British Columbia, Canada</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>McCaw</surname><given-names>James M</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>The University of Melbourne, AUSTRALIA</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: BP BD. Performed the experiments: KOK BD. Analyzed the data: KOK BD. Contributed reagents/materials/analysis tools: KOK SR BD. Wrote the paper: BP KOK BD SR.</p></fn><corresp id="cor001">* E-mail: <email>babak.p@ubc.ca</email></corresp></author-notes><pub-date pub-type="epub"><day>15</day><month>9</month><year>2015</year></pub-date><pub-date pub-type="collection"><year>2015</year></pub-date><volume>10</volume><issue>9</issue><elocation-id>e0137959</elocation-id><history><date date-type="received"><day>12</day><month>7</month><year>2013</year></date><date date-type="accepted"><day>24</day><month>8</month><year>2015</year></date></history><permissions><copyright-statement>© 2015 Kwok et al</copyright-statement><copyright-year>2015</copyright-year><copyright-holder>Kwok 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.0137959.pdf"></self-uri><abstract><p>Emerging and re-emerging infections such as SARS (2003) and pandemic H1N1 (2009) have caused concern for public health researchers and policy makers due to the increased burden of these diseases on health care systems. This concern has prompted the use of mathematical models to evaluate strategies to control disease spread, making these models invaluable tools to identify optimal intervention strategies. A particularly important quantity in infectious disease epidemiology is the basic reproduction number, R<sub>0.</sub> Estimation of this quantity is crucial for effective control responses in the early phase of an epidemic. In our previous study, an approach for estimating the basic reproduction number in real time was developed. This approach uses case notification data and the structure of potential transmission contacts to accurately estimate R<sub>0</sub> from the limited amount of information available at the early stage of an outbreak. Based on this approach, we extend the existing methodology; the most recent method features intra- and inter-age groups contact heterogeneity. Given the number of newly reported cases at the early stage of the outbreak, with parsimony assumptions on removal distribution and infectivity profile of the diseases, experiments to estimate real time R<sub>0</sub> under different levels of intra- and inter-group contact heterogeneity using two age groups are presented. We show that the new method converges more quickly to the actual value of R<sub>0</sub> than the previous one, in particular when there is high-level intra-group and inter-group contact heterogeneity. With the age specific contact patterns, number of newly reported cases, removal distribution, and information about the natural history of the 2009 pandemic influenza in Hong Kong, we also use the extended model to estimate R<sub>0</sub> and age-specific R<sub>0</sub>.</p></abstract><funding-group><funding-statement>BP would like to acknowledge the support of the Canadian Institutes of Health Research (Grant No. PTL-97126), the Michael Smith Foundation for Health Research (Grant No. CI-SSH-00931(06-1)), and the British Columbia Ministry of Health (Pandemic Preparedness Modeling Project). SR thanks the Medical Research Council (UK, Project Grant MR/J008761/1), Wellcome Trust (Project Award 093488/Z/10/Z) and the National Institute for Health Research (UK, Health Protection Research Unit). KOK would like to acknowledge the support by R01 TW008246-01 from Fogarty International Centre (FluScape); the Research Fund for the Control of Infectious Disease, Food and Health Bureau, Government of the Hong Kong SAR (RFCID Number: 110642 and PHE 21) and Research Grants Council, Hong Kong (Grant No. 776810). 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="1"></table-count><page-count count="12"></page-count></counts></article-meta></front><body><sec sec-type="intro" id="sec001"><title>Introduction</title><p>The emergence and global spread of a novel H1N1 influenza A virus (pH1N1) in early 2009 presented a number of new challenges to policy makers and epidemiologists [<xref rid="pone.0137959.ref001" ref-type="bibr">1</xref>]. In the early phases of the epidemic, with limited information on the number of infections, policy makers tried to mitigate disease spread as best they could. Mathematical modeling of infectious disease has made increasing contributions to public health decisions [<xref rid="pone.0137959.ref002" ref-type="bibr">2</xref>] [<xref rid="pone.0137959.ref003" ref-type="bibr">3</xref>] [<xref rid="pone.0137959.ref004" ref-type="bibr">4</xref>] with the aim of estimating the basic reproduction number, R<sub>0</sub> [<xref rid="pone.0137959.ref005" ref-type="bibr">5</xref>]. R<sub>0</sub> is defined as the number of secondary cases generated by one infectious individual during his/her entire infection period in a totally susceptible population [<xref rid="pone.0137959.ref006" ref-type="bibr">6</xref>]. When R<sub>0</sub> is greater than one, there is a likelihood that the infection will to continue to spread; however if R<sub>0</sub> is less than one, then the infection will eventually die out. For this reason, estimating the basic reproduction number is paramount for assessing the effectiveness of different intervention strategies.</p><p>Most classical mathematical models of infectious disease transmission, also known as compartmental models, assume homogeneous mixing among individuals. In other words, each individual has an equal chance of infecting others in a population. However, the estimates of R<sub>0</sub> from models with this parsimonious assumption may differ greatly from models featuring heterogeneous contacts patterns, resulting in different epidemiological outcomes under the same intervention policies [<xref rid="pone.0137959.ref007" ref-type="bibr">7</xref>]. Studies have suggested that many infectious diseases exhibit heterogeneity in transmission efficiency in certain individuals (for example “super-spreaders” in SARS (2003)), or in specific groups of individuals that appear to account for most of the transmission (such as school age children in H1N1 influenza A (2009)) [<xref rid="pone.0137959.ref008" ref-type="bibr">8</xref>] [<xref rid="pone.0137959.ref009" ref-type="bibr">9</xref>]. To ensure that this heterogeneity is modeled, contact network models are used to explicitly capture patterns of disease transmission. Contact network models can also be used to compare the intervention strategies for different infectious diseases when epidemiological data are limited, as was the case with SARS in 2003 and H1N1 pandemic influenza virus in 2009 [<xref rid="pone.0137959.ref005" ref-type="bibr">5</xref>] [<xref rid="pone.0137959.ref007" ref-type="bibr">7</xref>].</p><p>Heterogeneity in contact patterns among individuals in different age groups is an important factor in the transmission of infectious diseases. In this study, a new approach to network modeling is outlined to demonstrate how well the real time estimation of R<sub>0</sub> can be obtained in the early stages of a disease outbreak, using simulated data. Our approach will specifically take into account different levels of inter- and intra-age group heterogeneity. We also validate our analytical method using simulation data. Lastly, we use the tested methodology to estimate the basic reproduction number of pH1N1 in Hong Kong.</p><p>In the following section, we first introduce the homogeneous model (model A) without age-structure and then heterogeneous model (model B) that includes age structure. We then compare the performance of these two models using the simulation data in which inter- and intra-degree can be changed. We show that the second model converges very faster to real R<sub>0</sub> value when inter-degree is big. We finally use both models to estimate the R<sub>0</sub> value for the last pandemic data which occurred in Hong Kong.</p></sec><sec sec-type="materials|methods" id="sec002"><title>Methods</title><sec id="sec003"><title>a) Overview of a Network Approach to Estimate R<sub>0</sub></title><p>Contact network models characterize the contacts between individuals in a social setting that can potentially result in disease transmission. The contact patterns can vary between individuals in the population. We assume that there are <italic>N</italic> individuals in the population. Each individual is represented as a vertex, and an edge represents a pathway of possible infection transmission between two individuals. The number of contacts that an individual has is denoted by <italic>k</italic>., Also, the probability that a randomly chosen vertex has <italic>k</italic> contacts (<italic>k</italic> degrees) is denoted <italic>p</italic><sub><italic>k</italic></sub>. Useful information can be derived from the distribution of number of contacts, known as the degree distribution. For instance, the first moment of degree distribution is the average number of nearest neighbours to a randomly selected individual <italic>z</italic><sub>1</sub>, and the second moment is related to <italic>z</italic><sub>2</sub>, the number of second-degree neighbours (i.e., neighbours of neighbours) of a randomly chosen individual. The excess degree <italic>z</italic> is defined as the number of edges emanating from an infected individual, excluding the edge that was the source of the focal infection; it is simply the ratio <inline-formula id="pone.0137959.e001"><alternatives><graphic xlink:href="pone.0137959.e001.jpg" id="pone.0137959.e001g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M1"><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives></inline-formula> [<xref rid="pone.0137959.ref010" ref-type="bibr">10</xref>].</p><p>A continuous time approach to estimate ongoing basic reproduction number with three-state (susceptible, infected and removed) framework was previously described in [<xref rid="pone.0137959.ref011" ref-type="bibr">11</xref>]. Age structure was not explicitly included in this framework. The excess degree <inline-formula id="pone.0137959.e002"><alternatives><graphic xlink:href="pone.0137959.e002.jpg" id="pone.0137959.e002g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M2"><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives></inline-formula> was assumed to be the same to every individual. We defined <italic>λ</italic> as the transmission rate of an infected individual to a susceptible individual if the individual was not infected previously. Thus, <italic>T</italic>(<italic>t</italic>) = 1 − <italic>e</italic><sup>−<italic>λt</italic></sup> is the probability that an individual transmits the disease to one of its neighbours by time <italic>t</italic> (since acquiring the infection). In the same manner, <italic>λ</italic><sub><italic>r</italic></sub> was defined as the removal rate of an infected individual given that the individual was not removed from the infectious state, previously.</p><p>Thus <inline-formula id="pone.0137959.e003"><alternatives><graphic xlink:href="pone.0137959.e003.jpg" id="pone.0137959.e003g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M3"><mml:mrow><mml:mi>ψ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>−</mml:mo><mml:msub><mml:mi>λ</mml:mi><mml:mi>r</mml:mi></mml:msub><mml:mi>t</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:math></alternatives></inline-formula> and <inline-formula id="pone.0137959.e004"><alternatives><graphic xlink:href="pone.0137959.e004.jpg" id="pone.0137959.e004g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M4"><mml:mrow><mml:mi>ϕ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>λ</mml:mi><mml:mi>r</mml:mi></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>−</mml:mo><mml:msub><mml:mi>λ</mml:mi><mml:mi>r</mml:mi></mml:msub><mml:mi>t</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:math></alternatives></inline-formula> are probability function that the time for an individual to be removed at time <italic>t</italic> and its corresponding probability density function respectively. The basic reproduction number can be written as the product of the expected excess degree and the expected transmissibility <inline-formula id="pone.0137959.e005"><alternatives><graphic xlink:href="pone.0137959.e005.jpg" id="pone.0137959.e005g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M5"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mi>T</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></alternatives></inline-formula>.</p></sec><sec id="sec004"><title>b) Non-Age-Structured Model (Model A)</title><p>In this study, model A is simply the discretized version of the model in [<xref rid="pone.0137959.ref011" ref-type="bibr">11</xref>] when the approach is applied to daily count of incidence cases. A discrete-time approximation described in [<xref rid="pone.0137959.ref005" ref-type="bibr">5</xref>] is used throughout the derivation of the formula of R<sub>0</sub>. For ease of illustration of model A, we assume the discrete time step Δ<italic>t</italic> to be 1 day.</p><p>Given limited information about the natural history of newly emerging diseases, <italic>λ</italic> and <italic>λ</italic><sub><italic>r</italic></sub> were assumed to be constant. As the outbreak progresses, the number of susceptible <italic>S</italic>(<italic>t</italic>), infected <italic>I</italic>(<italic>t</italic>), and recovered <italic>R</italic>(<italic>t</italic>) individuals will vary simultaneously.</p><p>However, we cannot obtain these quantities directly, and the only available information during the epidemic is the number of newly reported cases at time t, denoted as <italic>J</italic>(<italic>t</italic>).</p><p>Using <italic>J</italic>(<italic>t</italic>). and <italic>ψ</italic>(<italic>t</italic>), the following equations can be obtained (discrete-value version of [<xref rid="pone.0137959.ref011" ref-type="bibr">11</xref>])
<disp-formula id="pone.0137959.e006"><alternatives><graphic xlink:href="pone.0137959.e006.jpg" id="pone.0137959.e006g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M6"><mml:mrow><mml:mi>R</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>t</mml:mi></mml:munderover><mml:mrow><mml:mo>{</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>ψ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>}</mml:mo><mml:mi>J</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>−</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow></mml:math></alternatives></disp-formula><disp-formula id="pone.0137959.e007"><alternatives><graphic xlink:href="pone.0137959.e007.jpg" id="pone.0137959.e007g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M7"><mml:mrow><mml:mi>I</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>t</mml:mi></mml:munderover><mml:mrow><mml:mi>ψ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>J</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>−</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow></mml:math></alternatives></disp-formula></p><p>The expected number of infectious and removed individuals who were infected by individuals already removed at time t, <italic>I</italic><sup><italic>r</italic></sup>(<italic>t</italic>), and <italic>R</italic><sup><italic>r</italic></sup>(<italic>t</italic>) respectively, can be estimated by <italic>ϕ</italic>(<italic>t</italic>), <italic>J</italic>(<italic>t</italic>) and <italic>ψ</italic>(<italic>t</italic>) as follows [<xref rid="pone.0137959.ref011" ref-type="bibr">11</xref>]:
<disp-formula id="pone.0137959.e008"><alternatives><graphic xlink:href="pone.0137959.e008.jpg" id="pone.0137959.e008g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M8"><mml:mrow><mml:msup><mml:mi>I</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:msup><mml:mi>R</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>t</mml:mi></mml:munderover><mml:mrow><mml:mi>J</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mfrac><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mi>J</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>ψ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mi>i</mml:mi><mml:mo>−</mml:mo><mml:mi>k</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle><mml:mi>ϕ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mi>J</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>ψ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle><mml:mi>ϕ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:mrow></mml:math></alternatives><label>(1)</label></disp-formula></p><p>At time <italic>t</italic>, the expected transmissibility of infectious individuals before they are removed can be calculated as
<disp-formula id="pone.0137959.e009"><alternatives><graphic xlink:href="pone.0137959.e009.jpg" id="pone.0137959.e009g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M9"><mml:mrow><mml:msup><mml:mi>T</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>t</mml:mi></mml:munderover><mml:mrow><mml:mi>T</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>ϕ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>−</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>t</mml:mi></mml:munderover><mml:mrow><mml:mi>ϕ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>−</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives><label>(2)</label></disp-formula>where <inline-formula id="pone.0137959.e010"><alternatives><graphic xlink:href="pone.0137959.e010.jpg" id="pone.0137959.e010g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M10"><mml:mrow><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>t</mml:mi></mml:munderover><mml:mrow><mml:mi>J</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow></mml:math></alternatives></inline-formula> is the cumulative number of new infections at time <italic>t</italic>.</p><p>There is an important relationship between <italic>T</italic><sup><italic>r</italic></sup>(<italic>t</italic>), <inline-formula id="pone.0137959.e011"><alternatives><graphic xlink:href="pone.0137959.e011.jpg" id="pone.0137959.e011g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M11"><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives></inline-formula>, <italic>I</italic><sup><italic>r</italic></sup>(<italic>t</italic>), <italic>R</italic><sup><italic>r</italic></sup>(<italic>t</italic>), and <italic>R</italic>(<italic>t</italic>) as in [<xref rid="pone.0137959.ref011" ref-type="bibr">11</xref>]:
<disp-formula id="pone.0137959.e012"><alternatives><graphic xlink:href="pone.0137959.e012.jpg" id="pone.0137959.e012g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M12"><mml:mrow><mml:msup><mml:mi>T</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msup><mml:mi>I</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msup><mml:mi>R</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>R</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives><label>(3)</label></disp-formula></p><p>By Eq (<xref rid="pone.0137959.e012" ref-type="disp-formula">3</xref>), <italic>λ</italic> can be estimated and we have <italic>T</italic>(<italic>t</italic>). Given that there is available information on the number of new infections up to time <italic>t</italic>, the basic reproduction number can be obtained at time <italic>t</italic> by the following formula
<disp-formula id="pone.0137959.e013"><alternatives><graphic xlink:href="pone.0137959.e013.jpg" id="pone.0137959.e013g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M13"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>t</mml:mi></mml:munderover><mml:mrow><mml:mi>ϕ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>T</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow></mml:math></alternatives><label>(4)</label></disp-formula></p></sec><sec id="sec005"><title>c) Age-Structured Model (Model B)</title><p>Model B is a modified version of model A, where age specific contact patterns are included. We defined <italic>k</italic> age groups with <italic>n</italic> individuals in the population. The excess degree of <italic>i</italic>th age group <italic>z</italic>(<italic>i</italic>) for <italic>i</italic> = 1, 2, …, <italic>k</italic> is calculated based on the degree distribution of <italic>n</italic>(<italic>i</italic>) individuals in <italic>i</italic>th age group. During the epidemic, we have information on the total number of new infections stratified by <italic>k</italic> age groups at time <italic>t</italic>. <italic>J</italic><sup><italic>B</italic></sup>(<italic>i</italic>, <italic>t</italic>) is defined as the total number of new infections in <italic>i</italic>th age group at time <italic>t</italic>. Thus, the excess degree of an infected population will vary across time and is defined as follows:
<disp-formula id="pone.0137959.e014"><alternatives><graphic xlink:href="pone.0137959.e014.jpg" id="pone.0137959.e014g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M14"><mml:mrow><mml:msub><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mrow><mml:mi>B</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mi>z</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:msup><mml:mi>J</mml:mi><mml:mi>B</mml:mi></mml:msup><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow><mml:mrow><mml:mi>J</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives></disp-formula></p><p>Eq (<xref rid="pone.0137959.e012" ref-type="disp-formula">3</xref>) then becomes
<disp-formula id="pone.0137959.e015"><alternatives><graphic xlink:href="pone.0137959.e015.jpg" id="pone.0137959.e015g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M15"><mml:mrow><mml:msup><mml:mi>T</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msup><mml:mi>I</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:msup><mml:mi>R</mml:mi><mml:mi>r</mml:mi></mml:msup><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>t</mml:mi></mml:munderover><mml:mrow><mml:mo>{</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>ψ</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>}</mml:mo><mml:mi>J</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>−</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle><mml:msub><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mrow><mml:mi>B</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>−</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives><label>(5)</label></disp-formula></p><p>Using Eq (<xref rid="pone.0137959.e015" ref-type="disp-formula">5</xref>) and <inline-formula id="pone.0137959.e016"><alternatives><graphic xlink:href="pone.0137959.e016.jpg" id="pone.0137959.e016g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M16"><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mi>n</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mi>z</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:mrow><mml:mi>n</mml:mi></mml:mfrac></mml:mrow></mml:math></alternatives></inline-formula> we estimate <italic>λ</italic>. The basic reproduction number from model B can then be obtained numerically by using Eq (<xref rid="pone.0137959.e013" ref-type="disp-formula">4</xref>) [<xref rid="pone.0137959.ref012" ref-type="bibr">12</xref>].</p></sec><sec id="sec006"><title>d) Simulation of Incidence Data with Known R<sub>0</sub> and Overall Degree Distribution</title><p>The time series data of the daily incidence of influenza is first simulated based on the binomial random network in [<xref rid="pone.0137959.ref005" ref-type="bibr">5</xref>] [<xref rid="pone.0137959.ref011" ref-type="bibr">11</xref>]. An overall degree distribution, a known R<sub>0</sub> or <italic>λ</italic><sub><italic>r</italic></sub> and the natural history of influenza are inputs of the network. We assume that the contact patterns of individuals would be similar to those seen in a typical day in the early phase of the epidemic.</p><p>To evaluate which model can rapidly provide a better estimate of R<sub>0</sub> in the early phase of an epidemic of influenza, the same overall degree distribution which is the same as simulated data and the number of daily reported cases are inputs of both model A (without age structure) [<xref rid="pone.0137959.ref005" ref-type="bibr">5</xref>] and the proposed model B (with age structure). Ongoing estimate of R<sub>0</sub> from each model will be provided. We also apply model B to the real data as an illustrative example with 2 age groups: school aged children (<20 years old) as age group 1 and adults (20 years old or above) as age group 2 to estimate R<sub>0</sub> for the 2009 outbreak of H1N1pdm in Hong Kong [<xref rid="pone.0137959.s001" ref-type="supplementary-material">S1 Data</xref>]. From the start of May in 2009, patients who were admitted to public hospitals with acute respiratory illness were tested by (RT)-PCR. Information from each individual tested was entered into an information system (E-flu) managed by the Hong Kong Hospital Authority–a statutory body responsible for managing Hong Kong's public hospitals and their services to the community. The E-flu database generates incidence data on children and adults by the date of symptoms onset. This database was gathered as a part of a public health response and therefore not subject to any ethics approval. Permission was granted to analyze these data anonymously for the purpose of this study.</p></sec></sec><sec sec-type="results" id="sec007"><title>Results</title><p>We used a time step in our model of 0.02 days. The two models were tested under different levels of intra- and inter-age group contact heterogeneity in different sizes of epidemics. For ease of illustration, only two age groups in the community are considered (adults and children), although the above equations could easily be generalized to situations where information about more refined age groups is available. With six different known values of R<sub>0</sub> and similar value for <italic>λ</italic><sub><italic>r</italic></sub>, we used computer simulations to generate incidence data for these two age groups under three levels of inter-group and within group heterogeneity in contact patterns. In these simulations we construct contact patterns such that individuals in each age group have uniform degree distributions with a given width.</p><p>Similar to assumptions in [<xref rid="pone.0137959.ref011" ref-type="bibr">11</xref>], we defined the value of <italic>λ</italic><sub><italic>r</italic></sub> as 0.25, an arbitrary value for generating synthetic incidence data. Six known values of R<sub>0</sub> (<xref rid="pone.0137959.g001" ref-type="fig">Fig 1A–1F</xref>) are 2.05, 2.58, 3.12 1.51, 1.52 and 1.53, respectively. In <xref rid="pone.0137959.g001" ref-type="fig">Fig 1A</xref>, we show that model B with age structure yields a better estimate of R<sub>0</sub> and converges at the true value (R<sub>0</sub> = 2.05) more quickly in the early phase of the epidemic (~day 9 to 15) over non-age-structured model A (day 18 to day 21). Between day 9 and day 15, the estimated R<sub>0</sub> values fluctuate ~ 0.4% around the true value. Between day 0 and day 15, less than 100 cases are required to estimate R<sub>0</sub>.</p><fig id="pone.0137959.g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0137959.g001</object-id><label>Fig 1</label><caption><title>Real-time estimation of R<sub>0</sub> using model A (black) and model B (green).</title><p>Both models are run under different hypothetical degree distributions and two known values of R<sub>0</sub> for the two age groups (shown in the sub-diagrams in the initial phase of the epidemic). The true values of R<sub>0</sub> (Fig 1A-1F) are 2.05, 2.58, 3.12 1.51, 1.52 and 1.53 respectively. The red line shows the epidemic curve and the dotted blue line is the true value of R<sub>0</sub>. The corresponding values of <italic>z</italic><sub><italic>2</italic></sub><italic>/z</italic><sub><italic>1</italic></sub> are (A) 19.1, (B) 24.1, (C) 29.1, (D) 14.0, (E)14.15, (F)14.3.</p></caption><graphic xlink:href="pone.0137959.g001"></graphic></fig><p>In densely populated area like Hong Kong, the between group and within-group variation in number of contacts is very high. To characterize this feature, we hypothesized 2 scenarios for 2 given values of R<sub>0</sub> as seen in <xref rid="pone.0137959.g001" ref-type="fig">Fig 1A–1C</xref> and <xref rid="pone.0137959.g001" ref-type="fig">Fig 1D–1F</xref>. Across <xref rid="pone.0137959.g001" ref-type="fig">Fig 1A–1C, we</xref> keep the degree distribution of individuals in the first age group unchanged (red bar chart in the sub-diagram) and shift the degree distribution of individuals in the second age group (blue bar chart in the sub-diagram) to the right, thereby increasing the inter-group heterogeneity of contacts. Model B shows a quick convergence to the true value of R<sub>0</sub>, while model A cannot accurately estimate R<sub>0</sub>, especially when inter-group heterogeneity in the contact network is high. This suggests that model B provides a better real-time estimate over model A whenever there is any level of inter-group contact heterogeneity.</p><p>We also test which model gives a more accurate estimate of R<sub>0</sub> with larger intra-group contact heterogeneity (a wider distribution in each age group). We run the simulations on each age group using three datasets where the contacts are uniformly distributed with a width of 3, 5, and 7 degrees respectively (<xref rid="pone.0137959.g001" ref-type="fig">Fig 1D–1F</xref>). For incidence data with 3 degrees on each group, model B provides more accurate estimates of R<sub>0</sub> between day 30 and day 45 than model A (as in <xref rid="pone.0137959.g001" ref-type="fig">Fig 1D</xref>). Under these conditions, about 50 to 100 cases are required for model B to converge on the true value of R<sub>0</sub>.</p><p>For incidence data with 5 degrees and 7 degrees on each group, model B also gives an accurate estimate of R<sub>0</sub> as represented in <xref rid="pone.0137959.g001" ref-type="fig">Fig 1E</xref> and <xref rid="pone.0137959.g001" ref-type="fig">Fig 1F</xref>. In <xref rid="pone.0137959.g001" ref-type="fig">Fig 1F</xref>, the distance between degree distributions of the two age groups is shortened compared with <xref rid="pone.0137959.g001" ref-type="fig">Fig 1D</xref>. In other words, the inter-group contact heterogeneity in the two age groups decreases across <xref rid="pone.0137959.g001" ref-type="fig">Fig 1D–1F</xref>. This explains why the differences between the two estimates from both model A and model B are not significant in <xref rid="pone.0137959.g001" ref-type="fig">Fig 1F</xref>. It also reveals that model A is a good choice if contacts are more uniformly distributed across different age groups.</p><p>We also applied model B to the hospital incidence data extracted from the E-flu system and the overall degree distribution from the contact survey to estimate R<sub>0</sub> for the 2009 outbreak of H1N1pdm in Hong Kong. In total, there were 28,338 laboratory-confirmed cases during the period between 30<sup>th</sup> April 2009 and 7<sup>th</sup> February 2010. The removal distribution can also be derived using the 24,873 laboratory-confirmed cases with complete information on symptomatic notification date and symptoms onset date (<xref rid="pone.0137959.g002" ref-type="fig">Fig 2A</xref>). The median number of days to be removed is 1 day. Using parameter estimates determined from the natural history of pandemic influenza in 2009 [<xref rid="pone.0137959.ref013" ref-type="bibr">13</xref>] [<xref rid="pone.0137959.ref014" ref-type="bibr">14</xref>] [<xref rid="pone.0137959.ref015" ref-type="bibr">15</xref>] [<xref rid="pone.0137959.ref016" ref-type="bibr">16</xref>] as shown in <xref rid="pone.0137959.t001" ref-type="table">Table 1</xref>, we re-defined our transmissibility function as follows,
<disp-formula id="pone.0137959.e017"><alternatives><graphic xlink:href="pone.0137959.e017.jpg" id="pone.0137959.e017g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M17"><mml:mrow><mml:mi>T</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mn>0</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>λ</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mi>t</mml:mi><mml:mo>−</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>l</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>λ</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msup></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mi>t</mml:mi><mml:mo><</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>l</mml:mi><mml:mo>≤</mml:mo><mml:mi>t</mml:mi><mml:mo><</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>l</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mi>t</mml:mi><mml:mo>≥</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>l</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>i</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math></alternatives></disp-formula>where Δ<italic>l</italic>, Δ<italic>s</italic> and Δ<italic>i</italic> represent the average duration of the latent phase of infection, average duration of the asymptomatic phase of infection, and the average duration of the symptomatic phase of infection respectively.</p><fig id="pone.0137959.g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0137959.g002</object-id><label>Fig 2</label><caption><title>Removal distribution and inverse cumulative function of number of close contacts.</title><p><bold>(A)</bold> Frequency distribution of the number of days required to remove infectious individuals after the date of symptoms onset, separated by age group (2–19 or 20 and above). (B) Inverse cumulative function of the number of close contacts, point estimates, and 95% confidence intervals for the variance in the number of close contacts in individuals aged 2–19 (children) and 20 or above (adults).</p></caption><graphic xlink:href="pone.0137959.g002"></graphic></fig><table-wrap id="pone.0137959.t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0137959.t001</object-id><label>Table 1</label><caption><title>Parameter values used and estimates in model A and model B.</title></caption><alternatives><graphic id="pone.0137959.t001g" xlink:href="pone.0137959.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></colgroup><thead><tr><th align="left" rowspan="1" colspan="1">Parameter values used in Model A and B</th><th align="center" rowspan="1" colspan="1">Value</th><th align="center" rowspan="1" colspan="1">Notes</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><bold>Average duration of the latent phase of infection (Δ<italic>l</italic>) (days).</bold></td><td align="center" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">Integer value within the 95% CI of the estimate (2.28, 3.12) from analysis of the 2009 influenza A (H1N1) pandemic in the province of Ontario, Canada [<xref rid="pone.0137959.ref013" ref-type="bibr">13</xref>].</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Average duration of asymptomatic phase of infection (Δ<italic>s</italic>) (days).</bold></td><td align="center" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">[<xref rid="pone.0137959.ref014" ref-type="bibr">14</xref>] [<xref rid="pone.0137959.ref015" ref-type="bibr">15</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Average duration of the symptomatic phase of infection (Δ<italic>i</italic>).</bold></td><td align="center" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">Integer value within the 95% CI of the estimate (2.06, 4.69) from analysis of the 2009 influenza A (H1N1) pandemic in the province of Ontario, Canada [<xref rid="pone.0137959.ref013" ref-type="bibr">13</xref>].</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Parameter estimates in Model A</bold></td><td align="center" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Overall excess degree in Hong Kong</bold></td><td align="center" rowspan="1" colspan="1">23.0</td><td align="left" rowspan="1" colspan="1">This value is derived from the degree distribution of 770 (<italic>n</italic>) subjects who provided either an exact or a range number of close contacts in the interview-led questionnaire survey in the Hong Kong population-based serological study in 2009 [<xref rid="pone.0137959.ref017" ref-type="bibr">17</xref>], [<xref rid="pone.0137959.ref016" ref-type="bibr">16</xref>].</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Basic reproduction number (R<sub>0</sub>)</bold></td><td align="center" rowspan="1" colspan="1">1.21</td><td align="left" rowspan="1" colspan="1">This value is extracted from the flat region in <xref rid="pone.0137959.g003" ref-type="fig">Fig 3</xref> from Day 82 to Day 102, after the first laboratory confirmed case on 30<sup>th</sup> April, 2009 (blue line in <xref rid="pone.0137959.g003" ref-type="fig">Fig 3</xref>).</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Parameter estimates in Model B</bold></td><td align="center" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Excess degree of children (less than 20 years old) in Hong Kong.(<italic>z</italic>(1))</bold></td><td align="center" rowspan="1" colspan="1">22.4</td><td align="left" rowspan="1" colspan="1">This value is derived from the degree distribution of 112 (<italic>n(1)</italic>) children who provided either an exact number of close contacts (n = 99) or a range number of close contacts (n = 13) on a pre-assigned random day in the interview-led questionnaire survey in the 2009 Hong Kong population-based serological study [<xref rid="pone.0137959.ref017" ref-type="bibr">17</xref>], [<xref rid="pone.0137959.ref016" ref-type="bibr">16</xref>].</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Excess degree of adults (20 years old or above) in Hong Kong. (<italic>z</italic>(2))</bold></td><td align="center" rowspan="1" colspan="1">18.9</td><td align="left" rowspan="1" colspan="1">This value is derived from the degree distribution of 658 (<italic>n(2)</italic>) adults who provided either an exact number of close contacts (n = 611) or a range number of close contacts (n = 47) for a pre-assigned random day in the interview-led questionnaire survey in the 2009 Hong Kong population-based serological study [<xref rid="pone.0137959.ref017" ref-type="bibr">17</xref>], [<xref rid="pone.0137959.ref016" ref-type="bibr">16</xref>].</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Overall excess degree in Hong Kong</bold></td><td align="center" rowspan="1" colspan="1">19.6</td><td align="left" rowspan="1" colspan="1">This value is derived from the weighted average of z(1) and z(2) and the proportion of the two age groups in the Hong Kong population in 2009 as per census and statistics department information.</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Basic reproduction number (R<sub>0</sub>): i) Overall; ii) Children; iii) Adults</bold></td><td align="center" rowspan="1" colspan="1">1.12; 1.28; 1.08</td><td align="left" rowspan="1" colspan="1">This value is extracted from the flat region in <xref rid="pone.0137959.g003" ref-type="fig">Fig 3</xref> from Day 82 to Day 102, after the first laboratory confirmed case on 30<sup>th</sup> April, 2009.</td></tr></tbody></table></alternatives></table-wrap><p>We derived the excess degree and degree distribution for each age group from an analysis of the number of close contacts, <italic>k</italic>, experienced by 770 participants aged 3 to 102 years taken from the interview-led social contact questionnaire in [<xref rid="pone.0137959.ref017" ref-type="bibr">17</xref>] [<xref rid="pone.0137959.ref018" ref-type="bibr">18</xref>]. The subjects in this study were required to record detailed contact histories for one pre-assigned day, providing information on how many contacts they had made, whom they had met, and how long they had met with each contact. Out of 770 subjects, 710 subjects recorded how many close contacts they had made. The rest of the 60 subjects provided a range in the number of close contacts they had experienced, and we defined their number of close contacts as the mid-point of this range. The mean number of close contacts for school-aged children and adults was 9.94 and 4.69 respectively. The degree distribution for individuals aged 2–19 is different from those aged 20 and above as shown in <xref rid="pone.0137959.g002" ref-type="fig">Fig 2B</xref>. In the sub-diagram of <xref rid="pone.0137959.g002" ref-type="fig">Fig 2B</xref>, the variance in the number of contacts of children and adults is 115.7 [90.4, 153.4] and 94.7 [85.3, 105.9] respectively. This shows that there is higher variation in the number of contacts in children than in adults.</p><p><xref rid="pone.0137959.g003" ref-type="fig">Fig 3</xref> reveals that while the daily incidence of laboratory-confirmed cases of pH1N1 among adults in Hong Kong increased monotonically from April until the main peak in September 2009, the incidence in children presented an anomaly. The incidence data for children has a minor peak at the early stage of the epidemic between June and July. In real time, at the time when the minor peak formed, and without prior knowledge of how the epidemic would unfold afterwards, this anomaly might have been interpreted as a stand-alone peak, with an inclining segment between June 1 and 21 and a declining segment between June 21 and July 1 as shown in the inset figure in <xref rid="pone.0137959.g003" ref-type="fig">Fig 3</xref>. Interestingly, using the inclining segment with model B, the estimate of R<sub>0</sub> converges very quickly to 1.2 in early June (black line in <xref rid="pone.0137959.g003" ref-type="fig">Fig 3</xref>). Since the declining segment destabilizes the early-stage R<sub>0</sub> estimator, the estimate after the minor peak fluctuates for a while until the epidemic curve begins to increase monotonically in July, when the estimator converges again to 1.12 between day 82 and day 102. The estimate of R<sub>0</sub> under model A converges to 1.21 in late August (Day 115), which is slightly higher than the estimate from model B. Estimates under both model A and B are similar to the estimate of the effective reproduction number between 1.1 and 1.2 that was proposed in summer of 2009 by [<xref rid="pone.0137959.ref018" ref-type="bibr">18</xref>]. The number of cases required for model B to converge to estimate R<sub>0</sub> of 1.2 and 1.12 up to Day 30 and 82 are 66 and ~2200, respectively, which are about 5 to 10% of the total number of laboratory confirmed cases in the whole epidemic. Using the excess degree of children and adults, the age specific R<sub>0</sub> can be estimated for both children and adults. We showed that R<sub>0</sub> for children is higher than that for adults (1.28 versus 1.08), which is consistent with other studies (see [<xref rid="pone.0137959.ref019" ref-type="bibr">19</xref>]). This implies that school-aged children likely contributed more infections than adults.</p><fig id="pone.0137959.g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0137959.g003</object-id><label>Fig 3</label><caption><title>Incidence of laboratory confirmed H1N1pdm cases during the 2009 pandemic influenza and real time estimation of R<sub>0</sub>.</title><p>Stacked bar chart shows the incidence of laboratory confirmed H1N1pdm cases in individuals aged 2–19 years and aged 20 or above. The line chart shows the real time estimation of R<sub>0</sub> using model B.</p></caption><graphic xlink:href="pone.0137959.g003"></graphic></fig></sec><sec sec-type="conclusions" id="sec008"><title>Discussion</title><p>In this paper, we examined the rate of convergence when estimating the basic reproduction number, R<sub>0</sub>, from case notification data during the early stage of an infectious disease outbreak. By explicitly taking in account the age structure of a population, this study improves on a network-based formalism that was introduced previously for a non-age-stratified population in [<xref rid="pone.0137959.ref011" ref-type="bibr">11</xref>]. Our proposed age-structured model (model B) provides a reliable estimate of R<sub>0</sub> in the early stage of an epidemic, even in the presence of noise or dominant stochastic fluctuations and potentially biased numbers of newly reported cases during the early stochastic regime. Using the simulated data, our analytical justification shows that network structures featuring inter-age group contact heterogeneity helps us provide an ongoing estimate of R<sub>0</sub> before the epidemic reaches the exponential growth as in <xref rid="pone.0137959.g001" ref-type="fig">Fig 1</xref>. Even for the real data, our proposed model B can still provide estimates of R<sub>0</sub> in the early stage of epidemics as of Day 30. Although the first peak deviates from the estimate of R<sub>0</sub>, it quickly converges again on the final estimate two and half months before the second peak in mid-September. The ongoing estimate of R<sub>0</sub> provides important information for policy makers to take immediate action to control further disease spread.</p><p>One of the important inputs to model B is the age-specific contact information. Thus, conducting a large-scale population-based survey on contact behaviour is suggested in every city. For densely populated cities like Hong Kong, contacts of individuals within and between different age groups may vary significantly. For this reason, a large scale population-based telephone survey in Hong Kong has recently been conducted to collect the contact histories of individuals at four different time points in a year, to examine the temporal changes of contacts, and generate unbiased estimates of contact between and among age groups [<xref rid="pone.0137959.ref020" ref-type="bibr">20</xref>].While model A provides a reasonable estimate for R<sub>0</sub> in the general population as a whole, model B–which takes age structure into account–portrays a more refined picture of patterns of disease spread across different age groups, thus providing an age-specific estimate of R<sub>0</sub>. To better prepare for re-emerging and emerging infectious diseases, the use of model B and further studies on the temporal change of contacts’ behaviour in major metropolitan areas may help refine these estimates even further, and provide more accurate insight into future outbreaks of infectious diseases in the early phase.</p><p>As we have shown in <xref rid="pone.0137959.g001" ref-type="fig">Fig 1A–1C</xref>, and <xref rid="pone.0137959.g001" ref-type="fig">Fig 1D–1F</xref>, estimates from model A are sensitive to the degree of contact heterogeneity while estimates from model B are not. For simplicity, we have illustrated different scenarios of degree distribution for only two age groups to test both model A and B. However, real populations have far more complex structures with very different contact patterns. Specifically, the degree distribution of individuals in each age group is neither uniform nor non-overlapping. The analytical/numerical framework presented in this paper can be readily applied to situations where these constraints are relaxed.</p><p>Our estimate of R<sub>0</sub> relies on assumptions of discrete time, constant infectivity, and removal functions. Also, there are some extrinsic factors, which may raise additional uncertainties. For example, we have not included the possibility of late reporting or underreporting of numbers of new cases in the early phases of an epidemic. The assumption that the number of reported cases accurately reflects the actual number of cases may not hold true in real situations since first, not all infected persons will seek medical attention; and second, there are strong biases for certain age groups (such as children) to seek medical attention even in a precautionary capacity. To our knowledge, this model is the first network-based modeling framework that includes heterogeneity in the estimation of R<sub>0</sub>. In this framework the confidence interval of R<sub>0,</sub> is not constructed because the stochastic variables such as <italic>I</italic><sup><italic>r</italic></sup>(<italic>t</italic>), <italic>R</italic><sup><italic>r</italic></sup>(<italic>t</italic>), <italic>T</italic><sup><italic>r</italic></sup>(<italic>t</italic>) used from Eqs (<xref rid="pone.0137959.e008" ref-type="disp-formula">1</xref>–<xref rid="pone.0137959.e015" ref-type="disp-formula">5</xref>) involves inherent uncertainties in arriving the estimator of R<sub>0</sub> in addition to the above extrinsic factors.</p><p>The proposed methodology in its current form, does not address an assessment of the uncertainty in the point estimate of R<sub>0</sub> at any given time; to achieve a more robust policy decision-support tool, the development of accompanying methodology to determine confidence intervals should be undertaken in future work.</p><p>Despite uncertainties about non-perfect observations collected from the field during an outbreak, this approach could still make an important contribution to the estimation of the basic reproduction number in the early phases of an epidemic, assuming the proportion of underreporting is not changing substantially. Improving case notification through real-time surveillance to obtain more accurate case data, combined with the proposed modeling framework, could provide an invaluable platform to inform policy makers by estimating a more accurate R<sub>0</sub> in real time, thus allowing policy makers to take immediate action to control epidemics while still in the early phases.</p></sec><sec sec-type="supplementary-material" id="sec009"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0137959.s001"><label>S1 Data</label><caption><title>Age-specific daily count of reported pH1N1 (2009) cases in Hong Kong from April 30, 2009 to February 7, 2010.</title><p>(CSV)</p></caption><media xlink:href="pone.0137959.s001.csv"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>B.P. would like to acknowledge the support of the Canadian Institutes of Health Research (Grant No. PTL-97126), the Michael Smith Foundation for Health Research (Grant No. CI-SSH-00931(06–1)), and the British Columbia Ministry of Health (Pandemic Preparedness Modeling Project). SR thanks the Medical Research Council (UK, Project Grant MR/J008761/1), Wellcome Trust (Project Award 093488/Z/10/Z) and the National Institute for Health Research (UK, Health Protection Research Unit). K.O.K. would like to acknowledge the support from Fogarty International Centre (FluScape, Grant No. R01 TW008246-01), Food and Health Bureau (Research Fund for the Control of Infectious Diseases), Government of the Hong Kong SAR (RFCID No. 110642 and PHE 21) and Research Grants Council, Hong Kong (Grant No. 776810). K.O.K. also thanks Prof. Ben Cowling for helpful discussion, Vivian Wei for graphics support and Marie Chi for administrative work.</p></ack><ref-list><title>References</title><ref id="pone.0137959.ref001"><label>1</label><mixed-citation publication-type="journal"><name><surname>Fraser</surname><given-names>C</given-names></name>, <name><surname>Donnelly</surname><given-names>CA</given-names></name>, <name><surname>Cauchemez</surname><given-names>S</given-names></name>, <name><surname>Hanage</surname><given-names>WP</given-names></name>, <name><surname>Van Kerkhove</surname><given-names>MD</given-names></name>, <name><surname>Hollingsworth</surname><given-names>TD</given-names></name>, <etal>et al</etal><article-title>Pandemic potential of a strain of influenza A (H1N1): early findings</article-title>. <source>Science</source>. <year>2009</year>;<volume>324</volume>(<issue>5934</issue>):<fpage>1557</fpage>–<lpage>61</lpage>. <pub-id pub-id-type="doi">10.1126/science.1176062</pub-id><pub-id pub-id-type="pmid">19433588</pub-id></mixed-citation></ref><ref id="pone.0137959.ref002"><label>2</label><mixed-citation publication-type="journal"><name><surname>Riley</surname><given-names>S</given-names></name>, <name><surname>Fraser</surname><given-names>C</given-names></name>, <name><surname>Donnelly</surname><given-names>CA</given-names></name>, <name><surname>Ghani</surname><given-names>AC</given-names></name>, <name><surname>Abu-Raddad</surname><given-names>LJ</given-names></name>, <name><surname>Hedley</surname><given-names>AJ</given-names></name>, <etal>et al</etal><article-title>Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions</article-title>. <source>Science</source>. <year>2003</year>;<volume>300</volume>(<issue>5627</issue>):<fpage>1961</fpage>–<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1126/science.1086478</pub-id>.<pub-id pub-id-type="pmid">12766206</pub-id></mixed-citation></ref><ref id="pone.0137959.ref003"><label>3</label><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>L-D</given-names></name>, <name><surname>Chen</surname><given-names>H-G</given-names></name>, <name><surname>Guo</surname><given-names>J-G</given-names></name>, <name><surname>Zeng</surname><given-names>X-J</given-names></name>, <name><surname>Hong</surname><given-names>X-L</given-names></name>, <name><surname>Xiong</surname><given-names>J-J</given-names></name>, <etal>et al</etal><article-title>A strategy to control transmission of Schistosoma japonicum in China</article-title>. <source>N Engl J Med</source>. <year>2009</year>;<volume>360</volume>(<issue>2</issue>):<fpage>121</fpage>–<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1056/NEJMoa0800135</pub-id>.<pub-id pub-id-type="pmid">19129526</pub-id></mixed-citation></ref><ref id="pone.0137959.ref004"><label>4</label><mixed-citation publication-type="journal"><name><surname>Wallinga</surname><given-names>J</given-names></name>, <name><surname>Van Boven</surname><given-names>M</given-names></name>, <name><surname>Lipsitch</surname><given-names>M</given-names></name>. <article-title>Supporting Information: Optimizing infectious disease interventions during an emerging epidemic</article-title>. <source>Proceedings of the National Academy of Sciences</source>. <year>2010</year>.</mixed-citation></ref><ref id="pone.0137959.ref005"><label>5</label><mixed-citation publication-type="journal"><name><surname>Pourbohloul</surname><given-names>B</given-names></name>, <name><surname>Ahued</surname><given-names>A</given-names></name>, <name><surname>Davoudi</surname><given-names>B</given-names></name>, <name><surname>Meza</surname><given-names>R</given-names></name>, <name><surname>Meyers</surname><given-names>LA</given-names></name>, <name><surname>Skowronski</surname><given-names>DM</given-names></name>, <etal>et al</etal><article-title>Initial human transmission dynamics of the pandemic (H1N1) 2009 virus in North America</article-title>. <source>Influenza Other Respi Viruses</source>. <year>2009</year>;<volume>3</volume>(<issue>5</issue>):<fpage>215</fpage>–<lpage>22</lpage>. .<pub-id pub-id-type="pmid">19702583</pub-id></mixed-citation></ref><ref id="pone.0137959.ref006"><label>6</label><mixed-citation publication-type="book"><name><surname>Anderson</surname><given-names>RM</given-names></name>, <name><surname>May</surname><given-names>RM</given-names></name>. <chapter-title>Infectious disease of humans: dynamics and control</chapter-title><publisher-name>Oxford Science Publications</publisher-name><year>1992</year>.</mixed-citation></ref><ref id="pone.0137959.ref007"><label>7</label><mixed-citation publication-type="journal"><name><surname>Meyers</surname><given-names>LA</given-names></name>, <name><surname>Pourbohloul</surname><given-names>B</given-names></name>, <name><surname>Newman</surname><given-names>ME</given-names></name>, <name><surname>Skowronski</surname><given-names>DM</given-names></name>, <name><surname>Brunham</surname><given-names>RC</given-names></name>. <article-title>Network theory and SARS: predicting outbreak diversity</article-title>. <source>J Theor Biol</source>. <year>2005</year>;<volume>232</volume>(<issue>1</issue>):<fpage>71</fpage>–<lpage>81</lpage>. .<pub-id pub-id-type="pmid">15498594</pub-id></mixed-citation></ref><ref id="pone.0137959.ref008"><label>8</label><mixed-citation publication-type="journal"><name><surname>Donnelly</surname><given-names>CA</given-names></name>, <name><surname>Ghani</surname><given-names>AC</given-names></name>, <name><surname>Leung</surname><given-names>GM</given-names></name>, <name><surname>Hedley</surname><given-names>AJ</given-names></name>, <name><surname>Fraser</surname><given-names>C</given-names></name>, <name><surname>Riley</surname><given-names>S</given-names></name>, <etal>et al</etal><article-title>Epidemiological determinants of spread of causal agent of severe acute respiratory syndrome in Hong Kong</article-title>. <source>Lancet</source>. <year>2003</year>;<volume>361</volume>(<issue>9371</issue>):<fpage>1761</fpage>–<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1016/S0140-6736(03)13410-1</pub-id>.<pub-id pub-id-type="pmid">12781533</pub-id></mixed-citation></ref><ref id="pone.0137959.ref009"><label>9</label><mixed-citation publication-type="journal"><name><surname>Wu</surname><given-names>JT</given-names></name>, <name><surname>Cowling</surname><given-names>BJ</given-names></name>, <name><surname>Lau</surname><given-names>EHY</given-names></name>, <name><surname>M Ip</surname><given-names>DK</given-names></name>, <name><surname>Ho</surname><given-names>LM</given-names></name>, <name><surname>Tsang</surname><given-names>T</given-names></name>, <etal>et al</etal><article-title>School closure and mitigation of pandemic (H1N1) 2009, Hong Kong</article-title>. <source>Emerging Infect Dis</source>. <year>2010</year>;<volume>16</volume>(<issue>3</issue>):<fpage>538</fpage>–<lpage>41</lpage>. <pub-id pub-id-type="doi">10.3201/eid1603.091216</pub-id><pub-id pub-id-type="pmid">20202441</pub-id></mixed-citation></ref><ref id="pone.0137959.ref010"><label>10</label><mixed-citation publication-type="journal"><name><surname>Noël</surname><given-names>PA</given-names></name>, <name><surname>Davoudi</surname><given-names>B</given-names></name>, <name><surname>Brunham</surname><given-names>RC</given-names></name>, <name><surname>Dube</surname><given-names>LJ</given-names></name>, <name><surname>Pourbohloul</surname><given-names>B</given-names></name>. <article-title>Time evolution of epidemic disease on finite and infinite networks</article-title>. <source>Physical review E, Statistical, nonlinear, and soft matter physics</source>. <year>2009</year>;<volume>79</volume>(<issue>2</issue> Pt 2):<fpage>026101</fpage> Epub 2009/04/28. .<pub-id pub-id-type="pmid">19391800</pub-id></mixed-citation></ref><ref id="pone.0137959.ref011"><label>11</label><mixed-citation publication-type="journal"><name><surname>Davoudi</surname><given-names>B</given-names></name>, <name><surname>Miller</surname><given-names>JC</given-names></name>, <name><surname>Meza</surname><given-names>R</given-names></name>, <name><surname>Meyers</surname><given-names>LA</given-names></name>, <name><surname>Earn</surname><given-names>DJD</given-names></name>, <name><surname>Pourbohloul</surname><given-names>B</given-names></name>. <article-title>Early Real-Time Estimation of the Basic Reproduction Number of Emerging Infectious Diseases</article-title>. <source>Physical Review X</source>. <year>2012</year>;<volume>2</volume>(<issue>3</issue>):<fpage>031005</fpage>.</mixed-citation></ref><ref id="pone.0137959.ref012"><label>12</label><mixed-citation publication-type="book"><name><surname>Press</surname><given-names>WH</given-names></name>, <name><surname>Teukolsky</surname><given-names>SA</given-names></name>, <name><surname>Vetterling</surname><given-names>WT</given-names></name>, <name><surname>Flannery</surname><given-names>BP</given-names></name>. <source>Numerical Recipes 3rd Edition: The Art of Scientific Computing</source>: <publisher-name>Cambridge University Press</publisher-name>; <year>2007</year>. 1256 p.</mixed-citation></ref><ref id="pone.0137959.ref013"><label>13</label><mixed-citation publication-type="journal"><name><surname>Tuite</surname><given-names>AR</given-names></name>, <name><surname>Greer</surname><given-names>AL</given-names></name>, <name><surname>Whelan</surname><given-names>M</given-names></name>, <name><surname>Winter</surname><given-names>AL</given-names></name>, <name><surname>Lee</surname><given-names>B</given-names></name>, <name><surname>Yan</surname><given-names>P</given-names></name>, <etal>et al</etal><article-title>Estimated epidemiologic parameters and morbidity associated with pandemic H1N1 influenza</article-title>. <source>CMAJ: Canadian Medical Association journal = journal de l'Association medicale canadienne</source>. <year>2010</year>;<volume>182</volume>(<issue>2</issue>):<fpage>131</fpage>–<lpage>6</lpage>. Epub 2009/12/05. <pub-id pub-id-type="doi">10.1503/cmaj.091807</pub-id><pub-id pub-id-type="pmid">19959592</pub-id></mixed-citation></ref><ref id="pone.0137959.ref014"><label>14</label><mixed-citation publication-type="journal"><name><surname>Ferguson</surname><given-names>NM</given-names></name>, <name><surname>Cummings</surname><given-names>DAT</given-names></name>, <name><surname>Cauchemez</surname><given-names>S</given-names></name>, <name><surname>Fraser</surname><given-names>C</given-names></name>, <name><surname>Riley</surname><given-names>S</given-names></name>, <name><surname>Meeyai</surname><given-names>A</given-names></name>, <etal>et al</etal><article-title>Strategies for containing an emerging influenza pandemic in Southeast Asia</article-title>. <source>Nature</source>. <year>2005</year>;<volume>437</volume>(<issue>7056</issue>):<fpage>209</fpage>–<lpage>14</lpage>. .<pub-id pub-id-type="pmid">16079797</pub-id></mixed-citation></ref><ref id="pone.0137959.ref015"><label>15</label><mixed-citation publication-type="journal"><name><surname>Longini</surname><given-names>IM</given-names><suffix>Jr.</suffix></name>, <name><surname>Nizam</surname><given-names>A</given-names></name>, <name><surname>Xu</surname><given-names>S</given-names></name>, <name><surname>Ungchusak</surname><given-names>K</given-names></name>, <name><surname>Hanshaoworakul</surname><given-names>W</given-names></name>, <name><surname>Cummings</surname><given-names>DA</given-names></name>, <etal>et al</etal><article-title>Containing pandemic influenza at the source</article-title>. <source>Science</source>. <year>2005</year>;<volume>309</volume>(<issue>5737</issue>):<fpage>1083</fpage>–<lpage>7</lpage>. Epub 2005/08/05. <pub-id pub-id-type="doi">10.1126/science.1115717</pub-id>.<pub-id pub-id-type="pmid">16079251</pub-id></mixed-citation></ref><ref id="pone.0137959.ref016"><label>16</label><mixed-citation publication-type="journal"><name><surname>Riley</surname><given-names>S</given-names></name>, <name><surname>Kwok</surname><given-names>KO</given-names></name>, <name><surname>Wu</surname><given-names>KM</given-names></name>, <name><surname>Ning</surname><given-names>DY</given-names></name>, <name><surname>Cowling</surname><given-names>BJ</given-names></name>, <name><surname>Wu</surname><given-names>JT</given-names></name>, <etal>et al</etal><article-title>Epidemiological characteristics of 2009 (H1N1) pandemic influenza based on paired sera from a longitudinal community cohort study</article-title>. <source>PLoS Med</source>. <year>2011</year>;<volume>8</volume>(<issue>6</issue>):<fpage>e1000442</fpage><pub-id pub-id-type="doi">10.1371/journal.pmed.1000442</pub-id><pub-id pub-id-type="pmid">21713000</pub-id></mixed-citation></ref><ref id="pone.0137959.ref017"><label>17</label><mixed-citation publication-type="journal"><name><surname>Kwok</surname><given-names>KO</given-names></name>, <name><surname>Cowling</surname><given-names>BJ</given-names></name>, <name><surname>Peiris</surname><given-names>JS</given-names></name>, <name><surname>Wei</surname><given-names>V</given-names></name>, <name><surname>Wu</surname><given-names>KM</given-names></name>, <name><surname>Ning</surname><given-names>DY</given-names></name>, <etal>et al</etal><article-title>Social contacts and the locations in which they occur as risk factors for influenza infection</article-title>. <source>Proceedings of the Royal Society B</source>. <year>2014</year>.</mixed-citation></ref><ref id="pone.0137959.ref018"><label>18</label><mixed-citation publication-type="journal"><name><surname>Cowling</surname><given-names>BJ</given-names></name>, <name><surname>Lau</surname><given-names>MS</given-names></name>, <name><surname>Ho</surname><given-names>LM</given-names></name>, <name><surname>Chuang</surname><given-names>SK</given-names></name>, <name><surname>Tsang</surname><given-names>T</given-names></name>, <name><surname>Liu</surname><given-names>SH</given-names></name>, <etal>et al</etal><article-title>The effective reproduction number of pandemic influenza: prospective estimation</article-title>. <source>Epidemiology</source>. <year>2010</year>;<volume>21</volume>(<issue>6</issue>):<fpage>842</fpage>–<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1097/EDE.0b013e3181f20977</pub-id><pub-id pub-id-type="pmid">20805752</pub-id></mixed-citation></ref><ref id="pone.0137959.ref019"><label>19</label><mixed-citation publication-type="journal"><name><surname>Glass</surname><given-names>K</given-names></name>, <name><surname>Kelly</surname><given-names>H</given-names></name>, <name><surname>Mercer</surname><given-names>GN</given-names></name>. <article-title>Pandemic influenza H1N1: reconciling serosurvey data with estimates of the reproduction number</article-title>. <source>Epidemiology</source>. <year>2012</year>;<volume>23</volume>(<issue>1</issue>):<fpage>86</fpage>–<lpage>94</lpage>. Epub 2011/11/18. <pub-id pub-id-type="doi">10.1097/EDE.0b013e31823a44a5</pub-id>.<pub-id pub-id-type="pmid">22089631</pub-id></mixed-citation></ref><ref id="pone.0137959.ref020"><label>20</label><mixed-citation publication-type="journal"><name><surname>Kwok</surname><given-names>KO</given-names></name>, <name><surname>Cowling</surname><given-names>BJ</given-names></name>, <name><surname>Wei</surname><given-names>V</given-names></name>, <name><surname>Riley</surname><given-names>S</given-names></name>, <name><surname>Read</surname><given-names>J</given-names></name>. <article-title>Temporal patterns of human social encounters and their potential exposure to infectious disease: A longitudinal study of Hong Kong residents</article-title>. <year>2014</year>.</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 001346 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 001346 | 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é=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |