Optimizing the Precision of Case Fatality Ratio Estimates Under the Surveillance Pyramid Approach
Identifieur interne : 000904 ( Pmc/Corpus ); précédent : 000903; suivant : 000905Optimizing the Precision of Case Fatality Ratio Estimates Under the Surveillance Pyramid Approach
Auteurs : Camille Pelat ; Neil M. Ferguson ; Peter J. White ; Carrie Reed ; Lyn Finelli ; Simon Cauchemez ; Christophe FraserSource :
- American Journal of Epidemiology [ 0002-9262 ] ; 2014.
Abstract
In the management of emerging infectious disease epidemics, precise and accurate estimation of severity indices, such as the probability of death after developing symptoms—the symptomatic case fatality ratio (sCFR)—is essential. Estimation of the sCFR may require merging data gathered through different surveillance systems and surveys. Since different surveillance strategies provide different levels of precision and accuracy, there is need for a theory to help investigators select the strategy that maximizes these properties. Here, we study the precision of sCFR estimators that combine data from several levels of the severity pyramid. We derive a formula for the standard error, which helps us find the estimator with the best precision given fixed resources. We further propose rules of thumb for guiding the choice of strategy: For example, should surveillance of a particular severity level be started? Which level should be preferred? We derive a formula for the optimal allocation of resources between chosen surveillance levels and provide a simple approximation that can be used in thinking more heuristically about planning surveillance. We illustrate these concepts with numerical examples corresponding to 3 influenza pandemic scenarios. Finally, we review the equally important issue of accuracy.
Url:
DOI: 10.1093/aje/kwu213
PubMed: 25255809
PubMed Central: 4240167
Links to Exploration step
PMC:4240167Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Optimizing the Precision of Case Fatality Ratio Estimates Under the Surveillance Pyramid Approach</title>
<author><name sortKey="Pelat, Camille" sort="Pelat, Camille" uniqKey="Pelat C" first="Camille" last="Pelat">Camille Pelat</name>
</author>
<author><name sortKey="Ferguson, Neil M" sort="Ferguson, Neil M" uniqKey="Ferguson N" first="Neil M." last="Ferguson">Neil M. Ferguson</name>
</author>
<author><name sortKey="White, Peter J" sort="White, Peter J" uniqKey="White P" first="Peter J." last="White">Peter J. White</name>
</author>
<author><name sortKey="Reed, Carrie" sort="Reed, Carrie" uniqKey="Reed C" first="Carrie" last="Reed">Carrie Reed</name>
</author>
<author><name sortKey="Finelli, Lyn" sort="Finelli, Lyn" uniqKey="Finelli L" first="Lyn" last="Finelli">Lyn Finelli</name>
</author>
<author><name sortKey="Cauchemez, Simon" sort="Cauchemez, Simon" uniqKey="Cauchemez S" first="Simon" last="Cauchemez">Simon Cauchemez</name>
</author>
<author><name sortKey="Fraser, Christophe" sort="Fraser, Christophe" uniqKey="Fraser C" first="Christophe" last="Fraser">Christophe Fraser</name>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">25255809</idno>
<idno type="pmc">4240167</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4240167</idno>
<idno type="RBID">PMC:4240167</idno>
<idno type="doi">10.1093/aje/kwu213</idno>
<date when="2014">2014</date>
<idno type="wicri:Area/Pmc/Corpus">000904</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000904</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Optimizing the Precision of Case Fatality Ratio Estimates Under the Surveillance Pyramid Approach</title>
<author><name sortKey="Pelat, Camille" sort="Pelat, Camille" uniqKey="Pelat C" first="Camille" last="Pelat">Camille Pelat</name>
</author>
<author><name sortKey="Ferguson, Neil M" sort="Ferguson, Neil M" uniqKey="Ferguson N" first="Neil M." last="Ferguson">Neil M. Ferguson</name>
</author>
<author><name sortKey="White, Peter J" sort="White, Peter J" uniqKey="White P" first="Peter J." last="White">Peter J. White</name>
</author>
<author><name sortKey="Reed, Carrie" sort="Reed, Carrie" uniqKey="Reed C" first="Carrie" last="Reed">Carrie Reed</name>
</author>
<author><name sortKey="Finelli, Lyn" sort="Finelli, Lyn" uniqKey="Finelli L" first="Lyn" last="Finelli">Lyn Finelli</name>
</author>
<author><name sortKey="Cauchemez, Simon" sort="Cauchemez, Simon" uniqKey="Cauchemez S" first="Simon" last="Cauchemez">Simon Cauchemez</name>
</author>
<author><name sortKey="Fraser, Christophe" sort="Fraser, Christophe" uniqKey="Fraser C" first="Christophe" last="Fraser">Christophe Fraser</name>
</author>
</analytic>
<series><title level="j">American Journal of Epidemiology</title>
<idno type="ISSN">0002-9262</idno>
<idno type="eISSN">1476-6256</idno>
<imprint><date when="2014">2014</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>In the management of emerging infectious disease epidemics, precise and accurate estimation of severity indices, such as the probability of death after developing symptoms—the symptomatic case fatality ratio (sCFR)—is essential. Estimation of the sCFR may require merging data gathered through different surveillance systems and surveys. Since different surveillance strategies provide different levels of precision and accuracy, there is need for a theory to help investigators select the strategy that maximizes these properties. Here, we study the precision of sCFR estimators that combine data from several levels of the severity pyramid. We derive a formula for the standard error, which helps us find the estimator with the best precision given fixed resources. We further propose rules of thumb for guiding the choice of strategy: For example, should surveillance of a particular severity level be started? Which level should be preferred? We derive a formula for the optimal allocation of resources between chosen surveillance levels and provide a simple approximation that can be used in thinking more heuristically about planning surveillance. We illustrate these concepts with numerical examples corresponding to 3 influenza pandemic scenarios. Finally, we review the equally important issue of accuracy.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Reed, C" uniqKey="Reed C">C Reed</name>
</author>
<author><name sortKey="Biggerstaff, M" uniqKey="Biggerstaff M">M Biggerstaff</name>
</author>
<author><name sortKey="Finelli, L" uniqKey="Finelli L">L Finelli</name>
</author>
</analytic>
</biblStruct>
<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>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Collins, Sd" uniqKey="Collins S">SD Collins</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Garske, T" uniqKey="Garske T">T Garske</name>
</author>
<author><name sortKey="Legrand, J" uniqKey="Legrand J">J Legrand</name>
</author>
<author><name sortKey="Donnelly, Ca" uniqKey="Donnelly C">CA Donnelly</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H Nishiura</name>
</author>
<author><name sortKey="Klinkenberg, D" uniqKey="Klinkenberg D">D Klinkenberg</name>
</author>
<author><name sortKey="Roberts, M" uniqKey="Roberts M">M Roberts</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mishra, Ac" uniqKey="Mishra A">AC Mishra</name>
</author>
<author><name sortKey="Chadha, Ms" uniqKey="Chadha M">MS Chadha</name>
</author>
<author><name sortKey="Choudhary, Ml" uniqKey="Choudhary M">ML Choudhary</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wilson, N" uniqKey="Wilson N">N Wilson</name>
</author>
<author><name sortKey="Baker, Mg" uniqKey="Baker M">MG Baker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Presanis, Am" uniqKey="Presanis A">AM Presanis</name>
</author>
<author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name>
</author>
<author><name sortKey="De Angelis, D" uniqKey="De Angelis D">D De Angelis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Presanis, Am" uniqKey="Presanis A">AM Presanis</name>
</author>
<author><name sortKey="De Angelis, D" uniqKey="De Angelis D">D De Angelis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Donaldson, Lj" uniqKey="Donaldson L">LJ Donaldson</name>
</author>
<author><name sortKey="Rutter, Pd" uniqKey="Rutter P">PD Rutter</name>
</author>
<author><name sortKey="Ellis, Bm" uniqKey="Ellis B">BM Ellis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brooks Pollock, E" uniqKey="Brooks Pollock E">E Brooks-Pollock</name>
</author>
<author><name sortKey="Tilston, N" uniqKey="Tilston N">N Tilston</name>
</author>
<author><name sortKey="Edmunds, Wj" uniqKey="Edmunds W">WJ Edmunds</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Steens, A" uniqKey="Steens A">A Steens</name>
</author>
<author><name sortKey="Waaijenborg, S" uniqKey="Waaijenborg S">S Waaijenborg</name>
</author>
<author><name sortKey="Teunis, Pf" uniqKey="Teunis P">PF Teunis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Presanis, Am" uniqKey="Presanis A">AM Presanis</name>
</author>
<author><name sortKey="Pebody, Rg" uniqKey="Pebody R">RG Pebody</name>
</author>
<author><name sortKey="Paterson, Bj" uniqKey="Paterson B">BJ Paterson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Armitage, P" uniqKey="Armitage P">P Armitage</name>
</author>
<author><name sortKey="Berry, G" uniqKey="Berry G">G Berry</name>
</author>
<author><name sortKey="Matthews, J" uniqKey="Matthews J">J Matthews</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>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wu, Jt" uniqKey="Wu J">JT Wu</name>
</author>
<author><name sortKey="Ho, A" uniqKey="Ho A">A Ho</name>
</author>
<author><name sortKey="Ma, Es" uniqKey="Ma E">ES Ma</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Atkinson, A" uniqKey="Atkinson A">A Atkinson</name>
</author>
<author><name sortKey="Donev, A" uniqKey="Donev A">A Donev</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Comin, A" uniqKey="Comin A">A Comin</name>
</author>
<author><name sortKey="Stegeman, A" uniqKey="Stegeman A">A Stegeman</name>
</author>
<author><name sortKey="Marangon, S" uniqKey="Marangon S">S Marangon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Alba, A" uniqKey="Alba A">A Alba</name>
</author>
<author><name sortKey="Casal, J" uniqKey="Casal J">J Casal</name>
</author>
<author><name sortKey="Napp, S" uniqKey="Napp S">S Napp</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Baptista, Fm" uniqKey="Baptista F">FM Baptista</name>
</author>
<author><name sortKey="Halasa, T" uniqKey="Halasa T">T Halasa</name>
</author>
<author><name sortKey="Alban, L" uniqKey="Alban L">L Alban</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hadorn, Dc" uniqKey="Hadorn D">DC Hadorn</name>
</author>
<author><name sortKey="Racloz, V" uniqKey="Racloz V">V Racloz</name>
</author>
<author><name sortKey="Schwermer, H" uniqKey="Schwermer H">H Schwermer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yamamoto, T" uniqKey="Yamamoto T">T Yamamoto</name>
</author>
<author><name sortKey="Tsutsui, T" uniqKey="Tsutsui T">T Tsutsui</name>
</author>
<author><name sortKey="Nishiguchi, A" uniqKey="Nishiguchi A">A Nishiguchi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Polgreen, Pm" uniqKey="Polgreen P">PM Polgreen</name>
</author>
<author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name>
</author>
<author><name sortKey="Segre, Am" uniqKey="Segre A">AM Segre</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Scarpino, Sv" uniqKey="Scarpino S">SV Scarpino</name>
</author>
<author><name sortKey="Dimitrov, Nb" uniqKey="Dimitrov N">NB Dimitrov</name>
</author>
<author><name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yang, P" uniqKey="Yang P">P Yang</name>
</author>
<author><name sortKey="Duan, W" uniqKey="Duan W">W Duan</name>
</author>
<author><name sortKey="Lv, M" uniqKey="Lv M">M Lv</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Byass, P" uniqKey="Byass P">P Byass</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ejima, K" uniqKey="Ejima K">K Ejima</name>
</author>
<author><name sortKey="Omori, R" uniqKey="Omori R">R Omori</name>
</author>
<author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Reich, Ng" uniqKey="Reich N">NG Reich</name>
</author>
<author><name sortKey="Lessler, J" uniqKey="Lessler J">J Lessler</name>
</author>
<author><name sortKey="Cummings, Da" uniqKey="Cummings D">DA Cummings</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ghani, Ac" uniqKey="Ghani A">AC Ghani</name>
</author>
<author><name sortKey="Donnelly, Ca" uniqKey="Donnelly C">CA Donnelly</name>
</author>
<author><name sortKey="Cox, Dr" uniqKey="Cox D">DR Cox</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wu, Jt" uniqKey="Wu J">JT Wu</name>
</author>
<author><name sortKey="Ma, Es" uniqKey="Ma E">ES Ma</name>
</author>
<author><name sortKey="Lee, Ck" uniqKey="Lee C">CK Lee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Payne, Am" uniqKey="Payne A">AM Payne</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">Am J Epidemiol</journal-id>
<journal-id journal-id-type="iso-abbrev">Am. J. Epidemiol</journal-id>
<journal-id journal-id-type="publisher-id">aje</journal-id>
<journal-id journal-id-type="hwp">amjepid</journal-id>
<journal-title-group><journal-title>American Journal of Epidemiology</journal-title>
</journal-title-group>
<issn pub-type="ppub">0002-9262</issn>
<issn pub-type="epub">1476-6256</issn>
<publisher><publisher-name>Oxford University Press</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">25255809</article-id>
<article-id pub-id-type="pmc">4240167</article-id>
<article-id pub-id-type="doi">10.1093/aje/kwu213</article-id>
<article-id pub-id-type="publisher-id">kwu213</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Practice of Epidemiology</subject>
</subj-group>
</article-categories>
<title-group><article-title>Optimizing the Precision of Case Fatality Ratio Estimates Under the Surveillance Pyramid Approach</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Pelat</surname>
<given-names>Camille</given-names>
</name>
<xref ref-type="corresp" rid="cor1">*</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ferguson</surname>
<given-names>Neil M.</given-names>
</name>
</contrib>
<contrib contrib-type="author"><name><surname>White</surname>
<given-names>Peter J.</given-names>
</name>
</contrib>
<contrib contrib-type="author"><name><surname>Reed</surname>
<given-names>Carrie</given-names>
</name>
</contrib>
<contrib contrib-type="author"><name><surname>Finelli</surname>
<given-names>Lyn</given-names>
</name>
</contrib>
<contrib contrib-type="author"><name><surname>Cauchemez</surname>
<given-names>Simon</given-names>
</name>
</contrib>
<contrib contrib-type="author"><name><surname>Fraser</surname>
<given-names>Christophe</given-names>
</name>
</contrib>
</contrib-group>
<author-notes><corresp id="cor1"><label>*</label>
Correspondence to Dr. Camille Pelat, Institut de Veille Sanitaire, Département des Maladies Infectieuses, 12 rue du Val d'Osne, 94415 Saint-Maurice Cedex, France (e-mail: <email>c.pelat@invs.sante.fr</email>
).</corresp>
<fn fn-type="con" id="AN1"><p>Abbreviations: sCFR, symptomatic case fatality ratio; SE, standard error.</p>
</fn>
</author-notes>
<pub-date pub-type="ppub"><day>15</day>
<month>11</month>
<year>2014</year>
</pub-date>
<pub-date pub-type="epub"><day>25</day>
<month>9</month>
<year>2014</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>25</day>
<month>9</month>
<year>2014</year>
</pub-date>
<pmc-comment> PMC Release delay is 0 months and 0 days and was based on the
. </pmc-comment>
<volume>180</volume>
<issue>10</issue>
<fpage>1036</fpage>
<lpage>1046</lpage>
<history><date date-type="received"><day>28</day>
<month>3</month>
<year>2013</year>
</date>
<date date-type="accepted"><day>17</day>
<month>7</month>
<year>2014</year>
</date>
</history>
<permissions><copyright-statement>© The Author 2014. Published by Oxford University Press on behalf of the Johns Hopkins Bloomberg School of Public Health.</copyright-statement>
<copyright-year>2014</copyright-year>
<license license-type="creative-commons" 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 (<uri xlink:type="simple" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</uri>
), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:type="simple" xlink:href="kwu213.pdf"></self-uri>
<abstract><p>In the management of emerging infectious disease epidemics, precise and accurate estimation of severity indices, such as the probability of death after developing symptoms—the symptomatic case fatality ratio (sCFR)—is essential. Estimation of the sCFR may require merging data gathered through different surveillance systems and surveys. Since different surveillance strategies provide different levels of precision and accuracy, there is need for a theory to help investigators select the strategy that maximizes these properties. Here, we study the precision of sCFR estimators that combine data from several levels of the severity pyramid. We derive a formula for the standard error, which helps us find the estimator with the best precision given fixed resources. We further propose rules of thumb for guiding the choice of strategy: For example, should surveillance of a particular severity level be started? Which level should be preferred? We derive a formula for the optimal allocation of resources between chosen surveillance levels and provide a simple approximation that can be used in thinking more heuristically about planning surveillance. We illustrate these concepts with numerical examples corresponding to 3 influenza pandemic scenarios. Finally, we review the equally important issue of accuracy.</p>
</abstract>
<kwd-group><kwd>case fatality ratio</kwd>
<kwd>emerging infectious diseases</kwd>
<kwd>influenza</kwd>
<kwd>pandemics</kwd>
<kwd>statistical planning</kwd>
<kwd>surveillance protocol</kwd>
</kwd-group>
</article-meta>
</front>
<body><p>During outbreaks of emerging infectious diseases, the symptomatic case fatality ratio (sCFR)—the probability of death following the development of symptoms—is a critical summary statistic characterizing disease severity, along with the probability of hospitalization given symptoms and the probability of symptoms given infection (<xref rid="KWU213C1" ref-type="bibr">1</xref>
, <xref rid="KWU213C2" ref-type="bibr">2</xref>
). Estimates of the sCFR influence the public health measures put in place to control an epidemic. Large-scale public health responses are expensive and socially disruptive; health authorities often face a difficult tradeoff between mitigation and the costs to society (<xref rid="KWU213C1" ref-type="bibr">1</xref>
). Accurate and precise estimates of the sCFR are essential for decision-making at the time response plans are being drawn up and as the plans are periodically revised during the course of the epidemic.</p>
<p>The 2009 A/H1N1 influenza pandemic illustrated well the challenges involved in rapid severity assessment. It was clear from early in the outbreak (<xref rid="KWU213C3" ref-type="bibr">3</xref>
) that the severity of the epidemic was substantially less than that of the 1918 A/H1N1 pandemic (for which the sCFR was approximately 0.02 (<xref rid="KWU213C4" ref-type="bibr">4</xref>
)), yet it proved more difficult to determine whether severity was <italic>intermediate</italic>
(sCFR ≈ 10<sup>−3</sup>
) or <italic>mild</italic>
(sCFR ≈ 10<sup>−4</sup>
) (<xref rid="KWU213C5" ref-type="bibr">5</xref>
–<xref rid="KWU213C13" ref-type="bibr">13</xref>
). Robust, precise estimates establishing that the pandemic had mild severity (sCFR = 4.5 × 10<sup>−4</sup>
, 95% credible interval: 2 × 10<sup>−4</sup>
, 9 × 10<sup>−4</sup>
) were published by July 2009 (<xref rid="KWU213C9" ref-type="bibr">9</xref>
).</p>
<p>In the context of an epidemic of <italic>mild</italic>
severity, sample sizes required to directly estimate the sCFR by recording fatalities occurring in a series of symptomatic cases may become prohibitively large. As a consequence, several authors have developed pyramidal approaches to estimating the sCFR as a product of conditional probabilities—for example, the probability of hospitalization given symptoms times the probability of death upon symptom-related hospitalization (<xref rid="KWU213C5" ref-type="bibr">5</xref>
, <xref rid="KWU213C10" ref-type="bibr">10</xref>
, <xref rid="KWU213C14" ref-type="bibr">14</xref>
). Strong assumptions underpin pyramidal approaches, particularly the assumption that deceased cases progress through the entire pyramid, which can limit their validity. However, these approaches present so many practical advantages that they were extensively used during the 2009 pandemic, and they can sometimes be the only available alternative.</p>
<p>Despite this fact, there is still little understanding of the statistical properties of pyramidal sCFR estimators. The standard error (SE), in particular, indicates when such estimators are more precise than direct estimation in a series of symptomatic cases. More generally, different pyramidal estimators, combining different surveillance systems and surveys, have different levels of precision. There is need for a theory to help in selecting the most precise estimator. Once the estimator has been selected, it is also unclear what the optimal allocation of resources between the different surveillance levels should be: Would it be better to conduct more outbreak investigations or to increase hospital-based surveillance? In an emerging infectious disease outbreak, where resources are finite, a clear strategy with which to efficiently reduce uncertainty around sCFR estimates is essential for management and planning.</p>
<p>In this article, we study the precision of pyramidal approaches to sCFR estimation and examine how the choice of the best approach depends upon outbreak characteristics. We propose rules of thumb for finding the most precise estimator and for optimizing resource allocation between surveillance levels. The general concepts are illustrated using 3 influenza pandemic scenarios with different levels of severity. We also discuss the issue of accuracy, which we propose may be best assessed on a case-by-case basis.</p>
<sec sec-type="methods" id="s2"><title>METHODS</title>
<p>Pyramidal approaches to sCFR estimation combine data from several levels of the severity pyramid, making assumptions about the course of the disease and about health-care organization (<xref rid="KWU213C5" ref-type="bibr">5</xref>
, <xref rid="KWU213C10" ref-type="bibr">10</xref>
, <xref rid="KWU213C14" ref-type="bibr">14</xref>
). For instance, in the severity pyramid of Figure <xref ref-type="fig" rid="KWU213F1">1</xref>
, symptomatic cases (<italic>S</italic>
) go through 2 severity levels before death (<italic>D</italic>
): medical attention (<italic>M</italic>
) and hospitalization (<italic>H</italic>
). Assuming that all cases go through <italic>M</italic>
before <italic>H</italic>
and go through <italic>H</italic>
before <italic>D</italic>
, the sCFR can be written as a product of progression probabilities: sCFR = <italic>P</italic>
(<italic>D</italic>
|<italic>H</italic>
) × <italic>P</italic>
(<italic>H</italic>
|<italic>M</italic>
) × <italic>P</italic>
(<italic>M</italic>
|<italic>S</italic>
) (see <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 1, available at http://aje.oxfordjournals.org/</ext-link>
). A natural estimate of the sCFR is the product of the progression probability estimates. According to the available data at each severity level, different pyramidal estimators are possible, as long as they are based on reasonable assumptions. Figure <xref ref-type="fig" rid="KWU213F1">1</xref>
illustrates 4 estimators of the sCFR.
<fig id="KWU213F1" position="float"><label>Figure 1.</label>
<caption><p>Classification of symptomatic cases of infectious disease into a severity pyramid comprising medical attention, hospitalization, and death (all related to symptomatic infection). A) Probabilities of progression to severity level <italic>B</italic>
for persons on level <italic>A</italic>
, <italic>P</italic>
(<italic>B</italic>
|<italic>A</italic>
), denoted <italic>p<sub>B</sub>
</italic>
<sub>|<italic>A</italic>
</sub>
; for example, the probability of medical attention following symptoms is <italic>p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
. Plain arrows indicate progression to the next level, and the dashed arrows indicate progression through intermediate level(s). B) Four estimators (numbered 1–4) of the symptomatic case fatality ratio, obtained either directly by following up symptomatic cases until death <inline-formula><mml:math id="ImEquation1"><mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
or by multiplying progression probabilities estimated from 2–3 surveillance levels. Outbreak investigations, household contact surveys, or symptomatic case series can be used to derive <inline-formula><mml:math id="ImEquation2"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>M</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
, <inline-formula><mml:math id="ImEquation3"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
, and <inline-formula><mml:math id="ImEquation4"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
; a surveillance system based on sentinel general practitioners can be used to derive <inline-formula><mml:math id="ImEquation5"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="ImEquation6"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
; and a hospital-based surveillance system can be used to derive <inline-formula><mml:math id="ImEquation7"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
.</p>
</caption>
<graphic xlink:href="kwu21301"></graphic>
</fig>
</p>
<p>More generally, we are interested in comparing <italic>K</italic>
pyramidal estimators, labeled <italic>k</italic>
= 1 … <italic>K</italic>
, with <italic>N<sub>k</sub>
</italic>
levels each. The sCFR in the <italic>k</italic>
th estimation strategy is <inline-formula><mml:math id="ImEquation8"><mml:mtext>sCFR</mml:mtext>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mo movablelimits="false">∏</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:msubsup>
<mml:mrow><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
, where <italic>p<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
denotes the probability of progression from level <italic>i</italic>
in strategy <italic>k</italic>
to the next level and is estimated in a sample of <italic>n<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
cases. The true value of the sCFR is independent of the estimation strategy, but the estimator <inline-formula><mml:math id="ImEquation9"><mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
<mml:msub><mml:mrow><mml:mover><mml:mrow><mml:mi mathvariant="normal">CFR</mml:mi>
</mml:mrow>
<mml:mo>ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mo movablelimits="false">∏</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:msubsup>
<mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
and its SE are not. Large SEs indicate low precision. A first-order approximation of the SE of pyramidal sCFR estimators is obtained with the delta method (<xref rid="KWU213C15" ref-type="bibr">15</xref>
) (<ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 2</ext-link>
):
<disp-formula id="KWU213M1"><label>(1)</label>
<mml:math id="DmEquation1"><mml:mtext>SE</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
<mml:msub><mml:mrow><mml:mover><mml:mrow><mml:mi mathvariant="normal">CFR</mml:mi>
</mml:mrow>
<mml:mo>ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>≈</mml:mo>
<mml:mtext>sCFR</mml:mtext>
<mml:msqrt><mml:munderover><mml:mo movablelimits="false">∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:munderover>
<mml:mrow><mml:mstyle><mml:mrow><mml:mfrac><mml:mn>1</mml:mn>
<mml:mrow><mml:mrow><mml:msub><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mfenced open="(" close=")"><mml:mrow><mml:mstyle><mml:mrow><mml:mfrac><mml:mn>1</mml:mn>
<mml:mrow><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mstyle>
</mml:mrow>
</mml:mfenced>
</mml:msqrt>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
</p>
<p>The core of our paper is to study the precision of sCFR pyramidal estimators to find 1) the most precise estimation strategy and 2) the allocation of resources between a strategy's surveillance levels that maximizes precision.</p>
<sec id="s2a"><title>Optimizing resource allocation in a pyramidal approach</title>
<p>Sample sizes are constrained by the available number of cases at each severity level but also by the finite amount of (monetary) resources available for surveillance. We denote as <italic>c<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
the cost of recruiting a case into the <italic>i</italic>
th surveillance level of strategy <italic>k</italic>
, so that <inline-formula><mml:math id="ImEquation10"><mml:msubsup><mml:mo movablelimits="false">∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:msubsup>
<mml:mrow><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow><mml:msub><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
is the cost of estimating sCFR with strategy <italic>k</italic>
and sample sizes {<italic>n<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
}. We model resource constraints by assuming that there is a fixed budget <italic>C</italic>
, and for each strategy <italic>k</italic>
we ask, “What is the optimal allocation of resources to the different surveillance levels?” This involves finding sample sizes {<italic>n<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
} that minimize the SE subject to the constraint of the fixed budget: <inline-formula><mml:math id="ImEquation11"><mml:msubsup><mml:mo movablelimits="false">∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:msubsup>
<mml:mrow><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow><mml:msub><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>C</mml:mi>
</mml:math>
</inline-formula>
. We show in <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 3</ext-link>
that this minimum is achieved for the following sample sizes:
<disp-formula id="KWU213M2"><label>(2)</label>
<mml:math id="DmEquation2"><mml:msubsup><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>∗</mml:mo>
</mml:msubsup>
<mml:mo>≈</mml:mo>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mi>C</mml:mi>
<mml:mrow><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mrow><mml:msqrt><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:msqrt>
</mml:mrow>
<mml:mrow><mml:munderover><mml:mo movablelimits="false">∑</mml:mo>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:munderover>
<mml:mrow><mml:msqrt><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:msqrt>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">∀</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo>,</mml:mo>
<mml:mspace width="thinmathspace"></mml:mspace>
<mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mstyle>
</mml:mstyle>
</mml:math>
</disp-formula>
</p>
<p>A first observation from equation 2 is that the sample size <inline-formula><mml:math id="ImEquation12"><mml:msubsup><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>∗</mml:mo>
</mml:msubsup>
</mml:math>
</inline-formula>
increases as <italic>p<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
decreases. This is obtained by studying the sign of <inline-formula><mml:math id="ImEquation13"><mml:mi mathvariant="normal">∂</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:msubsup><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>∗</mml:mo>
</mml:msubsup>
<mml:mo stretchy="false">)</mml:mo>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:mi mathvariant="normal">∂</mml:mi>
<mml:mrow><mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
: Being strictly positive on ]0;1[, it indicates that <inline-formula><mml:math id="ImEquation14"><mml:msubsup><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>∗</mml:mo>
</mml:msubsup>
</mml:math>
</inline-formula>
is a decreasing function of <italic>p<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
. Essentially, the rarer the events in each surveillance level, the more resources need to be assigned to it.</p>
<p>A second observation is that the optimal <italic>proportion of resources</italic>
allocated to each surveillance level <inline-formula><mml:math id="ImEquation15"><mml:mo stretchy="false">(</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>ρ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>∗</mml:mo>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mspace width="thickmathspace"></mml:mspace>
<mml:mo>×</mml:mo>
<mml:msubsup><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>∗</mml:mo>
</mml:msubsup>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:mi>C</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
is independent of the total budget <italic>C</italic>
. A very simple approximation of this optimal proportion is obtained when the probabilities <italic>p<sub>i,k</sub>
</italic>
are much smaller than 1 and all of the costs <italic>c<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
are equal (i.e., <italic>c</italic>
<sub><italic>i,k</italic>
</sub>
= <italic>c</italic>
, ∀<italic>i</italic>
):
<disp-formula id="KWU213M3"><label>(3)</label>
<mml:math id="DmEquation3"><mml:mstyle><mml:mrow><mml:mfrac><mml:mrow><mml:msubsup><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>∗</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>≈</mml:mo>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mrow><mml:mn>1</mml:mn>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:msqrt><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
<mml:mrow><mml:msubsup><mml:mo movablelimits="false">∑</mml:mo>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:msubsup>
<mml:mrow><mml:mn>1</mml:mn>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:msqrt><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">∀</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mspace width="thinmathspace"></mml:mspace>
<mml:mo>,</mml:mo>
<mml:mspace width="thinmathspace"></mml:mspace>
<mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mstyle>
</mml:mstyle>
</mml:math>
</disp-formula>
where <italic>n</italic>
is the total number of recruited cases (<italic>n</italic>
= <italic>C/c</italic>
). We suggest that this may serve as a useful heuristic for rapid surveillance setup: All else being equal, the proportion of cases that needs to be recruited at each surveillance level is proportional to the inverse square root of the probability of progression at that level.</p>
<p>We denote as <inline-formula><mml:math id="ImEquation16"><mml:mrow><mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
<mml:msubsup><mml:mrow><mml:mover><mml:mrow><mml:mi mathvariant="normal">CFR</mml:mi>
</mml:mrow>
<mml:mo>ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>∗</mml:mo>
</mml:msubsup>
</mml:math>
</inline-formula>
the estimator of strategy <italic>k</italic>
when precision is optimum. Its SE is obtained by replacing {<italic>n<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
} with <inline-formula><mml:math id="ImEquation17"><mml:mo fence="false" stretchy="false">{</mml:mo>
<mml:msubsup><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>∗</mml:mo>
</mml:msubsup>
<mml:mo fence="false" stretchy="false">}</mml:mo>
</mml:math>
</inline-formula>
in equation 1:
<disp-formula id="KWU213M4"><label>(4)</label>
<mml:math id="DmEquation4"><mml:mtext>SE</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
<mml:msubsup><mml:mrow><mml:mover><mml:mrow><mml:mi mathvariant="normal">CFR</mml:mi>
</mml:mrow>
<mml:mo>ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>∗</mml:mo>
</mml:msubsup>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>≈</mml:mo>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mrow><mml:mtext>sCFR</mml:mtext>
</mml:mrow>
<mml:mrow><mml:msqrt><mml:mi>C</mml:mi>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:munderover><mml:mo movablelimits="false">∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:munderover>
<mml:mrow><mml:msqrt><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mfenced open="(" close=")"><mml:mrow><mml:mstyle><mml:mrow><mml:mfrac><mml:mn>1</mml:mn>
<mml:mrow><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mstyle>
</mml:mrow>
</mml:mfenced>
</mml:msqrt>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mstyle>
</mml:math>
</disp-formula>
</p>
<p>Equation 4 can be inverted to calculate the budget necessary to reach a desired SE, σ, with strategy <italic>k</italic>
:
<disp-formula id="KWU213M5"><label>(5)</label>
<mml:math id="DmEquation5"><mml:mrow><mml:mi>C</mml:mi>
<mml:mo>≈</mml:mo>
<mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo>
<mml:mrow><mml:mfrac><mml:mrow><mml:mtext>sCFR</mml:mtext>
</mml:mrow>
<mml:mi>σ</mml:mi>
</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>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:munderover>
<mml:mrow><mml:msqrt><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mfrac><mml:mn>1</mml:mn>
<mml:mrow><mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>In some cases, it may be desirable to consider “recruitment needs” without reference to costs—for example, if it can be assumed that recruitment costs are equivalent between surveillance systems. Estimate precision is then a simple function of the total number of recruited cases <italic>n</italic>
:
<disp-formula id="KWU213M6"><label>(6)</label>
<mml:math id="DmEquation6"><mml:mtext>SE</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
<mml:msubsup><mml:mrow><mml:mover><mml:mrow><mml:mi mathvariant="normal">CFR</mml:mi>
</mml:mrow>
<mml:mo>ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>∗</mml:mo>
</mml:msubsup>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>≈</mml:mo>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mrow><mml:mtext>sCFR</mml:mtext>
</mml:mrow>
<mml:mrow><mml:msqrt><mml:mi>n</mml:mi>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:munderover><mml:mo movablelimits="false">∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:munderover>
<mml:mrow><mml:msqrt><mml:mstyle><mml:mrow><mml:mfrac><mml:mn>1</mml:mn>
<mml:mrow><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mstyle>
</mml:msqrt>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mstyle>
</mml:math>
</disp-formula>
and the overall number of cases to recruit to reach a targeted SE, σ, is
<disp-formula id="KWU213M7"><label>(7)</label>
<mml:math id="DmEquation7"><mml:mrow><mml:mi>n</mml:mi>
<mml:mo>≈</mml:mo>
<mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo>
<mml:mrow><mml:mfrac><mml:mrow><mml:mtext>sCFR</mml:mtext>
</mml:mrow>
<mml:mi>σ</mml:mi>
</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>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:munderover>
<mml:mrow><mml:msqrt><mml:mrow><mml:mfrac><mml:mn>1</mml:mn>
<mml:mrow><mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
In <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 3</ext-link>
, we study the best allocation of resources 1) made available partway during an outbreak (<ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Figure 1) and 2</ext-link>
) given that some surveillance systems have fixed sample sizes (e.g., when routine surveillance data are reused for sCFR estimation).</p>
</sec>
<sec id="s2b"><title>Choosing the most precise pyramidal estimator</title>
<p>While comparisons of SEs or necessary budgets can be carried out with equations 1, 4, and 5, we provide hereafter rules of thumb for rapid comparison of estimation strategies in the special case where all surveillance systems have equal recruitment costs and each estimator has a minimal SE thanks to appropriate allocation of resources.</p>
<p>We first study whether precision can be improved by adding a new surveillance level—that is, by splitting the progression probability <italic>p<sub>j</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
into 2 probabilities <italic>p′</italic>
and <italic>p</italic>
″ such that <italic>p<sub>j</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
= <italic>p</italic>
′ × <italic>p</italic>
″. In <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 4</ext-link>
, we show that this split increases precision only if <italic>p<sub>j</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
< 1/9 (approximately 0.11) and if <italic>p</italic>
′ and <italic>p</italic>
″ satisfy
<disp-formula id="KWU213M8"><label>(8)</label>
<mml:math id="DmEquation8"><mml:mspace width="-3pc"></mml:mspace>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mn>3</mml:mn>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:msqrt><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mn>10</mml:mn>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mn>9</mml:mn>
<mml:msup><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:msqrt>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
</mml:mrow>
<mml:mo><</mml:mo>
<mml:mo fence="false" stretchy="false">{</mml:mo>
<mml:msup><mml:mi>p</mml:mi>
<mml:mrow><mml:mi mathvariant="normal">′</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:msup><mml:mi>p</mml:mi>
<mml:mrow><mml:mi mathvariant="normal">′</mml:mi>
<mml:mi mathvariant="normal">′</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo fence="false" stretchy="false">}</mml:mo>
<mml:mo><</mml:mo>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mn>3</mml:mn>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msqrt><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mn>10</mml:mn>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mn>9</mml:mn>
<mml:msup><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:msqrt>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mstyle>
</mml:mstyle>
</mml:math>
</disp-formula>
</p>
<p>The smaller the value of <italic>p<sub>j</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
, the wider this interval, so the less important this second condition is. For example, when <italic>p<sub>j</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
is the 2009 sCFR, all splits where <italic>p</italic>
′ and <italic>p</italic>
″ are picked between 0.00025 and 0.999 allow a gain in precision. This is illustrated in <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Figure 2</ext-link>
.</p>
<p>We then study which pair <italic>p</italic>
′ × <italic>p</italic>
″ yields the most precise estimate. We show that the more similar <italic>p</italic>
′ is to <italic>p</italic>
″, the more precise the estimator, with maximum precision being reached when <inline-formula><mml:math id="ImEquation18"><mml:msup><mml:mi>p</mml:mi>
<mml:mrow><mml:mi mathvariant="normal">′</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>=</mml:mo>
<mml:msup><mml:mi>p</mml:mi>
<mml:mrow><mml:mi mathvariant="normal">′</mml:mi>
<mml:mi mathvariant="normal">′</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>=</mml:mo>
<mml:msqrt><mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msqrt>
</mml:math>
</inline-formula>
(see <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 4</ext-link>
). These decision rules are summarized with a decision tree in <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Figure 3</ext-link>
.</p>
</sec>
<sec id="s2c"><title>Optimization in the presence of uncertainty</title>
<p>It is easy to calculate optimal sample sizes when all <italic>p<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
are known, but in practice, of course, <italic>p<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
are unknown. Yet informed guesses, denoted <inline-formula><mml:math id="ImEquation19"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">~</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
, supported by the literature or by preliminary surveys, can be used to calculate the sCFR expected value, denoted <inline-formula><mml:math id="ImEquation20"><mml:mrow><mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
<mml:mrow><mml:mover><mml:mrow><mml:mi mathvariant="normal">CFR</mml:mi>
</mml:mrow>
<mml:mo>~</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mo movablelimits="false">∏</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:msubsup>
<mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">~</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
, and the sample sizes, denoted <inline-formula><mml:math id="ImEquation21"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>n</mml:mi>
<mml:mo stretchy="false">~</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
, that optimize estimator precision:
<disp-formula id="KWU213M9"><label>(9)</label>
<mml:math id="DmEquation9"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>n</mml:mi>
<mml:mo stretchy="false">~</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>≈</mml:mo>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mi>C</mml:mi>
<mml:mrow><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mstyle><mml:mrow><mml:mfrac><mml:mrow><mml:msqrt><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">~</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:msqrt>
</mml:mrow>
<mml:mrow><mml:munderover><mml:mo movablelimits="false">∑</mml:mo>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:munderover>
<mml:mrow><mml:msqrt><mml:mrow><mml:msub><mml:mi>c</mml:mi>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mrow><mml:mo>/</mml:mo>
</mml:mrow>
<mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">~</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mrow><mml:mspace width="thinmathspace"></mml:mspace>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:msqrt>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">∀</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mspace width="thinmathspace"></mml:mspace>
<mml:mo>…</mml:mo>
<mml:mspace width="thinmathspace"></mml:mspace>
<mml:mo>,</mml:mo>
<mml:mspace width="thinmathspace"></mml:mspace>
<mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mstyle>
</mml:mstyle>
</mml:math>
</disp-formula>
</p>
<p>Through sensitivity analysis, developed in <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 5</ext-link>
(<ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Figures 4–6</ext-link>
), we examine how deviations of <inline-formula><mml:math id="ImEquation22"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">~</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
from the true value <italic>p<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
affect the decided allocation of resources and, as a consequence, decrease estimator precision.</p>
<p>We also assess the robustness of the estimator choice to the initial uncertainty: Sensitivity techniques are used to construct 8 anticipation scenarios reflecting uncertainty about the severity of a particular outbreak as it starts. Robust performance ranking of sCFR estimators throughout all 8 scenarios is sought.</p>
</sec>
</sec>
<sec sec-type="results" id="s3"><title>RESULTS</title>
<p>The analytical results derived in the Methods section are illustrated below with numerical examples, using 3 influenza pandemic scenarios with different severity levels (Table <xref ref-type="table" rid="KWU213TB1">1</xref>
):severe (“1918-like”; sCFR = 2%), intermediate (“1957-like”; sCFR = 0.2%), and mild (“2009-like”; sCFR =0.025%). We study the 4 estimators presented in Figure <xref ref-type="fig" rid="KWU213F1">1</xref>
. Table <xref ref-type="table" rid="KWU213TB2">2</xref>
shows the key epidemiologic concepts illustrated in each example. For simplicity of illustration, we assume equal recruitment costs per case at all surveillance levels of all estimation strategies; consequently, the same budget allows for recruitment of the same number of cases. Other numerical examples, based on cost sets described <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 6</ext-link>
(<ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Tables 1 and 2</ext-link>
), are presented in <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendices 7 and 8</ext-link>
(<ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Tables 3–6</ext-link>
).
<table-wrap id="KWU213TB1" position="float"><label>Table 1.</label>
<caption><p>Severity Parameters Used for Simulating 3 Influenza Pandemic Scenarios With Different Levels of Severity</p>
</caption>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="left" span="1"></col>
<col align="char" char="." span="1"></col>
<col align="char" char="." span="1"></col>
<col align="char" char="." span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="2" colspan="1">Parameter</th>
<th align="center" rowspan="2" colspan="1">Description</th>
<th align="center" colspan="3" rowspan="1">Severity of Influenza Pandemic<hr></hr>
</th>
</tr>
<tr><th align="left" rowspan="1" colspan="1">Severe <break></break>
(1918-Like)</th>
<th align="center" rowspan="1" colspan="1">Intermediate (1957-Like)</th>
<th align="center" rowspan="1" colspan="1">Mild (2009-Like)</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1"><italic>p<sub>D</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
= sCFR</td>
<td rowspan="1" colspan="1">Probability of death following symptoms</td>
<td rowspan="1" colspan="1">0.02040 (<xref rid="KWU213C4" ref-type="bibr">4</xref>
)<sup>a</sup>
</td>
<td rowspan="1" colspan="1">0.00020 (<xref rid="KWU213C32" ref-type="bibr">32</xref>
)</td>
<td rowspan="1" colspan="1">0.00025 (<xref rid="KWU213C14" ref-type="bibr">14</xref>
)</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
</td>
<td rowspan="1" colspan="1">Probability of medical attention following symptoms</td>
<td rowspan="1" colspan="1">0.4000</td>
<td rowspan="1" colspan="1">0.2000</td>
<td rowspan="1" colspan="1">0.3500</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>p<sub>H</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
</td>
<td rowspan="1" colspan="1">Probability of hospitalization following medical attention</td>
<td rowspan="1" colspan="1">0.3500</td>
<td rowspan="1" colspan="1">0.1000</td>
<td rowspan="1" colspan="1">0.0157</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>p<sub>D</sub>
</italic>
<sub>|<italic>H</italic>
</sub>
</td>
<td rowspan="1" colspan="1">Probability of death following hospitalization</td>
<td rowspan="1" colspan="1">0.1457</td>
<td rowspan="1" colspan="1">0.1000</td>
<td rowspan="1" colspan="1">0.0455</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>p<sub>D</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
<italic>= p<sub>D</sub>
</italic>
<sub>|<italic>H</italic>
</sub>
<italic>× p<sub>H</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
</td>
<td rowspan="1" colspan="1">Probability of death following medical attention</td>
<td rowspan="1" colspan="1">0.0510</td>
<td rowspan="1" colspan="1">0.0100</td>
<td rowspan="1" colspan="1">0.0007</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
<italic>= p<sub>H</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
<italic>× p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
</td>
<td rowspan="1" colspan="1">Probability of hospitalization following symptoms</td>
<td rowspan="1" colspan="1">0.1400</td>
<td rowspan="1" colspan="1">0.0200</td>
<td rowspan="1" colspan="1">0.0055 (<xref rid="KWU213C14" ref-type="bibr">14</xref>
)</td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="TN1"><p>Abbreviation: sCFR, symptomatic case fatality ratio.</p>
</fn>
<fn id="TN2"><p><sup>a</sup>
For data found in the literature, source references are given in parentheses; the other figures are assumed.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="KWU213TB2" position="float"><label>Table 2.</label>
<caption><p>Key Concepts Illustrated in Numerical Examples Corresponding to 3 Influenza Pandemic Scenarios with Different Levels of Severity</p>
</caption>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="left" span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="1" colspan="1">Numerical Example</th>
<th align="center" rowspan="1" colspan="1">Methodological Concept(s)</th>
</tr>
</thead>
<tbody><tr><td rowspan="4" colspan="1">Precision of surveillance strategies' pyramidal <break></break>
estimators</td>
<td rowspan="1" colspan="1">The optimal allocation of resources for any pyramidal estimator can be calculated with equation 2.</td>
</tr>
<tr><td rowspan="1" colspan="1">The standard error of pyramidal sCFR estimators under optimal resource allocation can be calculated with equation 4.</td>
</tr>
<tr><td rowspan="1" colspan="1">The precision of sCFR estimators increases when progression probabilities inferior to 0.11<sup>a</sup>
are split in two.</td>
</tr>
<tr><td rowspan="1" colspan="1">The precision of sCFR estimators increases when multiplied progression probabilities are close to one another<sup>a</sup>
.</td>
</tr>
<tr><td rowspan="3" colspan="1">Necessary budget</td>
<td rowspan="1" colspan="1">The minimal necessary budget for any pyramidal estimator to reach a desired precision can be calculated with equation 5.</td>
</tr>
<tr><td rowspan="1" colspan="1">The minimal necessary budget for the single-level estimator increases dramatically as the sCFR decreases.</td>
</tr>
<tr><td rowspan="1" colspan="1">Using pyramidal estimators allows a great precision gain when the sCFR is small.</td>
</tr>
<tr><td rowspan="1" colspan="1">Robustness to initial uncertainty</td>
<td rowspan="1" colspan="1">See “Optimization in the presence of uncertainty” in the Methods section of the text.</td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="TN3"><p>Abbreviation: sCFR, symptomatic case fatality ratio.</p>
</fn>
<fn id="TN4"><p><sup>a</sup>
When recruitment costs at the concerned levels are equal.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</p>
<sec id="s3a"><title>Precision of pyramidal estimators</title>
<p>Table <xref ref-type="table" rid="KWU213TB3">3</xref>
presents the SE of the 4 sCFR estimators, obtained with equation 6, considering a budget allowing recruitment of 10,000 cases. The best possible precision is obtained for each estimator by appropriate allocation of resources between surveillance levels, calculated with equation 2. Figure <xref ref-type="fig" rid="KWU213F2">2</xref>
shows the 95% prediction intervals of the sCFR estimators, based on 10,000 recruited cases. Figure <xref ref-type="fig" rid="KWU213F3">3</xref>
presents SE ratios between the 3 pyramidal estimators and the single-level estimator. Precision usually increases with the number of surveillance levels in an estimation strategy, more so for small sCFRs. However, in our simulation, this was not systematic: In the severe scenario, a 2-level estimator <inline-formula><mml:math id="ImEquation23"><mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
was slightly more precise than the 3-level one <inline-formula><mml:math id="ImEquation24"><mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mspace width="thickmathspace"></mml:mspace>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>M</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
. This was predictable from the analytical results presented in the Methods section: The progression probability <italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
is above 0.11 (<italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
= 0.14), making it inefficient to further split it. The estimator <inline-formula><mml:math id="ImEquation25"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
was always more precise than <inline-formula><mml:math id="ImEquation26"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>M</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
, with the same number of surveillance levels. This was also predictable: Estimators that multiply close progression probabilities are the most precise, and in all 3 scenarios <italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
is closer to <italic>p<sub>D</sub>
</italic>
<sub>|<italic>H</italic>
</sub>
than <italic>p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
is to <italic>p<sub>D</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
(see Table <xref ref-type="table" rid="KWU213TB1">1</xref>
).
<table-wrap id="KWU213TB3" position="float"><label>Table 3.</label>
<caption><p>Minimal Standard Errors of Symptomatic Case Fatality Ratio Estimators (Based on 10,000 Recruited Cases) When Recruitment Costs Are Equal at All Surveillance Levels<sup>a</sup>
</p>
</caption>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="left" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" char="." span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" char="." span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" char="." span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="3" colspan="1">Estimator</th>
<th align="center" rowspan="3" colspan="1">Level<sup>b</sup>
</th>
<th align="center" rowspan="3" colspan="1">Event</th>
<th align="center" colspan="9" rowspan="1">Severity of Influenza Pandemic<hr></hr>
</th>
</tr>
<tr><th align="center" colspan="3" rowspan="1">Severe (1918-Like): sCFR = 2.04 × 10<sup>−2</sup>
<hr></hr>
</th>
<th align="center" colspan="3" rowspan="1">Intermediate (1957-Like): sCFR = 2 × 10<sup>−3</sup>
<hr></hr>
</th>
<th align="center" colspan="3" rowspan="1">Mild (2009-Like): sCFR = 2.5 × 10<sup>−4</sup>
<hr></hr>
</th>
</tr>
<tr><th align="left" rowspan="1" colspan="1">Optimal Sample Size, no.</th>
<th align="center" rowspan="1" colspan="1">Expected No. of Events</th>
<th align="center" rowspan="1" colspan="1">SE (×10<sup>−3</sup>
)</th>
<th align="center" rowspan="1" colspan="1">Optimal Sample Size, no.</th>
<th align="center" rowspan="1" colspan="1">Expected No. of Events</th>
<th align="center" rowspan="1" colspan="1">SE (×10<sup>−4</sup>
)</th>
<th align="center" rowspan="1" colspan="1">Optimal Sample Size, no.</th>
<th align="center" rowspan="1" colspan="1">Expected No. of Events</th>
<th align="center" rowspan="1" colspan="1">SE (×10<sup>−5</sup>
)</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1"><inline-formula><mml:math id="ImEquation27"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><italic>S</italic>
</td>
<td rowspan="1" colspan="1">Death</td>
<td rowspan="1" colspan="1">10,000</td>
<td rowspan="1" colspan="1">204</td>
<td rowspan="1" colspan="1">1.41</td>
<td rowspan="1" colspan="1">10,000</td>
<td rowspan="1" colspan="1">20</td>
<td rowspan="1" colspan="1">4.47</td>
<td rowspan="1" colspan="1">10,000</td>
<td rowspan="1" colspan="1">3</td>
<td rowspan="1" colspan="1">15.81</td>
</tr>
<tr><td rowspan="2" colspan="1"><inline-formula><mml:math id="ImEquation28"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>M</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><italic>S</italic>
</td>
<td rowspan="1" colspan="1">Medical attention</td>
<td rowspan="1" colspan="1">2,211</td>
<td rowspan="1" colspan="1">885</td>
<td rowspan="1" colspan="1">1.13</td>
<td rowspan="1" colspan="1">1,674</td>
<td rowspan="1" colspan="1">335</td>
<td rowspan="1" colspan="1">2.39</td>
<td rowspan="1" colspan="1">352</td>
<td rowspan="1" colspan="1">123</td>
<td rowspan="1" colspan="1">9.69</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>M</italic>
</td>
<td rowspan="1" colspan="1">Death</td>
<td rowspan="1" colspan="1">7,789</td>
<td rowspan="1" colspan="1">397</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">8,326</td>
<td rowspan="1" colspan="1">83</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">9,648</td>
<td rowspan="1" colspan="1">7</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="2" colspan="1"><inline-formula><mml:math id="ImEquation29"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><italic>S</italic>
</td>
<td rowspan="1" colspan="1">Hospitalization</td>
<td rowspan="1" colspan="1">5,058</td>
<td rowspan="1" colspan="1">708</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">7,000</td>
<td rowspan="1" colspan="1">140</td>
<td rowspan="1" colspan="1">2.00</td>
<td rowspan="1" colspan="1">7,459</td>
<td rowspan="1" colspan="1">41</td>
<td rowspan="1" colspan="1">4.51</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>H</italic>
</td>
<td rowspan="1" colspan="1">Death</td>
<td rowspan="1" colspan="1">4,942</td>
<td rowspan="1" colspan="1">720</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">3,000</td>
<td rowspan="1" colspan="1">300</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">2,541</td>
<td rowspan="1" colspan="1">116</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="3" colspan="1"><inline-formula><mml:math id="ImEquation30"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>M</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><italic>S</italic>
</td>
<td rowspan="1" colspan="1">Medical attention</td>
<td rowspan="1" colspan="1">2,445</td>
<td rowspan="1" colspan="1">978</td>
<td rowspan="1" colspan="1">1.02</td>
<td rowspan="1" colspan="1">2,500</td>
<td rowspan="1" colspan="1">500</td>
<td rowspan="1" colspan="1">1.60</td>
<td rowspan="1" colspan="1">983</td>
<td rowspan="1" colspan="1">344</td>
<td rowspan="1" colspan="1">3.47</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>M</italic>
</td>
<td rowspan="1" colspan="1">Hospitalization</td>
<td rowspan="1" colspan="1">2,721</td>
<td rowspan="1" colspan="1">952</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">3,750</td>
<td rowspan="1" colspan="1">375</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">5,713</td>
<td rowspan="1" colspan="1">90</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>H</italic>
</td>
<td rowspan="1" colspan="1">Death</td>
<td rowspan="1" colspan="1">4,834</td>
<td rowspan="1" colspan="1">704</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">3,750</td>
<td rowspan="1" colspan="1">375</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">3,304</td>
<td rowspan="1" colspan="1">150</td>
<td rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="TN5"><p>Abbreviations: sCFR, symptomatic case fatality ratio; SE, standard error.</p>
</fn>
<fn id="TN6"><p><sup>a</sup>
The optimal sample size at each surveillance level is obtained with equation 2, ensuring the minimum SE for each estimator. Expected numbers of events are calculated as sample size × <italic>p<sub>i</sub>
</italic>
<sub>|<italic>j</italic>
</sub>
.</p>
</fn>
<fn id="TN7"><p><sup>b</sup>
<italic>H</italic>
, hospitalized cases; <italic>M</italic>
, medically attended cases; <italic>S</italic>
, symptomatic cases.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<fig id="KWU213F2" position="float"><label>Figure 2.</label>
<caption><p>Ninety-five percent prediction intervals of symptomatic case fatality ratio (sCFR) estimators based on 10,000 recruited cases. Pyramidal sCFR estimators are obtained by multiplying estimates of <italic>p<sub>i</sub>
</italic>
<sub>|<italic>j</italic>
</sub>
, the probabilities of progressing to severity level <italic>i</italic>
for cases in level <italic>j</italic>
, which are based on 4 severity levels: symptomatic cases (<italic>S</italic>
), medically attended cases (<italic>M</italic>
), hospitalized cases (<italic>H</italic>
), and dead cases (<italic>D</italic>
). The 95% prediction intervals (T-shaped bars) are the intervals surrounding the true sCFRs (bullet points) in which future sCFR estimates will fall, with a probability of 95%, if the assumptions about progression probabilities are correct. Narrow prediction intervals indicate good precision. A) Severe (1918-like) influenza pandemic scenario; B) intermediate (1957-like) severity scenario; C) mild (2009-like) severity scenario. Maximal precision is ensured for each estimator through optimal resource allocation between surveillance levels. Recruitment costs are assumed to be equal at all surveillance levels of all estimation strategies.</p>
</caption>
<graphic xlink:href="kwu21302"></graphic>
</fig>
<fig id="KWU213F3" position="float"><label>Figure 3.</label>
<caption><p>Relative precision of symptomatic case fatality ratio (sCFR) pyramidal estimators versus the single-level estimator <inline-formula><mml:math id="ImEquation31"><mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
, given similar budgets. Relative precision is measured with the standard error (SE) ratio, obtained by dividing the minimal SE of each pyramidal estimator by that of the single-level estimator. Recruitment costs are assumed to be equal at all surveillance levels of all estimation strategies. It can be noted from equation 4 that the SE ratio of 2 estimators that have the same budget <italic>C</italic>
is independent of <italic>C</italic>
.</p>
</caption>
<graphic xlink:href="kwu21303"></graphic>
</fig>
</p>
</sec>
<sec id="s3b"><title>Necessary budget</title>
<p>In this second illustration, we compare the necessary budgets, in terms of sample size, for all estimators to achieve the same level of precision. Again, we assume that the best possible precision is obtained for each estimator by appropriate allocation of resources. In Table <xref ref-type="table" rid="KWU213TB4">4</xref>
, the number of cases needed to obtain a coefficient of variation (coefficient of variation = SE/sCFR) of 0.5 is calculated with equation 7. Figure <xref ref-type="fig" rid="KWU213F4">4</xref>
presents the necessary sample size of the single-level estimator for different targeted precisions. Figure <xref ref-type="fig" rid="KWU213F5">5</xref>
presents the relative sample sizes of the 3 pyramidal estimators. In short, choosing the optimal estimator substantially reduces recruitment efforts, particularly when the sCFR is low. This matters, because it is when the sCFR is low that such reduction is the most welcome, since necessary sample sizes for the single-level estimator are then prohibitive (Figure <xref ref-type="fig" rid="KWU213F4">4</xref>
). For example, compared with the severe scenario, the size of a case series would need to be 10- and 83-fold larger for the intermediate and mild severity scenarios, respectively, for a single-level sCFR estimator to have a similar coefficient of variation. However, those differences are reduced when using the optimal pyramidal estimator in each scenario: One would then need 2.67- and 8-fold more cases for the intermediate and mild severity scenarios, respectively.
<table-wrap id="KWU213TB4" position="float"><label>Table 4.</label>
<caption><p>Sample Sizes Necessary to Obtain a Coefficient of Variation of 0.5 for Symptomatic Case Fatality Ratio Estimators When Recruitment Costs Are Equal at All Surveillance Levels<sup>a</sup>
</p>
</caption>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="char" span="1"></col>
<col align="left" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
<col align="char" span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="3" colspan="1">Estimator</th>
<th align="center" rowspan="3" colspan="1">Level<sup>b</sup>
</th>
<th align="center" rowspan="3" colspan="1">Event</th>
<th align="center" colspan="9" rowspan="1">Severity of Influenza Pandemic<hr></hr>
</th>
</tr>
<tr><th align="center" colspan="3" rowspan="1">Severe (1918-Like): sCFR = 2.04 × 10<sup>−2</sup>
<hr></hr>
</th>
<th align="center" colspan="3" rowspan="1">Intermediate (1957-Like): sCFR = 2 × 10<sup>−3</sup>
<hr></hr>
</th>
<th align="center" colspan="3" rowspan="1">Mild (2009-Like): sCFR = 2.5 × 10<sup>−4</sup>
<hr></hr>
</th>
</tr>
<tr><th align="left" rowspan="1" colspan="1">Optimal Sample Size, no.</th>
<th align="center" rowspan="1" colspan="1">Expected No. of Events</th>
<th align="center" rowspan="1" colspan="1">Cumulated Sample Size, no.</th>
<th align="center" rowspan="1" colspan="1">Optimal Sample Size, no.</th>
<th align="center" rowspan="1" colspan="1">Expected No. of Events</th>
<th align="center" rowspan="1" colspan="1">Cumulated Sample Size, no.</th>
<th align="center" rowspan="1" colspan="1">Optimal Sample Size, no.</th>
<th align="center" rowspan="1" colspan="1">Expected No. of Events</th>
<th align="center" rowspan="1" colspan="1">Cumulated Sample Size, no.</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1"><inline-formula><mml:math id="ImEquation32"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><italic>S</italic>
</td>
<td rowspan="1" colspan="1">Death</td>
<td rowspan="1" colspan="1">192</td>
<td rowspan="1" colspan="1">4</td>
<td rowspan="1" colspan="1">192</td>
<td rowspan="1" colspan="1">1,996</td>
<td rowspan="1" colspan="1">4</td>
<td rowspan="1" colspan="1">1,996</td>
<td rowspan="1" colspan="1">15,996</td>
<td rowspan="1" colspan="1">4</td>
<td rowspan="1" colspan="1">15,996</td>
</tr>
<tr><td rowspan="2" colspan="1"><inline-formula><mml:math id="ImEquation33"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>M</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><italic>S</italic>
</td>
<td rowspan="1" colspan="1">Medical attention</td>
<td rowspan="1" colspan="1">27</td>
<td rowspan="1" colspan="1">11</td>
<td rowspan="1" colspan="1">123</td>
<td rowspan="1" colspan="1">96</td>
<td rowspan="1" colspan="1">19</td>
<td rowspan="1" colspan="1">572</td>
<td rowspan="1" colspan="1">211</td>
<td rowspan="1" colspan="1">74</td>
<td rowspan="1" colspan="1">6,011</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>M</italic>
</td>
<td rowspan="1" colspan="1">Death</td>
<td rowspan="1" colspan="1">96</td>
<td rowspan="1" colspan="1">5</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">476</td>
<td rowspan="1" colspan="1">5</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">5,800</td>
<td rowspan="1" colspan="1">4</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="2" colspan="1"><inline-formula><mml:math id="ImEquation34"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><italic>S</italic>
</td>
<td rowspan="1" colspan="1">Hospitalization</td>
<td rowspan="1" colspan="1">49</td>
<td rowspan="1" colspan="1">7</td>
<td rowspan="1" colspan="1">96</td>
<td rowspan="1" colspan="1">280</td>
<td rowspan="1" colspan="1">6</td>
<td rowspan="1" colspan="1">400</td>
<td rowspan="1" colspan="1">972</td>
<td rowspan="1" colspan="1">5</td>
<td rowspan="1" colspan="1">1,303</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>H</italic>
</td>
<td rowspan="1" colspan="1">Death</td>
<td rowspan="1" colspan="1">47</td>
<td rowspan="1" colspan="1">7</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">120</td>
<td rowspan="1" colspan="1">12</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">331</td>
<td rowspan="1" colspan="1">15</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="3" colspan="1"><inline-formula><mml:math id="ImEquation35"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>M</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td rowspan="1" colspan="1"><italic>S</italic>
</td>
<td rowspan="1" colspan="1">Medical attention</td>
<td rowspan="1" colspan="1">25</td>
<td rowspan="1" colspan="1">10</td>
<td rowspan="1" colspan="1">101</td>
<td rowspan="1" colspan="1">64</td>
<td rowspan="1" colspan="1">13</td>
<td rowspan="1" colspan="1">256</td>
<td rowspan="1" colspan="1">76</td>
<td rowspan="1" colspan="1">26</td>
<td rowspan="1" colspan="1">769</td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>M</italic>
</td>
<td rowspan="1" colspan="1">Hospitalization</td>
<td rowspan="1" colspan="1">27</td>
<td rowspan="1" colspan="1">10</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">96</td>
<td rowspan="1" colspan="1">10</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">439</td>
<td rowspan="1" colspan="1">7</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"><italic>H</italic>
</td>
<td rowspan="1" colspan="1">Death</td>
<td rowspan="1" colspan="1">49</td>
<td rowspan="1" colspan="1">7</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">96</td>
<td rowspan="1" colspan="1">10</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">254</td>
<td rowspan="1" colspan="1">12</td>
<td rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="TN8"><p>Abbreviation: sCFR, symptomatic case fatality ratio.</p>
</fn>
<fn id="TN9"><p><sup>a</sup>
Corresponding standard errors are 1.02 × 10<sup>−2</sup>
, 1.00 × 10<sup>−3</sup>
, and 1.25 × 10<sup>−4</sup>
for the severe, intermediate, and mild scenarios, respectively. Optimal sample sizes are obtained by optimally allocating the total number of recruited cases (cumulated sample size) between surveillance levels using equation 2. Expected numbers of events are calculated as sample size × <italic>p<sub>i</sub>
</italic>
<sub>|<italic>j</italic>
</sub>
.</p>
</fn>
<fn id="TN10"><p><sup>b</sup>
<italic>H</italic>
, hospitalized cases; <italic>M</italic>
, medically attended cases; <italic>S</italic>
, symptomatic cases.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<fig id="KWU213F4" position="float"><label>Figure 4.</label>
<caption><p>Sample size needed to reach different levels of precision with a single-level estimator of the symptomatic case fatality ratio. Low coefficients of variation indicate good precision. A single-level estimate may be obtained, for example, by counting fatalities in a symptomatic case series.</p>
</caption>
<graphic xlink:href="kwu21304"></graphic>
</fig>
<fig id="KWU213F5" position="float"><label>Figure 5.</label>
<caption><p>Relative necessary sample size of symptomatic case fatality ratio (sCFR) pyramidal estimators, in comparison with the single-level estimator <inline-formula><mml:math id="ImEquation36"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
. Sample size ratios divide the sample size needed to obtain a 0.5 coefficient of variation with pyramidal estimators by that of the single-level estimator. Maximal precision is ensured for each estimator through optimal resource allocation between surveillance levels. Recruitment costs are assumed to be equal at all surveillance levels of all strategies.</p>
</caption>
<graphic xlink:href="kwu21305"></graphic>
</fig>
</p>
</sec>
<sec id="s3c"><title>Robustness to initial uncertainty</title>
<p>True optimization requires knowledge of the probabilities we want to estimate, which is obviously circular. To be useful, optimization needs to be robust to the use of approximate values available early during an emerging infectious disease outbreak. We analyze how uncertainty in the preliminary estimates of <italic>p<sub>i</sub>
</italic>
<sub>,<italic>k</italic>
</sub>
may affect the optimal allocation of resources between surveillance levels. As an illustration, we study the estimator <inline-formula><mml:math id="ImEquation37"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>D</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>×</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>p</mml:mi>
<mml:mo stretchy="false">ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
: We assume that the probability of death following hospitalization (<italic>p<sub>D</sub>
</italic>
<sub>|<italic>H</italic>
</sub>
) is known but the probability of hospitalization following symptoms (<italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
) is initially uncertain. We allow the preliminary estimate of <italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
to vary between 0.0001 and 0.99. For each preliminary estimate of <italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
, we calculate the optimal allocation of resources between community-based surveillance systems and hospital-based surveillance systems. The more the preliminary estimate of <italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
deviates from the true value, the less optimal is the resource allocation and the less precise is the sCFR estimator. However, we find numerically that estimator precision is quite robust even to errors of several orders of magnitude in the preliminary estimate of <italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
(<ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Figure 4</ext-link>
).</p>
<p>For example, optimal precision for the mild scenario is obtained when 75% of cases are recruited in the community and 25% are recruited in hospitals. If the preliminary estimate of <italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
were 0.0003 (instead of 0.0055), the initial “best guess” of optimal resource allocation would be 92% of cases recruited in the community and 8% of cases recruited in the hospital, yielding an SE 1.21-fold larger than optimal. Even with an initial estimate of <italic>p<sub>H</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
of 0.1, giving us a best-guess resource allocation ratio of 40%–60%, the SE is only 1.23-fold larger than optimal. Similar results are observed for the intermediate and severe severity scenarios and with the other 2-level estimator.</p>
<p>This implies that optimization of resources between surveillance levels can be based on an informed guess about progression probabilities, since only if this guess is wrong by several orders of magnitude will precision decrease substantially. This also means that gains in precision obtained by optimal resource allocation with a particular pyramidal estimator are much less substantial than gains in precision obtained from a good choice of the estimator.</p>
<p>Given this conclusion, we analyze whether the optimal estimator can be identified in the presence of uncertainty at the start of an outbreak. For each studied pandemic, we combine the initial uncertainty bounds of <italic>p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
, <italic>p<sub>H</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
, and <italic>p<sub>D</sub>
</italic>
<sub>|<italic>H</italic>
</sub>
to obtain 8 anticipation scenarios. A consistent order of estimator precisions across all anticipation scenarios would strongly support the choice of the best one. In Figure <xref ref-type="fig" rid="KWU213F6">6</xref>
, we present this analysis for the mild severity scenario assuming uncertainty ranges of 0.2–0.5, 0.005–0.03, and 0.01–0.1 for <italic>p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
, <italic>p<sub>H</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
, and <italic>p<sub>D</sub>
</italic>
<sub>|<italic>H</italic>
</sub>
, respectively (for the other scenarios, see <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Figures 5 and 6</ext-link>
). For all 8 anticipation scenarios, the 3-level estimator is more precise than the others, with the single-level estimator always being the worst.
<fig id="KWU213F6" position="float"><label>Figure 6.</label>
<caption><p>Expected minimal standard error (SE) of symptomatic case fatality ratio (sCFR) estimators in the presence of uncertainty around severity parameters, in a mild (2009-like) influenza pandemic scenario. A recruitment capacity of 10,000 cases is assumed. Parameters <italic>p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
, <italic>p<sub>H</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
, and <italic>p<sub>D</sub>
</italic>
<sub>|<italic>H</italic>
</sub>
(the probabilities of medical attention upon symptoms, hospitalization upon medical attention, and death upon hospitalization, respectively) are assumed to be uncertain at the start of the pandemic. The true values of <italic>p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
, <italic>p<sub>H</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
, and <italic>p<sub>D</sub>
</italic>
<sub>|<italic>H</italic>
</sub>
are 0.3500, 0.0157, and 0.0455, respectively, with uncertainty bounds of 0.2–0.5, 0.005–0.03, and 0.01–0.1, respectively. The minimal SE of each sCFR estimator is plotted for the true pandemic scenario (bold line) and for the 8 anticipation scenarios constructed by combining the uncertainty bounds. Those minimal SEs are obtained by optimally allocating the 10,000 recruited cases between surveillance levels and are calculated using equation 6.</p>
</caption>
<graphic xlink:href="kwu21306"></graphic>
</fig>
</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4"><title>DISCUSSION</title>
<p>Producing precise and accurate estimates of case fatality ratios early during an emerging infectious disease outbreak is important for public health decision-making. We have shown here that statistical planning can help investigators find an estimation strategy that minimizes costs without sacrificing precision. We focused on the sCFR as the case fatality ratio of most interest; however, the methods presented here can be used to optimize the precision of any severity estimate that decomposes into a product of probabilities, like the infectious case fatality ratio (making use of serology data) (<xref rid="KWU213C16" ref-type="bibr">16</xref>
, <xref rid="KWU213C17" ref-type="bibr">17</xref>
).</p>
<p>We showed that the precision obtainable with a given budget depends greatly on the surveillance levels at which the estimation strategy is focused, and that pyramidal estimators with 3 surveillance levels can dramatically decrease the necessary budget in comparison with follow-up of a series of symptomatic cases. To a lesser extent, precision also depends on the allocation of resources between the surveillance levels of a given strategy.</p>
<p>As a result, we propose rules of thumb for finding the most precise estimator. First, precision generally increases with the number of surveillance levels, particularly when severity is low: Where feasible, surveillance levels should be further decomposed until the probability of progression to the next pyramid level is about 11% (at which point the gain in precision stops). Second, precision generally increases when the probabilities of progression multiplied together are similar. Finally, early identification of the optimal estimator is possible: Comparing the precision of all contemplated estimators over a plausible range of parameters guaranties a robust choice. Besides, as data are collected, severity parameters can be updated and surveillance strategies adapted.</p>
<p>Which estimator is ultimately preferred also depends on other criteria, such as the availability of data and the feasibility, timeliness, and quality of reporting, among others. Scaling up particular levels of a surveillance strategy might be difficult in real time; the limited precision gained by shuffling resources around surveillance levels must be contrasted with the associated inconvenience. Scaling up hospital-based surveillance or community case detection might be the most feasible option, the latter via outbreak investigation or Web or telephone cohorts. At the general practitioner level, one could imagine a system in which an aggregate number of cases is periodically reported, with outcome details on a proportion of these cases, and which could be scaled up or down depending on available data collection resources. When scaling up is done by enrolling more data providers, it is important to avoid introducing biases due to changes in the composition of the populations they serve (e.g., demographically). Finally, resource limitations might not be financial per se but might involve a limited number of staff; for example, suitably trained staff might be used to perform community contact tracing or might be sent to hospitals to extract data from patients’ case notes. This would be an example of reallocating resources between levels of a strategy.</p>
<p>Our study falls into the framework of <italic>optimal designs,</italic>
which aim at estimating parameters without bias and with minimum variance (<xref rid="KWU213C18" ref-type="bibr">18</xref>
); we have focused herein on the minimum variance issue. Few optimal design studies concern disease surveillance: Most have been simulation studies of surveillance protocols for animal diseases (<xref rid="KWU213C19" ref-type="bibr">19</xref>
–<xref rid="KWU213C23" ref-type="bibr">23</xref>
) and optimization studies of general practitioner recruitment for more precise estimation of influenza incidence in the community (<xref rid="KWU213C24" ref-type="bibr">24</xref>
–<xref rid="KWU213C27" ref-type="bibr">27</xref>
). More recently, Ejima et al. (<xref rid="KWU213C28" ref-type="bibr">28</xref>
) used simulations to assess the delay before precise estimation of the case fatality ratio in various severity scenarios. We believe that the simple examples presented here provide pedagogical insight into how statistical planning may help improve the precision of sCFR estimates. However, these simple analyses had limitations.</p>
<p>First, we studied only precision (minimum variance), not accuracy (absence of bias). In particular, pyramidal estimators can be biased if some underlying assumptions are not satisfied—for example, in Figure <xref ref-type="fig" rid="KWU213F1">1</xref>
, if <italic>D</italic>
⊂ <italic>H</italic>
⊂ <italic>M</italic>
⊂ <italic>S</italic>
is not true. This is the case if a large proportion of deaths occur outside the hospital or if symptomatic cases are hospitalized without previous medical attention. In the 2009 A/H1N1 pandemic, the reported proportion of deaths taking place outside of hospitals was 17% in the United States (<xref rid="KWU213C10" ref-type="bibr">10</xref>
), so it may be important to monitor death certification from primary care settings as well as hospitals, as was done in England (<xref rid="KWU213C11" ref-type="bibr">11</xref>
), and use them in a comprehensive evidence synthesis approach: For example, Presanis et al. (<xref rid="KWU213C10" ref-type="bibr">10</xref>
) estimated an overall sCFR by adding the sCFRs of hospitalized patients and nonhospitalized patients seen in general practice, both estimated with pyramidal estimators.</p>
<p>More generally, given how important severity estimates are early on in an outbreak and since, unfortunately, there is currently no real alternative to the pyramidal approach in a context of mild-to-intermediate sCFR, in the future it will be important to develop studies that can be conducted in parallel with pyramidal severity estimation to assess the validity and bias in the estimates. In particular, such studies should allow ascertainment of the likelihood of effectively going through each of the pyramids’ steps. For example, where a sample of symptomatic cases is used to estimate the probability of consulting a general practitioner upon development of symptoms, a subset can be used to also estimate the probability of hospitalization or death without a visit to a general practitioner. A random sample of death certificates can also be studied to identify influenza cases, reconstruct the patients’ health-care history, and determine the proportion that would not satisfy the pyramidal assumptions.</p>
<p>Biased pyramidal sCFR estimates can also result from reporting biases and inconsistent case ascertainment between surveillance levels. In this work, for simplicity's sake, we assumed that cases detected at any 1 level were representative and that the ascertainment of severity outcomes was complete. Reich et al. (<xref rid="KWU213C29" ref-type="bibr">29</xref>
) discussed in detail the impact of reporting rates on estimation of the case fatality ratio. As an illustration of accuracy issues in pyramidal approaches, Presanis et al. (<xref rid="KWU213C10" ref-type="bibr">10</xref>
) obtained 2 very different estimates of the 2009 A/H1N1 sCFR (0.048% and 0.007%) by using 2 pyramidal estimators, each one relying on a different set of reasonable assumptions. Each of these estimates had very tight credible intervals (i.e., good precision), and the 2 estimates were not consistent with each other.</p>
<p>This suggests a need for an analytical study of the accuracy of pyramidal estimators, which is outside the topic of the present paper. Garske et al. (<xref rid="KWU213C5" ref-type="bibr">5</xref>
) described some of the mechanisms resulting in biases in sCFR pyramidal estimates during the 2009 pandemic. Accuracy will likely need to be ascertained on a case-by-case basis, and it may depend on organizational and cultural norms surrounding health-care utilization.</p>
<p>Second, we estimated progression probabilities with simple ratios so as to find an analytical solution to our optimization problem. These estimates can be significantly biased when the incidence of infection increases quickly and the delay between severity outcomes is long<italic>.</italic>
Other authors have proposed estimation methods that account for right-truncation of severity events, which should be used in practice (<xref rid="KWU213C5" ref-type="bibr">5</xref>
, <xref rid="KWU213C6" ref-type="bibr">6</xref>
, <xref rid="KWU213C29" ref-type="bibr">29</xref>
, <xref rid="KWU213C30" ref-type="bibr">30</xref>
).</p>
<p>Third, in Table <xref ref-type="table" rid="KWU213TB4">4</xref>
we compared strategies in terms of budget, yet the most economical strategy might not be the timeliest. For example, to estimate a 1957-like sCFR with a coefficient of variation of 0.5, one must “wait” either until 4 deaths are reported within a case series of 1,996 symptomatic cases or until 10 deaths are reported within a 3-level estimation strategy of 256 cases overall. Thus, while the latter strategy is the most economical (256 cases to recruit vs. 1,996), one has to wait for more fatal events to be reported. If time to death is highly variable, the delay before precise sCFR estimation may increase in pyramidal approaches, offsetting their economic superiority.</p>
<p>Fourth, we assumed fixed costs per recruitment at all surveillance levels. In practice, these costs can vary during an outbreak. This limitation could be addressed by using nonlinear cost functions. In addition, surveys and surveillance systems used in pyramidal approaches to sCFR estimation often are not specifically created for that purpose but are “reused,” avoiding the delays and costs of setting up ad hoc studies. We took this possibility into account, and in <ext-link ext-link-type="uri" xlink:href="http://aje.oxfordjournals.org/lookup/suppl/doi:10.1093/aje/kwu213/-/DC1">Web Appendix 3</ext-link>
we provide an analysis of optimal resource allocation in situations where some data sets have a fixed size; in the case of reusing precollected data sets, their cumulative cost can be set to 0.</p>
<p>Fifth, the probabilities of progression along the severity pyramid have to be stable through time in order for our method to be feasible. This might not be the case if, for example, the symptomatic population is strongly encouraged to consult a general practitioner midway through an outbreak, making <italic>p<sub>M</sub>
</italic>
<sub>|<italic>S</italic>
</sub>
, <italic>p<sub>H</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
, and <italic>p<sub>D</sub>
</italic>
<sub>|<italic>M</italic>
</sub>
suddenly change.</p>
<p>Finally, symptomatic cases may not all be homogeneous with regard to risk of disease and death. For instance, for influenza, while the risk of infection is typically higher among children, the risks of disease and death conditional on infection generally increase with age, except for the very young (<xref rid="KWU213C31" ref-type="bibr">31</xref>
). Estimates of the sCFR can be stratified by risk group, provided that all surveillance levels used to generate estimates are similarly stratified, so that each stratum has its own independent pyramidal strategy. All of our considerations on precision and sample sizes then apply to each stratum.</p>
<p>To conclude, we have shown that statistical planning of surveillance strategies can help researchers achieve precise estimates of the sCFR with minimal costs through the judicious allocation of resources, even when those strategies are based on very approximate data early in an epidemic.</p>
</sec>
<sec sec-type="supplementary-material"><title>Supplementary Material</title>
<supplementary-material id="PMC_1" content-type="local-data"><caption><title>Web Material</title>
</caption>
<media mimetype="text" mime-subtype="html" xlink:href="supp_180_10_1036__index.html"></media>
<media xlink:role="associated-file" mimetype="application" mime-subtype="pdf" xlink:href="supp_kwu213_kwu213supp.pdf"></media>
</supplementary-material>
</sec>
</body>
<back><ack><title>ACKNOWLEDGMENTS</title>
<p>Author affiliations: MRC Centre for Outbreak Analysis and Modelling and National Institute for Health Research Health Protection Research Unit in Modelling Methodology at Imperial College London, Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom (Camille Pelat, Neil M. Ferguson, Peter J. White, Simon Cauchemez, Christophe Fraser); Modelling and Economics Unit, Public Health England, London, United Kingdom (Peter J. White); Influenza Division, National Center for Immunization and Respiratory Diseases, US Centers for Disease Control and Prevention, Atlanta, Georgia (Carrie Reed, Lyn Finelli); and Mathematical Modelling of Infectious Diseases Unit, Institut Pasteur, Paris, France (Simon Cauchemez).</p>
<p>S.C. and C.F. contributed equally to this article.</p>
<p>This work was supported by the methodology research program of the <funding-source id="cs1" rid="cn1">Medical Research Council</funding-source>
(grant <award-id id="cn1" rid="cs1">G0800596</award-id>
to C.F. and S.C.); the European Commission Seventh Framework Programme projects EMPERIE (European Management Platform for Emerging and Re-emerging Infectious Disease Entities; grant 223 498 to N.M.F.) and PREDEMICS (Preparedness, Prediction and Prevention of Emerging Zoonotic Viruses with Pandemic Potential using Multidisciplinary Approaches; grant <award-id id="cn2" rid="cs2">278433</award-id>
to S.C.); the Models of Infectious Disease Agent Study (MIDAS) network of the US <funding-source id="cs2" rid="cn2">National Institute of General Medical Sciences</funding-source>
(grant <award-id id="cn3" rid="cs3">5U54GM0884941-05</award-id>
to N.M.F. and S.C.); the United Kingdom <funding-source id="cs3" rid="cn3">National Institute for Health Research Health Protection Research Unit in Modelling Methodology at Imperial College London</funding-source>
, in partnership with <funding-source id="cs4" rid="cn4">Public Health England</funding-source>
(grant <award-id id="cn4" rid="cs4">HPRU-2012-10080</award-id>
to N.M.F., C.F., and P.J.W.); and the <funding-source id="cs5" rid="cn5">Integrative Biology of Emerging Infectious Diseases Laboratory of Excellence</funding-source>
at the Institut Pasteur (grant <award-id id="cn5" rid="cs5">ANR-10-LABX-62-IBEID</award-id>
to S.C.).</p>
<p>We thank Matthew Biggerstaff for his help with setting up the pandemic severity scenarios used in the presented simulations.</p>
<p>The funders played no role in study design, data collection and analysis, the decision to publish, or preparation of the manuscript. The views expressed are those of the authors and not necessarily those of the United Kingdom National Health Service, the United Kingdom National Institute for Health Research, the United Kingdom Department of Health, or Public Health England.</p>
<p>S.C. received consulting fees from Sanofi Pasteur MSD (Lyon, France) for a project on the modeling of varicella zoster virus. P.J.W. has received funds from Kendle Healthcare (London, United Kingdom) for providing an expert opinion on influenza and received funding from Elsevier (Reed Elsevier Group plc, Amsterdam, the Netherlands) in 2010 to attend an influenza-related meeting. All other authors have declared no competing interests.</p>
</ack>
<ref-list><title>REFERENCES</title>
<ref id="KWU213C1"><label>1</label>
<element-citation publication-type="book"><collab>World Health Organization</collab>
<article-title>Pandemic influenza preparedness and response: a WHO guidance document</article-title>
<year>2009</year>
<publisher-loc>Geneva, Switzerland</publisher-loc>
<publisher-name>World Health Organization</publisher-name>
<comment><uri xlink:type="simple" xlink:href="http://www.who.int/influenza/resources/documents/pandemic_guidance_04_2009/en/index.html">http://www.who.int/influenza/resources/documents/pandemic_guidance_04_2009/en/index.html</uri>
Accessed July 17, 2014</comment>
</element-citation>
</ref>
<ref id="KWU213C2"><label>2</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reed</surname>
<given-names>C</given-names>
</name>
<name><surname>Biggerstaff</surname>
<given-names>M</given-names>
</name>
<name><surname>Finelli</surname>
<given-names>L</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Novel framework for assessing epidemiologic effects of influenza epidemics and pandemics</article-title>
<source>Emerg Infect Dis</source>
<year>2013</year>
<volume>19</volume>
<issue>1</issue>
<fpage>85</fpage>
<lpage>91</lpage>
<pub-id pub-id-type="pmid">23260039</pub-id>
</element-citation>
</ref>
<ref id="KWU213C3"><label>3</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><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>
<etal></etal>
</person-group>
<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>1561</lpage>
<pub-id pub-id-type="pmid">19433588</pub-id>
</element-citation>
</ref>
<ref id="KWU213C4"><label>4</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Collins</surname>
<given-names>SD</given-names>
</name>
</person-group>
<article-title>Age and sex incidence of influenza and pneumonia morbidity and mortality in the epidemic of 1928–29 with comparative data for the epidemic of 1918–19: based on surveys of families in certain localities in the United States following the epidemics.</article-title>
<source>Public Health Rep</source>
<year>1931</year>
<volume>46</volume>
<issue>33</issue>
<fpage>1909</fpage>
<lpage>1937</lpage>
<pub-id pub-id-type="pmid">19315308</pub-id>
</element-citation>
</ref>
<ref id="KWU213C5"><label>5</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Garske</surname>
<given-names>T</given-names>
</name>
<name><surname>Legrand</surname>
<given-names>J</given-names>
</name>
<name><surname>Donnelly</surname>
<given-names>CA</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Assessing the severity of the novel influenza A/H1N1 pandemic</article-title>
<source>BMJ</source>
<year>2009</year>
<volume>339</volume>
<fpage>b2840</fpage>
<pub-id pub-id-type="pmid">19602714</pub-id>
</element-citation>
</ref>
<ref id="KWU213C6"><label>6</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nishiura</surname>
<given-names>H</given-names>
</name>
<name><surname>Klinkenberg</surname>
<given-names>D</given-names>
</name>
<name><surname>Roberts</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Early epidemiological assessment of the virulence of emerging infectious diseases: a case study of an influenza pandemic</article-title>
<source>PLoS One</source>
<year>2009</year>
<volume>4</volume>
<issue>8</issue>
<fpage>e6852</fpage>
<pub-id pub-id-type="pmid">19718434</pub-id>
</element-citation>
</ref>
<ref id="KWU213C7"><label>7</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mishra</surname>
<given-names>AC</given-names>
</name>
<name><surname>Chadha</surname>
<given-names>MS</given-names>
</name>
<name><surname>Choudhary</surname>
<given-names>ML</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Pandemic influenza (H1N1) 2009 is associated with severe disease in India</article-title>
<source>PLoS One</source>
<year>2010</year>
<volume>5</volume>
<issue>5</issue>
<fpage>e10540</fpage>
<pub-id pub-id-type="pmid">20479875</pub-id>
</element-citation>
</ref>
<ref id="KWU213C8"><label>8</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wilson</surname>
<given-names>N</given-names>
</name>
<name><surname>Baker</surname>
<given-names>MG</given-names>
</name>
</person-group>
<article-title>The emerging influenza pandemic: estimating the case fatality ratio</article-title>
<source>Euro Surveill</source>
<year>2009</year>
<volume>14</volume>
<issue>26</issue>
<comment>pii:19255</comment>
</element-citation>
</ref>
<ref id="KWU213C9"><label>9</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Presanis</surname>
<given-names>AM</given-names>
</name>
<name><surname>Lipsitch</surname>
<given-names>M</given-names>
</name>
<name><surname>De Angelis</surname>
<given-names>D</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The severity of pandemic H1N1 influenza in the United States, April–July 2009</article-title>
<source>PLoS Curr</source>
<year>2009</year>
<volume>1</volume>
<fpage>RRN1042</fpage>
<pub-id pub-id-type="pmid">20029614</pub-id>
</element-citation>
</ref>
<ref id="KWU213C10"><label>10</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Presanis</surname>
<given-names>AM</given-names>
</name>
<name><surname>De Angelis</surname>
<given-names>D</given-names>
</name>
<etal></etal>
</person-group>
<collab>The New York City Swine Flu Investigation Team</collab>
<article-title>The severity of pandemic H1N1 influenza in the United States, from April to July 2009: a Bayesian analysis</article-title>
<source>PLoS Med</source>
<year>2009</year>
<volume>6</volume>
<issue>12</issue>
<fpage>e1000207</fpage>
<pub-id pub-id-type="pmid">19997612</pub-id>
</element-citation>
</ref>
<ref id="KWU213C11"><label>11</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Donaldson</surname>
<given-names>LJ</given-names>
</name>
<name><surname>Rutter</surname>
<given-names>PD</given-names>
</name>
<name><surname>Ellis</surname>
<given-names>BM</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Mortality from pandemic A/H1N1 2009 influenza in England: public health surveillance study</article-title>
<source>BMJ</source>
<year>2009</year>
<volume>339</volume>
<fpage>b5213</fpage>
<pub-id pub-id-type="pmid">20007665</pub-id>
</element-citation>
</ref>
<ref id="KWU213C12"><label>12</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brooks-Pollock</surname>
<given-names>E</given-names>
</name>
<name><surname>Tilston</surname>
<given-names>N</given-names>
</name>
<name><surname>Edmunds</surname>
<given-names>WJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Using an online survey of healthcare-seeking behaviour to estimate the magnitude and severity of the 2009 H1N1v influenza epidemic in England</article-title>
<source>BMC Infect Dis</source>
<year>2011</year>
<volume>11</volume>
<fpage>68</fpage>
<pub-id pub-id-type="pmid">21410965</pub-id>
</element-citation>
</ref>
<ref id="KWU213C13"><label>13</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Steens</surname>
<given-names>A</given-names>
</name>
<name><surname>Waaijenborg</surname>
<given-names>S</given-names>
</name>
<name><surname>Teunis</surname>
<given-names>PF</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Age-dependent patterns of infection and severity explaining the low impact of 2009 influenza A (H1N1): evidence from serial serologic surveys in the Netherlands</article-title>
<source>Am J Epidemiol</source>
<year>2011</year>
<volume>174</volume>
<issue>11</issue>
<fpage>1307</fpage>
<lpage>1315</lpage>
<pub-id pub-id-type="pmid">22025354</pub-id>
</element-citation>
</ref>
<ref id="KWU213C14"><label>14</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Presanis</surname>
<given-names>AM</given-names>
</name>
<name><surname>Pebody</surname>
<given-names>RG</given-names>
</name>
<name><surname>Paterson</surname>
<given-names>BJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Changes in severity of 2009 pandemic A/H1N1 influenza in England: a Bayesian evidence synthesis</article-title>
<source>BMJ</source>
<year>2011</year>
<volume>343</volume>
<fpage>d5408</fpage>
<pub-id pub-id-type="pmid">21903689</pub-id>
</element-citation>
</ref>
<ref id="KWU213C15"><label>15</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Armitage</surname>
<given-names>P</given-names>
</name>
<name><surname>Berry</surname>
<given-names>G</given-names>
</name>
<name><surname>Matthews</surname>
<given-names>J</given-names>
</name>
</person-group>
<source>Statistical Methods in Medical Research</source>
<year>2002</year>
<publisher-loc>Oxford, United Kingdom</publisher-loc>
<publisher-name>Blackwell Scientific Publications</publisher-name>
</element-citation>
</ref>
<ref id="KWU213C16"><label>16</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><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>
<etal></etal>
</person-group>
<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="pmid">21713000</pub-id>
</element-citation>
</ref>
<ref id="KWU213C17"><label>17</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname>
<given-names>JT</given-names>
</name>
<name><surname>Ho</surname>
<given-names>A</given-names>
</name>
<name><surname>Ma</surname>
<given-names>ES</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Estimating infection attack rates and severity in real time during an influenza pandemic: analysis of serial cross-sectional serologic surveillance data</article-title>
<source>PLoS Med</source>
<year>2011</year>
<volume>8</volume>
<issue>10</issue>
<fpage>e1001103</fpage>
<pub-id pub-id-type="pmid">21990967</pub-id>
</element-citation>
</ref>
<ref id="KWU213C18"><label>18</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Atkinson</surname>
<given-names>A</given-names>
</name>
<name><surname>Donev</surname>
<given-names>A</given-names>
</name>
</person-group>
<source>Optimum Experimental Designs</source>
<year>1992</year>
<publisher-loc>Oxford, United Kingdom</publisher-loc>
<publisher-name>Clarendon Press</publisher-name>
</element-citation>
</ref>
<ref id="KWU213C19"><label>19</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Comin</surname>
<given-names>A</given-names>
</name>
<name><surname>Stegeman</surname>
<given-names>A</given-names>
</name>
<name><surname>Marangon</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Evaluating surveillance strategies for the early detection of low pathogenicity avian influenza infections</article-title>
<source>PLoS One</source>
<year>2012</year>
<volume>7</volume>
<issue>4</issue>
<fpage>e35956</fpage>
<pub-id pub-id-type="pmid">22545151</pub-id>
</element-citation>
</ref>
<ref id="KWU213C20"><label>20</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Alba</surname>
<given-names>A</given-names>
</name>
<name><surname>Casal</surname>
<given-names>J</given-names>
</name>
<name><surname>Napp</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Assessment of different surveillance systems for avian influenza in commercial poultry in Catalonia (North-Eastern Spain)</article-title>
<source>Prev Vet Med</source>
<year>2010</year>
<volume>97</volume>
<issue>2</issue>
<fpage>107</fpage>
<lpage>118</lpage>
<pub-id pub-id-type="pmid">20943281</pub-id>
</element-citation>
</ref>
<ref id="KWU213C21"><label>21</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Baptista</surname>
<given-names>FM</given-names>
</name>
<name><surname>Halasa</surname>
<given-names>T</given-names>
</name>
<name><surname>Alban</surname>
<given-names>L</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Modelling food safety and economic consequences of surveillance and control strategies for <italic>Salmonella</italic>
in pigs and pork</article-title>
<source>Epidemiol Infect</source>
<year>2011</year>
<volume>139</volume>
<issue>5</issue>
<fpage>754</fpage>
<lpage>764</lpage>
<pub-id pub-id-type="pmid">20653990</pub-id>
</element-citation>
</ref>
<ref id="KWU213C22"><label>22</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hadorn</surname>
<given-names>DC</given-names>
</name>
<name><surname>Racloz</surname>
<given-names>V</given-names>
</name>
<name><surname>Schwermer</surname>
<given-names>H</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Establishing a cost-effective national surveillance system for Bluetongue using scenario tree modelling</article-title>
<source>Vet Res</source>
<year>2009</year>
<volume>40</volume>
<issue>6</issue>
<fpage>57</fpage>
<pub-id pub-id-type="pmid">19607784</pub-id>
</element-citation>
</ref>
<ref id="KWU213C23"><label>23</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yamamoto</surname>
<given-names>T</given-names>
</name>
<name><surname>Tsutsui</surname>
<given-names>T</given-names>
</name>
<name><surname>Nishiguchi</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Evaluation of surveillance strategies for bovine brucellosis in Japan using a simulation model</article-title>
<source>Prev Vet Med</source>
<year>2008</year>
<volume>86</volume>
<issue>1-2</issue>
<fpage>57</fpage>
<lpage>74</lpage>
<pub-id pub-id-type="pmid">18440660</pub-id>
</element-citation>
</ref>
<ref id="KWU213C24"><label>24</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Polgreen</surname>
<given-names>PM</given-names>
</name>
<name><surname>Chen</surname>
<given-names>Z</given-names>
</name>
<name><surname>Segre</surname>
<given-names>AM</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Optimizing influenza sentinel surveillance at the state level</article-title>
<source>Am J Epidemiol</source>
<year>2009</year>
<volume>170</volume>
<issue>10</issue>
<fpage>1300</fpage>
<lpage>1306</lpage>
<pub-id pub-id-type="pmid">19822570</pub-id>
</element-citation>
</ref>
<ref id="KWU213C25"><label>25</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Scarpino</surname>
<given-names>SV</given-names>
</name>
<name><surname>Dimitrov</surname>
<given-names>NB</given-names>
</name>
<name><surname>Meyers</surname>
<given-names>LA</given-names>
</name>
</person-group>
<article-title>Optimizing provider recruitment for influenza surveillance networks</article-title>
<source>PLoS Comput Biol</source>
<year>2012</year>
<volume>8</volume>
<issue>4</issue>
<fpage>e1002472</fpage>
<pub-id pub-id-type="pmid">22511860</pub-id>
</element-citation>
</ref>
<ref id="KWU213C26"><label>26</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname>
<given-names>P</given-names>
</name>
<name><surname>Duan</surname>
<given-names>W</given-names>
</name>
<name><surname>Lv</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Review of an influenza surveillance system, Beijing, People's Republic of China</article-title>
<source>Emerg Infect Dis</source>
<year>2009</year>
<volume>15</volume>
<issue>10</issue>
<fpage>1603</fpage>
<lpage>1608</lpage>
<pub-id pub-id-type="pmid">19861053</pub-id>
</element-citation>
</ref>
<ref id="KWU213C27"><label>27</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Byass</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Empirical modelling of population sampling: lessons for designing sentinel surveillance</article-title>
<source>Public Health</source>
<year>2003</year>
<volume>117</volume>
<issue>1</issue>
<fpage>36</fpage>
<lpage>42</lpage>
<pub-id pub-id-type="pmid">12802903</pub-id>
</element-citation>
</ref>
<ref id="KWU213C28"><label>28</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ejima</surname>
<given-names>K</given-names>
</name>
<name><surname>Omori</surname>
<given-names>R</given-names>
</name>
<name><surname>Cowling</surname>
<given-names>BJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The time required to estimate the case fatality ratio of influenza using only the tip of an iceberg: joint estimation of the virulence and the transmission potential</article-title>
<source>Comput Math Methods Med</source>
<year>2012</year>
<volume>2012</volume>
<fpage>978901</fpage>
<pub-id pub-id-type="pmid">22649483</pub-id>
</element-citation>
</ref>
<ref id="KWU213C29"><label>29</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reich</surname>
<given-names>NG</given-names>
</name>
<name><surname>Lessler</surname>
<given-names>J</given-names>
</name>
<name><surname>Cummings</surname>
<given-names>DA</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Estimating absolute and relative case fatality ratios from infectious disease surveillance data</article-title>
<source>Biometrics</source>
<year>2012</year>
<volume>68</volume>
<issue>2</issue>
<fpage>598</fpage>
<lpage>606</lpage>
<pub-id pub-id-type="pmid">22276951</pub-id>
</element-citation>
</ref>
<ref id="KWU213C30"><label>30</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghani</surname>
<given-names>AC</given-names>
</name>
<name><surname>Donnelly</surname>
<given-names>CA</given-names>
</name>
<name><surname>Cox</surname>
<given-names>DR</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Methods for estimating the case fatality ratio for a novel, emerging infectious disease</article-title>
<source>Am J Epidemiol</source>
<year>2005</year>
<volume>162</volume>
<issue>5</issue>
<fpage>479</fpage>
<lpage>486</lpage>
<pub-id pub-id-type="pmid">16076827</pub-id>
</element-citation>
</ref>
<ref id="KWU213C31"><label>31</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname>
<given-names>JT</given-names>
</name>
<name><surname>Ma</surname>
<given-names>ES</given-names>
</name>
<name><surname>Lee</surname>
<given-names>CK</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The infection attack rate and severity of 2009 pandemic H1N1 influenza in Hong Kong</article-title>
<source>Clin Infect Dis</source>
<year>2010</year>
<volume>51</volume>
<issue>10</issue>
<fpage>1184</fpage>
<lpage>1191</lpage>
<pub-id pub-id-type="pmid">20964521</pub-id>
</element-citation>
</ref>
<ref id="KWU213C32"><label>32</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Payne</surname>
<given-names>AM</given-names>
</name>
</person-group>
<article-title>Some aspects of the epidemiology of the 1957 influenza pandemic</article-title>
<source>Proc R Soc Med</source>
<year>1958</year>
<volume>51</volume>
<issue>12</issue>
<fpage>1009</fpage>
<lpage>1015</lpage>
<pub-id pub-id-type="pmid">13623821</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/PandemieGrippaleV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000904 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000904 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= PandemieGrippaleV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:4240167 |texte= Optimizing the Precision of Case Fatality Ratio Estimates Under the Surveillance Pyramid Approach }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:25255809" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a PandemieGrippaleV1
![]() | This area was generated with Dilib version V0.6.34. | ![]() |