Deriving effective vaccine allocation strategies for pandemic influenza: Comparison of an agent-based simulation and a compartmental model
Identifieur interne : 000C50 ( Pmc/Corpus ); précédent : 000C49; suivant : 000C51Deriving effective vaccine allocation strategies for pandemic influenza: Comparison of an agent-based simulation and a compartmental model
Auteurs : Özden O. Dalg Ç ; Osman Y. Özalt N ; William A. Ciccotelli ; Fatih S. ErenaySource :
- PLoS ONE [ 1932-6203 ] ; 2017.
Abstract
Individuals are prioritized based on their risk profiles when allocating limited vaccine stocks during an influenza pandemic. Computationally expensive but realistic agent-based simulations and fast but stylized compartmental models are typically used to derive effective vaccine allocation strategies. A detailed comparison of these two approaches, however, is often omitted. We derive age-specific vaccine allocation strategies to mitigate a pandemic influenza outbreak in Seattle by applying derivative-free optimization to an agent-based simulation and also to a compartmental model. We compare the strategies derived by these two approaches under various infection aggressiveness and vaccine coverage scenarios. We observe that both approaches primarily vaccinate school children, however they may allocate the remaining vaccines in different ways. The vaccine allocation strategies derived by using the agent-based simulation are associated with up to 70% decrease in total cost and 34% reduction in the number of infections compared to the strategies derived by using the compartmental model. Nevertheless, the latter approach may still be competitive for very low and/or very high infection aggressiveness. Our results provide insights about potential differences between the vaccine allocation strategies derived by using agent-based simulations and those derived by using compartmental models.
Url:
DOI: 10.1371/journal.pone.0172261
PubMed: 28222123
PubMed Central: 5319753
Links to Exploration step
PMC:5319753Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Deriving effective vaccine allocation strategies for pandemic influenza: Comparison of an agent-based simulation and a compartmental model</title>
<author><name sortKey="Dalg C, Ozden O" sort="Dalg C, Ozden O" uniqKey="Dalg C O" first="Özden O." last="Dalg Ç">Özden O. Dalg Ç</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Management Sciences, University of Waterloo, Waterloo, Ontario, Canada</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ozalt N, Osman Y" sort="Ozalt N, Osman Y" uniqKey="Ozalt N O" first="Osman Y." last="Özalt N">Osman Y. Özalt N</name>
<affiliation><nlm:aff id="aff002"><addr-line>Edward P. Fitts Department of Industrial and Systems Engineering, North Carolina State University, Raleigh, North Carolina, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ciccotelli, William A" sort="Ciccotelli, William A" uniqKey="Ciccotelli W" first="William A." last="Ciccotelli">William A. Ciccotelli</name>
<affiliation><nlm:aff id="aff003"><addr-line>Department of Pathology and Molecular Medicine, McMaster University, Hamilton, Ontario, Canada</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff004"><addr-line>Grand River Hospital, Kitchener, Ontario, Canada</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Erenay, Fatih S" sort="Erenay, Fatih S" uniqKey="Erenay F" first="Fatih S." last="Erenay">Fatih S. Erenay</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Management Sciences, University of Waterloo, Waterloo, Ontario, Canada</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">28222123</idno>
<idno type="pmc">5319753</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5319753</idno>
<idno type="RBID">PMC:5319753</idno>
<idno type="doi">10.1371/journal.pone.0172261</idno>
<date when="2017">2017</date>
<idno type="wicri:Area/Pmc/Corpus">000C50</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000C50</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Deriving effective vaccine allocation strategies for pandemic influenza: Comparison of an agent-based simulation and a compartmental model</title>
<author><name sortKey="Dalg C, Ozden O" sort="Dalg C, Ozden O" uniqKey="Dalg C O" first="Özden O." last="Dalg Ç">Özden O. Dalg Ç</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Management Sciences, University of Waterloo, Waterloo, Ontario, Canada</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ozalt N, Osman Y" sort="Ozalt N, Osman Y" uniqKey="Ozalt N O" first="Osman Y." last="Özalt N">Osman Y. Özalt N</name>
<affiliation><nlm:aff id="aff002"><addr-line>Edward P. Fitts Department of Industrial and Systems Engineering, North Carolina State University, Raleigh, North Carolina, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ciccotelli, William A" sort="Ciccotelli, William A" uniqKey="Ciccotelli W" first="William A." last="Ciccotelli">William A. Ciccotelli</name>
<affiliation><nlm:aff id="aff003"><addr-line>Department of Pathology and Molecular Medicine, McMaster University, Hamilton, Ontario, Canada</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff004"><addr-line>Grand River Hospital, Kitchener, Ontario, Canada</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Erenay, Fatih S" sort="Erenay, Fatih S" uniqKey="Erenay F" first="Fatih S." last="Erenay">Fatih S. Erenay</name>
<affiliation><nlm:aff id="aff001"><addr-line>Department of Management Sciences, University of Waterloo, Waterloo, Ontario, Canada</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint><date when="2017">2017</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>Individuals are prioritized based on their risk profiles when allocating limited vaccine stocks during an influenza pandemic. Computationally expensive but realistic agent-based simulations and fast but stylized compartmental models are typically used to derive effective vaccine allocation strategies. A detailed comparison of these two approaches, however, is often omitted. We derive age-specific vaccine allocation strategies to mitigate a pandemic influenza outbreak in Seattle by applying derivative-free optimization to an agent-based simulation and also to a compartmental model. We compare the strategies derived by these two approaches under various infection aggressiveness and vaccine coverage scenarios. We observe that both approaches primarily vaccinate school children, however they may allocate the remaining vaccines in different ways. The vaccine allocation strategies derived by using the agent-based simulation are associated with up to 70% decrease in total cost and 34% reduction in the number of infections compared to the strategies derived by using the compartmental model. Nevertheless, the latter approach may still be competitive for very low and/or very high infection aggressiveness. Our results provide insights about potential differences between the vaccine allocation strategies derived by using agent-based simulations and those derived by using compartmental models.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Molinari, Nam" uniqKey="Molinari N">NAM Molinari</name>
</author>
<author><name sortKey="Ortega Sanchez, Ir" uniqKey="Ortega Sanchez I">IR Ortega-Sanchez</name>
</author>
<author><name sortKey="Messonnier, Ml" uniqKey="Messonnier M">ML Messonnier</name>
</author>
<author><name sortKey="Thompson, Ww" uniqKey="Thompson W">WW Thompson</name>
</author>
<author><name sortKey="Wortley, Pm" uniqKey="Wortley P">PM Wortley</name>
</author>
<author><name sortKey="Weintraub, E" uniqKey="Weintraub E">E Weintraub</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Parvin, Jd" uniqKey="Parvin J">JD Parvin</name>
</author>
<author><name sortKey="Moscona, A" uniqKey="Moscona A">A Moscona</name>
</author>
<author><name sortKey="Pan, W" uniqKey="Pan W">W Pan</name>
</author>
<author><name sortKey="Leider, J" uniqKey="Leider J">J Leider</name>
</author>
<author><name sortKey="Palese, P" uniqKey="Palese P">P Palese</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Taubenberger, Jk" uniqKey="Taubenberger J">JK Taubenberger</name>
</author>
<author><name sortKey="Morens, Dm" uniqKey="Morens D">DM Morens</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chao, Dl" uniqKey="Chao D">DL Chao</name>
</author>
<author><name sortKey="Halloran, Me" uniqKey="Halloran M">ME Halloran</name>
</author>
<author><name sortKey="Obenchain, Vj" uniqKey="Obenchain V">VJ Obenchain</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Talbot, Hk" uniqKey="Talbot H">HK Talbot</name>
</author>
<author><name sortKey="Zhu, Y" uniqKey="Zhu Y">Y Zhu</name>
</author>
<author><name sortKey="Chen, Q" uniqKey="Chen Q">Q Chen</name>
</author>
<author><name sortKey="Williams, Jv" uniqKey="Williams J">JV Williams</name>
</author>
<author><name sortKey="Thompson, Mg" uniqKey="Thompson M">MG Thompson</name>
</author>
<author><name sortKey="Griffin, Mr" uniqKey="Griffin M">MR Griffin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nichol, Kl" uniqKey="Nichol K">KL Nichol</name>
</author>
<author><name sortKey="Treanor, Jj" uniqKey="Treanor J">JJ Treanor</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hethcote, Hw" uniqKey="Hethcote H">HW Hethcote</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zaric, Gs" uniqKey="Zaric G">GS Zaric</name>
</author>
<author><name sortKey="Brandeau, Ml" uniqKey="Brandeau M">ML Brandeau</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dasbach, Ej" uniqKey="Dasbach E">EJ Dasbach</name>
</author>
<author><name sortKey="Elbasha, Eh" uniqKey="Elbasha E">EH Elbasha</name>
</author>
<author><name sortKey="Insinga, Rp" uniqKey="Insinga R">RP Insinga</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huang, G" uniqKey="Huang G">G Huang</name>
</author>
<author><name sortKey="Li, L" uniqKey="Li L">L Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Medlock, J" uniqKey="Medlock J">J Medlock</name>
</author>
<author><name sortKey="Galvani, Ap" uniqKey="Galvani A">AP Galvani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Medlock, J" uniqKey="Medlock J">J Medlock</name>
</author>
<author><name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name>
</author>
<author><name sortKey="Galvani, A" uniqKey="Galvani A">A Galvani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Eames, Kt" uniqKey="Eames K">KT Eames</name>
</author>
<author><name sortKey="Tilston, Nl" uniqKey="Tilston N">NL Tilston</name>
</author>
<author><name sortKey="Brooks Pollock, E" uniqKey="Brooks Pollock E">E Brooks-Pollock</name>
</author>
<author><name sortKey="Edmunds, Wj" uniqKey="Edmunds W">WJ Edmunds</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Araz, Om" uniqKey="Araz O">OM Araz</name>
</author>
<author><name sortKey="Galvani, A" uniqKey="Galvani A">A Galvani</name>
</author>
<author><name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Patel, R" uniqKey="Patel R">R Patel</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
<author><name sortKey="Halloran, Me" uniqKey="Halloran M">ME Halloran</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Das, Tk" uniqKey="Das T">TK Das</name>
</author>
<author><name sortKey="Savachkin, Aa" uniqKey="Savachkin A">AA Savachkin</name>
</author>
<author><name sortKey="Zhu, Y" uniqKey="Zhu Y">Y Zhu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Basta, Ne" uniqKey="Basta N">NE Basta</name>
</author>
<author><name sortKey="Chao, Dl" uniqKey="Chao D">DL Chao</name>
</author>
<author><name sortKey="Halloran, Me" uniqKey="Halloran M">ME Halloran</name>
</author>
<author><name sortKey="Matrajt, L" uniqKey="Matrajt L">L Matrajt</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lee, By" uniqKey="Lee B">BY Lee</name>
</author>
<author><name sortKey="Brown, St" uniqKey="Brown S">ST Brown</name>
</author>
<author><name sortKey="Korch, Gw" uniqKey="Korch G">GW Korch</name>
</author>
<author><name sortKey="Cooley, Pc" uniqKey="Cooley P">PC Cooley</name>
</author>
<author><name sortKey="Zimmerman, Rk" uniqKey="Zimmerman R">RK Zimmerman</name>
</author>
<author><name sortKey="Wheaton, Wd" uniqKey="Wheaton W">WD Wheaton</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Uribe Sanchez, A" uniqKey="Uribe Sanchez A">A Uribe-Sánchez</name>
</author>
<author><name sortKey="Savachkin, A" uniqKey="Savachkin A">A Savachkin</name>
</author>
<author><name sortKey="Santana, A" uniqKey="Santana A">A Santana</name>
</author>
<author><name sortKey="Prieto Santa, D" uniqKey="Prieto Santa D">D Prieto-Santa</name>
</author>
<author><name sortKey="Das, Tk" uniqKey="Das T">TK Das</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Andrad Ttir, S" uniqKey="Andrad Ttir S">S Andradóttir</name>
</author>
<author><name sortKey="Chiu, W" uniqKey="Chiu W">W Chiu</name>
</author>
<author><name sortKey="Goldsman, D" uniqKey="Goldsman D">D Goldsman</name>
</author>
<author><name sortKey="Lee, Ml" uniqKey="Lee M">ML Lee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dhamodharan, A" uniqKey="Dhamodharan A">A Dhamodharan</name>
</author>
<author><name sortKey="Proano, Ra" uniqKey="Proano R">RA Proano</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Feng, W" uniqKey="Feng W">W Feng</name>
</author>
<author><name sortKey="Kong, N" uniqKey="Kong N">N Kong</name>
</author>
<author><name sortKey="Wan, H" uniqKey="Wan H">H Wan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chhatwal, J" uniqKey="Chhatwal J">J Chhatwal</name>
</author>
<author><name sortKey="He, T" uniqKey="He T">T He</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rahmandad, H" uniqKey="Rahmandad H">H Rahmandad</name>
</author>
<author><name sortKey="Sterman, J" uniqKey="Sterman J">J Sterman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Germann, Tc" uniqKey="Germann T">TC Germann</name>
</author>
<author><name sortKey="Kadau, K" uniqKey="Kadau K">K Kadau</name>
</author>
<author><name sortKey="L, Im" uniqKey="L I">IM L</name>
</author>
<author><name sortKey="Macken, Ca" uniqKey="Macken C">CA Macken</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Anderson, Rm" uniqKey="Anderson R">RM Anderson</name>
</author>
<author><name sortKey="May, Rm" uniqKey="May R">RM May</name>
</author>
<author><name sortKey="Anderson, B" uniqKey="Anderson B">B Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mollison, D" uniqKey="Mollison D">D Mollison</name>
</author>
<author><name sortKey="Isham, V" uniqKey="Isham V">V Isham</name>
</author>
<author><name sortKey="Grenfell, B" uniqKey="Grenfell B">B Grenfell</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bansal, S" uniqKey="Bansal S">S Bansal</name>
</author>
<author><name sortKey="Grenfell, Bt" uniqKey="Grenfell B">BT Grenfell</name>
</author>
<author><name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lloyd, Al" uniqKey="Lloyd A">AL Lloyd</name>
</author>
<author><name sortKey="May, Rm" uniqKey="May R">RM May</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
<author><name sortKey="Cummings, Da" uniqKey="Cummings D">DA Cummings</name>
</author>
<author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name>
</author>
<author><name sortKey="Cajka, Jc" uniqKey="Cajka J">JC Cajka</name>
</author>
<author><name sortKey="Cooley, Pc" uniqKey="Cooley P">PC Cooley</name>
</author>
<author><name sortKey="Burke, Ds" uniqKey="Burke D">DS Burke</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Merler, S" uniqKey="Merler S">S Merler</name>
</author>
<author><name sortKey="Ajelli, M" uniqKey="Ajelli M">M Ajelli</name>
</author>
<author><name sortKey="Rizzo, C" uniqKey="Rizzo C">C Rizzo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hladish, T" uniqKey="Hladish T">T Hladish</name>
</author>
<author><name sortKey="Melamud, E" uniqKey="Melamud E">E Melamud</name>
</author>
<author><name sortKey="Barrera, La" uniqKey="Barrera L">LA Barrera</name>
</author>
<author><name sortKey="Galvani, A" uniqKey="Galvani A">A Galvani</name>
</author>
<author><name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Carrat, F" uniqKey="Carrat F">F Carrat</name>
</author>
<author><name sortKey="Vergu, E" uniqKey="Vergu E">E Vergu</name>
</author>
<author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
<author><name sortKey="Lemaitre, M" uniqKey="Lemaitre M">M Lemaitre</name>
</author>
<author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name>
</author>
<author><name sortKey="Leach, S" uniqKey="Leach S">S Leach</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Struchiner, Cj" uniqKey="Struchiner C">CJ Struchiner</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
<author><name sortKey="Halloran, Me" uniqKey="Halloran M">ME Halloran</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Atkinson, Ka" uniqKey="Atkinson K">KA Atkinson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Arias, E" uniqKey="Arias E">E Arias</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Basta, Ne" uniqKey="Basta N">NE Basta</name>
</author>
<author><name sortKey="Halloran, Me" uniqKey="Halloran M">ME Halloran</name>
</author>
<author><name sortKey="Matrajt, L" uniqKey="Matrajt L">L Matrajt</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Le Digabel, S" uniqKey="Le Digabel S">S Le Digabel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Heffernan, J" uniqKey="Heffernan J">J Heffernan</name>
</author>
<author><name sortKey="Smith, R" uniqKey="Smith R">R Smith</name>
</author>
<author><name sortKey="Wahl, L" uniqKey="Wahl L">L Wahl</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></biblStruct>
<biblStruct><analytic><author><name sortKey="Matrajt, L" uniqKey="Matrajt L">L Matrajt</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mylius, Sd" uniqKey="Mylius S">SD Mylius</name>
</author>
<author><name sortKey="Hagenaars, Tj" uniqKey="Hagenaars T">TJ Hagenaars</name>
</author>
<author><name sortKey="Lugner, Ak" uniqKey="Lugner A">AK Lugnér</name>
</author>
<author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name>
</author>
<author><name sortKey="Ng, S" uniqKey="Ng S">S Ng</name>
</author>
<author><name sortKey="Ma, Es" uniqKey="Ma E">ES Ma</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lazer, D" uniqKey="Lazer D">D Lazer</name>
</author>
<author><name sortKey="Kennedy, R" uniqKey="Kennedy R">R Kennedy</name>
</author>
<author><name sortKey="King, G" uniqKey="King G">G King</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pastor Satorras, R" uniqKey="Pastor Satorras R">R Pastor-Satorras</name>
</author>
<author><name sortKey="Castellano, C" uniqKey="Castellano C">C Castellano</name>
</author>
<author><name sortKey="Van Mieghem, P" uniqKey="Van Mieghem P">P Van Mieghem</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group><journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, CA USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">28222123</article-id>
<article-id pub-id-type="pmc">5319753</article-id>
<article-id pub-id-type="publisher-id">PONE-D-16-37663</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0172261</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject>
<subj-group><subject>Immunology</subject>
<subj-group><subject>Vaccination and Immunization</subject>
<subj-group><subject>Vaccines</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject>
<subj-group><subject>Immunology</subject>
<subj-group><subject>Vaccination and Immunization</subject>
<subj-group><subject>Vaccines</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject>
<subj-group><subject>Public and Occupational Health</subject>
<subj-group><subject>Preventive Medicine</subject>
<subj-group><subject>Vaccination and Immunization</subject>
<subj-group><subject>Vaccines</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Research and Analysis Methods</subject>
<subj-group><subject>Simulation and Modeling</subject>
<subj-group><subject>Agent-Based Modeling</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Computer and Information Sciences</subject>
<subj-group><subject>Systems Science</subject>
<subj-group><subject>Agent-Based Modeling</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Physical Sciences</subject>
<subj-group><subject>Mathematics</subject>
<subj-group><subject>Systems Science</subject>
<subj-group><subject>Agent-Based Modeling</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>People and Places</subject>
<subj-group><subject>Population Groupings</subject>
<subj-group><subject>Age Groups</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject>
<subj-group><subject>Immunology</subject>
<subj-group><subject>Vaccination and Immunization</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject>
<subj-group><subject>Immunology</subject>
<subj-group><subject>Vaccination and Immunization</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject>
<subj-group><subject>Public and Occupational Health</subject>
<subj-group><subject>Preventive Medicine</subject>
<subj-group><subject>Vaccination and Immunization</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Social Sciences</subject>
<subj-group><subject>Sociology</subject>
<subj-group><subject>Education</subject>
<subj-group><subject>Schools</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Research and Analysis Methods</subject>
<subj-group><subject>Simulation and Modeling</subject>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject>
<subj-group><subject>Epidemiology</subject>
<subj-group><subject>Infectious Disease Epidemiology</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject>
<subj-group><subject>Infectious Diseases</subject>
<subj-group><subject>Infectious Disease Epidemiology</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject>
<subj-group><subject>Infectious Diseases</subject>
<subj-group><subject>Viral Diseases</subject>
<subj-group><subject>Influenza</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>Deriving effective vaccine allocation strategies for pandemic influenza: Comparison of an agent-based simulation and a compartmental model</article-title>
<alt-title alt-title-type="running-head">Deriving effective vaccine allocation strategies for pandemic influenza</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Dalgıç</surname>
<given-names>Özden O.</given-names>
</name>
<xref ref-type="aff" rid="aff001"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Özaltın</surname>
<given-names>Osman Y.</given-names>
</name>
<xref ref-type="aff" rid="aff002"><sup>2</sup>
</xref>
<xref ref-type="corresp" rid="cor001">*</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ciccotelli</surname>
<given-names>William A.</given-names>
</name>
<xref ref-type="aff" rid="aff003"><sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff004"><sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Erenay</surname>
<given-names>Fatih S.</given-names>
</name>
<xref ref-type="aff" rid="aff001"><sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff001"><label>1</label>
<addr-line>Department of Management Sciences, University of Waterloo, Waterloo, Ontario, Canada</addr-line>
</aff>
<aff id="aff002"><label>2</label>
<addr-line>Edward P. Fitts Department of Industrial and Systems Engineering, North Carolina State University, Raleigh, North Carolina, United States of America</addr-line>
</aff>
<aff id="aff003"><label>3</label>
<addr-line>Department of Pathology and Molecular Medicine, McMaster University, Hamilton, Ontario, Canada</addr-line>
</aff>
<aff id="aff004"><label>4</label>
<addr-line>Grand River Hospital, Kitchener, Ontario, Canada</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Shaman</surname>
<given-names>Jeffrey</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Columbia University, UNITED STATES</addr-line>
</aff>
<author-notes><fn fn-type="COI-statement" id="coi001"><p><bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p><list list-type="simple"><list-item><p><bold>Conceptualization:</bold>
OD OO FE.</p>
</list-item>
<list-item><p><bold>Formal analysis:</bold>
OD.</p>
</list-item>
<list-item><p><bold>Funding acquisition:</bold>
OO FE.</p>
</list-item>
<list-item><p><bold>Methodology:</bold>
OD OO FE.</p>
</list-item>
<list-item><p><bold>Project administration:</bold>
OO FE.</p>
</list-item>
<list-item><p><bold>Software:</bold>
OD OO.</p>
</list-item>
<list-item><p><bold>Supervision:</bold>
OO FE WC.</p>
</list-item>
<list-item><p><bold>Writing – original draft:</bold>
OD OO FE.</p>
</list-item>
<list-item><p><bold>Writing – review & editing:</bold>
OD OO FE.</p>
</list-item>
</list>
</p>
</fn>
<corresp id="cor001">* E-mail: <email>oyozalti@ncsu.edu</email>
</corresp>
</author-notes>
<pub-date pub-type="collection"><year>2017</year>
</pub-date>
<pub-date pub-type="epub"><day>21</day>
<month>2</month>
<year>2017</year>
</pub-date>
<volume>12</volume>
<issue>2</issue>
<elocation-id>e0172261</elocation-id>
<history><date date-type="received"><day>20</day>
<month>9</month>
<year>2016</year>
</date>
<date date-type="accepted"><day>2</day>
<month>2</month>
<year>2017</year>
</date>
</history>
<permissions><copyright-statement>© 2017 Dalgıç et al</copyright-statement>
<copyright-year>2017</copyright-year>
<copyright-holder>Dalgıç et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>
, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="pone.0172261.pdf"></self-uri>
<abstract><p>Individuals are prioritized based on their risk profiles when allocating limited vaccine stocks during an influenza pandemic. Computationally expensive but realistic agent-based simulations and fast but stylized compartmental models are typically used to derive effective vaccine allocation strategies. A detailed comparison of these two approaches, however, is often omitted. We derive age-specific vaccine allocation strategies to mitigate a pandemic influenza outbreak in Seattle by applying derivative-free optimization to an agent-based simulation and also to a compartmental model. We compare the strategies derived by these two approaches under various infection aggressiveness and vaccine coverage scenarios. We observe that both approaches primarily vaccinate school children, however they may allocate the remaining vaccines in different ways. The vaccine allocation strategies derived by using the agent-based simulation are associated with up to 70% decrease in total cost and 34% reduction in the number of infections compared to the strategies derived by using the compartmental model. Nevertheless, the latter approach may still be competitive for very low and/or very high infection aggressiveness. Our results provide insights about potential differences between the vaccine allocation strategies derived by using agent-based simulations and those derived by using compartmental models.</p>
</abstract>
<funding-group><award-group id="award001"><funding-source><institution-wrap><institution-id institution-id-type="funder-id">http://dx.doi.org/10.13039/501100000038</institution-id>
<institution>Natural Sciences and Engineering Research Council of Canada</institution>
</institution-wrap>
</funding-source>
<award-id>DG 113788</award-id>
<principal-award-recipient><name><surname>Erenay</surname>
<given-names>Fatih S</given-names>
</name>
</principal-award-recipient>
</award-group>
<award-group id="award002"><funding-source><institution-wrap><institution-id institution-id-type="funder-id">http://dx.doi.org/10.13039/501100000038</institution-id>
<institution>Natural Sciences and Engineering Research Council of Canada</institution>
</institution-wrap>
</funding-source>
<award-id>DG 113790</award-id>
<principal-award-recipient><name><surname>Özaltın</surname>
<given-names>Osman Y</given-names>
</name>
</principal-award-recipient>
</award-group>
<funding-statement>This work was supported by the Natural Sciences and Engineering Research Council of Canada (DG 113788 and DG 113790). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts><fig-count count="3"></fig-count>
<table-count count="5"></table-count>
<page-count count="19"></page-count>
</counts>
<custom-meta-group><custom-meta id="data-availability"><meta-name>Data Availability</meta-name>
<meta-value>All relevant data are within the paper and its Supporting Information files.</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
<notes><title>Data Availability</title>
<p>All relevant data are within the paper and its Supporting Information files.</p>
</notes>
</front>
<body><sec sec-type="intro" id="sec001"><title>Introduction</title>
<p>Influenza is a highly contagious viral disease. Each year 5–15% of the world population is infected with influenza resulting in 3–5 million severe cases and 250,000–500,000 deaths [<xref rid="pone.0172261.ref001" ref-type="bibr">1</xref>
]. Together with pneumonia, influenza is the seventh leading cause of death in the U.S. [<xref rid="pone.0172261.ref002" ref-type="bibr">2</xref>
]. Considering the cost of medical treatments and working day losses, the annual burden of influenza epidemics to the U.S. economy was estimated as $87.1 billion in 2003 [<xref rid="pone.0172261.ref003" ref-type="bibr">3</xref>
].</p>
<p>An emerging virus that spreads globally may lead to an influenza pandemic [<xref rid="pone.0172261.ref004" ref-type="bibr">4</xref>
]. Unlike seasonal epidemics, pandemics occur abruptly and cause horrendous death tolls, e.g., the 1918 Spanish influenza pandemic infected around 500 million and killed approximately 50 million people worldwide [<xref rid="pone.0172261.ref005" ref-type="bibr">5</xref>
]. Among pandemic mitigation interventions (e.g., social distancing, public health measures, antiviral prophylaxis), vaccination provides the most efficient and durable response [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
, <xref rid="pone.0172261.ref007" ref-type="bibr">7</xref>
]. However, the supply of influenza vaccine is often scarce during a pandemic due to production limitations. Thus, the population is prioritized based on risk-factors related to influenza exposure and transmissibility when distributing the available doses [<xref rid="pone.0172261.ref008" ref-type="bibr">8</xref>
].</p>
<p>In the literature, compartmental models [<xref rid="pone.0172261.ref009" ref-type="bibr">9</xref>
–<xref rid="pone.0172261.ref016" ref-type="bibr">16</xref>
] and agent-based simulations [<xref rid="pone.0172261.ref017" ref-type="bibr">17</xref>
–<xref rid="pone.0172261.ref022" ref-type="bibr">22</xref>
] are frequently employed to make mitigation plans for influenza pandemics and evaluate the effectiveness of various public health interventions [<xref rid="pone.0172261.ref023" ref-type="bibr">23</xref>
–<xref rid="pone.0172261.ref025" ref-type="bibr">25</xref>
]. Although it is known that the infection propagation is different in these two approaches [<xref rid="pone.0172261.ref026" ref-type="bibr">26</xref>
], a detailed comparison of the strategies derived by using them is often omitted. On the one hand, compartmental models represent the number of individuals in each stage (or compartment) of the epidemic (e.g., susceptible, exposed, infected, recovered) by continuous-time state variables, and formulate the transitions among different compartments using differential equations. These models can rapidly evaluate many scenarios and intervention strategies, but they assume that individuals in each compartment mix uniformly and randomly with each other. Moreover, deterministic compartmental models do not consider the uncertainties in disease propagation (e.g., stochasticity in transmission events, incubation, and recovery periods) [<xref rid="pone.0172261.ref022" ref-type="bibr">22</xref>
]. Thus, such models may not accurately model infection dynamics, especially at the initial and final stages of a pandemic when few infectious individuals exist [<xref rid="pone.0172261.ref027" ref-type="bibr">27</xref>
]. Despite their simplifying assumptions, compartmental models have proved to be predictive [<xref rid="pone.0172261.ref028" ref-type="bibr">28</xref>
–<xref rid="pone.0172261.ref030" ref-type="bibr">30</xref>
], and they have been successfully extended to capture large-scale host heterogeneities. These extensions of the simple compartmental framework include age-specific contact patterns [<xref rid="pone.0172261.ref013" ref-type="bibr">13</xref>
] and heterogeneities induced by spatial structure [<xref rid="pone.0172261.ref031" ref-type="bibr">31</xref>
]. Agent-based simulations, on the other hand, consider uncertainties about the infection parameters, and they store individual-level information to model contact patterns in a population at the expense of increased computational burden [<xref rid="pone.0172261.ref032" ref-type="bibr">32</xref>
–<xref rid="pone.0172261.ref034" ref-type="bibr">34</xref>
].</p>
<p>Our goal in this paper is to identify possible scenarios under which performances of the effective age-specific vaccine allocation strategies derived by using compartmental models and agent-based simulations may differ significantly in practical settings. For this purpose, we consider an influenza pandemic in Seattle using a custom-built deterministic compartmental model and an agent-based simulation developed by Chao et al. (2010) [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]. The compartmental model is calibrated to closely approximate the results of the agent-based simulation under no vaccination. We apply mesh-adaptive optimization to derive effective age-specific vaccine allocation strategies based on four different objective functions. At each iteration of the optimization process, the performances of the newly generated vaccine allocation strategies are evaluated using the agent-based simulation in one set of experiments, and using the compartmental model in the other set of experiments. We perform sensitivity analysis to identify the differences between these two approaches under various vaccine coverage and infection aggressiveness scenarios.</p>
<p>We observe that the age-specific vaccine allocation strategies derived by using computationally expensive but more realistic agent-based simulation and those derived by using fast but more stylized compartmental model are different, although both models are calibrated to generate similar results under no vaccination. We use the agent-based simulation to evaluate the performances of strategies derived by the compartmental model. Our results show that the vaccine allocation strategies derived by the agent-based simulation are associated with up to 70% decrease in total cost and 34% reduction in the number of infections compared to the strategies derived by the compartmental model. Nevertheless, the latter approach may still be competitive for very low and/or very high infection aggressiveness scenarios.</p>
<p>It is clear that any two infectious disease spread models can differ from each other with respect to the assumptions in their design and parametrization. Furthermore, the degree to which each modeling approach allows for inclusion of heterogeneity and uncertainty varies. Therefore, the empirical comparison presented in this study is valid for the considered agent-based influenza pandemic simulation which is well-known and commonly used in the literature [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]. Our results, however, still provide important insights into the possible differences between the vaccine allocation strategies derived by using agent-based simulations and deterministic compartmental models.</p>
</sec>
<sec sec-type="materials|methods" id="sec002"><title>Materials and methods</title>
<p>We consider different levels of vaccine coverage and infection aggressiveness, and apply mesh-adaptive optimization [<xref rid="pone.0172261.ref035" ref-type="bibr">35</xref>
] to find effective strategies for allocating limited vaccine doses to different age-groups in the population with respect to four performance measures: total cost, number of deaths, number of infections, and years of life lost. We evaluate the performances of the trial vaccine allocation strategies at each iteration of the optimization algorithm using an agent-based simulation in one set of experiments, and a compartmental model in the other set of experiments.</p>
<sec id="sec003"><title>Agent-based simulation</title>
<p>We employ FluTe, an open-source and validated agent-based pandemic influenza simulation developed by Chao et al. (2010) [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]. FluTe’s contact network is composed of census tracts divided into communities of 500–3,000 individuals. Each community consists of randomly generated households of 1–7 individuals in one of the five age-groups: preschool children (0–4), school children (5–18), young adults (19–29), adults (30–64), and seniors (65 and over). Individuals can be members of multiple community-based mixing groups such as households, household clusters (composed of socially close households), neighborhoods and communities. The simulation has two time epochs for each day: day- and night-time. At night time, individuals can only make contacts within their community-based mixing groups, whereas they may contact other individuals in day-time if they share the same social mixing group, e.g., daycare, school, workplace.</p>
<p>In each time epoch, a contact for potential disease transmission between any two individuals sharing a mixing group is generated. During a contact between a susceptible and an infectious individual, influenza transmission may occur with a probability that depends on vaccine efficacy, virus load and symptoms of the infectious individual. Each infected person follows a predefined daily viral load profile representing the level of infectiousness on each day of the disease duration. Infected individuals may become symptomatic after an asymptomatic incubation period of one to three days [<xref rid="pone.0172261.ref036" ref-type="bibr">36</xref>
]. Symptomatic individuals are twice as infectious as the asymptomatic ones. Infected individuals recover and become immune after six days.</p>
<p>The vaccinated individuals have reduced likelihood of getting infected during a contact, becoming symptomatic when infected, and transmitting the disease [<xref rid="pone.0172261.ref037" ref-type="bibr">37</xref>
]. The vaccine efficacy reaches its maximum level in two weeks with exponential increments, and the maximum vaccine efficacy varies among the age groups. Due to incremental nature of vaccine efficacy, the timing of vaccine interventions affects success in containing influenza pandemics. FluTe allows administrating vaccines before (pre-vaccination) or after (reactive vaccination) the onset of the pandemic. We refer the reader to Chao et al. (2010) for further details about FluTe [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
].</p>
</sec>
<sec id="sec004"><title>Compartmental model</title>
<p>We develop a deterministic compartmental model that closely approximates the results of FluTe for Seattle under no vaccination. Similar to FluTe, we divide the population into five age groups, <italic>AG</italic>
= {preschool children (0–4), school children (5–18), young adults (19–29), adults (30–64), seniors (65+)}. Each age group <italic>i</italic>
∈ <italic>AG</italic>
includes vaccinated and unvaccinated individuals in five compartments: susceptible (S), exposed (E), infected (I), recovered (R), and dead (D). We denote the susceptible individuals in age group <italic>i</italic>
∈ <italic>AG</italic>
with vaccination status <italic>h</italic>
∈ <italic>H</italic>
= {(u)nvaccinated, (v)accinated} by <inline-formula id="pone.0172261.e001"><alternatives><graphic xlink:href="pone.0172261.e001.jpg" id="pone.0172261.e001g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M1"><mml:msubsup><mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
, for example. The exposed compartment (E) corresponds to the asymptomatic individuals in FluTe. In addition, the infected compartment (I) corresponds to the symptomatic individuals in FluTe. A proportion of asymptomatic individuals never develop symptoms in FluTe, therefore, we split the exposed compartment into two sub-compartments: those who eventually show disease symptoms (<italic>E</italic>
<sub>→</sub>
<italic>I</italic>
), and those who recover without showing symptoms (<italic>E</italic>
<sub>→</sub>
<italic>R</italic>
). The symptomatic individuals are twice as infectious as asymptomatic ones in both FluTe and the SEIR model.</p>
<p>The incidence rate of new infections in age group <italic>i</italic>
caused by infectious individuals in age group <italic>j</italic>
, denoted by λ<sub><italic>ij</italic>
</sub>
, is given by:
<disp-formula id="pone.0172261.e002"><alternatives><graphic xlink:href="pone.0172261.e002.jpg" id="pone.0172261.e002g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M2"><mml:mrow><mml:msub><mml:mo>λ</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:msub><mml:mo>Φ</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo>(</mml:mo>
<mml:msubsup><mml:mi>β</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mi>E</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mi>β</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mi>E</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:msub><mml:mi>N</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
</alternatives>
<label>(1)</label>
</disp-formula>
In <xref ref-type="disp-formula" rid="pone.0172261.e002">Eq (1)</xref>
, Φ<sub><italic>ij</italic>
</sub>
is the contact rate from age group <italic>i</italic>
to <italic>j</italic>
. Parameter <inline-formula id="pone.0172261.e003"><alternatives><graphic xlink:href="pone.0172261.e003.jpg" id="pone.0172261.e003g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M3"><mml:msubsup><mml:mi>β</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
(<inline-formula id="pone.0172261.e004"><alternatives><graphic xlink:href="pone.0172261.e004.jpg" id="pone.0172261.e004g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M4"><mml:msubsup><mml:mi>β</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
) denotes the transmission rate of unvaccinated (vaccinated) infectious individuals in age group <italic>j</italic>
given a single contact with a susceptible individual. Variables <inline-formula id="pone.0172261.e005"><alternatives><graphic xlink:href="pone.0172261.e005.jpg" id="pone.0172261.e005g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M5"><mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
(<inline-formula id="pone.0172261.e006"><alternatives><graphic xlink:href="pone.0172261.e006.jpg" id="pone.0172261.e006g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M6"><mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
) and <inline-formula id="pone.0172261.e007"><alternatives><graphic xlink:href="pone.0172261.e007.jpg" id="pone.0172261.e007g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M7"><mml:msubsup><mml:mi>E</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
(<inline-formula id="pone.0172261.e008"><alternatives><graphic xlink:href="pone.0172261.e008.jpg" id="pone.0172261.e008g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M8"><mml:msubsup><mml:mi>E</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
) represent the number of unvaccinated (vaccinated) infected and exposed individuals in age group <italic>j</italic>
, respectively. Note that <inline-formula id="pone.0172261.e009"><alternatives><graphic xlink:href="pone.0172261.e009.jpg" id="pone.0172261.e009g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M9"><mml:mrow><mml:msubsup><mml:mi>E</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>I</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>R</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
and <inline-formula id="pone.0172261.e010"><alternatives><graphic xlink:href="pone.0172261.e010.jpg" id="pone.0172261.e010g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M10"><mml:mrow><mml:msubsup><mml:mi>E</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>I</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>R</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
. Finally, <italic>N</italic>
<sub><italic>j</italic>
</sub>
is the size of age group <italic>j</italic>
, and ∑<sub><italic>j</italic>
∈<italic>AG</italic>
</sub>
<italic>N</italic>
<sub><italic>j</italic>
</sub>
= <italic>N</italic>
where <italic>N</italic>
denotes the total population size. The overall infection rate of individuals in age group <italic>i</italic>
is equal to λ<sub><italic>i</italic>
</sub>
= ∑<sub><italic>j</italic>
∈<italic>AG</italic>
</sub>
λ<sub><italic>ij</italic>
</sub>
. Note that vaccinated (<italic>v</italic>
) and unvaccinated (<italic>u</italic>
) compartments are interdependent because the infection rate λ<sub><italic>i</italic>
</sub>
depends both on the number of vaccinated and unvaccinated infectious individuals. <xref ref-type="fig" rid="pone.0172261.g001">Fig 1</xref>
depicts the transitions among the compartments, and the model equations are given by:
<disp-formula id="pone.0172261.e011"><alternatives><graphic xlink:href="pone.0172261.e011.jpg" id="pone.0172261.e011g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M11"><mml:mrow><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msubsup><mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo>=</mml:mo>
<mml:mo>-</mml:mo>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:msubsup><mml:mi>ϵ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub><mml:mo>λ</mml:mo>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msubsup><mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mspace width="2em"></mml:mspace>
<mml:mi>h</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>H</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:math>
</alternatives>
<label>(2a)</label>
</disp-formula>
<disp-formula id="pone.0172261.e012"><alternatives><graphic xlink:href="pone.0172261.e012.jpg" id="pone.0172261.e012g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M12"><mml:mrow><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>I</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo>=</mml:mo>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:msubsup><mml:mi>ϵ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub><mml:mo>λ</mml:mo>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msubsup><mml:mi>ω</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:msubsup><mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>-</mml:mo>
<mml:msubsup><mml:mi>τ</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>I</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mspace width="2em"></mml:mspace>
<mml:mi>h</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>H</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:math>
</alternatives>
<label>(2b)</label>
</disp-formula>
<disp-formula id="pone.0172261.e013"><alternatives><graphic xlink:href="pone.0172261.e013.jpg" id="pone.0172261.e013g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M13"><mml:mrow><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>R</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo>=</mml:mo>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:msubsup><mml:mi>ϵ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub><mml:mo>λ</mml:mo>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:msubsup><mml:mi>ω</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msubsup><mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>-</mml:mo>
<mml:msubsup><mml:mi>γ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>R</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mspace width="2em"></mml:mspace>
<mml:mi>h</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>H</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:math>
</alternatives>
<label>(2c)</label>
</disp-formula>
<disp-formula id="pone.0172261.e014"><alternatives><graphic xlink:href="pone.0172261.e014.jpg" id="pone.0172261.e014g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M14"><mml:mrow><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msubsup><mml:mi>I</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mi>τ</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>I</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>-</mml:mo>
<mml:mrow><mml:mo>(</mml:mo>
<mml:msubsup><mml:mi>ξ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mi>χ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mspace width="2em"></mml:mspace>
<mml:mi>h</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>H</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:math>
</alternatives>
<label>(2d)</label>
</disp-formula>
<disp-formula id="pone.0172261.e015"><alternatives><graphic xlink:href="pone.0172261.e015.jpg" id="pone.0172261.e015g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M15"><mml:mrow><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msubsup><mml:mi>R</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mi>ξ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mi>γ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:msubsup><mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>E</mml:mi>
<mml:mo>→</mml:mo>
</mml:msub>
<mml:mi>R</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mspace width="2em"></mml:mspace>
<mml:mi>h</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>H</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:math>
</alternatives>
<label>(2e)</label>
</disp-formula>
<disp-formula id="pone.0172261.e016"><alternatives><graphic xlink:href="pone.0172261.e016.jpg" id="pone.0172261.e016g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M16"><mml:mrow><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msub><mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo>=</mml:mo>
<mml:munder><mml:mo>∑</mml:mo>
<mml:mrow><mml:mi>h</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:munder>
<mml:msubsup><mml:mi>χ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mspace width="2em"></mml:mspace>
<mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:math>
</alternatives>
<label>(2f)</label>
</disp-formula>
</p>
<fig id="pone.0172261.g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0172261.g001</object-id>
<label>Fig 1</label>
<caption><title>Transition rates between compartments.</title>
</caption>
<graphic xlink:href="pone.0172261.g001"></graphic>
</fig>
<p>Parameter <inline-formula id="pone.0172261.e017"><alternatives><graphic xlink:href="pone.0172261.e017.jpg" id="pone.0172261.e017g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M17"><mml:mrow><mml:msubsup><mml:mi>ϵ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>∈</mml:mo>
<mml:mrow><mml:mo>[</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
denotes the reduction in the likelihood of becoming infected after vaccination, naturally <inline-formula id="pone.0172261.e018"><alternatives><graphic xlink:href="pone.0172261.e018.jpg" id="pone.0172261.e018g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M18"><mml:mrow><mml:msubsup><mml:mi>ϵ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
for unvaccinated individuals. Parameter <inline-formula id="pone.0172261.e019"><alternatives><graphic xlink:href="pone.0172261.e019.jpg" id="pone.0172261.e019g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M19"><mml:mrow><mml:msubsup><mml:mi>ω</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>∈</mml:mo>
<mml:mrow><mml:mo>[</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
denotes the proportion of exposed individuals who eventually become symptomatic. Parameters <inline-formula id="pone.0172261.e020"><alternatives><graphic xlink:href="pone.0172261.e020.jpg" id="pone.0172261.e020g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M20"><mml:msubsup><mml:mi>τ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
, <inline-formula id="pone.0172261.e021"><alternatives><graphic xlink:href="pone.0172261.e021.jpg" id="pone.0172261.e021g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M21"><mml:msubsup><mml:mi>γ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
, <inline-formula id="pone.0172261.e022"><alternatives><graphic xlink:href="pone.0172261.e022.jpg" id="pone.0172261.e022g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M22"><mml:msubsup><mml:mi>ξ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
, and <inline-formula id="pone.0172261.e023"><alternatives><graphic xlink:href="pone.0172261.e023.jpg" id="pone.0172261.e023g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M23"><mml:msubsup><mml:mi>χ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
denote the Exposed-to-Infected, Exposed-to-Recovered, Infected-to-Recovered, and Infected-to-Death transition rates, in that order. We set <inline-formula id="pone.0172261.e024"><alternatives><graphic xlink:href="pone.0172261.e024.jpg" id="pone.0172261.e024g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M24"><mml:mrow><mml:msubsup><mml:mi>γ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mn>1</mml:mn>
<mml:mrow><mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msubsup><mml:mi>ξ</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msubsup><mml:mi>τ</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
to ensure that the exposed individuals who recover without showing symptoms stay asymptomatic during the course of the disease. The <italic>D</italic>
compartment is included in the model for tracking the number of influenza-related deaths.</p>
<p>Let <italic>p</italic>
<sub><italic>i</italic>
</sub>
denote the vaccinated proportion of age group <italic>i</italic>
. Moreover, let <italic>b</italic>
<sub><italic>i</italic>
</sub>
denote the initial number of infecteds in age group <italic>i</italic>
. The boundary conditions ∀<italic>i</italic>
∈ <italic>AG</italic>
are given by <inline-formula id="pone.0172261.e025"><alternatives><graphic xlink:href="pone.0172261.e025.jpg" id="pone.0172261.e025g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M25"><mml:mrow><mml:msubsup><mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
, <inline-formula id="pone.0172261.e026"><alternatives><graphic xlink:href="pone.0172261.e026.jpg" id="pone.0172261.e026g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M26"><mml:mrow><mml:msubsup><mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow><mml:mo>(</mml:mo>
<mml:msub><mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
, <inline-formula id="pone.0172261.e027"><alternatives><graphic xlink:href="pone.0172261.e027.jpg" id="pone.0172261.e027g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M27"><mml:mrow><mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub><mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
, <inline-formula id="pone.0172261.e028"><alternatives><graphic xlink:href="pone.0172261.e028.jpg" id="pone.0172261.e028g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M28"><mml:mrow><mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub><mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
, <inline-formula id="pone.0172261.e029"><alternatives><graphic xlink:href="pone.0172261.e029.jpg" id="pone.0172261.e029g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M29"><mml:mrow><mml:msubsup><mml:mi>E</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mi>E</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msubsup><mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
. We solve the system of differential <xref ref-type="disp-formula" rid="pone.0172261.e011">Eq (2)</xref>
numerically using the fourth-order Runge-Kutta method [<xref rid="pone.0172261.ref038" ref-type="bibr">38</xref>
].</p>
</sec>
<sec id="sec005"><title>Formulating the optimization problem</title>
<p>The optimization aims to find an effective age-specific allocation of a given vaccine supply <italic>V</italic>
such that ∑<sub><italic>i</italic>
∈<italic>AG</italic>
</sub>
<italic>p</italic>
<sub><italic>i</italic>
</sub>
<italic>N</italic>
<sub><italic>i</italic>
</sub>
≤ <italic>V</italic>
. We consider four different performance measures:</p>
<list list-type="bullet"><list-item><p><italic>Total expected cost</italic>
is equal to the sum of vaccination, infection and mortality costs. Vaccination cost (<inline-formula id="pone.0172261.e030"><alternatives><graphic xlink:href="pone.0172261.e030.jpg" id="pone.0172261.e030g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M30"><mml:msubsup><mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>b</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
) includes the vaccine price, work time lost, and the cost of potential side effects. Infection cost, which is different for unvaccinated (<inline-formula id="pone.0172261.e031"><alternatives><graphic xlink:href="pone.0172261.e031.jpg" id="pone.0172261.e031g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M31"><mml:msubsup><mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
) and vaccinated (<inline-formula id="pone.0172261.e032"><alternatives><graphic xlink:href="pone.0172261.e032.jpg" id="pone.0172261.e032g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M32"><mml:msubsup><mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
) individuals, refers to the sum of medication, outpatient visits and hospitalization expenses. Mortality cost (<inline-formula id="pone.0172261.e033"><alternatives><graphic xlink:href="pone.0172261.e033.jpg" id="pone.0172261.e033g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M33"><mml:msubsup><mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
) stands for the terminal care expenses.</p>
</list-item>
<list-item><p><italic>Total number of infections</italic>
is equal to the number of individuals affected by the pandemic.</p>
</list-item>
<list-item><p><italic>Total number of deaths</italic>
is equal to the number of influenza-related deaths.</p>
</list-item>
<list-item><p><italic>Total years of life lost (YLL)</italic>
weighs each death with the expected remaining life time based on the U.S. life tables and the age distributions [<xref rid="pone.0172261.ref039" ref-type="bibr">39</xref>
, <xref rid="pone.0172261.ref040" ref-type="bibr">40</xref>
].</p>
</list-item>
</list>
<p>In particular, at any time during the course of the pandemic, the four performance measures are calculated as follows:
<disp-formula id="pone.0172261.e034"><alternatives><graphic xlink:href="pone.0172261.e034.jpg" id="pone.0172261.e034g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M34"><mml:mrow><mml:mtable columnalign="left"><mml:mtr columnalign="left"><mml:mtd columnalign="right"><mml:mrow><mml:mtext>Total cost </mml:mtext>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mi>T</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left"><mml:mrow><mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:munder><mml:mo>∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:msup><mml:mi>c</mml:mi>
<mml:mi>b</mml:mi>
</mml:msup>
<mml:msub><mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
<mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msubsup><mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
</mml:msubsup>
<mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msub><mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="right"><mml:mrow><mml:mtext>Total number of infections </mml:mtext>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mi>T</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left"><mml:mrow><mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:munder><mml:mo>∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="right"><mml:mrow><mml:mtext>Total number of deaths </mml:mtext>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mi>T</mml:mi>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left"><mml:mrow><mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:munder><mml:mo>∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msub><mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="right"><mml:mrow><mml:mtext>Total YLL </mml:mtext>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mi>T</mml:mi>
<mml:mi>Y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left"><mml:mrow><mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:munder><mml:mo>∑</mml:mo>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mo>∈</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:msub><mml:mi>Y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msub><mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
</alternatives>
</disp-formula>
where <italic>Y</italic>
<sub><italic>i</italic>
</sub>
denotes the YLL value of age group <italic>i</italic>
∈ <italic>AG</italic>
. Moreover, <inline-formula id="pone.0172261.e035"><alternatives><graphic xlink:href="pone.0172261.e035.jpg" id="pone.0172261.e035g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M35"><mml:msub><mml:mi>N</mml:mi>
<mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
(<inline-formula id="pone.0172261.e036"><alternatives><graphic xlink:href="pone.0172261.e036.jpg" id="pone.0172261.e036g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M36"><mml:msub><mml:mi>N</mml:mi>
<mml:msubsup><mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>v</mml:mi>
</mml:msubsup>
</mml:msub>
</mml:math>
</alternatives>
</inline-formula>
), and <inline-formula id="pone.0172261.e037"><alternatives><graphic xlink:href="pone.0172261.e037.jpg" id="pone.0172261.e037g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M37"><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msub><mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</alternatives>
</inline-formula>
represent the total number of infections among unvaccinated (vaccinated) individuals and influenza-related deaths in age group <italic>i</italic>
∈ <italic>AG</italic>
, respectively.</p>
</sec>
<sec id="sec006"><title>Calibration</title>
<p>We run FluTe using the population file for Seattle (around 560,000 residents) that is included in the software distribution package. Contact rates within mixing groups and infectious disease parameters of influenza are set based on the values estimated in Chao et al. (2010) so that attack rates are consistent with the 1957 Asian A(H2N2) and 2009 A(H1N1) influenza pandemics without vaccination [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]. We seed the model with 10 randomly generated infected people. The vaccinated individuals have 40% reduced probability of becoming infected, 40% reduced probability of becoming symptomatic given infection, and 67% reduced probability of transmitting infection. These values represent the effectiveness of a well-matched seasonal influenza vaccine [<xref rid="pone.0172261.ref041" ref-type="bibr">41</xref>
]. The vaccine is only 60% as effective in seniors as everyone else, since older people with weaker immune systems often have a lower immune response to influenza vaccine [<xref rid="pone.0172261.ref042" ref-type="bibr">42</xref>
].</p>
<p>The homogeneous mixing assumption of the compartmental model (referred to as SEIR model hereafter) results in faster and more diverse disease spread, whereas the infection follows a more tranquil pattern in FluTe, i.e., individuals can transmit the disease only to those in their contact list. For a fair comparison of the age-specific vaccine allocation strategies derived by using the SEIR model and FluTe, we calibrate the SEIR model so that the number of new infections for each day of the pandemic closely matches to the corresponding average outcome from FluTe over 100 replications under no vaccination. In particular, we calibrate the SEIR model by varying the contact rates (Φ<sub><italic>ij</italic>
</sub>
), transmission rates (<inline-formula id="pone.0172261.e038"><alternatives><graphic xlink:href="pone.0172261.e038.jpg" id="pone.0172261.e038g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M38"><mml:msubsup><mml:mi>β</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:math>
</alternatives>
</inline-formula>
), and initial number of infections (<italic>b</italic>
<sub><italic>i</italic>
</sub>
). As a goodness of fit measure, we use Pearson’s chi-square statistic (<italic>χ</italic>
<sup>2</sup>
measure), and employ the numerical optimization algorithm described in the following section to find parameters of the SEIR model that minimize the maximum <italic>χ</italic>
<sup>2</sup>
measure over all age groups.</p>
<p>We first perform the calibration for <italic>R</italic>
<sub>0</sub>
= 1.2. We then repeat the process for each <italic>R</italic>
<sub>0</sub>
value considered in our numerical analyses using the previous calibration results as the initial solution. We keep the contact rates the same as those found for <italic>R</italic>
<sub>0</sub>
= 1.2 (see <xref ref-type="supplementary-material" rid="pone.0172261.s001">S1 Table</xref>
) because initial tests show that further calibration of the contact rates for different <italic>R</italic>
<sub>0</sub>
values is not necessary to obtain good matches between the results of FluTe and the SEIR model. <xref ref-type="fig" rid="pone.0172261.g002">Fig 2a</xref>
depicts the cumulative number of infections in each day after the calibration process. The cumulative number of infections in different age groups are presented in <xref ref-type="fig" rid="pone.0172261.g002">Fig 2b–2f</xref>
. Furthermore, <xref ref-type="supplementary-material" rid="pone.0172261.s002">S2 Table</xref>
in the supplement shows the similarity between the age-specific attack rates in FluTe and in the SEIR model under no vaccination. The parameters of FluTe and the calibrated SEIR model are provided in <xref rid="pone.0172261.t001" ref-type="table">Table 1</xref>
. We do not present the parameters of FluTe related to network structure and virus load profile and refer the readers to Chao et al. (2010) for more details [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
].</p>
<fig id="pone.0172261.g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0172261.g002</object-id>
<label>Fig 2</label>
<caption><title>Cumulative number of infections in each age group of FluTe and the SEIR model after the calibration process for <italic>R</italic>
<sub>0</sub>
= 1.2 without vaccination.</title>
<p>(a) Total population. (b) Preschool children (0–4). (c) School children (5–18). (d) Young adults (19–29). (e) Adults (30–64). (f) Seniors (65+).</p>
</caption>
<graphic xlink:href="pone.0172261.g002"></graphic>
</fig>
<table-wrap id="pone.0172261.t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0172261.t001</object-id>
<label>Table 1</label>
<caption><title>Input parameters of FluTe and the SEIR model.</title>
</caption>
<alternatives><graphic id="pone.0172261.t001g" xlink:href="pone.0172261.t001"></graphic>
<table frame="box" rules="all" border="0"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<tbody><tr><td align="left" rowspan="1" colspan="1"></td>
<td align="center" rowspan="1" colspan="1">SEIR</td>
<td align="center" rowspan="1" colspan="1">FluTe</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="center" rowspan="1" colspan="1">Source</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Latent period duration (LPD)</td>
<td align="right" rowspan="1" colspan="1">1.9 days</td>
<td align="right" rowspan="1" colspan="1">1.9 days<xref ref-type="table-fn" rid="t001fn001"><sup>‡</sup>
</xref>
</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Total disease duration (TDD)</td>
<td align="right" rowspan="1" colspan="1">6 days</td>
<td align="right" rowspan="1" colspan="1">6 days</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Infectious period duration (IPD)</td>
<td align="right" rowspan="1" colspan="1">4.1 days</td>
<td align="right" rowspan="1" colspan="1">4.1 days</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">TDD−LPD</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Transition rate</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"> Exposed-to-Infectious(<italic>τ</italic>
<sup><italic>u</italic>
</sup>
)</td>
<td align="right" rowspan="1" colspan="1">0.52632</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(LPD)<sup>−1</sup>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"> Exposed-to-Recovered (<italic>γ</italic>
)</td>
<td align="right" rowspan="1" colspan="1">0.1667</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(TDD)<sup>−1</sup>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"> Infectious-to-Recovered (<italic>ξ</italic>
)</td>
<td align="right" rowspan="1" colspan="1">0.2439</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(IPD)<sup>−1</sup>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Likelihood</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"> of showing the symptoms (<italic>ω</italic>
<sup><italic>u</italic>
</sup>
)<xref ref-type="table-fn" rid="t001fn002"><sup>†</sup>
</xref>
</td>
<td align="right" rowspan="1" colspan="1">0.67</td>
<td align="right" rowspan="1" colspan="1">0.67</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="center" colspan="2" rowspan="1">Transmission rate/probability (<italic>β</italic>
<sup><italic>u</italic>
</sup>
)<xref ref-type="table-fn" rid="t001fn003">*</xref>
</td>
<td align="center" colspan="2" rowspan="1">Initial number of infections (<italic>b</italic>
)</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="center" rowspan="1" colspan="1">SEIR</td>
<td align="center" rowspan="1" colspan="1">FluTe</td>
<td align="center" rowspan="1" colspan="1">SEIR</td>
<td align="center" rowspan="1" colspan="1">FluTe</td>
<td align="right" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="7" colspan="1">Basic reproduction number (<italic>R</italic>
<sub>0</sub>
)</td>
<td align="center" rowspan="1" colspan="1">1.2</td>
<td align="right" rowspan="1" colspan="1">0.0044</td>
<td align="right" rowspan="1" colspan="1">0.2057</td>
<td align="right" rowspan="1" colspan="1">19.48</td>
<td align="right" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="7" colspan="1">Calibration and [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">1.4</td>
<td align="right" rowspan="1" colspan="1">0.0053</td>
<td align="right" rowspan="1" colspan="1">0.2435</td>
<td align="right" rowspan="1" colspan="1">22.07</td>
<td align="right" rowspan="1" colspan="1">10</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">1.6</td>
<td align="right" rowspan="1" colspan="1">0.0063</td>
<td align="right" rowspan="1" colspan="1">0.2812</td>
<td align="right" rowspan="1" colspan="1">7.88</td>
<td align="right" rowspan="1" colspan="1">10</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">1.8</td>
<td align="right" rowspan="1" colspan="1">0.0071</td>
<td align="right" rowspan="1" colspan="1">0.3190</td>
<td align="right" rowspan="1" colspan="1">5.01</td>
<td align="right" rowspan="1" colspan="1">10</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">2.0</td>
<td align="right" rowspan="1" colspan="1">0.0077</td>
<td align="right" rowspan="1" colspan="1">0.3567</td>
<td align="right" rowspan="1" colspan="1">5.01</td>
<td align="right" rowspan="1" colspan="1">10</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">2.2</td>
<td align="right" rowspan="1" colspan="1">0.0083</td>
<td align="right" rowspan="1" colspan="1">0.3944</td>
<td align="right" rowspan="1" colspan="1">5.01</td>
<td align="right" rowspan="1" colspan="1">10</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">2.4</td>
<td align="right" rowspan="1" colspan="1">0.0089</td>
<td align="right" rowspan="1" colspan="1">0.4322</td>
<td align="right" rowspan="1" colspan="1">5.01</td>
<td align="right" rowspan="1" colspan="1">10</td>
</tr>
<tr><td align="center" colspan="7" rowspan="1">Age group specific disease parameters</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="center" rowspan="1" colspan="1">0–4</td>
<td align="center" rowspan="1" colspan="1">5–18</td>
<td align="center" rowspan="1" colspan="1">19–29</td>
<td align="center" rowspan="1" colspan="1">30–64</td>
<td align="center" rowspan="1" colspan="1">65+</td>
<td align="right" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Population size</td>
<td align="right" rowspan="1" colspan="1">36,722</td>
<td align="right" rowspan="1" colspan="1">124,787</td>
<td align="right" rowspan="1" colspan="1">71,601</td>
<td align="right" rowspan="1" colspan="1">265,047</td>
<td align="right" rowspan="1" colspan="1">65,303</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Death rate (unvaccinated-<italic>χ</italic>
<sup><italic>u</italic>
</sup>
)</td>
<td align="right" rowspan="1" colspan="1">0.0049</td>
<td align="right" rowspan="1" colspan="1">0.002</td>
<td align="right" rowspan="1" colspan="1">0.0056</td>
<td align="right" rowspan="1" colspan="1">0.0038</td>
<td align="right" rowspan="1" colspan="1">0.001</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref013" ref-type="bibr">13</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Death rate (vaccinated-<italic>χ</italic>
<sup><italic>v</italic>
</sup>
)</td>
<td align="right" rowspan="1" colspan="1">0.0012</td>
<td align="right" rowspan="1" colspan="1">0.0005</td>
<td align="right" rowspan="1" colspan="1">0.0017</td>
<td align="right" rowspan="1" colspan="1">0.0011</td>
<td align="right" rowspan="1" colspan="1">0.0004</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref013" ref-type="bibr">13</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Vaccine efficacy in</td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"> becoming infected</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="right" rowspan="1" colspan="1">24%</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"> transmitting the disease</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="right" rowspan="1" colspan="1">24%</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"> showing the symptoms</td>
<td align="right" rowspan="1" colspan="1">67%</td>
<td align="right" rowspan="1" colspan="1">67%</td>
<td align="right" rowspan="1" colspan="1">67%</td>
<td align="right" rowspan="1" colspan="1">67%</td>
<td align="right" rowspan="1" colspan="1">40%</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
]</td>
</tr>
<tr><td align="center" colspan="7" rowspan="1">Performance related parameters</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="right" rowspan="1" colspan="1">0–4</td>
<td align="right" rowspan="1" colspan="1">5–18</td>
<td align="right" rowspan="1" colspan="1">19–29</td>
<td align="right" rowspan="1" colspan="1">30–64</td>
<td align="right" rowspan="1" colspan="1">65+</td>
<td align="right" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Infection cost (unvaccinated)($)</td>
<td align="right" rowspan="1" colspan="1">275.3</td>
<td align="right" rowspan="1" colspan="1">275.3</td>
<td align="right" rowspan="1" colspan="1">328.98</td>
<td align="right" rowspan="1" colspan="1">328.98</td>
<td align="right" rowspan="1" colspan="1">492.56</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref013" ref-type="bibr">13</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Infection cost (vaccinated)($)</td>
<td align="right" rowspan="1" colspan="1">231.58</td>
<td align="right" rowspan="1" colspan="1">231.58</td>
<td align="right" rowspan="1" colspan="1">264.71</td>
<td align="right" rowspan="1" colspan="1">264.71</td>
<td align="right" rowspan="1" colspan="1">404.54</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref013" ref-type="bibr">13</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Terminal care cost($)</td>
<td align="right" rowspan="1" colspan="1">3,435</td>
<td align="right" rowspan="1" colspan="1">3,435</td>
<td align="right" rowspan="1" colspan="1">7,605</td>
<td align="right" rowspan="1" colspan="1">7,605</td>
<td align="right" rowspan="1" colspan="1">8,309</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref013" ref-type="bibr">13</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Vaccination cost ($)</td>
<td align="right" rowspan="1" colspan="1">37.26</td>
<td align="right" rowspan="1" colspan="1">37.26</td>
<td align="right" rowspan="1" colspan="1">37.26</td>
<td align="right" rowspan="1" colspan="1">37.26</td>
<td align="right" rowspan="1" colspan="1">37.26</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref013" ref-type="bibr">13</xref>
]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">YLL (years)</td>
<td align="right" rowspan="1" colspan="1">79.42</td>
<td align="right" rowspan="1" colspan="1">70.71</td>
<td align="right" rowspan="1" colspan="1">57.92</td>
<td align="right" rowspan="1" colspan="1">36.45</td>
<td align="right" rowspan="1" colspan="1">13.37</td>
<td align="left" rowspan="1" colspan="1">[<xref rid="pone.0172261.ref039" ref-type="bibr">39</xref>
, <xref rid="pone.0172261.ref040" ref-type="bibr">40</xref>
]</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="t001fn001"><p><sup>‡</sup>
: Mean of a discrete distribution of 1, 2, or 3 days with probabilities 0.3, 0.5, and 0.2.</p>
</fn>
<fn id="t001fn002"><p><sup>†</sup>
: These numbers are for the unvaccinated individuals. For the vaccinated individuals, they should be multiplied by (1- vaccine efficacy in showing the symptoms).</p>
</fn>
<fn id="t001fn003"><p>*: These numbers are for the unvaccinated individuals. For the vaccinated individuals, they should be multiplied by (1- vaccine efficacy in transmitting the disease).</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec007"><title>Solution approach</title>
<p>The pandemic propagation is nonlinear because the incidence of new infections depends both on the current number of infectious and susceptible individuals. Moreover, the size and duration of outbreaks are uncertain in agent-based simulations like FluTe. All of these factors render traditional gradient-based optimization methods inapplicable. We therefore use a derivative-free approach, in particular, the mesh-adaptive direct search (MADS) algorithm as implemented in open-source software NOMAD [<xref rid="pone.0172261.ref043" ref-type="bibr">43</xref>
]. Starting from an initial solution, the MADS algorithm iteratively tries to improve the current best solution by generating trial points on a mesh, which is a discretization of the variable space. Each iteration is composed of two main steps: the search and the poll steps. The search step evaluates a number of trial mesh points. If an improved mesh point is found, then the next iteration is initiated with the new incumbent solution using a larger mesh size. Whenever the search step fails to generate an improved mesh point, then the poll step is invoked. The poll step explores the variable space near the current incumbent solution. If the poll step also fails to improve the current best solution, then the mesh size and poll size parameters are reduced in order to increase the search resolution. The MADS algorithm stops after a given number of iterations or when the mesh size reaches a precision limit. We refer the reader to Le Digabel (2011) for more information about the MADS algorithm and NOMAD [<xref rid="pone.0172261.ref044" ref-type="bibr">44</xref>
].</p>
<p>All numerical experiments are conducted using a PC with 48 cores (2.85 Ghz and 128 GB memory). We run both FluTe and the SEIR model up to one year. We generate 10 trial solutions in each iteration of the MADS algorithm. In one set of experiments, we use FluTe to evaluate the performance of trial solutions (FluTe+MADS), and in another set of experiments we use the SEIR model (SEIR+MADS). When using FluTe, we perform 24 replications to estimate the average performance of each trial solution. Based on our initial experiments, this sample size is sufficient to reduce the effects of sample variance on the results. We terminate the MADS algorithm after 1,000 trial solutions or when the mesh size is less than or equal to 10<sup>−13</sup>
. The SEIR model evaluates trial solutions much faster than FluTe. Therefore, when using the SEIR model, we terminate the MADS algorithm after 100,000 trial solutions or when the mesh size is less than or equal to 10<sup>−13</sup>
. We select <italic>p</italic>
<sub><italic>i</italic>
</sub>
= 0.5, ∀<italic>i</italic>
∈ <italic>AG</italic>
as the initial solution in the MADS algorithm, i.e., vaccinate 50% of the population in each age group. Although this solution may be infeasible under some vaccine coverage scenarios, it is still a proper initial solution because the MADS algorithm allows constraint violations in the intermediate iterations to diversify the search.</p>
</sec>
</sec>
<sec sec-type="results" id="sec008"><title>Results</title>
<p>We present the age-specific vaccine allocation strategies derived by FluTe+MADS and SEIR+MADS. We highlight the age groups prioritized by each approach and evaluate the relative performance of the proposed strategies under various scenarios for multiple objective functions.</p>
<p>The basic reproductive number (<italic>R</italic>
<sub>0</sub>
) represents the average number of infections generated by a typical infectious person in a completely susceptible population. If <italic>R</italic>
<sub>0</sub>
> 1, the infection spreads in the population. Otherwise, the infection eventually dies out without any intervention [<xref rid="pone.0172261.ref045" ref-type="bibr">45</xref>
]. Fraser et al. (2009) estimated the <italic>R</italic>
<sub>0</sub>
value of the 2009 H1N1 pandemic between 1.4 and 1.6 [<xref rid="pone.0172261.ref046" ref-type="bibr">46</xref>
]. Medlock and Galvani (2009) reported that 30% vaccine coverage can mitigate a pandemic like the 2009 H1N1 if there is no delay in response time (i.e., vaccination starts on the first day of the pandemic) [<xref rid="pone.0172261.ref013" ref-type="bibr">13</xref>
]. We therefore use <italic>R</italic>
<sub>0</sub>
= 1.6, 30% vaccine coverage, and no delay in response time as the base case in our experiments. Note that the result of Medlock and Galvani (2009) depends on vaccine efficacy. They assumed 80% vaccine efficacy against infection for 0 to 64-year-olds and 60% for those 65 and older. More than 30% coverage would be required for less effective vaccines (or a mismatched vaccine). Moreover, the whole process of producing a pandemic vaccine for a novel influenza virus takes four to six months [<xref rid="pone.0172261.ref047" ref-type="bibr">47</xref>
]. Therefore, a base case scenario with no response delay may seem unreasonable. However, we consider no other pandemic interventions such as quarantine, public health measures and antivirals, which may slow down the pandemic spread substantially. This effect may render the no response delay scenario more acceptable as a base case. In addition, analyzing the no response delay scenario is important to evaluate the maximum potential of a particular vaccination policy and coverage level to mitigate the pandemic.</p>
<p><xref rid="pone.0172261.t002" ref-type="table">Table 2</xref>
reports the age-specific vaccine allocation strategies derived by FluTe+MADS and SEIR+MADS in the base case. SEIR+MADS vaccinates only school children (5–18) to minimize the total cost (TC) objective. For other objectives, preschool (0–4) and school children are mainly vaccinated, and the remaining vaccines are allocated to young adults (19–29). FluTe+MADS vaccinates school children for the most part, and allocates the remaining vaccines to preschool children and young adults. Observe that SEIR+MADS uses fewer vaccine doses than FluTe+MADS for the TC objective, possibly because the effect of vaccination is more pronounced in the SEIR model as a result of the homogeneous mixing assumption.</p>
<table-wrap id="pone.0172261.t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0172261.t002</object-id>
<label>Table 2</label>
<caption><title>Vaccine allocation strategies obtained by FluTe+MADS and SEIR+MADS under the base-case scenario.</title>
</caption>
<alternatives><graphic id="pone.0172261.t002g" xlink:href="pone.0172261.t002"></graphic>
<table frame="box" rules="all" border="0"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="center" rowspan="2" colspan="1">Performance measure</th>
<th align="left" rowspan="2" colspan="1"></th>
<th align="center" colspan="5" rowspan="1">Vaccination fractions for each age group</th>
</tr>
<tr><th align="center" rowspan="1" colspan="1">0–4</th>
<th align="center" rowspan="1" colspan="1">5–18</th>
<th align="center" rowspan="1" colspan="1">19–29</th>
<th align="center" rowspan="1" colspan="1">30–64</th>
<th align="center" rowspan="1" colspan="1">65+</th>
</tr>
</thead>
<tbody><tr><td align="center" rowspan="2" colspan="1">TC</td>
<td align="center" rowspan="1" colspan="1">SEIR+MADS</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">87%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">FluTe+MADS</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">99%</td>
<td align="center" rowspan="1" colspan="1">34%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="2" colspan="1">TI</td>
<td align="center" rowspan="1" colspan="1">SEIR+MADS</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">FluTe+MADS</td>
<td align="center" rowspan="1" colspan="1">89%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
</tr>
<tr><td align="center" rowspan="2" colspan="1">TD</td>
<td align="center" rowspan="1" colspan="1">SEIR+MADS</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">FluTe+MADS</td>
<td align="center" rowspan="1" colspan="1">5%</td>
<td align="center" rowspan="1" colspan="1">97%</td>
<td align="center" rowspan="1" colspan="1">50%</td>
<td align="center" rowspan="1" colspan="1">3%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
</tr>
<tr><td align="center" rowspan="2" colspan="1">TY</td>
<td align="center" rowspan="1" colspan="1">SEIR+MADS</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">FluTe+MADS</td>
<td align="center" rowspan="1" colspan="1">21%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">50%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="t002fn001"><p><italic>R</italic>
<sub>0</sub>
= 1.6, 30% vaccine coverage, no delay in response time. <bold>TC:</bold>
Total cost, <bold>TI:</bold>
Total number of infections, <bold>TD:</bold>
Total number of deaths, <bold>TY:</bold>
Total YLL</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>In <xref ref-type="supplementary-material" rid="pone.0172261.s003">S1 Fig</xref>
, we compare the overall attack rates of the SEIR model and FluTe under the same vaccination policy: the available vaccine stocks are primarily allocated to school children, the remaining doses are first allocated to preschool children, and then to young adults. For each coverage level and <italic>R</italic>
<sub>0</sub>
value, the SEIR model has resulted in less overall attack rate than Flute, illustrating the more pronounced effect of vaccination in the former model. Furthermore, the critical vaccination fraction to stop the epidemic is lower in the SEIR model, and the amount of vaccine used relative to this threshold affects the optimal allocation. For example, if one can stop an epidemic by vaccinating only school children, there is no reason to vaccinate others.</p>
<p>We evaluate performances of the vaccine allocation strategies derived by FluTe+MADS and SEIR+MADS using FluTe with 100 replications. As seen in <xref rid="pone.0172261.t003" ref-type="table">Table 3</xref>
, the strategy from FluTe+MADS is significantly better than the one from SEIR+MADS for the TC objective. This is mainly due to the fact that the amount of vaccine used in the strategy from SEIR+MADS, while containing the disease effectively in the SEIR model, is insufficient to do so in FluTe.</p>
<table-wrap id="pone.0172261.t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0172261.t003</object-id>
<label>Table 3</label>
<caption><title>Objective values of recommended vaccine allocation strategies.</title>
</caption>
<alternatives><graphic id="pone.0172261.t003g" xlink:href="pone.0172261.t003"></graphic>
<table frame="box" rules="all" border="0"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="left" rowspan="2" colspan="1">Performance measure</th>
<th align="center" colspan="2" rowspan="1">Sample mean</th>
<th align="right" rowspan="2" colspan="1">Difference</th>
<th align="right" rowspan="2" colspan="1">p-value</th>
</tr>
<tr><th align="right" rowspan="1" colspan="1">SEIR+MADS</th>
<th align="right" rowspan="1" colspan="1">FluTe+MADS</th>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">TC ($M)</td>
<td align="char" char="." rowspan="1" colspan="1">23.2</td>
<td align="char" char="." rowspan="1" colspan="1">6.5</td>
<td align="char" char="." rowspan="1" colspan="1">16.7</td>
<td align="char" char="." rowspan="1" colspan="1"><0.001</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">TI (infections)</td>
<td align="char" char="." rowspan="1" colspan="1">1,170.3</td>
<td align="char" char="." rowspan="1" colspan="1">945.2</td>
<td align="char" char="." rowspan="1" colspan="1">225.1</td>
<td align="char" char="." rowspan="1" colspan="1">0.117</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">TD (deaths)</td>
<td align="char" char="." rowspan="1" colspan="1">10.7</td>
<td align="char" char="." rowspan="1" colspan="1">9.9</td>
<td align="char" char="." rowspan="1" colspan="1">0.8</td>
<td align="char" char="." rowspan="1" colspan="1">0.584</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">TY (life years lost)</td>
<td align="char" char="." rowspan="1" colspan="1">459.8</td>
<td align="char" char="." rowspan="1" colspan="1">367.5</td>
<td align="char" char="." rowspan="1" colspan="1">92.3</td>
<td align="char" char="." rowspan="1" colspan="1">0.094</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="t003fn001"><p><italic>R</italic>
<sub>0</sub>
= 1.6, 30% vaccine coverage, no delay in response time. The performance measures are calculated using FluTe with 100 replications.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<sec id="sec009"><title>Sensitivity analysis on R<sub>0</sub>
</title>
<p>We vary the basic reproduction number <italic>R</italic>
<sub>0</sub>
between 1.2 and 2.4 under 30% vaccine coverage with no delay in response time to analyze the sensitivity of the proposed vaccine allocation strategies to <italic>R</italic>
<sub>0</sub>
. For the total cost (TC) objective, SEIR+MADS increases the vaccinated proportion of school children from 47% to 100% as <italic>R</italic>
<sub>0</sub>
increases from 1.2 to 2.0 (see <xref ref-type="fig" rid="pone.0172261.g003">Fig 3a</xref>
). After covering school children, the remaining vaccines are allocated to preschool children when <italic>R</italic>
<sub>0</sub>
≥ 2.0. FluTe+MADS vaccinates 45% of school children as well as a small portion of preschool children and young adults when <italic>R</italic>
<sub>0</sub>
= 1.2 (see <xref ref-type="fig" rid="pone.0172261.g003">Fig 3b</xref>
). Furthermore, FluTe+MADS covers all school children when <italic>R</italic>
<sub>0</sub>
≥ 1.6. SEIR+MADS allocates 35% of the available vaccine when <italic>R</italic>
<sub>0</sub>
= 1.2 and all of the available vaccine when <italic>R</italic>
<sub>0</sub>
≥ 2.2. On the other hand, FluTe+MADS allocates 43% of the available vaccine when <italic>R</italic>
<sub>0</sub>
= 1.2 and all of the available vaccine when <italic>R</italic>
<sub>0</sub>
≥ 1.8. Intuitively, both methods use fewer vaccine doses for smaller <italic>R</italic>
<sub>0</sub>
values to keep the vaccination cost low.</p>
<fig id="pone.0172261.g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0172261.g003</object-id>
<label>Fig 3</label>
<caption><title>Vaccine allocation strategies derived by FluTe+MADS and SEIR+MADS under all objective functions for various R<sub>0</sub>
values (30% vaccine coverage, no delay in response time).</title>
<p>(a) SEIR+MADS with the TC objective. (b) FluTe+MADS with the TC objective. (c) SEIR+MADS with the TI objective. (d) FluTe+MADS with the TI objective. (e) SEIR+MADS with the TD objective. (f) FluTe+MADS with the TD objective. (g) SEIR+MADS with the TY objective. (h) FluTe+MADS with the TY objective.</p>
</caption>
<graphic xlink:href="pone.0172261.g003"></graphic>
</fig>
<p>For the total number of infections (TI) objective, SEIR+MADS recommends a similar strategy for each <italic>R</italic>
<sub>0</sub>
; vaccinate all preschool and school children as well as small portions of young adults and/or seniors (65+) (see <xref ref-type="fig" rid="pone.0172261.g003">Fig 3c</xref>
). However, vaccine allocation strategies from FluTe+MADS vary by <italic>R</italic>
<sub>0</sub>
(see <xref ref-type="fig" rid="pone.0172261.g003">Fig 3d</xref>
). These strategies cover all school children and allocate the remaining vaccines to preschool children and young adults when <italic>R</italic>
<sub>0</sub>
≥ 1.4. Note that both approaches use all available vaccine stocks, as vaccination cost is no longer a concern for the TI objective. This result is also valid for the total number of deaths (TD) and total YLL (TY) objectives.</p>
<p>For the TD and TY objectives, there is a trade off between reducing the number infections (by vaccinating school children) and reducing the casualties in high-mortality age groups (by vaccinating preschool children and young adults). SEIR+MADS again vaccinates all preschool and school children for all <italic>R</italic>
<sub>0</sub>
values (see <xref ref-type="fig" rid="pone.0172261.g003">Fig 3e and 3g</xref>
). FluTe+MADS covers all school children in addition to some proportions of preschool children and young adults when <italic>R</italic>
<sub>0</sub>
≤ 2.2 (see <xref ref-type="fig" rid="pone.0172261.g003">Fig 3f and 3h</xref>
). However, when <italic>R</italic>
<sub>0</sub>
= 2.4, FluTe+MADS vaccinates preschool children and young adults for the most part rather than school children. This unique case is due to high mortality rates of preschool children and young adults (see <xref rid="pone.0172261.t001" ref-type="table">Table 1</xref>
). The vaccine stocks become insufficient to contain the pandemic when <italic>R</italic>
<sub>0</sub>
= 2.4, therefore, the recommended strategy focuses on minimizing casualties in high-mortality age groups.</p>
<p>The vaccine allocation strategies from FluTe+MADS perform better than those from SEIR+MADS when <italic>R</italic>
<sub>0</sub>
is between 1.4 and 2.2, i.e., when the effective vaccine allocation becomes critical to control the pandemic (see <xref ref-type="supplementary-material" rid="pone.0172261.s004">S2 Fig</xref>
). However, SEIR+MADS performs better than FluTe+MADS or the difference is not significant when <italic>R</italic>
<sub>0</sub>
< 1.4 or <italic>R</italic>
<sub>0</sub>
> 2.2. Recall that we run the MADS algorithm for at most 1,000 trial solutions in FluTe+MADS, while we allow for 100,000 trial solutions in SEIR+MADS. Therefore, it is not unlikely that SEIR+MADS finds better strategies than FluTe+MADS, especially when a small vaccination level is enough to contain the infection (i.e., <italic>R</italic>
<sub>0</sub>
< 1.4). On the other hand, when the infection is very aggressive, i.e., <italic>R</italic>
<sub>0</sub>
> 2.2, 30% vaccine coverage in the base case is not enough for containing the pandemic, and therefore, there is limited room for improvement by optimizing the vaccine allocation. One exception to this pattern is for the TD and TY objectives where FluTe+MADS performs 7% and 9% better than SEIR+MADS, respectively, by prioritizing high-mortality age groups, i.e., preschool children and young adults when <italic>R</italic>
<sub>0</sub>
= 2.4.</p>
</sec>
<sec id="sec010"><title>Sensitivity analysis on vaccine coverage and response time</title>
<p>We set <italic>R</italic>
<sub>0</sub>
= 1.6 and analyze the sensitivity of the age-specific allocation strategies to vaccine coverage and response time. <xref rid="pone.0172261.t004" ref-type="table">Table 4</xref>
presents the strategies obtained under 20%, 30% and 40% vaccine coverage with no response delay (i.e., vaccination starts on the first day of the pandemic). SEIR+MADS vaccinates preschool and school children for the most part, and allocates the remaining vaccines to young adults for all objectives except the total cost (TC) objective for which 87% of school children is vaccinated under all coverage levels to reduce the vaccination cost. FluTe+MADS vaccinates school children, and allocates the remaining vaccines to preschool children and young adults.</p>
<table-wrap id="pone.0172261.t004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0172261.t004</object-id>
<label>Table 4</label>
<caption><title>Vaccine allocation strategies for different coverage scenarios.</title>
</caption>
<alternatives><graphic id="pone.0172261.t004g" xlink:href="pone.0172261.t004"></graphic>
<table frame="box" rules="all" border="0"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="right" rowspan="1" colspan="1"></th>
<th align="center" rowspan="1" colspan="1"></th>
<th align="center" colspan="5" rowspan="1">FluTe+MADS</th>
<th align="center" colspan="5" rowspan="1">SEIR+MADS</th>
</tr>
<tr><th align="right" rowspan="2" colspan="1"></th>
<th align="center" rowspan="2" colspan="1">Vaccine coverage</th>
<th align="center" colspan="5" rowspan="1">Vaccination fraction for each age group (<italic>p</italic>
<sub><italic>i</italic>
</sub>
)</th>
<th align="center" colspan="5" rowspan="1">Vaccination fraction for each age group (<italic>p</italic>
<sub><italic>i</italic>
</sub>
)</th>
</tr>
<tr><th align="center" rowspan="1" colspan="1">0–4</th>
<th align="center" rowspan="1" colspan="1">5–18</th>
<th align="center" rowspan="1" colspan="1">19–29</th>
<th align="center" rowspan="1" colspan="1">30–64</th>
<th align="center" rowspan="1" colspan="1">65+</th>
<th align="center" rowspan="1" colspan="1">0–4</th>
<th align="center" rowspan="1" colspan="1">5–18</th>
<th align="center" rowspan="1" colspan="1">19–29</th>
<th align="center" rowspan="1" colspan="1">30–64</th>
<th align="center" rowspan="1" colspan="1">65+</th>
</tr>
</thead>
<tbody><tr><td align="center" rowspan="3" colspan="1"><bold>TC</bold>
</td>
<td align="center" rowspan="1" colspan="1">20%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">88%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">87%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">30%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">99%</td>
<td align="center" rowspan="1" colspan="1">34%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">87%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">40%</td>
<td align="center" rowspan="1" colspan="1">34%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">22%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">3%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">87%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="3" colspan="1"><bold>TI</bold>
</td>
<td align="center" rowspan="1" colspan="1">20%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">87%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">90%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">30%</td>
<td align="center" rowspan="1" colspan="1">89%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">40%</td>
<td align="center" rowspan="1" colspan="1">64%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">99%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">84%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">3%</td>
</tr>
<tr><td align="center" rowspan="3" colspan="1"><bold>TD</bold>
</td>
<td align="center" rowspan="1" colspan="1">20%</td>
<td align="center" rowspan="1" colspan="1">15%</td>
<td align="center" rowspan="1" colspan="1">80%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">10%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">90%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">30%</td>
<td align="center" rowspan="1" colspan="1">5%</td>
<td align="center" rowspan="1" colspan="1">97%</td>
<td align="center" rowspan="1" colspan="1">50%</td>
<td align="center" rowspan="1" colspan="1">3%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">40%</td>
<td align="center" rowspan="1" colspan="1">6%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">96%</td>
<td align="center" rowspan="1" colspan="1">11%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">80%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="3" colspan="1"><bold>TY</bold>
</td>
<td align="center" rowspan="1" colspan="1">20%</td>
<td align="center" rowspan="1" colspan="1">34%</td>
<td align="center" rowspan="1" colspan="1">70%</td>
<td align="center" rowspan="1" colspan="1">6%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">12%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">90%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">30%</td>
<td align="center" rowspan="1" colspan="1">21%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">50%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">40%</td>
<td align="center" rowspan="1" colspan="1">64%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">84%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">95%</td>
<td align="center" rowspan="1" colspan="1">-</td>
<td align="center" rowspan="1" colspan="1">-</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="t004fn001"><p><italic>R</italic>
<sub>0</sub>
= 1.6, no response delay</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>The strategies from FluTe+MADS outperforms SEIR+MADS in all objective types when vaccine coverage is 30% (see <xref ref-type="supplementary-material" rid="pone.0172261.s005">S3 Fig</xref>
). However, strategies from SEIR+MADS for the total number of deaths (TD) and the total YLL (TY) objectives on average perform better than those derived by FluTe+MADS when the vaccine coverage is 20% or 40% (see <xref ref-type="supplementary-material" rid="pone.0172261.s005">S3 Fig</xref>
). If the vaccine coverage is 20%, there is a small set of effective vaccine allocation strategies that can contain the pandemic. On the other hand, if the vaccine coverage is 40%, the set of effective strategies is large with many local optimums. In both cases, a large number of trial solutions should be evaluated. FluTe+MADS may not find an effective strategy with 1,000 trial solutions when the vaccine coverage is 20%, and gets stuck at a local optimum when the vaccine coverage is 40%. Furthermore, the random noise around the objective values of each FluTe replication may make it difficult for the MADS algorithm to find an effective improvement direction.</p>
<p><xref rid="pone.0172261.t005" ref-type="table">Table 5</xref>
presents the allocation strategies for different response time scenarios with <italic>R</italic>
<sub>0</sub>
= 1.6 and 30% vaccine coverage. Prevaccination refers to the case where the vaccine is administered two weeks before the beginning of the pandemic so that it reaches its maximum effectiveness by the time the virus starts spreading. Others refer to the scenarios where the vaccination starts <italic>d</italic>
days after the onset of the outbreak (<italic>d</italic>
= 0, 10, 20, 30, 40, 60, 80, or 90). In these delayed response scenarios, the vaccine effectiveness will gradually increase and reach its maximum level in two weeks during the course of the pandemic both in the SEIR model and FluTe [<xref rid="pone.0172261.ref006" ref-type="bibr">6</xref>
].</p>
<table-wrap id="pone.0172261.t005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0172261.t005</object-id>
<label>Table 5</label>
<caption><title>Vaccine allocation strategies for different response time scenarios.</title>
</caption>
<alternatives><graphic id="pone.0172261.t005g" xlink:href="pone.0172261.t005"></graphic>
<table frame="box" rules="all" border="0"><colgroup span="1"><col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead><tr><th align="center" rowspan="1" colspan="1"></th>
<th align="center" rowspan="1" colspan="1"></th>
<th align="center" colspan="5" rowspan="1">SEIR+MADS</th>
<th align="center" colspan="5" rowspan="1">FluTe+MADS</th>
</tr>
<tr><th align="center" rowspan="2" colspan="1"></th>
<th align="center" rowspan="2" colspan="1">Response time (days)</th>
<th align="center" colspan="5" rowspan="1">Vaccination fraction for each age group (<italic>p</italic>
<sub><italic>i</italic>
</sub>
)</th>
<th align="center" colspan="5" rowspan="1">Vaccination fraction for each age group (<italic>p</italic>
<sub><italic>i</italic>
</sub>
)</th>
</tr>
<tr><th align="center" rowspan="1" colspan="1">0–4</th>
<th align="center" rowspan="1" colspan="1">5–18</th>
<th align="center" rowspan="1" colspan="1">19–29</th>
<th align="center" rowspan="1" colspan="1">30–64</th>
<th align="center" rowspan="1" colspan="1">65+</th>
<th align="center" rowspan="1" colspan="1">0–4</th>
<th align="center" rowspan="1" colspan="1">5–18</th>
<th align="center" rowspan="1" colspan="1">19–29</th>
<th align="center" rowspan="1" colspan="1">30–64</th>
<th align="center" rowspan="1" colspan="1">65+</th>
</tr>
</thead>
<tbody><tr><td align="center" rowspan="8" colspan="1"><bold>TC</bold>
</td>
<td align="center" rowspan="1" colspan="1">P</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">86%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">19%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">4%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">0</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">87%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">99%</td>
<td align="center" rowspan="1" colspan="1">34%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">10</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">89%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">77%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">15%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">5%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">20</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">91%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">46%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">31%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">4%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">40</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">96%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">96%</td>
<td align="center" rowspan="1" colspan="1">49%</td>
<td align="center" rowspan="1" colspan="1">5%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">60</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">23%</td>
<td align="center" rowspan="1" colspan="1">83%</td>
<td align="center" rowspan="1" colspan="1">71%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">8%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">80</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">48%</td>
<td align="center" rowspan="1" colspan="1">6%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">54%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">90</td>
<td align="center" rowspan="1" colspan="1">72%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">4%</td>
<td align="center" rowspan="1" colspan="1">5%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">55%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="8" colspan="1"><bold>TI</bold>
</td>
<td align="center" rowspan="1" colspan="1">P</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">16%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">43%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">0</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">89%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">10</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">67%</td>
<td align="center" rowspan="1" colspan="1">99%</td>
<td align="center" rowspan="1" colspan="1">16%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">20</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">7%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">39%</td>
<td align="center" rowspan="1" colspan="1">5%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">40</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">19%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">14%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">60</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">30%</td>
<td align="center" rowspan="1" colspan="1">93%</td>
<td align="center" rowspan="1" colspan="1">52%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">5%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">80</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">56%</td>
<td align="center" rowspan="1" colspan="1">65%</td>
<td align="center" rowspan="1" colspan="1">88%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">4%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">90</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">59%</td>
<td align="center" rowspan="1" colspan="1">40%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">36%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="8" colspan="1"><bold>TD</bold>
</td>
<td align="center" rowspan="1" colspan="1">P</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">13%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">17%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">0</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">5%</td>
<td align="center" rowspan="1" colspan="1">97%</td>
<td align="center" rowspan="1" colspan="1">50%</td>
<td align="center" rowspan="1" colspan="1">3%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">10</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">60%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">20</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">14%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">27%</td>
<td align="center" rowspan="1" colspan="1">6%</td>
<td align="center" rowspan="1" colspan="1">6%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">40</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">60%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">60</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">5%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">66%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">23%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">8%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">80</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">22%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">90</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">87%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">24%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="8" colspan="1"><bold>TY</bold>
</td>
<td align="center" rowspan="1" colspan="1">P</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">71%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">10%</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">6%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">0</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">21%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">50%</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">10</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">93%</td>
<td align="center" rowspan="1" colspan="1">98%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">4%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">20</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">36%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">37%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">4%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">40</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">1%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">59%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">60</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">1%</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">80</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">47%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="center" rowspan="1" colspan="1">90</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">9%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">100%</td>
<td align="center" rowspan="1" colspan="1">–</td>
<td align="center" rowspan="1" colspan="1">2%</td>
<td align="center" rowspan="1" colspan="1">–</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="t005fn001"><p><italic>R</italic>
<sub>0</sub>
= 1.6, 30% vaccine coverage. <bold>P</bold>
: Prevaccination.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>In <xref rid="pone.0172261.t005" ref-type="table">Table 5</xref>
, SEIR+MADS vaccinates preschool children and school children for all objectives except for the TC objective. More interestingly, the recommended allocation strategy is rather insensitive to response delay for the TI, TD and TY objectives. This result is due to the fact that the impact of vaccination is so pronounced in the SEIR model that the infection can still be contained by vaccinating high transmitters even after 90-day response delay.</p>
<p>The vaccine allocation strategies from FluTe+MADS, however, vary significantly with response time, although vaccinating the school children is still the main focus in most cases. For the total number of deaths (TD) objective, FluTe+MADS increases the vaccinated proportion of age groups who have high mortality risk as the response is further delayed. When the response delay is 80 days, for example, the vaccinated proportion of school children drops to 0% and young adults are fully covered. This result complies with the literature which suggests prioritizing vaccination of high transmitters (e.g., school children) earlier and prioritizing vaccination of those who have high mortality rate (e.g., young adults, preschool children) later in the pandemic [<xref rid="pone.0172261.ref014" ref-type="bibr">14</xref>
, <xref rid="pone.0172261.ref048" ref-type="bibr">48</xref>
, <xref rid="pone.0172261.ref049" ref-type="bibr">49</xref>
].</p>
<p>The strategies from FluTe+MADS generally outperform those from SEIR+MADS for different response time scenarios (see <xref ref-type="supplementary-material" rid="pone.0172261.s005">S3 Fig</xref>
). The performance difference mainly decreases as the response is further delayed, except for the TD objective, where significant differences are obtained by focusing on high-risk age groups even when the response delay is more than 60 days.</p>
</sec>
</sec>
<sec sec-type="conclusions" id="sec011"><title>Conclusions</title>
<p>Health policy makers commonly use agent-based simulations and compartmental models when evaluating and designing effective vaccine allocation strategies. Our study shines a light on the possible differences between the strategies obtained by these two approaches. In particular, we observe that age-specific vaccine allocation strategies derived by using a computationally expensive but realistic agent-based simulation and those derived by using a fast but stylized compartmental model may be different from each other. They, however, both recommend vaccinating school children for the most part, which complies with the literature [<xref rid="pone.0172261.ref019" ref-type="bibr">19</xref>
, <xref rid="pone.0172261.ref049" ref-type="bibr">49</xref>
].</p>
<p>The age-specific vaccine allocation strategies derived by using the agent-based simulation significantly outperforms those derived using the compartmental model especially for moderate levels of basic reproduction number (i.e., 1.4 ≤ <italic>R</italic>
<sub>0</sub>
≤ 2.2) when vaccine stocks are not very scarce. In other cases, either it is rather easy to control the pandemic (e.g., when <italic>R</italic>
<sub>0</sub>
= 1.2), so there are several strategies that can effectively control the infection, or vaccination is not sufficient to control the pandemic. In such extreme scenarios, the performance gap between the two approaches is small because there is limited room for improvement.</p>
<p>We also note the following two observations. First, SEIR+MADS tries to limit the number of infected individuals for all <italic>R</italic>
<sub>0</sub>
levels using almost the same vaccine allocation strategy for all objective functions under all scenarios. On the other hand, strategies from FluTe+MADS vary significantly for different objective functions and parameter scenarios seeking a balance between the number of infections and influenza-related deaths. Second, the impact of vaccination is very significant in SEIR+MADS possibly due to the homogeneous mixing assumption. Therefore, once evaluated by FluTe, the performance of the strategies from SEIR+MADS appears to be less effective compared those derived by FluTe+MADS. Similar observations are reported by other studies that evaluate the performances of strategies derived by compartmental models using agent-based simulations [<xref rid="pone.0172261.ref020" ref-type="bibr">20</xref>
].</p>
<p>Our observations summarized above are based on the comparison of a specific agent-based simulation and a compartmental model. However, FluTe is commonly used and well-received in the literature [<xref rid="pone.0172261.ref050" ref-type="bibr">50</xref>
–<xref rid="pone.0172261.ref052" ref-type="bibr">52</xref>
]. Therefore, although they are not directly generalizable, our observations are likely to hold for other agent-based simulations and compartmental models whose assumptions are similar to those analyzed in this study. The differences between the age-specific vaccine allocation strategies derived by SEIR+MADS and FluTe+MADS are possibly due to the combined effect of considering the heterogeneity in contact patterns and stochasticity in disease progression. However, we left measuring the individual effect of heterogeneity and stochasticity for future research as such an analysis requires a more simplified methodological setting.</p>
<p>There are a few limitations of our analysis. First, we examine age-specific vaccine allocation but do not consider other important risk factors, e.g., chronic medical conditions and pregnancy. Moreover, we use a deterministic numerical optimization algorithm, which does not incorporate the variance in the simulation replications when updating the search direction. Using a faster agent-based simulation, the number of replications can be increased to reduce the variance. Alternatively, ranking and selection methods can be applied to find the proper number of replications for each vaccine allocation strategy, which is left for future research.</p>
</sec>
<sec sec-type="supplementary-material" id="sec012"><title>Supporting information</title>
<supplementary-material content-type="local-data" id="pone.0172261.s001"><label>S1 Table</label>
<caption><title>The daily contact rate of age group i and age group j (Φ<sub>ij</sub>
) used in the SEIR model after the calibration.</title>
<p>We assume that the total number of daily contacts between age groups <italic>i</italic>
and <italic>j</italic>
is symmetric. To calculate the daily contact rate of someone in age group <italic>i</italic>
with age group <italic>j</italic>
, we divided the total number of daily contacts between <italic>i</italic>
and <italic>j</italic>
by the population of age group <italic>i</italic>
. So, the age-specific contact rates are asymmetric because the population of age groups <italic>i</italic>
and <italic>j</italic>
are different.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0172261.s001.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0172261.s002"><label>S2 Table</label>
<caption><title>Age specific influenza attack rates in Asian A(H2N2) pandemic and simulations of metropolitan Seattle using FluTe and the SEIR model for R<sub>0</sub>
= 1.6 without vaccination.</title>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0172261.s002.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0172261.s003"><label>S1 Fig</label>
<caption><title>Overall attack rates for different vaccine coverage levels and R<sub>0</sub>
values.</title>
<p>(a) FluTe. (b) SEIR model. Both Flute and the SEIR model evaluate the same vaccine allocation policy for all vaccine coverage levels. In this policy, available vaccine stocks are primarily allocated to school children, the remaining doses are first allocated to preschool children, and then to young adults.</p>
<p>(EPS)</p>
</caption>
<media xlink:href="pone.0172261.s003.eps"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0172261.s004"><label>S2 Fig</label>
<caption><title>Objective values of vaccine allocation strategies for different R<sub>0</sub>
values (30% vaccine coverage, no delay in response time).</title>
<p>(a) Total cost. (b) Number of infections. (c) Mortality. (d) Years of life lost.</p>
<p>(EPS)</p>
</caption>
<media xlink:href="pone.0172261.s004.eps"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0172261.s005"><label>S3 Fig</label>
<caption><title>Objective values of vaccine allocation strategies for different coverage scenarios (R<sub>0</sub>
= 1.6, no delay in response time) and response time scenarios (R<sub>0</sub>
= 1.6, 30% vaccine coverage).</title>
<p>(a) Total cost. (b) Number of infections. (c) Number of deaths. (d) Years of life lost.</p>
<p>(EPS)</p>
</caption>
<media xlink:href="pone.0172261.s005.eps"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ref-list><title>References</title>
<ref id="pone.0172261.ref001"><label>1</label>
<mixed-citation publication-type="journal"><collab>CDC</collab>
. <article-title>Estimates of Deaths Associated with Seasonal Influenza—United States, 1976–2007</article-title>
. <source>Morbidity and Mortality Weekly Report</source>
. <year>2010</year>
;<volume>59</volume>
(<issue>33</issue>
):<fpage>1057</fpage>
–<lpage>1062</lpage>
.<pub-id pub-id-type="pmid">20798667</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref002"><label>2</label>
<mixed-citation publication-type="other">CDC. Leading causes of death; 2014. Available from: <ext-link ext-link-type="uri" xlink:href="http://www.cdc.gov/nchs/fastats/leading-causes-of-death.htm">http://www.cdc.gov/nchs/fastats/leading-causes-of-death.htm</ext-link>
.</mixed-citation>
</ref>
<ref id="pone.0172261.ref003"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Molinari</surname>
<given-names>NAM</given-names>
</name>
, <name><surname>Ortega-Sanchez</surname>
<given-names>IR</given-names>
</name>
, <name><surname>Messonnier</surname>
<given-names>ML</given-names>
</name>
, <name><surname>Thompson</surname>
<given-names>WW</given-names>
</name>
, <name><surname>Wortley</surname>
<given-names>PM</given-names>
</name>
, <name><surname>Weintraub</surname>
<given-names>E</given-names>
</name>
, <etal>et al</etal>
<article-title>The annual impact of seasonal influenza in the U.S.: measuring disease burden and costs</article-title>
. <source>Vaccine</source>
. <year>2007</year>
;<volume>25</volume>
(<issue>27</issue>
):<fpage>5086</fpage>
–<lpage>5096</lpage>
. <pub-id pub-id-type="doi">10.1016/j.vaccine.2007.03.046</pub-id>
<pub-id pub-id-type="pmid">17544181</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref004"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Parvin</surname>
<given-names>JD</given-names>
</name>
, <name><surname>Moscona</surname>
<given-names>A</given-names>
</name>
, <name><surname>Pan</surname>
<given-names>W</given-names>
</name>
, <name><surname>Leider</surname>
<given-names>J</given-names>
</name>
, <name><surname>Palese</surname>
<given-names>P</given-names>
</name>
. <article-title>Measurement of the mutation rates of animal viruses: influenza A virus and poliovirus type 1</article-title>
. <source>Journal of Virology</source>
. <year>1986</year>
;<volume>59</volume>
(<issue>2</issue>
):<fpage>377</fpage>
–<lpage>383</lpage>
. <pub-id pub-id-type="pmid">3016304</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref005"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Taubenberger</surname>
<given-names>JK</given-names>
</name>
, <name><surname>Morens</surname>
<given-names>DM</given-names>
</name>
. <article-title>1918 Influenza: the mother of all pandemics</article-title>
. <source>Emerging Infectious Diseases</source>
. <year>2006</year>
;<volume>12</volume>
(<issue>1</issue>
):<fpage>15</fpage>
–<lpage>22</lpage>
. <pub-id pub-id-type="doi">10.3201/eid1201.050979</pub-id>
<pub-id pub-id-type="pmid">16494711</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref006"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Chao</surname>
<given-names>DL</given-names>
</name>
, <name><surname>Halloran</surname>
<given-names>ME</given-names>
</name>
, <name><surname>Obenchain</surname>
<given-names>VJ</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
. <article-title>FluTE, a publicly available stochastic influenza epidemic simulation model</article-title>
. <source>PLoS Computational Biology</source>
. <year>2010</year>
;<volume>6</volume>
(<issue>1</issue>
):<fpage>e1000656</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1000656</pub-id>
<pub-id pub-id-type="pmid">20126529</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref007"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Talbot</surname>
<given-names>HK</given-names>
</name>
, <name><surname>Zhu</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>Q</given-names>
</name>
, <name><surname>Williams</surname>
<given-names>JV</given-names>
</name>
, <name><surname>Thompson</surname>
<given-names>MG</given-names>
</name>
, <name><surname>Griffin</surname>
<given-names>MR</given-names>
</name>
. <article-title>Effectiveness of influenza vaccine for preventing laboratory-confirmed influenza hospitalizations in adults, 2011–2012 influenza season</article-title>
. <source>Clinical Infectious Diseases</source>
. <year>2013</year>
;<volume>56</volume>
(<issue>12</issue>
):<fpage>1774</fpage>
–<lpage>1777</lpage>
. <pub-id pub-id-type="doi">10.1093/cid/cit124</pub-id>
<pub-id pub-id-type="pmid">23449269</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref008"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>Nichol</surname>
<given-names>KL</given-names>
</name>
, <name><surname>Treanor</surname>
<given-names>JJ</given-names>
</name>
. <article-title>Vaccines for seasonal and pandemic influenza</article-title>
. <source>Journal of Infectious Diseases</source>
. <year>2006</year>
;<volume>194</volume>
(<issue>Supplement 2</issue>
):<fpage>S111</fpage>
–<lpage>S118</lpage>
. <pub-id pub-id-type="doi">10.1086/507544</pub-id>
<pub-id pub-id-type="pmid">17163383</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref009"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Hethcote</surname>
<given-names>HW</given-names>
</name>
. <article-title>The mathematics of infectious diseases</article-title>
. <source>SIAM Review</source>
. <year>2000</year>
;<volume>42</volume>
(<issue>4</issue>
):<fpage>599</fpage>
–<lpage>653</lpage>
. <pub-id pub-id-type="doi">10.1137/S0036144500371907</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref010"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Zaric</surname>
<given-names>GS</given-names>
</name>
, <name><surname>Brandeau</surname>
<given-names>ML</given-names>
</name>
. <article-title>Resource allocation for epidemic control over short time horizons</article-title>
. <source>Mathematical Biosciences</source>
. <year>2001</year>
;<volume>171</volume>
(<issue>1</issue>
):<fpage>33</fpage>
–<lpage>58</lpage>
. <pub-id pub-id-type="doi">10.1016/S0025-5564(01)00050-5</pub-id>
<pub-id pub-id-type="pmid">11325383</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref011"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Dasbach</surname>
<given-names>EJ</given-names>
</name>
, <name><surname>Elbasha</surname>
<given-names>EH</given-names>
</name>
, <name><surname>Insinga</surname>
<given-names>RP</given-names>
</name>
. <article-title>Mathematical models for predicting the epidemiologic and economic impact of vaccination against human papillomavirus infection and disease</article-title>
. <source>Epidemiologic Reviews</source>
. <year>2006</year>
;<volume>28</volume>
(<issue>1</issue>
):<fpage>88</fpage>
–<lpage>100</lpage>
. <pub-id pub-id-type="doi">10.1093/epirev/mxj006</pub-id>
<pub-id pub-id-type="pmid">16740585</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref012"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Huang</surname>
<given-names>G</given-names>
</name>
, <name><surname>Li</surname>
<given-names>L</given-names>
</name>
. <article-title>A mathematical model of infectious diseases</article-title>
. <source>Annals of Operations Research</source>
. <year>2009</year>
;<volume>168</volume>
(<issue>1</issue>
):<fpage>41</fpage>
–<lpage>80</lpage>
. <pub-id pub-id-type="doi">10.1007/s10479-008-0373-7</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref013"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Medlock</surname>
<given-names>J</given-names>
</name>
, <name><surname>Galvani</surname>
<given-names>AP</given-names>
</name>
. <article-title>Optimizing influenza vaccine distribution</article-title>
. <source>Science</source>
. <year>2009</year>
;<volume>325</volume>
(<issue>5948</issue>
):<fpage>1705</fpage>
–<lpage>1708</lpage>
. <pub-id pub-id-type="doi">10.1126/science.1175570</pub-id>
<pub-id pub-id-type="pmid">19696313</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref014"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Medlock</surname>
<given-names>J</given-names>
</name>
, <name><surname>Meyers</surname>
<given-names>LA</given-names>
</name>
, <name><surname>Galvani</surname>
<given-names>A</given-names>
</name>
. <article-title>Optimizing allocation for a delayed influenza vaccination campaign</article-title>
. <source>PLoS Currents Influenza</source>
. <year>2009</year>
;. <pub-id pub-id-type="doi">10.1371/currents.RRN1134</pub-id>
<pub-id pub-id-type="pmid">20033093</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref015"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Eames</surname>
<given-names>KT</given-names>
</name>
, <name><surname>Tilston</surname>
<given-names>NL</given-names>
</name>
, <name><surname>Brooks-Pollock</surname>
<given-names>E</given-names>
</name>
, <name><surname>Edmunds</surname>
<given-names>WJ</given-names>
</name>
. <article-title>Measured dynamic social contact patterns explain the spread of H1N1v influenza</article-title>
. <source>PLoS Computational Biology</source>
. <year>2012</year>
;<volume>8</volume>
(<issue>3</issue>
):<fpage>e1002425</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002425</pub-id>
<pub-id pub-id-type="pmid">22412366</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref016"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Araz</surname>
<given-names>OM</given-names>
</name>
, <name><surname>Galvani</surname>
<given-names>A</given-names>
</name>
, <name><surname>Meyers</surname>
<given-names>LA</given-names>
</name>
. <article-title>Geographic prioritization of distributing pandemic influenza vaccines</article-title>
. <source>Health Care Management Science</source>
. <year>2012</year>
;<volume>15</volume>
(<issue>3</issue>
):<fpage>175</fpage>
–<lpage>187</lpage>
. <pub-id pub-id-type="doi">10.1007/s10729-012-9199-6</pub-id>
<pub-id pub-id-type="pmid">22618029</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref017"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Patel</surname>
<given-names>R</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
, <name><surname>Halloran</surname>
<given-names>ME</given-names>
</name>
. <article-title>Finding optimal vaccination strategies for pandemic influenza using genetic algorithms</article-title>
. <source>Journal of Theoretical Biology</source>
. <year>2005</year>
;<volume>234</volume>
(<issue>2</issue>
):<fpage>201</fpage>
–<lpage>212</lpage>
. <pub-id pub-id-type="doi">10.1016/j.jtbi.2004.11.032</pub-id>
<pub-id pub-id-type="pmid">15757679</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref018"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Das</surname>
<given-names>TK</given-names>
</name>
, <name><surname>Savachkin</surname>
<given-names>AA</given-names>
</name>
, <name><surname>Zhu</surname>
<given-names>Y</given-names>
</name>
. <article-title>A large scale simulation model of pandemic influenza outbreaks for development of dynamic mitigation strategies</article-title>
. <source>IIE Transactions</source>
. <year>2008</year>
;<volume>40</volume>
(<issue>9</issue>
):<fpage>893</fpage>
–<lpage>905</lpage>
. <pub-id pub-id-type="doi">10.1080/07408170802165856</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref019"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Basta</surname>
<given-names>NE</given-names>
</name>
, <name><surname>Chao</surname>
<given-names>DL</given-names>
</name>
, <name><surname>Halloran</surname>
<given-names>ME</given-names>
</name>
, <name><surname>Matrajt</surname>
<given-names>L</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
. <article-title>Strategies for pandemic and seasonal influenza vaccination of schoolchildren in the United States</article-title>
. <source>American Journal of Epidemiology</source>
. <year>2009</year>
;<volume>170</volume>
(<issue>6</issue>
):<fpage>679</fpage>
–<lpage>686</lpage>
. <pub-id pub-id-type="doi">10.1093/aje/kwp237</pub-id>
<pub-id pub-id-type="pmid">19679750</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref020"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Lee</surname>
<given-names>BY</given-names>
</name>
, <name><surname>Brown</surname>
<given-names>ST</given-names>
</name>
, <name><surname>Korch</surname>
<given-names>GW</given-names>
</name>
, <name><surname>Cooley</surname>
<given-names>PC</given-names>
</name>
, <name><surname>Zimmerman</surname>
<given-names>RK</given-names>
</name>
, <name><surname>Wheaton</surname>
<given-names>WD</given-names>
</name>
, <etal>et al</etal>
<article-title>A computer simulation of vaccine prioritization, allocation, and rationing during the 2009 H1N1 influenza pandemic</article-title>
. <source>Vaccine</source>
. <year>2010</year>
;<volume>28</volume>
(<issue>31</issue>
):<fpage>4875</fpage>
–<lpage>4879</lpage>
. <pub-id pub-id-type="doi">10.1016/j.vaccine.2010.05.002</pub-id>
<pub-id pub-id-type="pmid">20483192</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref021"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>Uribe-Sánchez</surname>
<given-names>A</given-names>
</name>
, <name><surname>Savachkin</surname>
<given-names>A</given-names>
</name>
, <name><surname>Santana</surname>
<given-names>A</given-names>
</name>
, <name><surname>Prieto-Santa</surname>
<given-names>D</given-names>
</name>
, <name><surname>Das</surname>
<given-names>TK</given-names>
</name>
. <article-title>A predictive decision-aid methodology for dynamic mitigation of influenza pandemics</article-title>
. <source>OR Spectrum</source>
. <year>2011</year>
;<volume>33</volume>
(<issue>3</issue>
):<fpage>751</fpage>
–<lpage>786</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0172261.ref022"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Andradóttir</surname>
<given-names>S</given-names>
</name>
, <name><surname>Chiu</surname>
<given-names>W</given-names>
</name>
, <name><surname>Goldsman</surname>
<given-names>D</given-names>
</name>
, <name><surname>Lee</surname>
<given-names>ML</given-names>
</name>
. <article-title>Simulation of influenza propagation: Model development, parameter estimation, and mitigation strategies</article-title>
. <source>IIE Transactions on Healthcare Systems Engineering</source>
. <year>2014</year>
;<volume>4</volume>
(<issue>1</issue>
):<fpage>27</fpage>
–<lpage>48</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0172261.ref023"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>Dhamodharan</surname>
<given-names>A</given-names>
</name>
, <name><surname>Proano</surname>
<given-names>RA</given-names>
</name>
. <article-title>Determining the optimal vaccine vial size in developing countries: a Monte Carlo simulation approach</article-title>
. <source>Health Care Management Science</source>
. <year>2012</year>
;<volume>15</volume>
(<issue>3</issue>
):<fpage>188</fpage>
–<lpage>196</lpage>
. <pub-id pub-id-type="doi">10.1007/s10729-012-9200-4</pub-id>
<pub-id pub-id-type="pmid">22528136</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref024"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>Feng</surname>
<given-names>W</given-names>
</name>
, <name><surname>Kong</surname>
<given-names>N</given-names>
</name>
, <name><surname>Wan</surname>
<given-names>H</given-names>
</name>
. <article-title>A simulation study of cadaveric liver allocation with a single-score patient prioritization formula</article-title>
. <source>Journal of Simulation</source>
. <year>2013</year>
;<volume>7</volume>
(<issue>2</issue>
):<fpage>109</fpage>
–<lpage>125</lpage>
. <pub-id pub-id-type="doi">10.1057/jos.2012.21</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref025"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Chhatwal</surname>
<given-names>J</given-names>
</name>
, <name><surname>He</surname>
<given-names>T</given-names>
</name>
. <article-title>Economic evaluations with agent-based modelling: an introduction</article-title>
. <source>Pharmacoeconomics</source>
. <year>2015</year>
;<volume>33</volume>
(<issue>5</issue>
):<fpage>423</fpage>
–<lpage>433</lpage>
. <pub-id pub-id-type="doi">10.1007/s40273-015-0254-2</pub-id>
<pub-id pub-id-type="pmid">25609398</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref026"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Rahmandad</surname>
<given-names>H</given-names>
</name>
, <name><surname>Sterman</surname>
<given-names>J</given-names>
</name>
. <article-title>Heterogeneity and Network Structure in the Dynamics of Diffusion: Comapring Agent-Based and Differential Equation Models</article-title>
. <source>Management Science</source>
. <year>2008</year>
;<volume>54</volume>
(<issue>5</issue>
):<fpage>998</fpage>
–<lpage>1014</lpage>
. <pub-id pub-id-type="doi">10.1287/mnsc.1070.0787</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref027"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Germann</surname>
<given-names>TC</given-names>
</name>
, <name><surname>Kadau</surname>
<given-names>K</given-names>
</name>
, <name><surname>L</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
, <name><surname>Macken</surname>
<given-names>CA</given-names>
</name>
. <article-title>Mitigation strategies for pandemic influenza in the United States</article-title>
. <source>Proceedings of the National Academy of Sciences</source>
. <year>2006</year>
;<volume>103</volume>
(<issue>15</issue>
):<fpage>5935</fpage>
–<lpage>5640</lpage>
. <pub-id pub-id-type="doi">10.1073/pnas.0601266103</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref028"><label>28</label>
<mixed-citation publication-type="book"><name><surname>Anderson</surname>
<given-names>RM</given-names>
</name>
, <name><surname>May</surname>
<given-names>RM</given-names>
</name>
, <name><surname>Anderson</surname>
<given-names>B</given-names>
</name>
. <source>Infectious diseases of humans: dynamics and control</source>
. <publisher-loc>Oxford</publisher-loc>
: <publisher-name>Oxford University Press</publisher-name>
; <year>1991</year>
.</mixed-citation>
</ref>
<ref id="pone.0172261.ref029"><label>29</label>
<mixed-citation publication-type="journal"><name><surname>Mollison</surname>
<given-names>D</given-names>
</name>
, <name><surname>Isham</surname>
<given-names>V</given-names>
</name>
, <name><surname>Grenfell</surname>
<given-names>B</given-names>
</name>
. <article-title>Epidemics: models and data</article-title>
. <source>Journal of the Royal Statistical Society Series A (Statistics in Society)</source>
. <year>1994</year>
; p. <fpage>115</fpage>
–<lpage>149</lpage>
. <pub-id pub-id-type="doi">10.2307/2983509</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref030"><label>30</label>
<mixed-citation publication-type="journal"><name><surname>Bansal</surname>
<given-names>S</given-names>
</name>
, <name><surname>Grenfell</surname>
<given-names>BT</given-names>
</name>
, <name><surname>Meyers</surname>
<given-names>LA</given-names>
</name>
. <article-title>When individual behaviour matters: homogeneous and network models in epidemiology</article-title>
. <source>Journal of the Royal Society Interface</source>
. <year>2007</year>
;<volume>4</volume>
(<issue>16</issue>
):<fpage>879</fpage>
–<lpage>891</lpage>
. <pub-id pub-id-type="doi">10.1098/rsif.2007.1100</pub-id>
<pub-id pub-id-type="pmid">17640863</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref031"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>Lloyd</surname>
<given-names>AL</given-names>
</name>
, <name><surname>May</surname>
<given-names>RM</given-names>
</name>
. <article-title>Spatial heterogeneity in epidemic models</article-title>
. <source>Journal of theoretical biology</source>
. <year>1996</year>
;<volume>179</volume>
(<issue>1</issue>
):<fpage>1</fpage>
–<lpage>11</lpage>
. <pub-id pub-id-type="doi">10.1006/jtbi.1996.0042</pub-id>
<pub-id pub-id-type="pmid">8733427</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref032"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
, <name><surname>Cummings</surname>
<given-names>DA</given-names>
</name>
, <name><surname>Fraser</surname>
<given-names>C</given-names>
</name>
, <name><surname>Cajka</surname>
<given-names>JC</given-names>
</name>
, <name><surname>Cooley</surname>
<given-names>PC</given-names>
</name>
, <name><surname>Burke</surname>
<given-names>DS</given-names>
</name>
. <article-title>Strategies for mitigating an influenza pandemic</article-title>
. <source>Nature</source>
. <year>2006</year>
;<volume>442</volume>
(<issue>7101</issue>
):<fpage>448</fpage>
–<lpage>452</lpage>
. <pub-id pub-id-type="doi">10.1038/nature04795</pub-id>
<pub-id pub-id-type="pmid">16642006</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref033"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>Merler</surname>
<given-names>S</given-names>
</name>
, <name><surname>Ajelli</surname>
<given-names>M</given-names>
</name>
, <name><surname>Rizzo</surname>
<given-names>C</given-names>
</name>
. <article-title>Age-prioritized use of antivirals during an influenza pandemic</article-title>
. <source>BMC Infectious Diseases</source>
. <year>2009</year>
;<volume>9</volume>
:<fpage>117</fpage>
–<lpage>117</lpage>
. <pub-id pub-id-type="doi">10.1186/1471-2334-9-117</pub-id>
<pub-id pub-id-type="pmid">19638194</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref034"><label>34</label>
<mixed-citation publication-type="journal"><name><surname>Hladish</surname>
<given-names>T</given-names>
</name>
, <name><surname>Melamud</surname>
<given-names>E</given-names>
</name>
, <name><surname>Barrera</surname>
<given-names>LA</given-names>
</name>
, <name><surname>Galvani</surname>
<given-names>A</given-names>
</name>
, <name><surname>Meyers</surname>
<given-names>LA</given-names>
</name>
. <article-title>EpiFire: An open source C++ library and application for contact network epidemiology</article-title>
. <source>BMC Bioinformatics</source>
. <year>2012</year>
;<volume>13</volume>
(<issue>76</issue>
):<fpage>1</fpage>
–<lpage>12</lpage>
. <pub-id pub-id-type="doi">10.1186/1471-2105-13-76</pub-id>
<pub-id pub-id-type="pmid">22214541</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref035"><label>35</label>
<mixed-citation publication-type="other">Audet C, Le Digabel S, Tribes C. NOMAD user guide; 2009. Available from: <ext-link ext-link-type="uri" xlink:href="https://www.gerad.ca/nomad/Downloads/user_guide.pdf">https://www.gerad.ca/nomad/Downloads/user_guide.pdf</ext-link>
.</mixed-citation>
</ref>
<ref id="pone.0172261.ref036"><label>36</label>
<mixed-citation publication-type="journal"><name><surname>Carrat</surname>
<given-names>F</given-names>
</name>
, <name><surname>Vergu</surname>
<given-names>E</given-names>
</name>
, <name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
, <name><surname>Lemaitre</surname>
<given-names>M</given-names>
</name>
, <name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <name><surname>Leach</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
<article-title>Time lines of infection and disease in human influenza: a review of volunteer challenge studies</article-title>
. <source>American Journal of Epidemiology</source>
. <year>2008</year>
;<volume>167</volume>
(<issue>7</issue>
):<fpage>775</fpage>
–<lpage>785</lpage>
. <pub-id pub-id-type="doi">10.1093/aje/kwm375</pub-id>
<pub-id pub-id-type="pmid">18230677</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref037"><label>37</label>
<mixed-citation publication-type="book"><name><surname>Struchiner</surname>
<given-names>CJ</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
, <name><surname>Halloran</surname>
<given-names>ME</given-names>
</name>
. <source>Design and analysis of vaccine studies</source>
. <publisher-loc>New York</publisher-loc>
: <publisher-name>Springer</publisher-name>
; <year>2010</year>
.</mixed-citation>
</ref>
<ref id="pone.0172261.ref038"><label>38</label>
<mixed-citation publication-type="book"><name><surname>Atkinson</surname>
<given-names>KA</given-names>
</name>
. <source>An introduction to numerical analysis (2nd ed.)</source>
. <publisher-loc>New York</publisher-loc>
: <publisher-name>John Wiley & Sons</publisher-name>
; <year>1989</year>
.</mixed-citation>
</ref>
<ref id="pone.0172261.ref039"><label>39</label>
<mixed-citation publication-type="journal"><name><surname>Arias</surname>
<given-names>E</given-names>
</name>
. <article-title>United States life tables, 2010</article-title>
. <source>National Vital Statistics Reports</source>
. <year>2014</year>
;<volume>63</volume>
(<issue>7</issue>
):<fpage>1</fpage>
–<lpage>63</lpage>
. <pub-id pub-id-type="pmid">25383611</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref040"><label>40</label>
<mixed-citation publication-type="other">United States Census Bureau. Population estimates; 2014.</mixed-citation>
</ref>
<ref id="pone.0172261.ref041"><label>41</label>
<mixed-citation publication-type="journal"><name><surname>Basta</surname>
<given-names>NE</given-names>
</name>
, <name><surname>Halloran</surname>
<given-names>ME</given-names>
</name>
, <name><surname>Matrajt</surname>
<given-names>L</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
. <article-title>Estimating influenza vaccine efficacy from challenge and community-based study data</article-title>
. <source>American Journal of Epidemiology</source>
. <year>2008</year>
;<volume>168</volume>
(<issue>12</issue>
):<fpage>1343</fpage>
–<lpage>1352</lpage>
. <pub-id pub-id-type="doi">10.1093/aje/kwn259</pub-id>
<pub-id pub-id-type="pmid">18974084</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref042"><label>42</label>
<mixed-citation publication-type="other">CDC. Vaccine effectiveness—How well does the flu vaccine work?; 2016.</mixed-citation>
</ref>
<ref id="pone.0172261.ref043"><label>43</label>
<mixed-citation publication-type="other">Abramson MA, Audet C, Couture G, Dennis, Jr JE, Le Digabel S, Tribes C. The NOMAD project; 2009.</mixed-citation>
</ref>
<ref id="pone.0172261.ref044"><label>44</label>
<mixed-citation publication-type="journal"><name><surname>Le Digabel</surname>
<given-names>S</given-names>
</name>
. <article-title>Algorithm 909: NOMAD: Nonlinear Optimization with the MADS algorithm</article-title>
. <source>ACM Transactions on Mathematical Software</source>
. <year>2011</year>
;<volume>37</volume>
(<issue>4</issue>
):<fpage>1</fpage>
–<lpage>15</lpage>
. <pub-id pub-id-type="doi">10.1145/1916461.1916468</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref045"><label>45</label>
<mixed-citation publication-type="journal"><name><surname>Heffernan</surname>
<given-names>J</given-names>
</name>
, <name><surname>Smith</surname>
<given-names>R</given-names>
</name>
, <name><surname>Wahl</surname>
<given-names>L</given-names>
</name>
. <article-title>Perspectives on the basic reproductive ratio</article-title>
. <source>Journal of the Royal Society Interface</source>
. <year>2005</year>
;<volume>2</volume>
(<issue>4</issue>
):<fpage>281</fpage>
–<lpage>293</lpage>
. <pub-id pub-id-type="doi">10.1098/rsif.2005.0042</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref046"><label>46</label>
<mixed-citation publication-type="journal"><name><surname>Fraser</surname>
<given-names>C</given-names>
</name>
, <name><surname>Donnelly</surname>
<given-names>CA</given-names>
</name>
, <name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
<article-title>Pandemic potential of a strain of influenza A (H1N1): early findings</article-title>
. <source>Science</source>
. <year>2009</year>
;<volume>324</volume>
(<issue>5934</issue>
):<fpage>1557</fpage>
–<lpage>1561</lpage>
. <pub-id pub-id-type="doi">10.1126/science.1176062</pub-id>
<pub-id pub-id-type="pmid">19433588</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref047"><label>47</label>
<mixed-citation publication-type="other">WHO. Pandemic influenza vaccine manufacturing process and timeline; 2009. Available from: <ext-link ext-link-type="uri" xlink:href="http://www.who.int/csr/disease/swineflu/notes/h1n1_vaccine_20090806/en/">www.who.int/csr/disease/swineflu/notes/h1n1_vaccine_20090806/en/</ext-link>
.</mixed-citation>
</ref>
<ref id="pone.0172261.ref048"><label>48</label>
<mixed-citation publication-type="journal"><name><surname>Matrajt</surname>
<given-names>L</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
. <article-title>Optimizing vaccine allocation at different points in time during an epidemic</article-title>
. <source>PloS ONE</source>
. <year>2010</year>
;<volume>5</volume>
(<issue>11</issue>
):<fpage>e13767</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0013767</pub-id>
<pub-id pub-id-type="pmid">21085686</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref049"><label>49</label>
<mixed-citation publication-type="journal"><name><surname>Mylius</surname>
<given-names>SD</given-names>
</name>
, <name><surname>Hagenaars</surname>
<given-names>TJ</given-names>
</name>
, <name><surname>Lugnér</surname>
<given-names>AK</given-names>
</name>
, <name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
. <article-title>Optimal allocation of pandemic influenza vaccine depends on age, risk and timing</article-title>
. <source>Vaccine</source>
. <year>2008</year>
;<volume>26</volume>
(<issue>29</issue>
):<fpage>3742</fpage>
–<lpage>3749</lpage>
. <pub-id pub-id-type="doi">10.1016/j.vaccine.2008.04.043</pub-id>
<pub-id pub-id-type="pmid">18524428</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref050"><label>50</label>
<mixed-citation publication-type="journal"><name><surname>Cowling</surname>
<given-names>BJ</given-names>
</name>
, <name><surname>Ng</surname>
<given-names>S</given-names>
</name>
, <name><surname>Ma</surname>
<given-names>ES</given-names>
</name>
, <etal>et al</etal>
<article-title>Protective efficacy of seasonal influenza vaccination against seasonal and pandemic influenza virus infection during 2009 in Hong Kong</article-title>
. <source>Clinical Infectious Diseases</source>
. <year>2010</year>
;<volume>51</volume>
(<issue>12</issue>
):<fpage>1370</fpage>
–<lpage>1379</lpage>
. <pub-id pub-id-type="doi">10.1086/657311</pub-id>
<pub-id pub-id-type="pmid">21067351</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref051"><label>51</label>
<mixed-citation publication-type="journal"><name><surname>Lazer</surname>
<given-names>D</given-names>
</name>
, <name><surname>Kennedy</surname>
<given-names>R</given-names>
</name>
, <name><surname>King</surname>
<given-names>G</given-names>
</name>
, <name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
. <article-title>The parable of Google flu: traps in big data analysis</article-title>
. <source>Science</source>
. <year>2014</year>
;<volume>343</volume>
(<issue>6176</issue>
):<fpage>1203</fpage>
–<lpage>1205</lpage>
. <pub-id pub-id-type="doi">10.1126/science.1248506</pub-id>
<pub-id pub-id-type="pmid">24626916</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0172261.ref052"><label>52</label>
<mixed-citation publication-type="journal"><name><surname>Pastor-Satorras</surname>
<given-names>R</given-names>
</name>
, <name><surname>Castellano</surname>
<given-names>C</given-names>
</name>
, <name><surname>Van Mieghem</surname>
<given-names>P</given-names>
</name>
, <name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
. <article-title>Epidemic processes in complex networks</article-title>
. <source>Reviews of Modern Physics</source>
. <year>2015</year>
;<volume>87</volume>
(<issue>3</issue>
):<fpage>925</fpage>
–<lpage>979</lpage>
. <pub-id pub-id-type="doi">10.1103/RevModPhys.87.925</pub-id>
</mixed-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 000C50 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000C50 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= PandemieGrippaleV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:5319753 |texte= Deriving effective vaccine allocation strategies for pandemic influenza: Comparison of an agent-based simulation and a compartmental model }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:28222123" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a PandemieGrippaleV1
![]() | This area was generated with Dilib version V0.6.34. | ![]() |