New coronavirus outbreak. Lessons learned from the severe acute respiratory syndrome epidemic
Identifieur interne : 000C50 ( Pmc/Corpus ); précédent : 000C49; suivant : 000C51New coronavirus outbreak. Lessons learned from the severe acute respiratory syndrome epidemic
Auteurs : E. Álvarez ; J. Donado-Campos ; F. MorillaSource :
- Epidemiology and Infection [ 0950-2688 ] ; 2015.
Abstract
System dynamics approach offers great potential for addressing how intervention policies can affect the spread of emerging infectious diseases in complex and highly networked systems. Here, we develop a model that explains the severe acute respiratory syndrome coronavirus (SARS-CoV) epidemic that occurred in Hong Kong in 2003. The dynamic model developed with system dynamics methodology included 23 variables (five states, four flows, eight auxiliary variables, six parameters), five differential equations and 12 algebraic equations. The parameters were optimized following an iterative process of simulation to fit the real data from the epidemics. Univariate and multivariate sensitivity analyses were performed to determine the reliability of the model. In addition, we discuss how further testing using this model can inform community interventions to reduce the risk in current and future outbreaks, such as the recently Middle East respiratory syndrome coronavirus (MERS-CoV) epidemic.
Url:
DOI: 10.1017/S095026881400377X
PubMed: 25591619
PubMed Central: 7108648
Links to Exploration step
PMC:7108648Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">New coronavirus outbreak. Lessons learned from the severe acute respiratory syndrome epidemic</title>
<author><name sortKey="Alvarez, E" sort="Alvarez, E" uniqKey="Alvarez E" first="E." last="Álvarez">E. Álvarez</name>
<affiliation><nlm:aff id="aff1"><addr-line>Centro de Biología Molecular Severo Ochoa</addr-line>
,<addr-line>Consejo Superior de Investigaciones Científicas</addr-line>
,<institution>Universidad Autónoma de Madrid (CSIC-UAM)</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Donado Campos, J" sort="Donado Campos, J" uniqKey="Donado Campos J" first="J." last="Donado-Campos">J. Donado-Campos</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Preventive Medicine and Public Health</addr-line>
,<addr-line>School of Medicine</addr-line>
,<institution>Universidad Autónoma de Madrid/IdiPAZ</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Morilla, F" sort="Morilla, F" uniqKey="Morilla F" first="F." last="Morilla">F. Morilla</name>
<affiliation><nlm:aff id="aff3"><addr-line>Departamento de Informática y Automática,ETSI Informática</addr-line>
,<institution>Universidad Nacional de Educación a Distancia UNED</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">25591619</idno>
<idno type="pmc">7108648</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC7108648</idno>
<idno type="RBID">PMC:7108648</idno>
<idno type="doi">10.1017/S095026881400377X</idno>
<date when="2015">2015</date>
<idno type="wicri:Area/Pmc/Corpus">000C50</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000C50</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">New coronavirus outbreak. Lessons learned from the severe acute respiratory syndrome epidemic</title>
<author><name sortKey="Alvarez, E" sort="Alvarez, E" uniqKey="Alvarez E" first="E." last="Álvarez">E. Álvarez</name>
<affiliation><nlm:aff id="aff1"><addr-line>Centro de Biología Molecular Severo Ochoa</addr-line>
,<addr-line>Consejo Superior de Investigaciones Científicas</addr-line>
,<institution>Universidad Autónoma de Madrid (CSIC-UAM)</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Donado Campos, J" sort="Donado Campos, J" uniqKey="Donado Campos J" first="J." last="Donado-Campos">J. Donado-Campos</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Preventive Medicine and Public Health</addr-line>
,<addr-line>School of Medicine</addr-line>
,<institution>Universidad Autónoma de Madrid/IdiPAZ</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Morilla, F" sort="Morilla, F" uniqKey="Morilla F" first="F." last="Morilla">F. Morilla</name>
<affiliation><nlm:aff id="aff3"><addr-line>Departamento de Informática y Automática,ETSI Informática</addr-line>
,<institution>Universidad Nacional de Educación a Distancia UNED</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Epidemiology and Infection</title>
<idno type="ISSN">0950-2688</idno>
<idno type="eISSN">1469-4409</idno>
<imprint><date when="2015">2015</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><title>SUMMARY</title>
<p>System dynamics approach offers great potential for addressing how intervention policies can affect the spread of emerging infectious diseases in complex and highly networked systems. Here, we develop a model that explains the severe acute respiratory syndrome coronavirus (SARS-CoV) epidemic that occurred in Hong Kong in 2003. The dynamic model developed with system dynamics methodology included 23 variables (five states, four flows, eight auxiliary variables, six parameters), five differential equations and 12 algebraic equations. The parameters were optimized following an iterative process of simulation to fit the real data from the epidemics. Univariate and multivariate sensitivity analyses were performed to determine the reliability of the model. In addition, we discuss how further testing using this model can inform community interventions to reduce the risk in current and future outbreaks, such as the recently Middle East respiratory syndrome coronavirus (MERS-CoV) epidemic.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="De Groot, Rj" uniqKey="De Groot R">RJ de Groot</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Zaki, Am" uniqKey="Zaki A">AM Zaki</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Assiri, A" uniqKey="Assiri A">A Assiri</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mailles, A" uniqKey="Mailles A">A Mailles</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Drosten, C" uniqKey="Drosten C">C Drosten</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Marra, Ma" uniqKey="Marra M">MA Marra</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Peiris, Js" uniqKey="Peiris J">JS Peiris</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Annan, A" uniqKey="Annan A">A Annan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cotten, M" uniqKey="Cotten M">M Cotten</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, W" uniqKey="Li W">W Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chau, Ph" uniqKey="Chau P">PH Chau</name>
</author>
<author><name sortKey="Yip, Ps" uniqKey="Yip P">PS Yip</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sterman, J" uniqKey="Sterman J">J Sterman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Homer, Jb" uniqKey="Homer J">JB Homer</name>
</author>
<author><name sortKey="Hirsch, Gb" uniqKey="Hirsch G">GB Hirsch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jones, Ap" uniqKey="Jones A">AP Jones</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Koopman, Js" uniqKey="Koopman J">JS Koopman</name>
</author>
<author><name sortKey="Jacquez, G" uniqKey="Jacquez G">G Jacquez</name>
</author>
<author><name sortKey="Chick, Se" uniqKey="Chick S">SE Chick</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="Omony, J" uniqKey="Omony J">J Omony</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Pitzer, Ve" uniqKey="Pitzer V">VE Pitzer</name>
</author>
<author><name sortKey="Leung, Gm" uniqKey="Leung G">GM Leung</name>
</author>
<author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcbryde, Es" uniqKey="Mcbryde E">ES McBryde</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hirose, H" uniqKey="Hirose H">H Hirose</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fine, Pe" uniqKey="Fine P">PE Fine</name>
</author>
<author><name sortKey="Clarkson, Ja" uniqKey="Clarkson J">JA Clarkson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Anderson, Rm" uniqKey="Anderson R">RM Anderson</name>
</author>
<author><name sortKey="Grenfell, Bt" uniqKey="Grenfell B">BT Grenfell</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="Edmunds, Wj" uniqKey="Edmunds W">WJ Edmunds</name>
</author>
<author><name sortKey="O Callaghan, Cj" uniqKey="O Callaghan C">CJ O'Callaghan</name>
</author>
<author><name sortKey="Nokes, Dj" uniqKey="Nokes D">DJ Nokes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mossong, J" uniqKey="Mossong J">J Mossong</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
<author><name sortKey="Teunis, P" uniqKey="Teunis P">P Teunis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Royston, G" uniqKey="Royston G">G Royston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Morilla, F" uniqKey="Morilla F">F Morilla</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hirsch, G" uniqKey="Hirsch G">G Hirsch</name>
</author>
<author><name sortKey="Miller, S" uniqKey="Miller S">S Miller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Luginbuhl, Wh" uniqKey="Luginbuhl W">WH Luginbuhl</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Roberts, C" uniqKey="Roberts C">C Roberts</name>
</author>
<author><name sortKey="Dangerfield, B" uniqKey="Dangerfield B">B Dangerfield</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dangerfield, B" uniqKey="Dangerfield B">B Dangerfield</name>
</author>
<author><name sortKey="Fang, Y" uniqKey="Fang Y">Y Fang</name>
</author>
<author><name sortKey="Roberts, C" uniqKey="Roberts C">C Roberts</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Halder, N" uniqKey="Halder N">N Halder</name>
</author>
<author><name sortKey="Kelso, Jk" uniqKey="Kelso J">JK Kelso</name>
</author>
<author><name sortKey="Milne, Gj" uniqKey="Milne G">GJ Milne</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pitman, Rj" uniqKey="Pitman R">RJ Pitman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Grassly, Nc" uniqKey="Grassly N">NC Grassly</name>
</author>
<author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hui, Ds" uniqKey="Hui D">DS Hui</name>
</author>
<author><name sortKey="Memish, Za" uniqKey="Memish Z">ZA Memish</name>
</author>
<author><name sortKey="Zumla, A" uniqKey="Zumla A">A Zumla</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">Epidemiol Infect</journal-id>
<journal-id journal-id-type="iso-abbrev">Epidemiol. Infect</journal-id>
<journal-id journal-id-type="publisher-id">HYG</journal-id>
<journal-title-group><journal-title>Epidemiology and Infection</journal-title>
</journal-title-group>
<issn pub-type="ppub">0950-2688</issn>
<issn pub-type="epub">1469-4409</issn>
<publisher><publisher-name>Cambridge University Press</publisher-name>
<publisher-loc>Cambridge, UK</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">25591619</article-id>
<article-id pub-id-type="pmc">7108648</article-id>
<article-id pub-id-type="doi">10.1017/S095026881400377X</article-id>
<article-id pub-id-type="pii">S095026881400377X</article-id>
<article-id pub-id-type="publisher-id">00377</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Original Papers</subject>
</subj-group>
<subj-group subj-group-type="section"><subject>Coronavirus</subject>
</subj-group>
</article-categories>
<title-group><article-title>New coronavirus outbreak. Lessons learned from the severe acute respiratory syndrome epidemic</article-title>
<alt-title alt-title-type="left-running">E. Álvarez and others</alt-title>
<alt-title alt-title-type="right-running">SARS-CoV outbreak modelling</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>ÁLVAREZ</surname>
<given-names>E.</given-names>
</name>
<xref ref-type="aff" rid="aff1">1</xref>
<xref ref-type="corresp" rid="cor1">*</xref>
</contrib>
<contrib contrib-type="author"><name><surname>DONADO-CAMPOS</surname>
<given-names>J.</given-names>
</name>
<xref ref-type="aff" rid="aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>MORILLA</surname>
<given-names>F.</given-names>
</name>
<xref ref-type="aff" rid="aff3">3</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>Centro de Biología Molecular Severo Ochoa</addr-line>
,<addr-line>Consejo Superior de Investigaciones Científicas</addr-line>
,<institution>Universidad Autónoma de Madrid (CSIC-UAM)</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>Department of Preventive Medicine and Public Health</addr-line>
,<addr-line>School of Medicine</addr-line>
,<institution>Universidad Autónoma de Madrid/IdiPAZ</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>Departamento de Informática y Automática,ETSI Informática</addr-line>
,<institution>Universidad Nacional de Educación a Distancia UNED</institution>
,<addr-line>Madrid</addr-line>
,<country>Spain</country>
</aff>
<author-notes><corresp id="cor1"><label>*</label>
<addr-line>Author for correspondence: Dr E. Álvarez</addr-line>
, <addr-line>Centro de Biología Molecular Severo Ochoa</addr-line>
. <addr-line>Consejo Superior de Investigaciones Científicas</addr-line>
. <institution>Universidad Autónoma de Madrid (CSIC-UAM)</institution>
. <addr-line>C/ Nicolás Cabrera 1</addr-line>
, <addr-line>28049 Madrid</addr-line>
, <country>Spain</country>
. (Email: <email>ealvarezmba@gmail.com</email>
)</corresp>
</author-notes>
<pub-date pub-type="collection"><month>10</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="epub"><day>16</day>
<month>1</month>
<year>2015</year>
</pub-date>
<volume>143</volume>
<issue>13</issue>
<fpage>2882</fpage>
<lpage>2893</lpage>
<history><date date-type="received"><day>17</day>
<month>5</month>
<year>2014</year>
</date>
<date date-type="rev-recd"><day>06</day>
<month>11</month>
<year>2014</year>
</date>
<date date-type="accepted"><day>12</day>
<month>12</month>
<year>2014</year>
</date>
</history>
<permissions><copyright-statement>© Cambridge University Press 2015</copyright-statement>
<copyright-year>2015</copyright-year>
<copyright-holder>Cambridge University Press</copyright-holder>
<license license-type="open-access"><license-p>This article is made available via the PMC Open Access Subset for unrestricted re-use and analyses in any form or by any means subject to acknowledgement of the original source. These permissions are granted for the duration of the COVID-19 pandemic or until permissions are revoked in writing. Upon expiration of these permissions, PMC is granted a perpetual license to make this article available via PMC and Europe PMC, consistent with existing copyright protections.</license-p>
</license>
</permissions>
<self-uri xlink:title="pdf" xlink:href="S095026881400377Xa.pdf"></self-uri>
<abstract abstract-type="normal"><title>SUMMARY</title>
<p>System dynamics approach offers great potential for addressing how intervention policies can affect the spread of emerging infectious diseases in complex and highly networked systems. Here, we develop a model that explains the severe acute respiratory syndrome coronavirus (SARS-CoV) epidemic that occurred in Hong Kong in 2003. The dynamic model developed with system dynamics methodology included 23 variables (five states, four flows, eight auxiliary variables, six parameters), five differential equations and 12 algebraic equations. The parameters were optimized following an iterative process of simulation to fit the real data from the epidemics. Univariate and multivariate sensitivity analyses were performed to determine the reliability of the model. In addition, we discuss how further testing using this model can inform community interventions to reduce the risk in current and future outbreaks, such as the recently Middle East respiratory syndrome coronavirus (MERS-CoV) epidemic.</p>
</abstract>
<kwd-group><title>Key words</title>
<kwd>Emerging infections</kwd>
<kwd>epidemics</kwd>
<kwd>MERS-CoV</kwd>
<kwd>modelling</kwd>
<kwd>SARS-CoV</kwd>
</kwd-group>
<counts><fig-count count="7"></fig-count>
<table-count count="4"></table-count>
<ref-count count="44"></ref-count>
<page-count count="12"></page-count>
</counts>
</article-meta>
</front>
<body><sec sec-type="intro" id="sec1"><title>INTRODUCTION</title>
<p>Middle East respiratory syndrome (MERS), is a respiratory illness caused by a novel coronavirus (CoV) [<xref rid="ref1" ref-type="bibr">1</xref>
]. The disease was reported for the first time in Saudi Arabia in June 2012 and spread to several countries in Africa, Asia, Americas and Europe [<xref rid="ref2" ref-type="bibr">2</xref>
, <xref rid="ref3" ref-type="bibr">3</xref>
]. The capability of human-to-human transmission has been observed in at least four hospital settings [<xref rid="ref4" ref-type="bibr">4</xref>
–<xref rid="ref7" ref-type="bibr">7</xref>
]. Significantly, MERS-CoV shares certain similarities with the severe acute respiratory syndrome (SARS)-CoV that produced a global epidemic with more than 8000 human cases in 2002–2003 [<xref rid="ref8" ref-type="bibr">8</xref>
, <xref rid="ref9" ref-type="bibr">9</xref>
]. First, a number of patients infected with both viruses developed an acute respiratory disease that in some cases resulted in death [<xref rid="ref7" ref-type="bibr">7</xref>
, <xref rid="ref9" ref-type="bibr">9</xref>
]. In this sense, MERS-CoV appears to be highly pathogenic with an estimated case-fatality rate of around 50%, although this might be an overestimation as many infected patients may not have sought hospital assistance [<xref rid="ref2" ref-type="bibr">2</xref>
]. Second, both MERS-CoV and SARS-CoV, belong to the genus <italic>Betacoronavirus</italic>
and are closely related to coronaviruses isolated from bats [<xref rid="ref10" ref-type="bibr">10</xref>
–<xref rid="ref12" ref-type="bibr">12</xref>
]. This strongly suggests that MERS-CoV and SARS-CoV may have been transmitted from bats to humans through intermediate species (e.g. camels for MERS-CoV and civet cats for SARS-CoV). Third, since they are new emerging viruses, there are no effective vaccines or antiviral treatments. The SARS epidemic is a clear example of how a networked health system can respond to a new threat to human health.</p>
<p>One of the best-characterized outbreaks during the SARS epidemic was in Hong Kong in 2003 where there were 1755 confirmed cases with 299 deaths (WHO; <uri xlink:type="simple" xlink:href="http://www.who.int/csr/sars/country/en/index.html">http://www.who.int/csr/sars/country/en/index.html</uri>
). The outbreak began in mid-February caused by an infected person who travelled from Guangdong to Hong Kong [<xref rid="ref13" ref-type="bibr">13</xref>
]. An important fact in the generation of a model of the outbreak is that the Hong Kong health authorities quickly implemented contagion control procedures [<xref rid="ref14" ref-type="bibr">14</xref>
]. In general, two interventions were introduced to prevent the spread of SARS-CoV. The first was implementation of quarantine measures to isolate healthy people who had been in contact with infected people and therefore potentially in contact with the virus; isolating those that could be infected and asymptomatic during the incubation period and isolating and treating patients who had developed the disease. The other intervention was the application of protective measures by healthy people who were in contact with infected people to avoid becoming infected, such as respiratory protection for healthcare workers and daily disinfection of the environment of affected rooms [<xref rid="ref14" ref-type="bibr">14</xref>
]. These control interventions were implemented progressively in the Hong Kong Special Administrative Region from mid-March to late April [<xref rid="ref15" ref-type="bibr">15</xref>
]. The application of these procedures allowed the rapid control of the outbreak in the subsequent months. In this sense, it is estimated that the epidemic in Hong Kong ended in late June.</p>
<p>System dynamics has been proved to be a powerful instrument for analysing social, economic, ecological and biological systems [<xref rid="ref16" ref-type="bibr">16</xref>
]. In addition, disease epidemiology has been studied using this approach, whereby system dynamics offers the practical application of concepts by computerized models that allow the systematic test of different scenarios and alternative policies [<xref rid="ref17" ref-type="bibr">17</xref>
–<xref rid="ref19" ref-type="bibr">19</xref>
]. In this work, we performed a modelling of the Hong Kong SARS-CoV outbreak using system dynamics. The developed model contains five states, four flows, eight auxiliary variables and six parameters that interact through five differential and 12 algebraic equations. The parameters of the model were optimized following an iterative process of simulation to obtain a model that largely fits the data available to the epidemic. Moreover, the credibility of our model and its parameters are supported by both univariate and multivariate sensitivity analyses. The model reproduces how the implementation of control measures was effective in preventing the spread of infection to the rest of the population. Basically, these measures result in a sustained reduction in the frequency of contacts. At present, the application of similar measures for infection containment can help to prevent the spread of new emerging epidemics, such as the outbreak caused by MERS-CoV.</p>
</sec>
<sec sec-type="methods" id="sec2"><title>METHODS</title>
<sec id="sec2-1"><title>Data sources and collection</title>
<p>Data on the total population of Hong Kong in 2003 was obtained from statistics of the Census and Statistic Department, The Government of the Hong Kong Special Administrative Region (<uri xlink:type="simple" xlink:href="http://www.censtatd.gov.hk/home/">http://www.censtatd.gov.hk/home/</uri>
). Data on cumulative cases, deaths and recoveries during the SARS epidemics in Hong Kong was obtained from Global Alert and Response databases of the WHO (<uri xlink:type="simple" xlink:href="http://www.who.int/csr/sars/country/en/index.html">http://www.who.int/csr/sars/country/en/index.html</uri>
).</p>
</sec>
<sec id="sec2-2"><title>Modelling</title>
<p>The model was developed following the four-step sequence proposed by system dynamics methodology [<xref rid="ref16" ref-type="bibr">16</xref>
]. First, the real data from the Hong Kong SARS epidemics (<xref ref-type="fig" rid="fig01">Fig. 1</xref>
) together with other evidences and our professional experience were used to create a mental modelling of the reality of the outbreak. Second, the model structure that is able to explain the evolution of the epidemics was represented as a Forrester diagram (<xref ref-type="fig" rid="fig02">Fig. 2</xref>
). Third, the outbreak was mathematically modelled as a continuous dynamic process represented by a set of differential and algebraic equations (<xref rid="tab01 tab02 tab03" ref-type="table">Tables 1–3</xref>
). Finally, the model was optimized to fit the real data from <xref ref-type="fig" rid="fig01">Figure 1</xref>
.
<fig id="fig01" orientation="portrait" position="float"><label>Fig. 1.</label>
<caption><p>Real data from the Hong Kong SARS outbreak. (<italic>a</italic>
) Daily reported new cases, deaths, and recovered. (<italic>b</italic>
) Cumulative cases, dead, and recovered.</p>
</caption>
<graphic xlink:href="S095026881400377X_fig1"></graphic>
</fig>
<fig id="fig02" orientation="portrait" position="float"><label>Fig. 2.</label>
<caption><p>Graphical description of the model using Vensim software. Stock variables are represented inside boxes; flow variables are designed as arrows with a valve; the rest of the variables are auxiliary (inside circles) or parameters (in bold); single lines with arrow are used as connectors to specify that the destination variable is affected by the variable of origin.</p>
</caption>
<graphic xlink:href="S095026881400377X_fig2"></graphic>
</fig>
<table-wrap id="tab01" orientation="portrait" position="float"><label>Table 1.</label>
<caption><p>Differential equations for the five model stocks</p>
</caption>
<alternatives><graphic xlink:href="S095026881400377X_tab1"></graphic>
<table frame="hsides" rules="groups"><col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<thead><tr><th align="left" rowspan="1" colspan="1">Eqn no.</th>
<th align="left" rowspan="1" colspan="1">Stock</th>
<th align="left" rowspan="1" colspan="1">Equation</th>
<th align="left" rowspan="1" colspan="1">Initial value</th>
<th align="left" rowspan="1" colspan="1">Units</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">(1)</td>
<td rowspan="1" colspan="1"><sc>susceptible</sc>
</td>
<td rowspan="1" colspan="1">∂ <sc>susceptible</sc>
(<italic>t</italic>
)/∂<italic>t</italic>
= −Incidence(<italic>t</italic>
)</td>
<td rowspan="1" colspan="1">6·5 × 10<sup>6</sup>
<xref ref-type="table-fn" rid="tfn1_1">*</xref>
</td>
<td rowspan="1" colspan="1">People</td>
</tr>
<tr><td rowspan="1" colspan="1">(2)</td>
<td rowspan="1" colspan="1"><sc>latent</sc>
</td>
<td rowspan="1" colspan="1">∂ <sc>latent</sc>
(<italic>t</italic>
))/∂<italic>t</italic>
= Incidence(<italic>t</italic>
) − Sick per day(<italic>t</italic>
)</td>
<td rowspan="1" colspan="1">0</td>
<td rowspan="1" colspan="1">People</td>
</tr>
<tr><td rowspan="1" colspan="1">(3)</td>
<td rowspan="1" colspan="1"><sc>infected</sc>
</td>
<td rowspan="1" colspan="1">∂ <sc>infected</sc>
(<italic>t</italic>
)/∂<italic>t</italic>
= Sick per day(<italic>t</italic>
) − Daily recovered(<italic>t</italic>
) − Daily deaths(<italic>t</italic>
)</td>
<td rowspan="1" colspan="1">1<xref ref-type="table-fn" rid="tfn1_2">†</xref>
</td>
<td rowspan="1" colspan="1">People</td>
</tr>
<tr><td rowspan="1" colspan="1">(4)</td>
<td rowspan="1" colspan="1"><sc>recovered</sc>
</td>
<td rowspan="1" colspan="1">∂ <sc>recovered</sc>
(<italic>t</italic>
)/∂<italic>t</italic>
= Daily recovered(<italic>t</italic>
)</td>
<td rowspan="1" colspan="1">0</td>
<td rowspan="1" colspan="1">People</td>
</tr>
<tr><td rowspan="1" colspan="1">(5)</td>
<td rowspan="1" colspan="1"><sc>dead</sc>
</td>
<td rowspan="1" colspan="1">∂ <sc>dead</sc>
(<italic>t</italic>
)/∂<italic>t</italic>
= Daily deaths(<italic>t</italic>
)</td>
<td rowspan="1" colspan="1">0</td>
<td rowspan="1" colspan="1">People</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="tfn1_1"><label>*</label>
<p>Hong Kong's population in 2003 was about 6·5 million people who were susceptible to infection with the SARS virus.</p>
</fn>
<fn id="tfn1_2"><label>†</label>
<p>It is considered that the Hong Kong outbreak originated from an infected person from the province of Guangdong [<xref rid="ref13" ref-type="bibr">13</xref>
].</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="tab02" orientation="portrait" position="float"><label>Table 2.</label>
<caption><p>Algebraic equations for the four model flows</p>
</caption>
<alternatives><graphic xlink:href="S095026881400377X_tab2"></graphic>
<table frame="hsides" rules="groups"><col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<thead><tr><th align="left" rowspan="1" colspan="1">Eqn no.</th>
<th align="left" rowspan="1" colspan="1">Flow</th>
<th align="left" rowspan="1" colspan="1">Equation</th>
<th align="left" rowspan="1" colspan="1">Units</th>
<th rowspan="1" colspan="1"></th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">(1)</td>
<td rowspan="1" colspan="1">Incidence</td>
<td rowspan="1" colspan="1">Incidence(<italic>t</italic>
) = <sc>susceptible</sc>
(<italic>t</italic>
) × Prevalence(<italic>t</italic>
) × Contagion rate(<italic>t</italic>
)</td>
<td colspan="2" align="left" rowspan="1">People day<sup>−1</sup>
</td>
</tr>
<tr><td rowspan="1" colspan="1">(2)</td>
<td rowspan="1" colspan="1">Sick per day</td>
<td rowspan="1" colspan="1">Sick per day(<italic>t</italic>
) = <sc>latent</sc>
(<italic>t</italic>
)/Incubation period</td>
<td colspan="2" align="left" rowspan="1">People day<sup>−1</sup>
</td>
</tr>
<tr><td rowspan="1" colspan="1">(3)</td>
<td rowspan="1" colspan="1">Daily recovered</td>
<td rowspan="1" colspan="1">Daily recovered(<italic>t</italic>
) = <sc>infected</sc>
(<italic>t</italic>
) × Recovery rate(<italic>t</italic>
)</td>
<td colspan="2" align="left" rowspan="1">People day<sup>−1</sup>
</td>
</tr>
<tr><td rowspan="1" colspan="1">(4)</td>
<td rowspan="1" colspan="1">Daily deaths</td>
<td rowspan="1" colspan="1">Daily deaths(<italic>t</italic>
) = <sc>infected</sc>
(<italic>t</italic>
) × Case fatality/Disease duration</td>
<td colspan="2" align="left" rowspan="1">People day<sup>−1</sup>
</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<table-wrap id="tab03" orientation="portrait" position="float"><label>Table 3.</label>
<caption><p>Algebraic equations for the eight model auxiliary variables</p>
</caption>
<alternatives><graphic xlink:href="S095026881400377X_tab3"></graphic>
<table frame="hsides" rules="groups"><col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<thead><tr><th align="left" rowspan="1" colspan="1">Eqn no.</th>
<th colspan="1" align="left" rowspan="1">Auxiliary variable</th>
<th colspan="1" align="left" rowspan="1">Equation</th>
<th colspan="1" align="left" rowspan="1">Units</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">(1)</td>
<td rowspan="1" colspan="1">Prevalence</td>
<td colspan="1" align="left" rowspan="1">Prevalence(<italic>t</italic>
) = <sc>infected</sc>
(<italic>t</italic>
)/Population(<italic>t</italic>
)</td>
<td colspan="1" align="left" rowspan="1">Dimensionless</td>
</tr>
<tr><td rowspan="1" colspan="1">(2)</td>
<td rowspan="1" colspan="1">Recovery rate</td>
<td colspan="1" align="left" rowspan="1">Recovery rate(<italic>t</italic>
) = (1 − Case fatality)/(Disease duration)</td>
<td colspan="1" align="left" rowspan="1">Day<sup>−1</sup>
</td>
</tr>
<tr><td rowspan="1" colspan="1">(3)</td>
<td rowspan="1" colspan="1">Contagion rate</td>
<td colspan="1" align="left" rowspan="1">Contagion rate(<italic>t</italic>
) = Frequency of contacts(t) × Infectivity</td>
<td colspan="1" align="left" rowspan="1">Day<sup>−1</sup>
</td>
</tr>
<tr><td rowspan="1" colspan="1">(4)</td>
<td rowspan="1" colspan="1">Population</td>
<td colspan="1" align="left" rowspan="1">Population(<italic>t</italic>
) = <sc>infected</sc>
(<italic>t</italic>
) + <sc>latent</sc>
(<italic>t</italic>
) + <sc>recovered</sc>
(<italic>t</italic>
) + <sc>susceptible</sc>
(<italic>t</italic>
)</td>
<td colspan="1" align="left" rowspan="1">People</td>
</tr>
<tr><td rowspan="1" colspan="1">(5)</td>
<td rowspan="1" colspan="1">Cumulative cases</td>
<td colspan="1" align="left" rowspan="1">Cumulative cases(<italic>t)</italic>
= <sc>dead</sc>
(<italic>t</italic>
) + <sc>recovered</sc>
(<italic>t</italic>
) + <sc>infected</sc>
(<italic>t</italic>
)</td>
<td colspan="1" align="left" rowspan="1">People</td>
</tr>
<tr><td rowspan="1" colspan="1">(6)</td>
<td rowspan="1" colspan="1">Cumulative attack rate</td>
<td colspan="1" align="left" rowspan="1">Cumulative attack rate (<italic>t</italic>
) = Cumulative cases(<italic>t</italic>
)/Population(<italic>t</italic>
)</td>
<td colspan="1" align="left" rowspan="1">Dimensionless</td>
</tr>
<tr><td rowspan="1" colspan="1">(7)</td>
<td rowspan="1" colspan="1">Frequency of contacts</td>
<td colspan="1" align="left" rowspan="1">Frequency of contacts(<italic>t</italic>
) = Daily contacts /(1 + (Cumulative attack rate/Threshold of cumulative attack rate)<sup>3</sup>
)</td>
<td colspan="1" align="left" rowspan="1">Day<sup>−1</sup>
</td>
</tr>
<tr><td rowspan="1" colspan="1">(8)</td>
<td rowspan="1" colspan="1">Basic reproductive number (<italic>R</italic>
<sub>0</sub>
)</td>
<td colspan="1" align="left" rowspan="1">Basic reproductive number(<italic>t</italic>
) = Contagion rate(<italic>t</italic>
)/Recovery rate(<italic>t</italic>
)</td>
<td colspan="1" align="left" rowspan="1">Dimensionless</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</p>
<p>A dynamic compartmental model provides a framework for the study of transport between different compartments of a system, i.e. well known epidemiological compartmental models [<xref rid="ref20" ref-type="bibr">20</xref>
]. To explain how the protective measures taken by the government of Hong Kong allowed the rapid control of the epidemic we consider a dynamic model with five compartments (states) and four transitions (physical flows) between them. This model is based on two assumptions. First, the individuals are classified into five subgroups (<sc>susceptible, latent, infected, recovered, dead</sc>
). Although, the last subgroup is not strictly needed, it is used to keep an account of those dead. Second, every day there is a different number of people: who are infected without symptoms of the disease (incidence); who develop signs and symptoms of the illness (sick per day); who recover from the disease (daily recovered), and who die (daily deaths) as consequence of the outbreak.</p>
<p>The model structure built with Vensim software (Ventana Systems Inc., USA) (<xref ref-type="fig" rid="fig02">Fig. 2</xref>
) contains the five states and the four flows mentioned above and also eight auxiliary variables (inside circles) and six parameters (in bold). These elements are linked by the physical flows (double line with arrow) and by the information transmissions (single line with arrow) according to the mathematical model represented by the set of differential and algebraic equations (<xref rid="tab01 tab02 tab03" ref-type="table">Tables 1–3</xref>
).</p>
<p>The five differential equations in <xref rid="tab01" ref-type="table">Table 1</xref>
establish the mass balance (inflows minus outflows) in the compartments, and as such, they describe the changes in the number of people (stocks) in the five subgroups. The four algebraic equations of <xref rid="tab02" ref-type="table">Table 2</xref>
express that the physical transitions depend directly from the stock in the compartment of origin and indirectly from other stocks and parameters through the corresponding auxiliary variables. These equations involve three stocks (<sc>susceptible, latent, infected</sc>
), three auxiliary variables (prevalence, contagion rate, recovery rate) and three parameters (incubation period, case fatality, disease duration). For instance, the incidence flow is proportional to the number of susceptible people, the prevalence and the contagion rate. The contagion rate or transmission coefficient theoretically depends on the number of contacts per unit time and the probability of effective contact, i.e. the probability that a contact between an infectious source and susceptible host results in a successful transfer whereby the susceptible host becomes infected.</p>
<p>The first three algebraic equations of <xref rid="tab03" ref-type="table">Table 3</xref>
express the auxiliary variables (prevalence, contagion rate, recovery rate) used in the equations of <xref rid="tab02" ref-type="table">Table 2</xref>
. In turn, these depend on other auxiliary variables, stocks and parameters. For instance, prevalence is defined as the ratio between number of infected people and the total population [<xref rid="tab03" ref-type="table">Table 3</xref>
, equation (1)]. The recovery rate depends on the illness duration and case fatality [<xref rid="tab03" ref-type="table">Table 3</xref>
, equation (2)]. Equation (3) in <xref rid="tab03" ref-type="table">Table 3</xref>
expresses the contagion rate which is directly dependent on the auxiliary variable ‘frequency of contacts’ and the parameter ‘infectivity’.</p>
<p>One of the key variables in the model is frequency of contacts because it tries to reproduce the control measures carried out by the Hong Kong Government to control the SARS outbreak. We assumed that: (1) the control measures gave, as consequence, a marked reduction of daily contacts in Hong Kong; (2) the control measures were based on the cumulative attack rate, measured as the ratio between cumulative cases and total population [<xref rid="tab03" ref-type="table">Table 3</xref>
, equation (6)]. Therefore, we decided to use the frequency of contacts as the daily contacts modulated by the repression Hill function (<xref ref-type="fig" rid="fig03">Fig. 3</xref>
), in accordance with the equation (7) of <xref rid="tab03" ref-type="table">Table 3</xref>
. Despite the complexity level of biological systems several cases have been modelled using the Hill function, in order to simulate repressor activities of enzymatic reactions and the regulation mechanisms of several transcription factors [<xref rid="ref21" ref-type="bibr">21</xref>
]. In our model this repressor function allows the establishment of the relationship between frequency of contacts and the cumulative attack rate. Note that the cumulative attack rate is used as an active repressor, so half-maximal repression occurs when the cumulative attack rate is equal to the threshold, and almost total repression occurs when this cumulative attack rate is double the threshold.
<fig id="fig03" orientation="portrait" position="float"><label>Fig. 3.</label>
<caption><p>Normalized representation of repression Hill function. The solid line corresponds to a value of <italic>n</italic>
= 3 [as used in equation (7) of <xref rid="tab03" ref-type="table">Table 3</xref>
]. The dotted line corresponds to the threshold repression function (when <italic>n</italic>
tends to infinity). The representation is normalized to the threshold of cumulative attack rate on the abscissa axis and the value of the discontinuity (daily contacts) on the ordinate axis.</p>
</caption>
<graphic xlink:href="S095026881400377X_fig3"></graphic>
</fig>
</p>
<p>Equation (8) in <xref rid="tab03" ref-type="table">Table 3</xref>
is used to report on the basic reproductive number (<italic>R</italic>
<sub>0</sub>
), which is defined as the expected number of secondary infectious cases generated by an average infectious case in an entirely susceptible population [<xref rid="ref20" ref-type="bibr">20</xref>
]. In our model this number is calculated as the ratio between the contagion rate and the recovery rate. If <italic>R</italic>
<sub>0</sub>
< 1, then the infected individuals in the total population fail to replace themselves, and the disease does not spread. However, if <italic>R</italic>
<sub>0</sub>
> 1, the number of cases generally increases over time and the disease spreads.</p>
</sec>
<sec id="sec2-3"><title>Parameterization</title>
<p>One of the most challenging issues in system dynamics modelling is to establish the value of the model parameters. The parameters can be estimated taking advantage of the known information available in the literature and can be optimized by means of an iterative process of simulation. Using this approach, we set values for the six parameters of our model which are summarized in <xref rid="tab04" ref-type="table">Table 4</xref>
.
<table-wrap id="tab04" orientation="portrait" position="float"><label>Table 4.</label>
<caption><p>Values for the six model parameters</p>
</caption>
<alternatives><graphic xlink:href="S095026881400377X_tab4"></graphic>
<table frame="hsides" rules="groups"><col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<thead><tr><th align="left" rowspan="1" colspan="1">Parameter</th>
<th align="left" rowspan="1" colspan="1">Value</th>
<th align="left" rowspan="1" colspan="1">Units</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">Infectivity</td>
<td rowspan="1" colspan="1">0·022</td>
<td rowspan="1" colspan="1">Dimensionless</td>
</tr>
<tr><td rowspan="1" colspan="1">Incubation period</td>
<td rowspan="1" colspan="1">5·3</td>
<td rowspan="1" colspan="1">Day</td>
</tr>
<tr><td rowspan="1" colspan="1">Disease duration</td>
<td rowspan="1" colspan="1">26</td>
<td rowspan="1" colspan="1">Day</td>
</tr>
<tr><td rowspan="1" colspan="1">Case fatality</td>
<td rowspan="1" colspan="1">0·17</td>
<td rowspan="1" colspan="1">Dimensionless</td>
</tr>
<tr><td rowspan="1" colspan="1">Daily contacts</td>
<td rowspan="1" colspan="1">16·8</td>
<td rowspan="1" colspan="1">Day<sup>−1</sup>
</td>
</tr>
<tr><td rowspan="1" colspan="1">Threshold of cumulative attack rate</td>
<td rowspan="1" colspan="1">7·8 × 10<sup>−5</sup>
</td>
<td rowspan="1" colspan="1">Dimensionless</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</p>
<p>Infectivity expresses the ability of the pathogen to penetrate, survive and multiply in the host and it is measured through the secondary attack rate which is defined as the probability that infection occurs in susceptible persons within a reasonable incubation period following known contact with an infectious person or an infectious source. Epidemiological studies in Singapore showed that 80% of individuals infected with SARS-CoV did not cause secondary infections, suggesting low infectivity [<xref rid="ref22" ref-type="bibr">22</xref>
]. After the optimization process, the final value set for infectivity was 0·022, which is in agreement with a previous work examining the probability of transmission for SARS-CoV [<xref rid="ref23" ref-type="bibr">23</xref>
].</p>
<p>The incubation period of the disease, during which time the individual is asymptomatic appears to be between 2 and 10 days with an average value of 5·3 days [<xref rid="ref24" ref-type="bibr">24</xref>
]. Similarly, although the duration of the disease is variable depending on the case, the average of the most severe and the mildest cases is 26 days [<xref rid="ref24" ref-type="bibr">24</xref>
]. During this time period the virus can be transmitted to other people.</p>
<p>Several studies have shown that the case fatality of SARS-CoV was also variable in the 2003 outbreak depending on different factors. It has been estimated that the mean case fatality for SARS in Hong Kong was around 0·17 [<xref rid="ref25" ref-type="bibr">25</xref>
].</p>
<p>Social contact patterns have been shown to be highly relevant on the transmission dynamics of respiratory infections such as measles, rubella and influenza [<xref rid="ref26" ref-type="bibr">26</xref>
–<xref rid="ref28" ref-type="bibr">28</xref>
]. Quantifying this parameter is critical for estimating the impact of such infections, for designing and targeting preventive interventions, and for modelling their impact. In our model the daily contacts parameter, indicating the mean number of contacts of a person in 1 day, was set to 16·8 after the optimization process, which is in agreement with previous studies quantifying these social mixing patterns [<xref rid="ref29" ref-type="bibr">29</xref>
, <xref rid="ref30" ref-type="bibr">30</xref>
].</p>
<p>The parameter ‘threshold of cumulative attack rate’ is critical in our model since it allows setting the value of the cumulative attack rate at which control measures were established by the Hong Kong authorities. This parameter was finally set at 7·8 × 10<sup>−5</sup>
. Therefore, according to the Hill function of <xref ref-type="fig" rid="fig03">Figure 3</xref>
, for any number of the cumulative attack rate below the 78 cases by million, the frequency of contacts will be greater than 8·4 (half of the daily contacts), and in the opposite case the frequency of contacts will be markedly reduced.</p>
<p>As explained before, our dynamic model was subjected to successive rounds of simulation and optimization in order to fine-tune the parameters of the model. In all these simulations we assumed that the Hong Kong epidemic originated from a traveller from Guangdong in China [<xref rid="ref13" ref-type="bibr">13</xref>
], and at the beginning of the outbreak the entire population of Hong Kong was susceptible to SARS-CoV infection. These assumptions are represented by the initial values in the five stocks (<xref rid="tab01" ref-type="table">Table 1</xref>
). Moreover, based on the reported data, the simulation period was set to 146 days with a time step of 1 day.</p>
</sec>
<sec id="sec2-4"><title>Sensitivity analysis</title>
<p>Parameters of system dynamics models are subject to uncertainty. Sensitivity analyses were conducted to provide insight into how uncertainty in the parameters affects the model outputs and which parameters tend to drive these variations. In this task it is essential to define the probabilistic distribution patterns of the model parameters, which are shown in <xref rid="tab05" ref-type="table">Table 5</xref>
, together with the references that support these patterns.
<table-wrap id="tab05" orientation="portrait" position="float"><label>Table 5.</label>
<caption><p>Probabilistic distribution patterns of the model parameters</p>
</caption>
<alternatives><graphic xlink:href="S095026881400377X_tab5"></graphic>
<table frame="hsides" rules="groups"><col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<col align="left" width="1*" span="1"></col>
<thead><tr><th align="left" rowspan="1" colspan="1">Parameter</th>
<th align="left" rowspan="1" colspan="1">Distribution</th>
<th align="left" rowspan="1" colspan="1">Range</th>
<th align="left" rowspan="1" colspan="1">Distr. values<xref ref-type="table-fn" rid="tfn5_1">*</xref>
</th>
<th align="left" rowspan="1" colspan="1">Reference</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">Infectivity</td>
<td rowspan="1" colspan="1">Gamma</td>
<td rowspan="1" colspan="1">(0·001–0·2)</td>
<td rowspan="1" colspan="1">1·1/0/0·025</td>
<td rowspan="1" colspan="1">[<xref rid="ref23" ref-type="bibr">23</xref>
]</td>
</tr>
<tr><td rowspan="1" colspan="1">Incubation period</td>
<td rowspan="1" colspan="1">Gamma</td>
<td rowspan="1" colspan="1">(1–30)</td>
<td rowspan="1" colspan="1">1·4/1/3·85</td>
<td rowspan="1" colspan="1">[<xref rid="ref24" ref-type="bibr">24</xref>
]</td>
</tr>
<tr><td rowspan="1" colspan="1">Disease duration</td>
<td rowspan="1" colspan="1">Gamma</td>
<td rowspan="1" colspan="1">(1–70)</td>
<td rowspan="1" colspan="1">5·6/0/4·55</td>
<td rowspan="1" colspan="1">[<xref rid="ref24" ref-type="bibr">24</xref>
]</td>
</tr>
<tr><td rowspan="1" colspan="1">Case fatality</td>
<td rowspan="1" colspan="1">Normal</td>
<td rowspan="1" colspan="1">(0·14–0·21)</td>
<td rowspan="1" colspan="1">0·171/0·01</td>
<td rowspan="1" colspan="1">[<xref rid="ref25" ref-type="bibr">25</xref>
]</td>
</tr>
<tr><td rowspan="1" colspan="1">Daily contacts</td>
<td rowspan="1" colspan="1">Normal</td>
<td rowspan="1" colspan="1">(1–40)</td>
<td rowspan="1" colspan="1">16·8/8·5</td>
<td rowspan="1" colspan="1">[<xref rid="ref29" ref-type="bibr">29</xref>
]</td>
</tr>
<tr><td rowspan="1" colspan="1">Threshold of cumulative attack rate</td>
<td rowspan="1" colspan="1">Normal</td>
<td rowspan="1" colspan="1">(5·0 × 10<sup>−5</sup>
–1·1 × 10<sup>−4</sup>
)</td>
<td rowspan="1" colspan="1">7·8 × 10<sup>−5</sup>
/1 × 10<sup>−5</sup>
</td>
<td rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="tfn5_1"><label>*</label>
<p>The distribution values for normal distributions are mean and standard deviation, whereas for gamma distributions they are order/shift/stretch.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</p>
<p>The most influential parameters were estimated by univariate analyses, in which changes in the model output were studied after disturbance in each parameter value independently. In complex models, univariate sensitivity analysis can be insufficient for a comprehensive study of the model. Simultaneous fluctuations in the value of more than one parameter may create an unexpected output change due to nonlinear relationships in different model components. Thus, to test the influence of simultaneous changes in the model parameters, the univariate analyses were followed by Monte-Carlo multivariate sensitivity analysis, in which the values of the six parameters were altered at the same time.</p>
</sec>
<sec id="sec2-5"><title>Software</title>
<p>Modelling process, simulations and sensitivity analyses were performed using Vensim DSS software v. 5.7a (Ventana Systems).</p>
</sec>
</sec>
<sec sec-type="results" id="sec3"><title>RESULTS</title>
<p><xref ref-type="fig" rid="fig04">Figure 4</xref>
shows a graphic comparison between the simulation results using the parameters of <xref rid="tab04" ref-type="table">Table 4</xref>
and the real data of the Hong Kong outbreak. We compared only the variables of the model for which records were found in the public databases, which are the same six variables shown in <xref ref-type="fig" rid="fig01">Figure 1</xref>
.
<fig id="fig04" orientation="portrait" position="float"><label>Fig. 4.</label>
<caption><p>Simulation of the model. For simulation, the time period was set to 146 days with a time step of 1 day. The graphs represent the output of the simulation (solid lines) and the data reported for the SARS outbreak (dotted lines). (<italic>a</italic>
) Sick per day, (<italic>b</italic>
) cumulative cases, (<italic>c</italic>
) recovered, (<italic>d</italic>
) dead, (<italic>e</italic>
) daily recovered, and (<italic>f</italic>
) daily deaths.</p>
</caption>
<graphic xlink:href="S095026881400377X_fig4"></graphic>
</fig>
</p>
<p>The simulation output for the variable ‘sick per day’ fit the data reported by the Hong Kong authorities (<xref ref-type="fig" rid="fig04">Fig. 4<italic>a</italic>
</xref>
), suggesting that the model was able to reproduce the epidemic curve. We observed that the number of new cases per day obtained in the simulation grew during the first 46 days. From this time, the number of new infections gradually fell to values <1 at later stages of the outbreak.</p>
<p>As a consequence, the auxiliary variable that stores the cumulative SARS cases showed a characteristic sigmoidal growth (<xref ref-type="fig" rid="fig04">Fig. 4<italic>b</italic>
</xref>
), consistent with the real data. We observed that the number of SARS cases grew from one at the beginning of the epidemic to around 1800 at the final stage of the outbreak, similar to the 1755 cases reported by the authorities. Moreover, we observed a high fit between the model predictions and the real data.</p>
<p>As expected, the number of deaths and recovered people in the simulation also grew with a sigmoidal shape to reach values similar to those found in the public databases (<xref ref-type="fig" rid="fig04">Fig. 4<italic>c</italic>
, <italic>d</italic>
</xref>
). However, we observed a partial fit of these stock variables to the data reported during the outbreak. These slight mismatches were also observed in the flow variables for recovered and daily deaths (<xref ref-type="fig" rid="fig04">Fig. 4<italic>e</italic>
, <italic>f</italic>
</xref>
) that may be due to delays in reporting of cases by the authorities.</p>
<p>In short, looking the simulation results of <xref ref-type="fig" rid="fig03">Figure 3</xref>
we can conclude that our model is able to reproduce largely the most important indicators of the SARS epidemic that occurred in Hong Kong in 2003.</p>
<p>Focusing on the evolution of the basic reproduction number (<xref ref-type="fig" rid="fig05">Fig. 5</xref>
), we note that during the early stages of the epidemic, <italic>R</italic>
<sub>0</sub>
>1, which is consistent with the observed disease spread. Moreover, after day 30, <italic>R</italic>
<sub>0</sub>
starts to drop to values <1, probably due to the implementation of containment measures by the Hong Kong authorities after the issuance of the first global alert against SARS on 12 March 2003. These results are consistent with a previous report showing the basic reproductive numbers for different SARS epidemic curves, which supports the notion that our model is able to largely replicate the disease outbreak in Hong Kong [<xref rid="ref31" ref-type="bibr">31</xref>
].
<fig id="fig05" orientation="portrait" position="float"><label>Fig. 5.</label>
<caption><p>Basic reproductive number obtained from the simulation. <italic>R</italic>
<sub>0</sub>
was calculated as the ratio between the contagion rate and the recovery rate.</p>
</caption>
<graphic xlink:href="S095026881400377X_fig5"></graphic>
</fig>
</p>
<p>The results of the univariate sensitivity analysis are shown in <xref ref-type="fig" rid="fig06">Figure 6</xref>
. We focused our attention on the epidemic wave although the analysis is possible in other variables, as shown in the multivariate analysis of <xref ref-type="fig" rid="fig07">Figure 7</xref>
. Variations in case fatality, threshold of cumulative attack rate and disease duration induced little changes in the epidemic curve (<xref ref-type="fig" rid="fig06">Fig. 6–<italic>c</italic>
</xref>
), while more extensive alterations in the epidemic wave were observed after changes in infectivity, daily contacts and incubation period (<xref ref-type="fig" rid="fig06">Fig. 6<italic>d–f</italic>
</xref>
). Variations in the case-fatality parameter do not alter the output of the variable sick per day (<xref ref-type="fig" rid="fig06">Fig. 6<italic>a</italic>
</xref>
), although other variables from the model such as recovered and dead are highly impacted (data not shown). Small changes in the shape and the maximum of the curve are observed after modification of the parameters ‘threshold of cumulative attack rate’ and ‘disease duration’ (<xref ref-type="fig" rid="fig06">Fig. 6<italic>b</italic>
, <italic>c</italic>
</xref>
). By contrast, changes in infectivity, daily contacts and to a lesser extent in incubation period significantly alter both the position and height of the maximum of the epidemic curve (<xref ref-type="fig" rid="fig06">Fig. 6<italic>d</italic>
–<italic>f</italic>
</xref>
).
<fig id="fig06" orientation="portrait" position="float"><label>Fig. 6.</label>
<caption><p>Sensitivity of the model to changes in parameters. Univariate analysis with 1000 simulations was performed for each parameter: (<italic>a</italic>
) case fatality, (<italic>b</italic>
) threshold of cumulative attack rate, (<italic>c</italic>
) disease duration, (<italic>d</italic>
) infectivity, (<italic>e</italic>
) daily contacts, and (<italic>f</italic>
) incubation period. The parameters were varied considering ranges and distributions shown in <xref rid="tab05" ref-type="table">Table 5</xref>
. The solid black line represents the simulation output and the grey area represents the 95% confidence bounds.</p>
</caption>
<graphic xlink:href="S095026881400377X_fig6"></graphic>
</fig>
<fig id="fig07" orientation="portrait" position="float"><label>Fig. 7.</label>
<caption><p>Multivariate sensitivity analysis of the model. A Monte-Carlo sensitivity analysis with 1000 simulations was performed and the results for (<italic>a</italic>
) sick per day, (<italic>b</italic>
) infected, (<italic>c</italic>
) recovered, and (<italic>d</italic>
) dead are shown. The values of the six parameters were altered at the same time following the distributions described in <xref rid="tab05" ref-type="table">Table 5</xref>
. The solid black line represents the simulation output and the grey area represents the 95% confidence bounds.</p>
</caption>
<graphic xlink:href="S095026881400377X_fig7"></graphic>
</fig>
</p>
<p>The results of the multivariate sensitivity analysis are shown in <xref ref-type="fig" rid="fig07">Figure 7</xref>
, we analysed the output of four variables: sick per day, infected, recovered, and dead (<xref ref-type="fig" rid="fig07">Fig. 7<italic>a</italic>
–<italic>d</italic>
</xref>
, respectively). Variations in the model parameters clearly change the shape of the epidemic curve, altering both the position and the height of the maximum of this variable (<xref ref-type="fig" rid="fig07">Fig. 7<italic>a</italic>
</xref>
). Similarly, the output of the variables infected and recovered is highly impacted by the changes in parameters carried out in the multivariate sensitivity analysis (<xref ref-type="fig" rid="fig07">Fig. 7<italic>b</italic>
, <italic>c</italic>
</xref>
). The alterations of variable outputs are clearly exemplified by the variable ‘dead’ (<xref ref-type="fig" rid="fig07">Fig. 7<italic>d</italic>
</xref>
). Certain changes in the parameters of the model can explain an increase in the total number deaths during the epidemic rising from about 300 to about 700. The observed variations in model output when the value of the parameters is changed support the idea that this model might be able to simulate different scenarios and epidemic conditions.</p>
</sec>
<sec sec-type="discussion" id="sec4"><title>DISCUSSION</title>
<p>System dynamics modelling has been successfully applied to study complex public health issues such as the design of optimal policies in healthcare [<xref rid="ref32" ref-type="bibr">32</xref>
], the impact of public health intervention in different situations [<xref rid="ref33" ref-type="bibr">33</xref>
, <xref rid="ref34" ref-type="bibr">34</xref>
], and to study disease epidemiology [<xref rid="ref35" ref-type="bibr">35</xref>
–<xref rid="ref37" ref-type="bibr">37</xref>
]. In the latter, system dynamics technology has become a powerful tool to understand and predict the impact of infectious diseases. Epidemiological models can help health authorities to make recommendations regarding intervention to fight the spread of directly transmissible pathogens, especially when empirical data is limited. In this sense, mathematical models have been previously used to advise health policies against diseases such as pandemic influenza [<xref rid="ref38" ref-type="bibr">38</xref>
, <xref rid="ref39" ref-type="bibr">39</xref>
] and SARS [<xref rid="ref31" ref-type="bibr">31</xref>
, <xref rid="ref40" ref-type="bibr">40</xref>
, <xref rid="ref41" ref-type="bibr">41</xref>
].</p>
<p>Several models studying SARS transmission and interventions have been published. These are detailed hybrid stochastic and compartmental models that successfully explain the behaviour of the epidemic [<xref rid="ref14" ref-type="bibr">14</xref>
, <xref rid="ref31" ref-type="bibr">31</xref>
, <xref rid="ref42" ref-type="bibr">42</xref>
]. Here, trying to follow Ockham's razor principle, we have built a simpler deterministic model, which is also able to reproduce the behaviour of the epidemic based on its natural history and the intervention measures taken in Hong Kong. The use of a less complicated model could be helpful in understanding the disease epidemic and also facilitating its reuse under other conditions.</p>
<p>The epidemiological models depend on the consistency of the chosen parameters. Therefore, the accurate quantification of these parameters is critical to estimate the path of a disease, to predict the impact of possible interventions, and to inform planning and decision making. Here, we have combined reported information from the SARS epidemic with an iterative optimization process to set the final values for the model parameters. Under these conditions, the model output fits the epidemic curve observed in the Hong Kong SARS-CoV outbreak (<xref ref-type="fig" rid="fig04">Fig. 4</xref>
).</p>
<p>Of the factors that influence the dynamics of infectious diseases, the person-to-person contact pattern has been shown to be essential in disease spread [<xref rid="ref43" ref-type="bibr">43</xref>
]. Our model takes this essential factor into account through the auxiliary variable ‘frequency of contacts’, which is dependent on the auxiliary variable ‘cumulative attack rate’ and the parameters ‘daily contacts’ and ‘threshold of cumulative attack rate’ (<xref ref-type="fig" rid="fig02">Fig. 2</xref>
). A previous work showed that mixing patterns and contact characteristics were remarkably similar across the different European countries analysed in that study even though the average number of contacts recorded differed. Interestingly, the authors suggests that the results may well be applicable to other countries with similar social structures, and that the initial epidemic phase of an emerging infection in susceptible populations, such as SARS was in 2003, is likely to be very similar [<xref rid="ref30" ref-type="bibr">30</xref>
]. During the SARS outbreak, health authorities, hospitals, and the overall population progressively implemented quarantine and protection policies to prevent the transmission of the disease [<xref rid="ref15" ref-type="bibr">15</xref>
]. With this in mind, we made the assumption that the intervention of the health authorities caused a decrease in the frequency of contacts, which in turn led a decrease in the rate of contagion. To mimic this event in our model, the number of daily contacts is regulated by means of a repression Hill function, which has been used to simulate repressor activities in complex biological systems [<xref rid="ref21" ref-type="bibr">21</xref>
].</p>
<p>The basic reproductive number is a key epidemiological variable that characterizes the potential of a disease to spread. Several works have estimated that prior to the first global alert the basic reproductive number for SARS was >1, correlating with an exponential increase in the number of cases. However, the implementation of effective control measures, such as quarantine, isolation, and strict hygiene practice in hospitals led a sudden decrease in <italic>R</italic>
<sub>0</sub>
[<xref rid="ref14" ref-type="bibr">14</xref>
, <xref rid="ref31" ref-type="bibr">31</xref>
]. The fact that, in our model, <italic>R</italic>
<sub>0</sub>
drops to values <1 around the date of the global alert reinforces the idea that a decrease in the frequency of contacts is able to effectively simulate the effects of the control measures established in the first stages of the epidemic. It is important to note that estimations of <italic>R</italic>
<sub>0</sub>
in previous publications are considerably lower (around 2–4) than ours, which is almost 12 at the beginning of the outbreak [<xref rid="ref14" ref-type="bibr">14</xref>
, <xref rid="ref31" ref-type="bibr">31</xref>
, <xref rid="ref42" ref-type="bibr">42</xref>
]. Nevertheless, the estimated value of <italic>R</italic>
<sub>0</sub>
also differs in these works and the credible intervals surrounding these deterministic estimations were wide, reaching superior values of almost 8. This high variability can be explained in part by the superspreading events that occurred in SARS epidemics. Superspreading is an unusual situation, in which a single individual directly infects a large number of other people that has a large influence on the early course of the epidemic [<xref rid="ref42" ref-type="bibr">42</xref>
]. Interestingly, the fact that our model does not explicitly account these events could partially explain the very high estimated value of <italic>R</italic>
<sub>0</sub>
at the beginning of the outbreak.</p>
<p>The reliability of the model parameters is supported by univariate and multivariate sensitivity analyses (<xref ref-type="fig" rid="fig06">Figs 6</xref>
and <xref ref-type="fig" rid="fig07">7</xref>
). Furthermore, sensitivity analysis is a powerful tool to analyse the influence of certain decision making in the evolution of the epidemic. Taken together, these findings strongly suggest that our model, together with other system dynamics models can be used by epidemiologists to investigate the likely consequences of future re-emergences of SARS-CoV based on analysis of the previous known epidemics. In addition, by adapting the key parameters of these models or with a little change in the model structures, they can be used to face emerging outbreaks of infectious diseases, such as the recent MERS-CoV epidemic. In this regard, there are several similarities and differences which should be taken into account when using this model. Both SARS-CoV and MERS-CoV may cause severe respiratory failure, extrapulmonary features such as diarrhoea and also mild or asymptomatic cases. In contrast with SARS, MERS has lower human-to-human transmission potential, affects predominantly older people with more comorbid illness and has a higher case-fatality rate [<xref rid="ref44" ref-type="bibr">44</xref>
]. These factors would affect the output of the model variables (e.g. epidemic curve, cumulative cases, number of dead, etc.).</p>
</sec>
</body>
<back><ack id="ack"><title>ACKNOWLEDGEMENTS</title>
<p>This research was supported by an intramural contract from Universidad Autónoma de Madrid awarded to E. Álvarez. The Institutional Grant awarded to the Centro de Biología Molecular ‘Severo Ochoa’ by the Fundación Ramón Areces is also acknowledged.</p>
<sec id="nts1" sec-type="COI-statement"><title>DECLARATION OF INTEREST</title>
<p>None.</p>
</sec>
</ack>
<ref-list id="reflist1"><title>REFERENCES</title>
<ref id="ref1"><label>1.</label>
<mixed-citation publication-type="journal"><name><surname>de Groot</surname>
<given-names>RJ</given-names>
</name>
, <etal></etal>
<article-title>Middle East respiratory syndrome coronavirus (MERS-CoV): announcement of the Coronavirus Study Group</article-title>
. <source>Journal of Virology</source>
<year>2013</year>
; <volume>87</volume>
: <fpage>7790</fpage>
–<lpage>7792</lpage>
.<pub-id pub-id-type="pmid">23678167</pub-id>
</mixed-citation>
</ref>
<ref id="ref2"><label>2.</label>
<mixed-citation publication-type="journal"><collab>CDC</collab>
. <article-title>Update: severe respiratory illness associated with Middle East Respiratory Syndrome Coronavirus (MERS-CoV) – worldwide, 2012–2013</article-title>
. <source>Morbidity and Mortality Weekly Report</source>
<year>2013</year>
; <volume>62</volume>
: <fpage>480</fpage>
–<lpage>483</lpage>
.<pub-id pub-id-type="pmid">23760190</pub-id>
</mixed-citation>
</ref>
<ref id="ref3"><label>3.</label>
<mixed-citation publication-type="journal"><name><surname>Zaki</surname>
<given-names>AM</given-names>
</name>
, <etal></etal>
<article-title>Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia</article-title>
. <source>New England Journal of Medicine</source>
<year>2012</year>
; <volume>367</volume>
: <fpage>1814</fpage>
–<lpage>1820</lpage>
.<pub-id pub-id-type="pmid">23075143</pub-id>
</mixed-citation>
</ref>
<ref id="ref4"><label>4.</label>
<mixed-citation publication-type="journal">Anon. <article-title>Evidence of person-to-person transmission within a family cluster of novel coronavirus infections, United Kingdom, February 2013</article-title>
. <source>Eurosurveillance</source>
<year>2013</year>
; <volume>18</volume>
: <fpage>20427</fpage>
.<pub-id pub-id-type="pmid">23517868</pub-id>
</mixed-citation>
</ref>
<ref id="ref5"><label>5.</label>
<mixed-citation publication-type="journal"><name><surname>Assiri</surname>
<given-names>A</given-names>
</name>
, <etal></etal>
<article-title>Hospital outbreak of Middle East respiratory syndrome coronavirus</article-title>
. <source>New England Journal of Medicine</source>
<year>2013</year>
; <volume>369</volume>
: <fpage>407</fpage>
–<lpage>416</lpage>
.<pub-id pub-id-type="pmid">23782161</pub-id>
</mixed-citation>
</ref>
<ref id="ref6"><label>6.</label>
<mixed-citation publication-type="journal"><name><surname>Mailles</surname>
<given-names>A</given-names>
</name>
, <etal></etal>
<article-title>First cases of Middle East Respiratory Syndrome Coronavirus (MERS-CoV) infections in France, investigations and implications for the prevention of human-to-human transmission, France, May 2013</article-title>
. <source>Eurosurveillance</source>
<year>2013</year>
; <volume>18</volume>
.</mixed-citation>
</ref>
<ref id="ref7"><label>7.</label>
<mixed-citation publication-type="journal"><name><surname>Drosten</surname>
<given-names>C</given-names>
</name>
, <etal></etal>
<article-title>Clinical features and virological analysis of a case of Middle East respiratory syndrome coronavirus infection</article-title>
. <source>Lancet Infectious Diseases</source>
<year>2013</year>
; <volume>13</volume>
: <fpage>745</fpage>
–<lpage>751</lpage>
.<pub-id pub-id-type="pmid">23782859</pub-id>
</mixed-citation>
</ref>
<ref id="ref8"><label>8.</label>
<mixed-citation publication-type="journal"><name><surname>Marra</surname>
<given-names>MA</given-names>
</name>
, <etal></etal>
<article-title>The genome sequence of the SARS-associated coronavirus</article-title>
. <source>Science</source>
<year>2003</year>
; <volume>300</volume>
: <fpage>1399</fpage>
–<lpage>1404</lpage>
.<pub-id pub-id-type="pmid">12730501</pub-id>
</mixed-citation>
</ref>
<ref id="ref9"><label>9.</label>
<mixed-citation publication-type="journal"><name><surname>Peiris</surname>
<given-names>JS</given-names>
</name>
, <etal></etal>
<article-title>Coronavirus as a possible cause of severe acute respiratory syndrome</article-title>
. <source>Lancet</source>
<year>2003</year>
; <volume>361</volume>
: <fpage>1319</fpage>
–<lpage>1325</lpage>
.<pub-id pub-id-type="pmid">12711465</pub-id>
</mixed-citation>
</ref>
<ref id="ref10"><label>10.</label>
<mixed-citation publication-type="journal"><name><surname>Annan</surname>
<given-names>A</given-names>
</name>
, <etal></etal>
<article-title>Human betacoronavirus 2c EMC/2012-related viruses in bats, Ghana and Europe</article-title>
. <source>Emerging Infectious Diseases</source>
<year>2013</year>
; <volume>19</volume>
: <fpage>456</fpage>
–<lpage>459</lpage>
.<pub-id pub-id-type="pmid">23622767</pub-id>
</mixed-citation>
</ref>
<ref id="ref11"><label>11.</label>
<mixed-citation publication-type="journal"><name><surname>Cotten</surname>
<given-names>M</given-names>
</name>
, <etal></etal>
<article-title>Full-genome deep sequencing and phylogenetic analysis of novel human betacoronavirus</article-title>
. <source>Emerging Infectious Diseases</source>
<year>2013</year>
; <volume>19</volume>
: <fpage>736</fpage>
–<lpage>742B</lpage>
.<pub-id pub-id-type="pmid">23693015</pub-id>
</mixed-citation>
</ref>
<ref id="ref12"><label>12.</label>
<mixed-citation publication-type="journal"><name><surname>Li</surname>
<given-names>W</given-names>
</name>
, <etal></etal>
<article-title>Bats are natural reservoirs of SARS-like coronaviruses</article-title>
. <source>Science</source>
<year>2005</year>
; <volume>310</volume>
: <fpage>676</fpage>
–<lpage>679</lpage>
.<pub-id pub-id-type="pmid">16195424</pub-id>
</mixed-citation>
</ref>
<ref id="ref13"><label>13.</label>
<mixed-citation publication-type="journal"><collab>CDC</collab>
. <article-title>Update: outbreak of severe acute respiratory syndrome–worldwide, 2003</article-title>
. <source>Morbidity and Mortality Weekly Report</source>
<year>2003</year>
; <volume>52</volume>
: <fpage>241</fpage>
–<lpage>246</lpage>
, 248.<pub-id pub-id-type="pmid">12680518</pub-id>
</mixed-citation>
</ref>
<ref id="ref14"><label>14.</label>
<mixed-citation publication-type="journal"><name><surname>Riley</surname>
<given-names>S</given-names>
</name>
, <etal></etal>
<article-title>Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions</article-title>
. <source>Science</source>
<year>2003</year>
; <volume>300</volume>
: <fpage>1961</fpage>
–<lpage>1966</lpage>
.<pub-id pub-id-type="pmid">12766206</pub-id>
</mixed-citation>
</ref>
<ref id="ref15"><label>15.</label>
<mixed-citation publication-type="journal"><name><surname>Chau</surname>
<given-names>PH</given-names>
</name>
, <name><surname>Yip</surname>
<given-names>PS</given-names>
</name>
. <article-title>Monitoring the severe acute respiratory syndrome epidemic and assessing effectiveness of interventions in Hong Kong Special Administrative Region</article-title>
. <source>Journal of Epidemiology and Community Health</source>
<year>2003</year>
; <volume>57</volume>
: <fpage>766</fpage>
–<lpage>769</lpage>
.<pub-id pub-id-type="pmid">14573569</pub-id>
</mixed-citation>
</ref>
<ref id="ref16"><label>16.</label>
<mixed-citation publication-type="book"><name><surname>Sterman</surname>
<given-names>J</given-names>
</name>
. <source>Business Dynamics: Systems Thinking and Modeling for a Complex World</source>
. <publisher-loc>Boston, USA</publisher-loc>
: <publisher-name>Irwin/McGraw-Hill</publisher-name>
, <year>2000</year>
.</mixed-citation>
</ref>
<ref id="ref17"><label>17.</label>
<mixed-citation publication-type="journal"><name><surname>Homer</surname>
<given-names>JB</given-names>
</name>
, <name><surname>Hirsch</surname>
<given-names>GB</given-names>
</name>
. <article-title>System dynamics modeling for public health: background and opportunities</article-title>
. <source>American Journal of Public Health</source>
<year>2006</year>
; <volume>96</volume>
: <fpage>452</fpage>
–<lpage>458</lpage>
.<pub-id pub-id-type="pmid">16449591</pub-id>
</mixed-citation>
</ref>
<ref id="ref18"><label>18.</label>
<mixed-citation publication-type="journal"><name><surname>Jones</surname>
<given-names>AP</given-names>
</name>
, <etal></etal>
<article-title>Understanding diabetes population dynamics through simulation modeling and experimentation</article-title>
. <source>American Journal of Public Health</source>
<year>2006</year>
; <volume>96</volume>
: <fpage>488</fpage>
–<lpage>494</lpage>
.<pub-id pub-id-type="pmid">16449587</pub-id>
</mixed-citation>
</ref>
<ref id="ref19"><label>19.</label>
<mixed-citation publication-type="journal"><name><surname>Koopman</surname>
<given-names>JS</given-names>
</name>
, <name><surname>Jacquez</surname>
<given-names>G</given-names>
</name>
, <name><surname>Chick</surname>
<given-names>SE</given-names>
</name>
. <article-title>New data and tools for integrating discrete and continuous population modeling strategies</article-title>
. <source>Annals of the New York Academy of Sciences</source>
<year>2001</year>
; <volume>954</volume>
: <fpage>268</fpage>
–<lpage>294</lpage>
.<pub-id pub-id-type="pmid">11797861</pub-id>
</mixed-citation>
</ref>
<ref id="ref20"><label>20.</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>
. <source>Infectious Diseases of Humans: Dynamics and Control</source>
. <publisher-loc>New York</publisher-loc>
: <publisher-name>Oxford Science Publication</publisher-name>
, <year>1991</year>
.</mixed-citation>
</ref>
<ref id="ref21"><label>21.</label>
<mixed-citation publication-type="journal"><name><surname>Omony</surname>
<given-names>J</given-names>
</name>
, <etal></etal>
<article-title>Modeling and analysis of the dynamic behavior of the XlnR regulon in Aspergillus niger</article-title>
. <source>BMC Systems Biology</source>
<year>2011</year>
; <volume>5</volume>
(<issue>Suppl. 1</issue>
): <fpage>S14</fpage>
.</mixed-citation>
</ref>
<ref id="ref22"><label>22.</label>
<mixed-citation publication-type="journal"><collab>CDC</collab>
. <article-title>Severe acute respiratory syndrome–Singapore, 2003</article-title>
. <source>Morbidity and Mortality Weekly Report</source>
<year>2003</year>
; <volume>52</volume>
: <fpage>405</fpage>
–<lpage>411</lpage>
.<pub-id pub-id-type="pmid">12807088</pub-id>
</mixed-citation>
</ref>
<ref id="ref23"><label>23.</label>
<mixed-citation publication-type="journal"><name><surname>Pitzer</surname>
<given-names>VE</given-names>
</name>
, <name><surname>Leung</surname>
<given-names>GM</given-names>
</name>
, <name><surname>Lipsitch</surname>
<given-names>M</given-names>
</name>
. <article-title>Estimating variability in the transmission of severe acute respiratory syndrome to household contacts in Hong Kong, China</article-title>
. <source>American Journal of Epidemiology</source>
<year>2007</year>
; <volume>166</volume>
: <fpage>355</fpage>
–<lpage>363</lpage>
.<pub-id pub-id-type="pmid">17493952</pub-id>
</mixed-citation>
</ref>
<ref id="ref24"><label>24.</label>
<mixed-citation publication-type="journal"><name><surname>McBryde</surname>
<given-names>ES</given-names>
</name>
, <etal></etal>
<article-title>Bayesian modelling of an epidemic of severe acute respiratory syndrome</article-title>
. <source>Bulletin of Mathematical Biology</source>
<year>2006</year>
; <volume>68</volume>
: <fpage>889</fpage>
–<lpage>917</lpage>
.<pub-id pub-id-type="pmid">16802088</pub-id>
</mixed-citation>
</ref>
<ref id="ref25"><label>25.</label>
<mixed-citation publication-type="journal"><name><surname>Hirose</surname>
<given-names>H</given-names>
</name>
. <article-title>The mixed trunsored model with applications to SARS</article-title>
. <source>Mathematics and Computers in Simulation</source>
<year>2007</year>
; <volume>74</volume>
: <fpage>443</fpage>
–<lpage>453</lpage>
.</mixed-citation>
</ref>
<ref id="ref26"><label>26.</label>
<mixed-citation publication-type="journal"><name><surname>Fine</surname>
<given-names>PE</given-names>
</name>
, <name><surname>Clarkson</surname>
<given-names>JA</given-names>
</name>
. <article-title>Measles in England and Wales–I: An analysis of factors underlying seasonal patterns</article-title>
. <source>International Journal of Epidemiology</source>
<year>1982</year>
; <volume>11</volume>
: <fpage>5</fpage>
–<lpage>14</lpage>
.<pub-id pub-id-type="pmid">7085179</pub-id>
</mixed-citation>
</ref>
<ref id="ref27"><label>27.</label>
<mixed-citation publication-type="journal"><name><surname>Anderson</surname>
<given-names>RM</given-names>
</name>
, <name><surname>Grenfell</surname>
<given-names>BT</given-names>
</name>
. <article-title>Quantitative investigations of different vaccination policies for the control of congenital rubella syndrome (CRS) in the United Kingdom</article-title>
. <source>Journal of Hygiene</source>
<year>1986</year>
; <volume>96</volume>
: <fpage>305</fpage>
–<lpage>333</lpage>
.<pub-id pub-id-type="pmid">3701044</pub-id>
</mixed-citation>
</ref>
<ref id="ref28"><label>28.</label>
<mixed-citation publication-type="journal"><name><surname>Anderson</surname>
<given-names>RM</given-names>
</name>
, <name><surname>May</surname>
<given-names>RM</given-names>
</name>
. <article-title>Age-related changes in the rate of disease transmission: implications for the design of vaccination programmes</article-title>
. <source>Journal of Hygiene</source>
<year>1985</year>
; <volume>94</volume>
: <fpage>365</fpage>
–<lpage>436</lpage>
.<pub-id pub-id-type="pmid">4008922</pub-id>
</mixed-citation>
</ref>
<ref id="ref29"><label>29.</label>
<mixed-citation publication-type="journal"><name><surname>Edmunds</surname>
<given-names>WJ</given-names>
</name>
, <name><surname>O'Callaghan</surname>
<given-names>CJ</given-names>
</name>
, <name><surname>Nokes</surname>
<given-names>DJ</given-names>
</name>
. <article-title>Who mixes with whom? A method to determine the contact patterns of adults that may lead to the spread of airborne infections</article-title>
. <source>Proceedings of the Royal Society of London, Series B: Biological Sciences</source>
<year>1997</year>
; <volume>264</volume>
: <fpage>949</fpage>
–<lpage>957</lpage>
.<pub-id pub-id-type="pmid">9263464</pub-id>
</mixed-citation>
</ref>
<ref id="ref30"><label>30.</label>
<mixed-citation publication-type="journal"><name><surname>Mossong</surname>
<given-names>J</given-names>
</name>
, <etal></etal>
<article-title>Social contacts and mixing patterns relevant to the spread of infectious diseases</article-title>
. <source>PLoS Medicine</source>
<year>2008</year>
; <volume>5</volume>
: <elocation-id>e74</elocation-id>
.<pub-id pub-id-type="pmid">18366252</pub-id>
</mixed-citation>
</ref>
<ref id="ref31"><label>31.</label>
<mixed-citation publication-type="journal"><name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
, <name><surname>Teunis</surname>
<given-names>P</given-names>
</name>
. <article-title>Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures</article-title>
. <source>American Journal of Epidemiology</source>
<year>2004</year>
; <volume>160</volume>
: <fpage>509</fpage>
–<lpage>516</lpage>
.<pub-id pub-id-type="pmid">15353409</pub-id>
</mixed-citation>
</ref>
<ref id="ref32"><label>32.</label>
<mixed-citation publication-type="journal"><name><surname>Royston</surname>
<given-names>G</given-names>
</name>
, <etal></etal>
<article-title>Using system dynamics to help develop and implement policies and programmes in health care in England</article-title>
. <source>System Dynamics Review</source>
<year>1999</year>
; <volume>15</volume>
: <fpage>293</fpage>
–<lpage>313</lpage>
.</mixed-citation>
</ref>
<ref id="ref33"><label>33.</label>
<mixed-citation publication-type="journal"><name><surname>Morilla</surname>
<given-names>F</given-names>
</name>
, <etal></etal>
<article-title>Work disability related to musculoskeletal pain: a system dynamics approach</article-title>
. <source>Journal of Musculoskeletal Pain</source>
<year>2014</year>
; <volume>22</volume>
: <fpage>51</fpage>
–<lpage>61</lpage>
.</mixed-citation>
</ref>
<ref id="ref34"><label>34.</label>
<mixed-citation publication-type="journal"><name><surname>Hirsch</surname>
<given-names>G</given-names>
</name>
, <name><surname>Miller</surname>
<given-names>S</given-names>
</name>
. <article-title>Evaluating HMO policies with a computer simulation model</article-title>
. <source>Medical Care</source>
<year>1974</year>
; <volume>12</volume>
: <fpage>668</fpage>
–<lpage>681</lpage>
.<pub-id pub-id-type="pmid">4604449</pub-id>
</mixed-citation>
</ref>
<ref id="ref35"><label>35.</label>
<mixed-citation publication-type="journal"><name><surname>Luginbuhl</surname>
<given-names>WH</given-names>
</name>
, <etal></etal>
<article-title>Prevention and rehabilitation as a means of cost containment: the example of myocardial infarction</article-title>
. <source>Journal of Public Health Policy</source>
<year>1981</year>
; <volume>2</volume>
: <fpage>103</fpage>
–<lpage>115</lpage>
.<pub-id pub-id-type="pmid">6788801</pub-id>
</mixed-citation>
</ref>
<ref id="ref36"><label>36.</label>
<mixed-citation publication-type="journal"><name><surname>Roberts</surname>
<given-names>C</given-names>
</name>
, <name><surname>Dangerfield</surname>
<given-names>B</given-names>
</name>
. <article-title>Modelling the epidemiological consequences of HIV infection and AIDS: a contribution from operational research</article-title>
. <source>Journal of the Operational Research Society</source>
<year>1990</year>
; <volume>41</volume>
: <fpage>273</fpage>
–<lpage>289</lpage>
</mixed-citation>
</ref>
<ref id="ref37"><label>37.</label>
<mixed-citation publication-type="journal"><name><surname>Dangerfield</surname>
<given-names>B</given-names>
</name>
, <name><surname>Fang</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Roberts</surname>
<given-names>C</given-names>
</name>
. <article-title>Model based scenarios for the epidemiology of HIV/AIDS: the consequences of highly active antiretroviral therapy</article-title>
. <source>System Dynamics Review</source>
<year>2001</year>
; <volume>17</volume>
: <fpage>119</fpage>
–<lpage>150</lpage>
.</mixed-citation>
</ref>
<ref id="ref38"><label>38.</label>
<mixed-citation publication-type="journal"><name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr.</suffix>
</name>
, <etal></etal>
<article-title>Containing pandemic influenza at the source</article-title>
. <source>Science</source>
<year>2005</year>
; <volume>309</volume>
: <fpage>1083</fpage>
–<lpage>1087</lpage>
.<pub-id pub-id-type="pmid">16079251</pub-id>
</mixed-citation>
</ref>
<ref id="ref39"><label>39.</label>
<mixed-citation publication-type="journal"><name><surname>Halder</surname>
<given-names>N</given-names>
</name>
, <name><surname>Kelso</surname>
<given-names>JK</given-names>
</name>
, <name><surname>Milne</surname>
<given-names>GJ</given-names>
</name>
. <article-title>Cost-effective strategies for mitigating a future influenza pandemic with H1N1 2009 characteristics</article-title>
. <source>PLoS ONE</source>
<year>2011</year>
; <volume>6</volume>
: <elocation-id>e22087</elocation-id>
.<pub-id pub-id-type="pmid">21760957</pub-id>
</mixed-citation>
</ref>
<ref id="ref40"><label>40.</label>
<mixed-citation publication-type="journal"><name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <etal></etal>
<article-title>Real-time estimates in early detection of SARS</article-title>
. <source>Emerging Infectious Diseases</source>
<year>2006</year>
; <volume>12</volume>
: <fpage>110</fpage>
–<lpage>113</lpage>
.<pub-id pub-id-type="pmid">16494726</pub-id>
</mixed-citation>
</ref>
<ref id="ref41"><label>41.</label>
<mixed-citation publication-type="journal"><name><surname>Pitman</surname>
<given-names>RJ</given-names>
</name>
, <etal></etal>
<article-title>Entry screening for severe acute respiratory syndrome (SARS) or influenza: policy evaluation</article-title>
. <source>British Medical Journal (Clinical Research Edition)</source>
<year>2005</year>
; <volume>331</volume>
: <fpage>1242</fpage>
–<lpage>1243</lpage>
.</mixed-citation>
</ref>
<ref id="ref42"><label>42.</label>
<mixed-citation publication-type="journal"><name><surname>Lipsitch</surname>
<given-names>M</given-names>
</name>
, <etal></etal>
<article-title>Transmission dynamics and control of severe acute respiratory syndrome</article-title>
. <source>Science</source>
<year>2003</year>
; <volume>300</volume>
: <fpage>1966</fpage>
–<lpage>1970</lpage>
.<pub-id pub-id-type="pmid">12766207</pub-id>
</mixed-citation>
</ref>
<ref id="ref43"><label>43.</label>
<mixed-citation publication-type="journal"><name><surname>Grassly</surname>
<given-names>NC</given-names>
</name>
, <name><surname>Fraser</surname>
<given-names>C</given-names>
</name>
. <article-title>Mathematical models of infectious disease transmission</article-title>
. <source>Nature Reviews Microbiology</source>
<year>2008</year>
; <volume>6</volume>
: <fpage>477</fpage>
–<lpage>487</lpage>
.<pub-id pub-id-type="pmid">18533288</pub-id>
</mixed-citation>
</ref>
<ref id="ref44"><label>44.</label>
<mixed-citation publication-type="journal"><name><surname>Hui</surname>
<given-names>DS</given-names>
</name>
, <name><surname>Memish</surname>
<given-names>ZA</given-names>
</name>
, <name><surname>Zumla</surname>
<given-names>A</given-names>
</name>
. <article-title>Severe acute respiratory syndrome vs</article-title>
. <article-title>the Middle East respiratory syndrome</article-title>
. <source>Current Opinion in Pulmonary Medicine</source>
<year>2014</year>
; <volume>20</volume>
: <fpage>233</fpage>
–<lpage>241</lpage>
.</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000C50 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000C50 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= MersV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:7108648 |texte= New coronavirus outbreak. Lessons learned from the severe acute respiratory syndrome epidemic }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:25591619" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a MersV1
![]() | This area was generated with Dilib version V0.6.33. | ![]() |