Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Enhancement of Collective Immunity in Tokyo Metropolitan Area by Selective Vaccination against an Emerging Influenza Pandemic</title>
<author><name sortKey="Saito, Masaya M" sort="Saito, Masaya M" uniqKey="Saito M" first="Masaya M." last="Saito">Masaya M. Saito</name>
<affiliation><nlm:aff id="aff1"><addr-line>Institute of Statistical Mathematics, Tachikawa, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Imoto, Seiya" sort="Imoto, Seiya" uniqKey="Imoto S" first="Seiya" last="Imoto">Seiya Imoto</name>
<affiliation><nlm:aff id="aff2"><addr-line>Human Genome Center, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Yamaguchi, Rui" sort="Yamaguchi, Rui" uniqKey="Yamaguchi R" first="Rui" last="Yamaguchi">Rui Yamaguchi</name>
<affiliation><nlm:aff id="aff2"><addr-line>Human Genome Center, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Tsubokura, Masaharu" sort="Tsubokura, Masaharu" uniqKey="Tsubokura M" first="Masaharu" last="Tsubokura">Masaharu Tsubokura</name>
<affiliation><nlm:aff id="aff3"><addr-line>Division of Social Communication Systems for Advanced Clinical Research, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kami, Masahiro" sort="Kami, Masahiro" uniqKey="Kami M" first="Masahiro" last="Kami">Masahiro Kami</name>
<affiliation><nlm:aff id="aff3"><addr-line>Division of Social Communication Systems for Advanced Clinical Research, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Nakada, Haruka" sort="Nakada, Haruka" uniqKey="Nakada H" first="Haruka" last="Nakada">Haruka Nakada</name>
<affiliation><nlm:aff id="aff3"><addr-line>Division of Social Communication Systems for Advanced Clinical Research, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Sato, Hiroki" sort="Sato, Hiroki" uniqKey="Sato H" first="Hiroki" last="Sato">Hiroki Sato</name>
<affiliation><nlm:aff id="aff4"><addr-line>Department of Preventive Medicine and Public Health, National Defense Medical College, Tokorozawa, Saitama, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Miyano, Satoru" sort="Miyano, Satoru" uniqKey="Miyano S" first="Satoru" last="Miyano">Satoru Miyano</name>
<affiliation><nlm:aff id="aff2"><addr-line>Human Genome Center, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Higuchi, Tomoyuki" sort="Higuchi, Tomoyuki" uniqKey="Higuchi T" first="Tomoyuki" last="Higuchi">Tomoyuki Higuchi</name>
<affiliation><nlm:aff id="aff1"><addr-line>Institute of Statistical Mathematics, Tachikawa, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">24058445</idno>
<idno type="pmc">3776821</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3776821</idno>
<idno type="RBID">PMC:3776821</idno>
<idno type="doi">10.1371/journal.pone.0072866</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">000D72</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000D72</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Enhancement of Collective Immunity in Tokyo Metropolitan Area by Selective Vaccination against an Emerging Influenza Pandemic</title>
<author><name sortKey="Saito, Masaya M" sort="Saito, Masaya M" uniqKey="Saito M" first="Masaya M." last="Saito">Masaya M. Saito</name>
<affiliation><nlm:aff id="aff1"><addr-line>Institute of Statistical Mathematics, Tachikawa, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Imoto, Seiya" sort="Imoto, Seiya" uniqKey="Imoto S" first="Seiya" last="Imoto">Seiya Imoto</name>
<affiliation><nlm:aff id="aff2"><addr-line>Human Genome Center, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Yamaguchi, Rui" sort="Yamaguchi, Rui" uniqKey="Yamaguchi R" first="Rui" last="Yamaguchi">Rui Yamaguchi</name>
<affiliation><nlm:aff id="aff2"><addr-line>Human Genome Center, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Tsubokura, Masaharu" sort="Tsubokura, Masaharu" uniqKey="Tsubokura M" first="Masaharu" last="Tsubokura">Masaharu Tsubokura</name>
<affiliation><nlm:aff id="aff3"><addr-line>Division of Social Communication Systems for Advanced Clinical Research, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kami, Masahiro" sort="Kami, Masahiro" uniqKey="Kami M" first="Masahiro" last="Kami">Masahiro Kami</name>
<affiliation><nlm:aff id="aff3"><addr-line>Division of Social Communication Systems for Advanced Clinical Research, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Nakada, Haruka" sort="Nakada, Haruka" uniqKey="Nakada H" first="Haruka" last="Nakada">Haruka Nakada</name>
<affiliation><nlm:aff id="aff3"><addr-line>Division of Social Communication Systems for Advanced Clinical Research, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Sato, Hiroki" sort="Sato, Hiroki" uniqKey="Sato H" first="Hiroki" last="Sato">Hiroki Sato</name>
<affiliation><nlm:aff id="aff4"><addr-line>Department of Preventive Medicine and Public Health, National Defense Medical College, Tokorozawa, Saitama, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Miyano, Satoru" sort="Miyano, Satoru" uniqKey="Miyano S" first="Satoru" last="Miyano">Satoru Miyano</name>
<affiliation><nlm:aff id="aff2"><addr-line>Human Genome Center, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Higuchi, Tomoyuki" sort="Higuchi, Tomoyuki" uniqKey="Higuchi T" first="Tomoyuki" last="Higuchi">Tomoyuki Higuchi</name>
<affiliation><nlm:aff id="aff1"><addr-line>Institute of Statistical Mathematics, Tachikawa, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint><date when="2013">2013</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>Vaccination is a preventive measure against influenza that does not require placing restrictions on social activities. However, since the stockpile of vaccine that can be prepared before the arrival of an emerging pandemic strain is generally quite limited, one has to select priority target groups to which the first stockpile is distributed. In this paper, we study a simulation-based priority target selection method with the goal of enhancing the collective immunity of the whole population. To model the region in which the disease spreads, we consider an urban area composed of suburbs and central areas connected by a single commuter train line. Human activity is modelled following an agent-based approach. The degree to which collective immunity is enhanced is judged by the attack rate in unvaccinated people. The simulation results show that if students and office workers are given exclusive priority in the first three months, the attack rate can be reduced from <inline-formula><inline-graphic xlink:href="pone.0072866.e001.jpg"></inline-graphic>
</inline-formula>
in the baseline case down to 1–2%. In contrast, random vaccination only slightly reduces the attack rate. It should be noted that giving preference to active social groups does not mean sacrificing those at high risk, which corresponds to the elderly in our simulation model. Compared with the random administration of vaccine to all social groups, this design successfully reduces the attack rate across all age groups.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Kerkhove, Mdv" uniqKey="Kerkhove M">MDV Kerkhove</name>
</author>
<author><name sortKey="Vandemaele, Kah" uniqKey="Vandemaele K">KAH Vandemaele</name>
</author>
<author><name sortKey="Shinde, V" uniqKey="Shinde V">V Shinde</name>
</author>
<author><name sortKey="Jaramillo Gutierrez, G" uniqKey="Jaramillo Gutierrez G">G Jaramillo-Gutierrez</name>
</author>
<author><name sortKey="Koukounari, A" uniqKey="Koukounari A">A Koukounari</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sha, L" uniqKey="Sha L">L Sha</name>
</author>
<author><name sortKey="Zhu, R" uniqKey="Zhu R">R Zhu</name>
</author>
<author><name sortKey="Cao, L" uniqKey="Cao L">L Cao</name>
</author>
<author><name sortKey="Yuan, Y" uniqKey="Yuan Y">Y Yuan</name>
</author>
<author><name sortKey="Li, Y" uniqKey="Li Y">Y Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
<author><name sortKey="Cummings, Dat" uniqKey="Cummings D">DAT Cummings</name>
</author>
<author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name>
</author>
<author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name>
</author>
<author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
<author><name sortKey="Nizam, A" uniqKey="Nizam A">A Nizam</name>
</author>
<author><name sortKey="Xu, S" uniqKey="Xu S">S Xu</name>
</author>
<author><name sortKey="Ungchusak, K" uniqKey="Ungchusak K">K Ungchusak</name>
</author>
<author><name sortKey="Hanshaoworakul, W" uniqKey="Hanshaoworakul W">W Hanshaoworakul</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
<author><name sortKey="Cummings, Dat" uniqKey="Cummings D">DAT Cummings</name>
</author>
<author><name sortKey="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>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sato, H" uniqKey="Sato H">H Sato</name>
</author>
<author><name sortKey="Nakada, H" uniqKey="Nakada H">H Nakada</name>
</author>
<author><name sortKey="Yamaguchi, R" uniqKey="Yamaguchi R">R Yamaguchi</name>
</author>
<author><name sortKey="Imoto, S" uniqKey="Imoto S">S Imoto</name>
</author>
<author><name sortKey="Miyano, S" uniqKey="Miyano S">S Miyano</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Kermack, Wo" uniqKey="Kermack W">WO Kermack</name>
</author>
<author><name sortKey="Mckendrick, Ag" uniqKey="Mckendrick A">AG McKendrick</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dawood, F" uniqKey="Dawood F">F Dawood</name>
</author>
<author><name sortKey="Jain, S" uniqKey="Jain S">S Jain</name>
</author>
<author><name sortKey="Finelli, L" uniqKey="Finelli L">L Finelli</name>
</author>
<author><name sortKey="Shaw, M" uniqKey="Shaw M">M Shaw</name>
</author>
<author><name sortKey="Lindstrom, S" uniqKey="Lindstrom S">S Lindstrom</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="Ira, M" uniqKey="Ira M">M Ira</name>
</author>
<author><name sortKey="Longini, J" uniqKey="Longini J">J Longini</name>
</author>
<author><name sortKey="Macken, Ca" uniqKey="Macken C">CA Macken</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cooley, P" uniqKey="Cooley P">P Cooley</name>
</author>
<author><name sortKey="Brown, S" uniqKey="Brown S">S Brown</name>
</author>
<author><name sortKey="Cajka, J" uniqKey="Cajka J">J Cajka</name>
</author>
<author><name sortKey="Chasteen, B" uniqKey="Chasteen B">B Chasteen</name>
</author>
<author><name sortKey="Ganapathi, L" uniqKey="Ganapathi L">L Ganapathi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Feigin, R" uniqKey="Feigin R">R Feigin</name>
</author>
<author><name sortKey="Baker, C" uniqKey="Baker C">C Baker</name>
</author>
<author><name sortKey="Herwaldt, L" uniqKey="Herwaldt L">L Herwaldt</name>
</author>
<author><name sortKey="Lampe, R" uniqKey="Lampe R">R Lampe</name>
</author>
<author><name sortKey="Mason, E" uniqKey="Mason E">E Mason</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dick, E" uniqKey="Dick E">E Dick</name>
</author>
<author><name sortKey="Jennings, L" uniqKey="Jennings L">L Jennings</name>
</author>
<author><name sortKey="Mink, K" uniqKey="Mink K">K Mink</name>
</author>
<author><name sortKey="Wartgow, C" uniqKey="Wartgow C">C Wartgow</name>
</author>
<author><name sortKey="Inhorn, S" uniqKey="Inhorn S">S Inhorn</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Duguid, J" uniqKey="Duguid J">J Duguid</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Jefferson, T" uniqKey="Jefferson T">T Jefferson</name>
</author>
<author><name sortKey="Rivetti, A" uniqKey="Rivetti A">A Rivetti</name>
</author>
<author><name sortKey="Harnden, A" uniqKey="Harnden A">A Harnden</name>
</author>
<author><name sortKey="Pietrantonj, Cd" uniqKey="Pietrantonj C">CD Pietrantonj</name>
</author>
<author><name sortKey="Demicheli, V" uniqKey="Demicheli V">V Demicheli</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jefferson, T" uniqKey="Jefferson T">T Jefferson</name>
</author>
<author><name sortKey="Pietrantonj, Cd" uniqKey="Pietrantonj C">CD Pietrantonj</name>
</author>
<author><name sortKey="Rivetti, A" uniqKey="Rivetti A">A Rivetti</name>
</author>
<author><name sortKey="Bawazeer, Ga" uniqKey="Bawazeer G">GA Bawazeer</name>
</author>
<author><name sortKey="Al Ansary, La" uniqKey="Al Ansary L">LA Al-Ansary</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jefferson, T" uniqKey="Jefferson T">T Jefferson</name>
</author>
<author><name sortKey="Pietrantonj, Cd" uniqKey="Pietrantonj C">CD Pietrantonj</name>
</author>
<author><name sortKey="Al Ansary, L" uniqKey="Al Ansary L">L Al-Ansary</name>
</author>
<author><name sortKey="Ferroni, E" uniqKey="Ferroni E">E Ferroni</name>
</author>
<author><name sortKey="Thorning, S" uniqKey="Thorning S">S Thorning</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Geoffery, Aw" uniqKey="Geoffery A">AW Geoffery</name>
</author>
<author><name sortKey="Perter, Gs" uniqKey="Perter G">GS Perter</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Inaida, S" uniqKey="Inaida S">S Inaida</name>
</author>
<author><name sortKey="Yasui, Y" uniqKey="Yasui Y">Y Yasui</name>
</author>
<author><name sortKey="Tada, Y" uniqKey="Tada Y">Y Tada</name>
</author>
<author><name sortKey="Taniguchi, K" uniqKey="Taniguchi K">K Taniguchi</name>
</author>
<author><name sortKey="Okabe, N" uniqKey="Okabe N">N Okabe</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, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">24058445</article-id>
<article-id pub-id-type="pmc">3776821</article-id>
<article-id pub-id-type="publisher-id">PONE-D-12-17684</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0072866</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Enhancement of Collective Immunity in Tokyo Metropolitan Area by Selective Vaccination against an Emerging Influenza Pandemic</article-title>
<alt-title alt-title-type="running-head">Enhancement of Collective Immunity for Influenza</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Saito</surname>
<given-names>Masaya M.</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Imoto</surname>
<given-names>Seiya</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Yamaguchi</surname>
<given-names>Rui</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Tsubokura</surname>
<given-names>Masaharu</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Kami</surname>
<given-names>Masahiro</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Nakada</surname>
<given-names>Haruka</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Sato</surname>
<given-names>Hiroki</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Miyano</surname>
<given-names>Satoru</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Higuchi</surname>
<given-names>Tomoyuki</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>Institute of Statistical Mathematics, Tachikawa, Tokyo, Japan</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>Human Genome Center, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>Division of Social Communication Systems for Advanced Clinical Research, Institute of Medical Science, University of Tokyo, Minato-ku, Tokyo, Japan</addr-line>
</aff>
<aff id="aff4"><label>4</label>
<addr-line>Department of Preventive Medicine and Public Health, National Defense Medical College, Tokorozawa, Saitama, Japan</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Boni</surname>
<given-names>Maciej F.</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>University of Oxford, Viet Nam</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>saitohm@ism.ac.jp</email>
</corresp>
<fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p>Conceived and designed simulations: SI RY TH SM. Performed the simulations: MMS. Wrote the paper: MMS SI RY. Contributed the development of the simulator: MMS. Gave medical advice on the simulation setting: MT MK HN HS.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><year>2013</year>
</pub-date>
<pub-date pub-type="epub"><day>18</day>
<month>9</month>
<year>2013</year>
</pub-date>
<volume>8</volume>
<issue>9</issue>
<elocation-id>e72866</elocation-id>
<history><date date-type="received"><day>21</day>
<month>6</month>
<year>2012</year>
</date>
<date date-type="accepted"><day>20</day>
<month>7</month>
<year>2013</year>
</date>
</history>
<permissions><copyright-statement>© 2013 Saito et al</copyright-statement>
<copyright-year>2013</copyright-year>
<copyright-holder>Saito et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p>
</license>
</permissions>
<abstract><p>Vaccination is a preventive measure against influenza that does not require placing restrictions on social activities. However, since the stockpile of vaccine that can be prepared before the arrival of an emerging pandemic strain is generally quite limited, one has to select priority target groups to which the first stockpile is distributed. In this paper, we study a simulation-based priority target selection method with the goal of enhancing the collective immunity of the whole population. To model the region in which the disease spreads, we consider an urban area composed of suburbs and central areas connected by a single commuter train line. Human activity is modelled following an agent-based approach. The degree to which collective immunity is enhanced is judged by the attack rate in unvaccinated people. The simulation results show that if students and office workers are given exclusive priority in the first three months, the attack rate can be reduced from <inline-formula><inline-graphic xlink:href="pone.0072866.e001.jpg"></inline-graphic>
</inline-formula>
in the baseline case down to 1–2%. In contrast, random vaccination only slightly reduces the attack rate. It should be noted that giving preference to active social groups does not mean sacrificing those at high risk, which corresponds to the elderly in our simulation model. Compared with the random administration of vaccine to all social groups, this design successfully reduces the attack rate across all age groups.</p>
</abstract>
<funding-group><funding-statement>The authors have no support or funding to report.</funding-statement>
</funding-group>
<counts><page-count count="9"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>New influenza strains continuously appear via mutations and/or reassortment. Once a new strain acquires stable human-human transmissibility, a pandemic can easily arise since few people will have immunity against that strain. Recently, we experienced the 2009 pandemic. Fortunately, the 2009 flu was low-pathogenic strain, and its social impact was similar to that of seasonal influenza <xref rid="pone.0072866-Kerkhove1" ref-type="bibr">[1]</xref>
, although the symptoms were more severe in young children <xref rid="pone.0072866-Sha1" ref-type="bibr">[2]</xref>
. However, there is a highly pathogenic strain of avian flu A (subtype H5N1), and sporadic cases of human infection with this strain began to be reported in the late 1990s <xref rid="pone.0072866-World1" ref-type="bibr">[3]</xref>
. If this strain were to acquire stable human-human transmissibility, the resulting pandemic would have dire consequences.</p>
<p>Pioneering studies of the application of agent-based simulation to influenza prevention were performed for Southeast Asia <xref rid="pone.0072866-Ferguson1" ref-type="bibr">[4]</xref>
, <xref rid="pone.0072866-Longini1" ref-type="bibr">[5]</xref>
and the US and UK <xref rid="pone.0072866-Ferguson2" ref-type="bibr">[6]</xref>
. This research studied the feasibility of containment, that is, preventing the spread from one small population to another within a targeted region/country, as a function of the global reproduction number. We here mention similarity and difference between our model and the ones that were used in these studies. In Ferguson et al., a household for each individual, a single group to which each individual belonged, and the community to which the individual belonged were included in the calculation of the force of infection (FOI). The contribution from the community was modelled by a gravity model, and the corresponding FOI term decreased with the distance from infectious individuals. Instead of the effect of the community, Longini et al. allowed individuals to belong to multiple groups. In our simulation model, the concept of belonging to a group was implemented by the individuals' daily schedule, and the disease was transmitted inside these temporary groups. This approach requires the time step of the simulation to be small (1 minute in our simulation compared to 6 hours in Ferguson et al. and 24 hours in Longini et al.). We took this computationally expensive approach in order to incorporate the effect of individuals moving around the target city via trains.</p>
<p>The present study uses large-scale multi-agent simulations to evaluate the effectiveness of vaccinations in suppressing an influenza epidemic in an urban area. During epidemics inside a city, trains (or other public transport systems) play an important role in the increase in the number of infected people. Trains not only transport people to their workplaces or schools, but also gather them in high-density groups. One impractical approach for suppressing the spread of epidemics is to suspend train services and/or implement furloughs to prevent the concentration of people in places where transmission would be high. However, since these measures are obviously unacceptable due to decreased quality of life and economic impacts, other measures are desired that can be applied without halting the activity of the city. The administration of vaccines to enhance the collective immunity is one means of suppressing the spread of influenza. Since we cannot assume a long period of time from the discovery of the strain to its arrival in any particular region (e.g., the 2009 pandemic strain arrived in Japan two weeks after discovery) and since border quarantine inspections cannot fully detect infected people <xref rid="pone.0072866-Sato1" ref-type="bibr">[7]</xref>
, the administration of vaccine to the residents of the region will be simultaneous with the spread of the influenza. We were interested in determining if the degree of collective immunity could be enhanced if vaccines were distributed preferentially to some selected groups rather than equally to all individuals. To investigate this, we carried out simulations in which different groups received priority for vaccination.</p>
</sec>
<sec sec-type="methods" id="s2"><title>Method</title>
<sec id="s2a"><title>Modelling of human activity</title>
<p>We employed a multi-agent simulation approach to describe groups of people who are able to contact each other. Individuals move around different kinds of places in an urban area according to a role-dependent schedule. The disease (influenza) is stochastically transmitted in these local groups, which are formed as a consequence of the actions of individuals.</p>
<p>Our simulator contains several towns in an urban area that are connected by a commuter train line, along which influenza can spread. Each town consists of several types of places and residents. The places include workplaces, schools, homes, and public places, such as shops and parks. Residents are classified into employees, students, and domiciliaries, where domiciliaries referes to individuals who spend the majority of their day in their home (e.g. stay-at-home caretakers, some elderly persons). The members of the first two classes have a specific place where they visit every day and relatively many opportunities to contact with other people. The members of the third class (domiciliaries) spend most of the day in their households and have fewer contacts. The types of places and residents are related by the schedules followed by individuals. For example, if an individual has the role of “employee,” then their schedule requires that they go toward a particular “workplace” at a randomly selected time each morning (for the other elements in schedules, see Section S1.4 of <xref ref-type="supplementary-material" rid="pone.0072866.s001">File S1</xref>
). As a consequence of the evaluation of these schedules for the residents, temporary groups are formed, such as “workplace,” “schools,” and so on. It is assumed at any moment that the disease is transmitted only inside such groups.</p>
<p>In this study, the model was constructed to render a simplified Tokyo metropolitan area. The urban area has five towns connected by a commuter train line. A schematic illustration of the model city and the simulation is shown in <xref ref-type="fig" rid="pone-0072866-g001">Fig. 1</xref>
. For the numbers of residents and schools in towns A–E (<xref ref-type="table" rid="pone-0072866-t001">Table 1</xref>
), we used data for the year 2005 as published by the Tokyo metropolitan government <xref rid="pone.0072866-Tokyo1" ref-type="bibr">[8]</xref>
. However, the numbers of workplaces are arbitrarily given so as to be concentrated at the central areas (towns D and E) and their scale (the number of employees) is assumed to follow a Pareto law (Fig. S3 in <xref ref-type="supplementary-material" rid="pone.0072866.s001">File S1</xref>
). The bottom three panels of <xref ref-type="fig" rid="pone-0072866-g001">Fig. 1</xref>
show how epidemic spread over workplaces. It takes a fortnight or longer for most of workplaces to encounter infectious persons (middle panel) and an epidemic peaks at around the 75-th day (right panel; here a Pareto low distribution of workplaces can be seen). We shall do two different scales of simulations. In full-scale simulations, the numbers of places and residents are the same as <xref ref-type="table" rid="pone-0072866-t001">Table 1</xref>
(the number of agents is 1,246,732; see Fig. S2 in <xref ref-type="supplementary-material" rid="pone.0072866.s001">File S1</xref>
for age-specific population of each role). In the scaled-down simulations, these numbers are reduced by a factor of one-tenth, namely, the number of agents is 124,673.</p>
<fig id="pone-0072866-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.g001</object-id>
<label>Figure 1</label>
<caption><title>Schematic illustration of the simulator and the model city.</title>
<p>Individual activities are explicitly described and infectious transmissions are stochastically evaluated at local places according to the number density of infectious visitors. The bottom three panels show how influenza spreads over many workplaces in a typical baseline simulation run, where the other places are blank due to a small number of patients there.</p>
</caption>
<graphic xlink:href="pone.0072866.g001"></graphic>
</fig>
<table-wrap id="pone-0072866-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.t001</object-id>
<label>Table 1</label>
<caption><title>Places and Population in the model city.</title>
</caption>
<alternatives><graphic id="pone-0072866-t001-1" xlink:href="pone.0072866.t001"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">town</td>
<td align="left" rowspan="1" colspan="1">school</td>
<td align="left" rowspan="1" colspan="1">workplace</td>
<td align="left" rowspan="1" colspan="1">park</td>
<td align="left" rowspan="1" colspan="1">population</td>
<td align="left" rowspan="1" colspan="1">shop</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">A</td>
<td align="left" rowspan="1" colspan="1">70</td>
<td align="left" rowspan="1" colspan="1">100</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">571,641</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">B</td>
<td align="left" rowspan="1" colspan="1">20</td>
<td align="left" rowspan="1" colspan="1">100</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">176,866</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">C</td>
<td align="left" rowspan="1" colspan="1">12</td>
<td align="left" rowspan="1" colspan="1">100</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">138,684</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">D</td>
<td align="left" rowspan="1" colspan="1">29</td>
<td align="left" rowspan="1" colspan="1">2000</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">314,861</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">E</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">2000</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">44,680</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</sec>
<sec id="s2b"><title>Modelling of Transmission and Vaccine</title>
<p>In our simulation, the health state and the progression of the disease are described in a manner similar to the traditional Susceptible-Exposed-Infectious-Recovered (SEIR) model <xref rid="pone.0072866-Kermack1" ref-type="bibr">[9]</xref>
except that here the disease is stochastically transmitted in local groups. The health status of any individual is described by one of four states: <italic>s</italic>
(susceptible), <italic>e</italic>
(exposed; infected but not yet infectious), <italic>i</italic>
(infected; infected and infectious), and <italic>r</italic>
(removed; either recovered or dead).</p>
<p>For the stochastic transition of an individual in state <inline-formula><inline-graphic xlink:href="pone.0072866.e002.jpg"></inline-graphic>
</inline-formula>
to state <inline-formula><inline-graphic xlink:href="pone.0072866.e003.jpg"></inline-graphic>
</inline-formula>
, the probability per unit time of the transition is proportional to the density of individuals in state <inline-formula><inline-graphic xlink:href="pone.0072866.e004.jpg"></inline-graphic>
</inline-formula>
in the same place and the transmission efficiency assigned to that place. The transmission efficiencies of the various types of places are shown in terms of the basic reproduction number in <xref ref-type="table" rid="pone-0072866-t002">Table 2</xref>
. Further transitions of each individual's health status to state <inline-formula><inline-graphic xlink:href="pone.0072866.e005.jpg"></inline-graphic>
</inline-formula>
and state <inline-formula><inline-graphic xlink:href="pone.0072866.e006.jpg"></inline-graphic>
</inline-formula>
are made with constant probabilities. The inverses of these probabilities are the mean latent and infectious periods and are assigned typical values <xref rid="pone.0072866-Dawood1" ref-type="bibr">[10]</xref>
of 3.5 days and 3 days, respectively.</p>
<table-wrap id="pone-0072866-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.t002</object-id>
<label>Table 2</label>
<caption><title>Transmission efficiencies of each place type in terms of basic reproduction number.</title>
</caption>
<alternatives><graphic id="pone-0072866-t002-2" xlink:href="pone.0072866.t002"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">place</td>
<td align="left" rowspan="1" colspan="1">train</td>
<td align="left" rowspan="1" colspan="1">school</td>
<td align="left" rowspan="1" colspan="1">workplace</td>
<td align="left" rowspan="1" colspan="1">household</td>
<td align="left" rowspan="1" colspan="1">park</td>
<td align="left" rowspan="1" colspan="1">shop</td>
<td align="left" rowspan="1" colspan="1">AR</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">case A</td>
<td align="left" rowspan="1" colspan="1">1.5</td>
<td align="left" rowspan="1" colspan="1">2.0</td>
<td align="left" rowspan="1" colspan="1">2.0</td>
<td align="left" rowspan="1" colspan="1">2.0</td>
<td align="left" rowspan="1" colspan="1">1.0</td>
<td align="left" rowspan="1" colspan="1">0.6</td>
<td align="left" rowspan="1" colspan="1">0.30</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">case B</td>
<td align="left" rowspan="1" colspan="1">1.0</td>
<td align="left" rowspan="1" colspan="1">2.5</td>
<td align="left" rowspan="1" colspan="1">2.5</td>
<td align="left" rowspan="1" colspan="1">2.0</td>
<td align="left" rowspan="1" colspan="1">1.0</td>
<td align="left" rowspan="1" colspan="1">0.6</td>
<td align="left" rowspan="1" colspan="1">0.29</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">case C</td>
<td align="left" rowspan="1" colspan="1">1.0</td>
<td align="left" rowspan="1" colspan="1">2.0</td>
<td align="left" rowspan="1" colspan="1">2.0</td>
<td align="left" rowspan="1" colspan="1">4.0</td>
<td align="left" rowspan="1" colspan="1">1.0</td>
<td align="left" rowspan="1" colspan="1">0.6</td>
<td align="left" rowspan="1" colspan="1">0.28</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>The effect of the vaccine is controlled by two parameters: <inline-formula><inline-graphic xlink:href="pone.0072866.e007.jpg"></inline-graphic>
</inline-formula>
, which indicates whether a vaccine enhances the immunity of a dosed person, and <inline-formula><inline-graphic xlink:href="pone.0072866.e008.jpg"></inline-graphic>
</inline-formula>
, the degree to which their immunity is enhanced. If a vaccine is given, the immunity of the target person is enhanced with probability <inline-formula><inline-graphic xlink:href="pone.0072866.e009.jpg"></inline-graphic>
</inline-formula>
, and the degree of immunity is enhanced and the symptoms relaxed by a factor of <inline-formula><inline-graphic xlink:href="pone.0072866.e010.jpg"></inline-graphic>
</inline-formula>
. Suppose that <inline-formula><inline-graphic xlink:href="pone.0072866.e011.jpg"></inline-graphic>
</inline-formula>
is the probability that an unvaccinated susceptible person transitions into an exposed state. If the person is vaccinated and the vaccine is effective, then the probability is reduced to <inline-formula><inline-graphic xlink:href="pone.0072866.e012.jpg"></inline-graphic>
</inline-formula>
and the average latent period <inline-formula><inline-graphic xlink:href="pone.0072866.e013.jpg"></inline-graphic>
</inline-formula>
is shortened to <inline-formula><inline-graphic xlink:href="pone.0072866.e014.jpg"></inline-graphic>
</inline-formula>
. The parameter <inline-formula><inline-graphic xlink:href="pone.0072866.e015.jpg"></inline-graphic>
</inline-formula>
takes a value between 0 (perfect protection) to 1 (useless vaccine). For further details of the implementation, see Section S1.2 of <xref ref-type="supplementary-material" rid="pone.0072866.s001">File S1</xref>
(Fig. S1 in <xref ref-type="supplementary-material" rid="pone.0072866.s001">File S1</xref>
).</p>
<p>In general, it is not known in which place, households or workplaces/schools, there is a higher likelihood of transmission. We thus considered three representative cases while keeping the resultant global attack rate (AR) around 30%: (A) household dominant, (B) workplace/school dominant, and (C) all three places being equal. In order to inform relative importance, we show the number of transmission in each kind of places in <xref ref-type="table" rid="pone-0072866-t003">Table 3</xref>
. Previous studies <xref rid="pone.0072866-Ferguson1" ref-type="bibr">[4]</xref>
, <xref rid="pone.0072866-Longini1" ref-type="bibr">[5]</xref>
, <xref rid="pone.0072866-Germann1" ref-type="bibr">[11]</xref>
, <xref rid="pone.0072866-Cooley1" ref-type="bibr">[12]</xref>
devoted to agent-based influenza simulations assigned 33% to the AR; this was derived from a retrospective analysis of the 1957–1958 pandemic. It would be useful to compare the value we used for the transmissibility with the values used in previous works. Because we used a different transmission model, however, we need to first examine the corresponding model parameters. Let us consider a situation in which a person remains stationary and receives <inline-formula><inline-graphic xlink:href="pone.0072866.e016.jpg"></inline-graphic>
</inline-formula>
visitors, <inline-formula><inline-graphic xlink:href="pone.0072866.e017.jpg"></inline-graphic>
</inline-formula>
of whom are infectious. In these studies, the transmission process is governed by two parameters, the transmission probability <inline-formula><inline-graphic xlink:href="pone.0072866.e018.jpg"></inline-graphic>
</inline-formula>
per contact and the number <inline-formula><inline-graphic xlink:href="pone.0072866.e019.jpg"></inline-graphic>
</inline-formula>
of persons with whom a person has contact on a given day. The probability <inline-formula><inline-graphic xlink:href="pone.0072866.e020.jpg"></inline-graphic>
</inline-formula>
with which the person transfers to the exposed state is given by<disp-formula id="pone.0072866.e021"><graphic xlink:href="pone.0072866.e021.jpg" position="anchor" orientation="portrait"></graphic>
</disp-formula>
In our model, the counterpart of <inline-formula><inline-graphic xlink:href="pone.0072866.e022.jpg"></inline-graphic>
</inline-formula>
is given by<disp-formula id="pone.0072866.e023"><graphic xlink:href="pone.0072866.e023.jpg" position="anchor" orientation="portrait"></graphic>
</disp-formula>
where <inline-formula><inline-graphic xlink:href="pone.0072866.e024.jpg"></inline-graphic>
</inline-formula>
is the site-specific transmission efficiency in terms of the reproduction number and <inline-formula><inline-graphic xlink:href="pone.0072866.e025.jpg"></inline-graphic>
</inline-formula>
is the time spent visiting that location. Comparing the two forms for <inline-formula><inline-graphic xlink:href="pone.0072866.e026.jpg"></inline-graphic>
</inline-formula>
, we obtain the formula <inline-formula><inline-graphic xlink:href="pone.0072866.e027.jpg"></inline-graphic>
</inline-formula>
to relate the parameters of previous studies with ours. <xref ref-type="table" rid="pone-0072866-t004">Table 4</xref>
compares the transmission efficiencies of case A in <xref ref-type="table" rid="pone-0072866-t002">Table 2</xref>
with those of Cooley et al. <xref rid="pone.0072866-Cooley1" ref-type="bibr">[12]</xref>
, which studied the importance of transmission in subways in New York City. Our values for transmission in workplaces and households are comparable to the estimations of Cooley et al.; their values for commuter and noncommuter passengers differ greatly, and it is difficult to say anything beyond that our value is between these two.</p>
<table-wrap id="pone-0072866-t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.t003</object-id>
<label>Table 3</label>
<caption><title>Number of transmissions in each place type in baseline simulations.</title>
</caption>
<alternatives><graphic id="pone-0072866-t003-3" xlink:href="pone.0072866.t003"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">place</td>
<td align="left" rowspan="1" colspan="1">train</td>
<td align="left" rowspan="1" colspan="1">school</td>
<td align="left" rowspan="1" colspan="1">workplace</td>
<td align="left" rowspan="1" colspan="1">household</td>
<td align="left" rowspan="1" colspan="1">park</td>
<td align="left" rowspan="1" colspan="1">shop</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">case A</td>
<td align="left" rowspan="1" colspan="1">13498</td>
<td align="left" rowspan="1" colspan="1">2558</td>
<td align="left" rowspan="1" colspan="1">12623</td>
<td align="left" rowspan="1" colspan="1">10537</td>
<td align="left" rowspan="1" colspan="1">85</td>
<td align="left" rowspan="1" colspan="1">413</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">case B</td>
<td align="left" rowspan="1" colspan="1">8736</td>
<td align="left" rowspan="1" colspan="1">2816</td>
<td align="left" rowspan="1" colspan="1">15330</td>
<td align="left" rowspan="1" colspan="1">10052</td>
<td align="left" rowspan="1" colspan="1">51</td>
<td align="left" rowspan="1" colspan="1">395</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">case C</td>
<td align="left" rowspan="1" colspan="1">7777</td>
<td align="left" rowspan="1" colspan="1">2145</td>
<td align="left" rowspan="1" colspan="1">10597</td>
<td align="left" rowspan="1" colspan="1">12663</td>
<td align="left" rowspan="1" colspan="1">93</td>
<td align="left" rowspan="1" colspan="1">363</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<table-wrap id="pone-0072866-t004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.t004</object-id>
<label>Table 4</label>
<caption><title>Comparison of our transmission efficiencies with those of Cooley et al. <xref rid="pone.0072866-Cooley1" ref-type="bibr">[12]</xref>
.</title>
</caption>
<alternatives><graphic id="pone-0072866-t004-4" xlink:href="pone.0072866.t004"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">place</td>
<td align="left" rowspan="1" colspan="1">train 1</td>
<td align="left" rowspan="1" colspan="1">train 2</td>
<td align="left" rowspan="1" colspan="1">workplace</td>
<td align="left" rowspan="1" colspan="1">household</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1"><italic>P</italic>
(this paper)</td>
<td align="left" rowspan="1" colspan="1">0.1667</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0.3333</td>
<td align="left" rowspan="1" colspan="1">0.3333</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><inline-formula><inline-graphic xlink:href="pone.0072866.e028.jpg"></inline-graphic>
</inline-formula>
</td>
<td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">2</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><italic>τ</italic>
</td>
<td align="left" rowspan="1" colspan="1">2 hr</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">12 hr</td>
<td align="left" rowspan="1" colspan="1">12 hr</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><italic>P</italic>
(Cooley et al.)</td>
<td align="left" rowspan="1" colspan="1">1.9481</td>
<td align="left" rowspan="1" colspan="1">0.0324</td>
<td align="left" rowspan="1" colspan="1">0.2116</td>
<td align="left" rowspan="1" colspan="1">0.3688</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><italic>x</italic>
</td>
<td align="left" rowspan="1" colspan="1">0.0575</td>
<td align="left" rowspan="1" colspan="1">0.0048</td>
<td align="left" rowspan="1" colspan="1">0.0575</td>
<td align="left" rowspan="1" colspan="1">0.4</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><italic>c</italic>
</td>
<td align="left" rowspan="1" colspan="1">33.88</td>
<td align="left" rowspan="1" colspan="1">6.75</td>
<td align="left" rowspan="1" colspan="1">3.68</td>
<td align="left" rowspan="1" colspan="1">0.922</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt101"><label></label>
<p>Cooley et al. assigned different values of <inline-formula><inline-graphic xlink:href="pone.0072866.e029.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pone.0072866.e030.jpg"></inline-graphic>
</inline-formula>
to commuter passengers (train 1) and noncommuter ones (train 2). The value of <inline-formula><inline-graphic xlink:href="pone.0072866.e031.jpg"></inline-graphic>
</inline-formula>
assigned to adult-adult contact in households is used here, while different values were assigned to adult-children and child-child contacts in the original work.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>In a crowded train, one passenger may approach very close to another, and there is then a higher probability than usual that the disease will be transmitted. In order to incorporate this effect, we allowed the transmission efficiency on trains to increase with the mean person-person distance, which is calculated from the number of passengers. Based on studies of epidemics <xref rid="pone.0072866-Feigin1" ref-type="bibr">[13]</xref>
and aerosol dynamics <xref rid="pone.0072866-Dick1" ref-type="bibr">[14]</xref>
,<xref rid="pone.0072866-Duguid1" ref-type="bibr">[15]</xref>
, droplet transmission is very effective at distances of up to 1 meter, but may still be possible up to a distance of 3 meters <xref rid="pone.0072866-Siegel1" ref-type="bibr">[16]</xref>
. Based on this, we employed a transmission efficiency function that has a cutoff at 3 meters. The reproduction number in trains thus takes the value of 1.0 at this cutoff and increases up to 1.93, as shown in <xref ref-type="table" rid="pone-0072866-t002">Table 2</xref>
(see <xref ref-type="supplementary-material" rid="pone.0072866.s001">File S1</xref>
, Section S1.3 for further details).</p>
</sec>
<sec id="s2c"><title>Baseline epidemic simulation</title>
<p>As a baseline for the later evaluation of vaccination programs, we calculated the pattern in which influenza spread in the model urban area when there was no intervention (See Table S1 in <xref ref-type="supplementary-material" rid="pone.0072866.s001">File S1</xref>
for a complete list of parameters to configure simulation). This area was an abstraction of the towns along a commuter rail in the Tokyo metropolitan area. We assumed the initial number of infected people was 600 to ensure an outbreak would arise, that is, we did not consider cases in which the introduced strains became extinct. <xref ref-type="fig" rid="pone-0072866-g002">Figure 2</xref>
shows the spread pattern for a typical simulation run of case A. The epidemic peak occurs on approximately the 75-th day (left panel). We modelled that office workers and students remained at home on nights and weekends. This effect was reflected in the small daily but large weekly oscillations in the number of patients. The proportion of residents who are eventually infected is called the illness <italic>attack rate</italic>
(AR). The cumulative number of infected people, shown in the middle panel, indicates that the AR is <inline-formula><inline-graphic xlink:href="pone.0072866.e032.jpg"></inline-graphic>
</inline-formula>
in the baseline simulation.</p>
<fig id="pone-0072866-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.g002</object-id>
<label>Figure 2</label>
<caption><title>Baseline (without vaccination) simulation run.</title>
<p>Variation of the numbers of residents in states <inline-formula><inline-graphic xlink:href="pone.0072866.e033.jpg"></inline-graphic>
</inline-formula>
is shown for cases A, B and C of transmissiblity.</p>
</caption>
<graphic xlink:href="pone.0072866.g002"></graphic>
</fig>
</sec>
<sec id="s2d"><title>Calibration of vaccine strength</title>
<p>Prepared vaccines vary greatly in effectiveness depending on the target virus. Here, we calibrate the effect of the vaccine in simulations to reproduce an efficacy value estimated by retrospective studies for real seasonal influenza <xref rid="pone.0072866-Jefferson1" ref-type="bibr">[17]</xref>
, <xref rid="pone.0072866-Jefferson2" ref-type="bibr">[18]</xref>
, <xref rid="pone.0072866-Jefferson3" ref-type="bibr">[19]</xref>
. The efficacy is defined in terms of the ARs of the vaccinated and unvaccinated groups <xref rid="pone.0072866-Geoffery1" ref-type="bibr">[20]</xref>
,<disp-formula id="pone.0072866.e034"><graphic xlink:href="pone.0072866.e034.jpg" position="anchor" orientation="portrait"></graphic>
</disp-formula>
In real investigations, the ARV and ARU are calculated for different populations that satisfy the definitions, whereas in the simulations these can be calculated from two simulation runs, with vaccination and without vaccination, respectively, for a chosen small population that does not affect the collective immunity of the entire city (we set this to be 0.4% of the city residents).</p>
<p><xref ref-type="table" rid="pone-0072866-t005">Table 5</xref>
shows the efficacy calculated in this manner for different settings of vaccination strength. In this table, <inline-formula><inline-graphic xlink:href="pone.0072866.e035.jpg"></inline-graphic>
</inline-formula>
denotes the value of <inline-formula><inline-graphic xlink:href="pone.0072866.e036.jpg"></inline-graphic>
</inline-formula>
for people aged <inline-formula><inline-graphic xlink:href="pone.0072866.e037.jpg"></inline-graphic>
</inline-formula>
to <inline-formula><inline-graphic xlink:href="pone.0072866.e038.jpg"></inline-graphic>
</inline-formula>
years. In this experiment, the age-based values <inline-formula><inline-graphic xlink:href="pone.0072866.e039.jpg"></inline-graphic>
</inline-formula>
of the probability <inline-formula><inline-graphic xlink:href="pone.0072866.e040.jpg"></inline-graphic>
</inline-formula>
were chosen to be <inline-formula><inline-graphic xlink:href="pone.0072866.e041.jpg"></inline-graphic>
</inline-formula>
, <inline-formula><inline-graphic xlink:href="pone.0072866.e042.jpg"></inline-graphic>
</inline-formula>
, and <inline-formula><inline-graphic xlink:href="pone.0072866.e043.jpg"></inline-graphic>
</inline-formula>
. Based on studies <xref rid="pone.0072866-Jefferson1" ref-type="bibr">[17]</xref>
, <xref rid="pone.0072866-Jefferson2" ref-type="bibr">[18]</xref>
, <xref rid="pone.0072866-Jefferson3" ref-type="bibr">[19]</xref>
that state that the efficacy is approximately 60% in adults, we set <inline-formula><inline-graphic xlink:href="pone.0072866.e044.jpg"></inline-graphic>
</inline-formula>
, <inline-formula><inline-graphic xlink:href="pone.0072866.e045.jpg"></inline-graphic>
</inline-formula>
, and <inline-formula><inline-graphic xlink:href="pone.0072866.e046.jpg"></inline-graphic>
</inline-formula>
.</p>
<table-wrap id="pone-0072866-t005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.t005</object-id>
<label>Table 5</label>
<caption><title>Values of vaccination strength and resultant efficacies.</title>
</caption>
<alternatives><graphic id="pone-0072866-t005-5" xlink:href="pone.0072866.t005"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1"><italic>ν</italic>
<sub>1–5</sub>
</td>
<td align="left" rowspan="1" colspan="1"><italic>ν</italic>
<sub>5–65</sub>
</td>
<td align="left" rowspan="1" colspan="1"><italic>ν</italic>
<sub>65–</sub>
</td>
<td align="left" rowspan="1" colspan="1">efficacy</td>
<td align="left" rowspan="1" colspan="1">illUnc</td>
<td align="left" rowspan="1" colspan="1">illVac</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">0.2</td>
<td align="left" rowspan="1" colspan="1">0.1</td>
<td align="left" rowspan="1" colspan="1">0.2</td>
<td align="left" rowspan="1" colspan="1">0.7264</td>
<td align="left" rowspan="1" colspan="1">0.5338</td>
<td align="left" rowspan="1" colspan="1">0.1460</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">0.4</td>
<td align="left" rowspan="1" colspan="1">0.2</td>
<td align="left" rowspan="1" colspan="1">0.4</td>
<td align="left" rowspan="1" colspan="1">0.6073</td>
<td align="left" rowspan="1" colspan="1">0.5338</td>
<td align="left" rowspan="1" colspan="1">0.2096</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">0.8</td>
<td align="left" rowspan="1" colspan="1">0.4</td>
<td align="left" rowspan="1" colspan="1">0.8</td>
<td align="left" rowspan="1" colspan="1">0.4533</td>
<td align="left" rowspan="1" colspan="1">0.5338</td>
<td align="left" rowspan="1" colspan="1">0.2918</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">0.9</td>
<td align="left" rowspan="1" colspan="1">0.5</td>
<td align="left" rowspan="1" colspan="1">0.9</td>
<td align="left" rowspan="1" colspan="1">0.3844</td>
<td align="left" rowspan="1" colspan="1">0.5338</td>
<td align="left" rowspan="1" colspan="1">0.3286</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt102"><label></label>
<p>Attack rates in the unvaccinated and the vaccinated populations are denoted by illUnc and illVac, respectively.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
</sec>
<sec id="s3"><title>Results</title>
<p>We studied the performance of vaccination programs for the enhancement of collective immunity. Each program was defined by the target group and the coverage of the population of the urban area. The performance of each program was evaluated by the reduction in the AR of the unvaccinated population. For target groups, we considered arbitrarily chosen residents (Arb), employees (Emp), students (Stu), the combined group of students and employees (StuEmp), domiciliaries (Dom), train passengers (Tr), and the subset of StuEmp who use trains to get to their offices or schools (TrStuEmp). For each target, we carried out simulation runs in which 16%, 32%, and 48% of the chosen group were vaccinated. We assumed that vaccines could be distributed to 48% of the residents in three months; we also considered half and twice this reference rate. In all three cases, we assumed for simplicity that the vaccine stockpile was uniformly distributed. This reference speed was chosen according to the 2009 vaccination schedule in Japan <xref rid="pone.0072866-Ministry1" ref-type="bibr">[21]</xref>
. In this schedule, domestic production covers 45% of the Japanese population in six months, from October to March. The priority groups included pregnant women, individuals with diseases, children, and elderly persons. The rest of the population was covered by the remainder of the domestic vaccine, and supplemented by imported vaccine. We note here the discrepancy between the plan mentioned above and the estimation <xref rid="pone.0072866-Ministry2" ref-type="bibr">[22]</xref>
of the actual number of vaccines. Looking at <xref ref-type="table" rid="pone-0072866-t006">Table 6</xref>
, we see that the planned rate and the actual rate are similar until the middle of December, after which the actual rate is less than the planned rate. According to reports of a learned discussion <xref rid="pone.0072866-Ministry3" ref-type="bibr">[23]</xref>
, 2009 flu may influence people hesitate to be vaccinated and so there was an oversupply after January.</p>
<table-wrap id="pone-0072866-t006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.t006</object-id>
<label>Table 6</label>
<caption><title>Schedule of domestic vaccine supply for Japan in the 2009/2010 season and estimation of actual consumed courses.</title>
</caption>
<alternatives><graphic id="pone-0072866-t006-6" xlink:href="pone.0072866.t006"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Month</td>
<td colspan="2" align="left" rowspan="1">2009–10</td>
<td colspan="2" align="left" rowspan="1">11</td>
<td colspan="3" align="left" rowspan="1">12</td>
<td colspan="2" align="left" rowspan="1">2010–01</td>
<td colspan="2" align="left" rowspan="1">2</td>
<td colspan="2" align="left" rowspan="1">3</td>
<td align="left" rowspan="1" colspan="1">Total</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Period</td>
<td align="left" rowspan="1" colspan="1">I</td>
<td align="left" rowspan="1" colspan="1">II</td>
<td align="left" rowspan="1" colspan="1">I</td>
<td align="left" rowspan="1" colspan="1">II</td>
<td align="left" rowspan="1" colspan="1">I</td>
<td align="left" rowspan="1" colspan="1">II</td>
<td align="left" rowspan="1" colspan="1">III</td>
<td align="left" rowspan="1" colspan="1">I</td>
<td align="left" rowspan="1" colspan="1">II</td>
<td align="left" rowspan="1" colspan="1">I</td>
<td align="left" rowspan="1" colspan="1">II</td>
<td align="left" rowspan="1" colspan="1">I</td>
<td align="left" rowspan="1" colspan="1">II</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">Plan(×10<sup>4</sup>
)</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">118</td>
<td align="left" rowspan="1" colspan="1">134</td>
<td align="left" rowspan="1" colspan="1">364</td>
<td align="left" rowspan="1" colspan="1">452</td>
<td align="left" rowspan="1" colspan="1">572</td>
<td align="left" rowspan="1" colspan="1">515</td>
<td align="left" rowspan="1" colspan="1">459</td>
<td align="left" rowspan="1" colspan="1">659</td>
<td align="left" rowspan="1" colspan="1">649</td>
<td align="left" rowspan="1" colspan="1">557</td>
<td align="left" rowspan="1" colspan="1">349</td>
<td align="left" rowspan="1" colspan="1">560</td>
<td align="left" rowspan="1" colspan="1">5,388</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Actual(×10<sup>4</sup>
)</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">158</td>
<td align="left" rowspan="1" colspan="1">303</td>
<td align="left" rowspan="1" colspan="1">142</td>
<td align="left" rowspan="1" colspan="1">418</td>
<td align="left" rowspan="1" colspan="1">627</td>
<td align="left" rowspan="1" colspan="1">227</td>
<td align="left" rowspan="1" colspan="1">162</td>
<td align="left" rowspan="1" colspan="1">170</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">2,207</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p><xref ref-type="fig" rid="pone-0072866-g003">Figure 3</xref>
shows the AR in the unvaccinated population as a function of the number of vaccinated people. In the top panels, the denominator of the AR is the unvaccinated population (hereinafter this rate is called the unvaccinated AR), and the denominator in the bottom panels is the population of the entire city (hereinafter this rate is called the total AR). For each different vaccination program, we carried out eight runs, each with a different random initial distribution of infected individuals. Administration of the vaccine to employees was the most effective way to reduce the AR, followed by vaccinating students and domiciliaries. This is reasonable since office works are the most active residents. Targeting office workers was advantageous for all the different cases of transmission efficiencies, including case C, which assumed that households have higher transmission rates than do workplaces. If office workers, corresponding to 32% of the entire population, are vaccinated, the total AR is reduced to 2%. Observing the number of infected individuals over time, we find that transmissions are suppressed and become sporadic. When other groups are targeted instead, there is some reduction in transmission but not to this extent.</p>
<fig id="pone-0072866-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.g003</object-id>
<label>Figure 3</label>
<caption><title>Dependence of the illness attack rate (AR) on the target selection.</title>
<p>The abscissa denotes the proportion of vaccinated people and the ordinate denotes the AR in the unvaccinated population (left panel) or of all residents (right panel). For each target, the coverage is chosen to be 16%, 32%, or 48% in terms of whole the residents, and the speed of administration of vaccine is chosen so that 16% of the residents are distributed in 30 days.</p>
</caption>
<graphic xlink:href="pone.0072866.g003"></graphic>
</fig>
<p>Vaccination of domiciliaries produced a very limited effect. The total AR was almost the same as that of the baseline. Note that the unvaccinated AR was larger than the baseline AR, and it even increased with increased vaccine coverage due to changes in the denominator. As the number of vaccinated domiciliaries increases, the proportion of office workers in the unvaccinated group increases. This increases the unvaccinated AR.</p>
<p>Targeting students had a slightly better effect, but it did not bring benefits to unvaccinated individuals. When students making up 16% of the city residents were vaccinated, the total AR was reduced to 20–25%, depending on the efficiency of transmission. The transmission rate in the different types of location provides information about the extent of benefit to unvaccinated individuals, as shown in <xref ref-type="fig" rid="pone-0072866-g004">Fig. 4</xref>
. Vaccinating students somewhat reduced transmissions in schools, comparable to the effect in offices when office workers were targeted. The number of transmissions in other locations, however, remained near the baseline level.</p>
<fig id="pone-0072866-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.g004</object-id>
<label>Figure 4</label>
<caption><title>Number of transmissions in each place type in simulations with vaccination.</title>
</caption>
<graphic xlink:href="pone.0072866.g004"></graphic>
</fig>
<p>Vaccinating arbitrarily chosen people provided an intermediate effect. Unlike with students, the AR can be further reduced by extending the coverage. Targeting random train passengers for vaccination was slightly better than a fully random program. This is because the proportion of office workers is larger in the train passenger group than in city as a whole.</p>
<p>Prioritizing active members, as suggested by the simulation results, is contrary to conventional approaches, which prioritize high-risk individuals (e.g., the elderly and people with disease) followed by medical practitioners. One may wonder if prioritizing office workers would disadvantage high-risk individuals, but as we show in <xref ref-type="table" rid="pone-0072866-t007">Table 7</xref>
, prioritizing active members results in a smaller attack rate of the elderly than does random vaccination.</p>
<table-wrap id="pone-0072866-t007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.t007</object-id>
<label>Table 7</label>
<caption><title>Age-based AR for each vaccination target in case A.</title>
</caption>
<alternatives><graphic id="pone-0072866-t007-7" xlink:href="pone.0072866.t007"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Target</td>
<td align="left" rowspan="1" colspan="1">Age 0–5</td>
<td align="left" rowspan="1" colspan="1">Age 5–65</td>
<td align="left" rowspan="1" colspan="1">Age 65–</td>
</tr>
</thead>
<tbody><tr><td colspan="4" align="left" rowspan="1">coverage 16%</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Arbitrary</td>
<td align="left" rowspan="1" colspan="1">0.205 (0.014)</td>
<td align="left" rowspan="1" colspan="1">0.240 (0.018)</td>
<td align="left" rowspan="1" colspan="1">0.076 (0.006)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Domiciliary</td>
<td align="left" rowspan="1" colspan="1">0.301 (0.007)</td>
<td align="left" rowspan="1" colspan="1">0.347 (0.006)</td>
<td align="left" rowspan="1" colspan="1">0.106 (0.005)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Student</td>
<td align="left" rowspan="1" colspan="1">0.142 (0.008)</td>
<td align="left" rowspan="1" colspan="1">0.276 (0.007)</td>
<td align="left" rowspan="1" colspan="1">0.083 (0.004)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Employee</td>
<td align="left" rowspan="1" colspan="1">0.100 (0.037)</td>
<td align="left" rowspan="1" colspan="1">0.086 (0.031)</td>
<td align="left" rowspan="1" colspan="1">0.029 (0.011)</td>
</tr>
<tr><td colspan="4" align="left" rowspan="1">coverage 32%</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Arbitrary</td>
<td align="left" rowspan="1" colspan="1">0.103 (0.037)</td>
<td align="left" rowspan="1" colspan="1">0.118 (0.043)</td>
<td align="left" rowspan="1" colspan="1">0.036 (0.015)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Domiciliary</td>
<td align="left" rowspan="1" colspan="1">0.295 (0.014)</td>
<td align="left" rowspan="1" colspan="1">0.340 (0.006)</td>
<td align="left" rowspan="1" colspan="1">0.094 (0.003)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Student</td>
<td align="left" rowspan="1" colspan="1">0.124 (0.007)</td>
<td align="left" rowspan="1" colspan="1">0.268 (0.008)</td>
<td align="left" rowspan="1" colspan="1">0.080 (0.003)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Employee</td>
<td align="left" rowspan="1" colspan="1">0.029 (0.015)</td>
<td align="left" rowspan="1" colspan="1">0.020 (0.009)</td>
<td align="left" rowspan="1" colspan="1">0.007 (0.003)</td>
</tr>
<tr><td colspan="4" align="left" rowspan="1">coverage 48%</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Arbitrary</td>
<td align="left" rowspan="1" colspan="1">0.079 (0.032)</td>
<td align="left" rowspan="1" colspan="1">0.087 (0.035)</td>
<td align="left" rowspan="1" colspan="1">0.025 (0.010)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Domiciliary</td>
<td align="left" rowspan="1" colspan="1">0.290 (0.008)</td>
<td align="left" rowspan="1" colspan="1">0.335 (0.004)</td>
<td align="left" rowspan="1" colspan="1">0.089 (0.003)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Student</td>
<td align="left" rowspan="1" colspan="1">0.146 (0.011)</td>
<td align="left" rowspan="1" colspan="1">0.282 (0.012)</td>
<td align="left" rowspan="1" colspan="1">0.085 (0.006)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Employee</td>
<td align="left" rowspan="1" colspan="1">0.038 (0.013)</td>
<td align="left" rowspan="1" colspan="1">0.028 (0.010)</td>
<td align="left" rowspan="1" colspan="1">0.009 (0.003)</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt103"><label></label>
<p><sup>*</sup>
All data are given as average and standard deviation (in paren.) of 64 simulation runs.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>We have shown the results of vaccinating a single group. We demonstrated that, with the same distribution speed (96% of residents are covered in 180 days), targeting offices workers for vaccination resulted in the greatest reduction in the AR, followed by targeting students and then domiciliaries. Through this same demonstration, we also evaluated the degree to which the collective immunity would be degraded due to a delay in the design or manufacture of a vaccine. <xref ref-type="fig" rid="pone-0072866-g005">Figure 5</xref>
shows the number of patients over time for three different transmission efficiencies and for the 0-th, 30-th, and 60-th day as the start date of the administration of vaccine. Each panel corresponds to 64 simulation runs (8 different choices of initial patients and 8 different random seeds). The differences among the three cases is small. A small epidemic occurred (AR is 0.3–0.4%) and about 90 days were needed for extinction of the transmission chains, even if the administration of vaccine was begun at the start of the simulation. If the initiation of vaccination was delayed, the scale of the outbreak increased: AR was 1–2% and the extinction was delayed until the 150-th day. This reduction of the AR to a small value despite a delay indicates that vaccination is sufficient to protect against low-pathogenic strains, whereas aggressive intervention (e.g., school closure and traffic blockades) should be combined with vaccination against highly pathogenic strains, since vaccination alone will not lead to early extinction.</p>
<fig id="pone-0072866-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0072866.g005</object-id>
<label>Figure 5</label>
<caption><title>The dependence of <inline-formula><inline-graphic xlink:href="pone.0072866.e047.jpg"></inline-graphic>
</inline-formula>
on the onset of vaccination.</title>
<p>To show the effect of the selection of initial patients and the stochasticity of transmission, 64 runs (8 different choices of initial patients and 8 different random seeds) are carried out in each case. The administration covers all residents in 180 days and vaccinations are carried out in the order of TrStuEmp, the rest of StuEmp, and finally Dom.</p>
</caption>
<graphic xlink:href="pone.0072866.g005"></graphic>
</fig>
</sec>
<sec id="s4"><title>Discussion</title>
<p>We have showed that the AR can be reduced to <inline-formula><inline-graphic xlink:href="pone.0072866.e048.jpg"></inline-graphic>
</inline-formula>
or less if students and employees are intensively vaccinated in the first 90 days. This is the result of an intervention program that relies solely on vaccinations, and the AR can be further reduced by individual protection efforts (e.g., wearing masks and avoiding crowded places). If the encountered virus is not highly pathogenic, this value is acceptable. In this case, the goal of intervention is to avoid an excess of patients going to medical practitioners. However, early extinction of transmission chains is required in highly pathogenic cases, and vaccination alone is not sufficient. To achieve early extinction by using only the collective immunity induced by vaccinations, administration of the vaccine would need to be carried out at least three times as fast as the typical speed.</p>
<p>It would be worthwhile to interpret our simulation results along with the situation of the 2009 pandemic (2009pdm) in Japan. Inaida et al. <xref rid="pone.0072866-Inaida1" ref-type="bibr">[24]</xref>
studied the 2009 pandemic geographically and temporally in four metropolitan areas of Japan, analysing time courses obtained from about 5,000 sentinel observation points. The 2009pdm spread similarly in four metropolitan areas and the growth of the number of patients is relatively rapid compared to the surrounding rural areas. However, a correlation between the AR and the rate of large families was recognized in Tokyo and Nagoya, but not in Osaka and Fukuoka. Different transmission models of households for central and suburb areas may need to be introduced into the simulation, since such large families mainly exist in suburban areas.</p>
<p>Our simulations indicate that a two-month delay in vaccination retains an effect on the enhancement of collective immunity to a certain extent. Retracing the timeline of the 2009pdm, we know that this delay is allowable. The analysis of Inadia et al. <xref rid="pone.0072866-Inaida1" ref-type="bibr">[24]</xref>
showed that 2009pdm spread sporadically in Japan during its early stages. The pandemic virus landed during the end of April to the middle of May. The number of reported cases per sentinel site exceeded unity, which is used as a working definition of the onset of an epidemic, in the beginning of August and exhibited further slow growth for more than 6 weeks. The distribution of vaccines began in October. Therefore, the delay of vaccine distribution from the onset of the epidemic was around one month. Nonetheless, the sporadic transmission was perhaps due to the onset in summer, which is unlike seasonal influenza.</p>
<p>The vaccination plan <xref rid="pone.0072866-Ministry1" ref-type="bibr">[21]</xref>
in 2009/2010 season gave a higher priority to certain groups including pregnant women, individuals with underlying diseases, children, and elderly persons. The rest of the population was intended to be covered by the remainder of the domestic vaccine, and supplemented by imported vaccine. The estimated numbers of courses distributed to respective groups <xref rid="pone.0072866-Ministry2" ref-type="bibr">[22]</xref>
in November and December indicate that the distribution of vaccines in the early stage was actually focused on the elderly and high-risk persons as had been planned. In November, the number of distributed courses was 920,000 to people with underlying disease, 1,335,000 to people without underlying diseases and 1,011,000 to the age group above 65, respectively. In December, these numbers were 1,468,000, 3,956,000 and 2,224,000, respectively, where distribution to small children (≤8 years) contributed to a large number of people without underlying disease receiving the vaccine in December. The distribution to the age group of 9–65 is almost the same in these two months (73 and 82.6 thousands, respectively). From this document, we also know that the planned rate and the actual rate are similar until the middle of December, after which the actual rate is less than the planned rate. According to reports of a learned discussion <xref rid="pone.0072866-Ministry3" ref-type="bibr">[23]</xref>
, 2009 flu may have influenced people to avoid vaccination, resulting in an oversupply after January.</p>
</sec>
<sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0072866.s001"><label>File S1</label>
<caption><p><bold>Enhancement of collective immunity by selective vaccination against an emerging influenza pandemic.</bold>
Contains Figure S1, Figure S2, Figure S3, and Table S1. <bold>Figure S1. Pseudo code of a single step of the simulation.</bold>
All information on the simulated urban area is contained in the structure instance city, which is located in shared memory. Iterations and branches are in a Fortran-like code, but the structure name and its field are split by a dot. Simulations are carried out in parallel based on OpenMP, and iterations marked by !OMP DO are adequately split by the compiler and carried out in parallel. Calculation of the transition probability places{key = v}.pr is implemented to conform to the diagram of Eq. (2). <bold>Figure S2. Age-specific distribution of the population of Tokyo in 2005.</bold>
We sampled from this distribution to obtain the ages of individuals, and their roles were assigned according to their ages. The proportions of roles in the simulation are represented by different colors (blue: students, red: employees, and yellow: domiciliaries). <bold>Figure S3. Distribution of corporation sizes.</bold>
The rank in the corporation size versus the number of employees. <bold>Table S1. List of parameters configuring simulation.</bold>
</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0072866.s001.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ref-list><title>References</title>
<ref id="pone.0072866-Kerkhove1"><label>1</label>
<mixed-citation publication-type="journal"><name><surname>Kerkhove</surname>
<given-names>MDV</given-names>
</name>
, <name><surname>Vandemaele</surname>
<given-names>KAH</given-names>
</name>
, <name><surname>Shinde</surname>
<given-names>V</given-names>
</name>
, <name><surname>Jaramillo-Gutierrez</surname>
<given-names>G</given-names>
</name>
, <name><surname>Koukounari</surname>
<given-names>A</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Risk factors for severe outcomes following 2010 influenza A (H1N1)</article-title>
. <source>PLoS Med</source>
<volume>8</volume>
.</mixed-citation>
</ref>
<ref id="pone.0072866-Sha1"><label>2</label>
<mixed-citation publication-type="journal"><name><surname>Sha</surname>
<given-names>L</given-names>
</name>
, <name><surname>Zhu</surname>
<given-names>R</given-names>
</name>
, <name><surname>Cao</surname>
<given-names>L</given-names>
</name>
, <name><surname>Yuan</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Li</surname>
<given-names>Y</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Analysis of clinical manifestations of hospitalized children infected with seasonal influenza A virus and 2009 novel influenza A (H1N1) virus in Beijing</article-title>
. <source>Chinese Journal of Pediatrics</source>
<volume>49</volume>
: <fpage>539</fpage>
–<lpage>544</lpage>
.<pub-id pub-id-type="pmid">22088186</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-World1"><label>3</label>
<mixed-citation publication-type="other">World Health Organization (2011) H5N1 avian influenza: timeline of major events. Global Alert and Response (GAR). Available: <ext-link ext-link-type="uri" xlink:href="http://web.archive.org/web/20110903044611/http://www.who.int/csr/disease/avian-influenza/ai-t">http://web.archive.org/web/20110903044611/http://www.who.int/csr/disease/avian-influenza/ai-t</ext-link>
Original source disappeared. Accessed: 15 Aug 2013.</mixed-citation>
</ref>
<ref id="pone.0072866-Ferguson1"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
, <name><surname>Cummings</surname>
<given-names>DAT</given-names>
</name>
, <name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <name><surname>Fraser</surname>
<given-names>C</given-names>
</name>
, <name><surname>Riley</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
(<year>2005</year>
) <article-title>Strategies for containing an emerging influenza pandemic in southeast asia</article-title>
. <source>Nature</source>
<volume>437</volume>
: <fpage>209</fpage>
–<lpage>214</lpage>
.<pub-id pub-id-type="pmid">16079797</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-Longini1"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
, <name><surname>Nizam</surname>
<given-names>A</given-names>
</name>
, <name><surname>Xu</surname>
<given-names>S</given-names>
</name>
, <name><surname>Ungchusak</surname>
<given-names>K</given-names>
</name>
, <name><surname>Hanshaoworakul</surname>
<given-names>W</given-names>
</name>
, <etal>et al</etal>
(<year>2005</year>
) <article-title>Containing pandemic influenza at the source</article-title>
. <source>Science</source>
<volume>309</volume>
: <fpage>1084</fpage>
–<lpage>1087</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0072866-Ferguson2"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
, <name><surname>Cummings</surname>
<given-names>DAT</given-names>
</name>
, <name><surname>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>
(<year>2006</year>
) <article-title>Strategies for mitigating an influenza pandemic</article-title>
. <source>Nature</source>
<volume>442</volume>
: <fpage>448</fpage>
–<lpage>452</lpage>
.<pub-id pub-id-type="pmid">16642006</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-Sato1"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Sato</surname>
<given-names>H</given-names>
</name>
, <name><surname>Nakada</surname>
<given-names>H</given-names>
</name>
, <name><surname>Yamaguchi</surname>
<given-names>R</given-names>
</name>
, <name><surname>Imoto</surname>
<given-names>S</given-names>
</name>
, <name><surname>Miyano</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>When should we intervene to control the 2009 influenza A (H1N1) pandemic?</article-title>
<source>Euro Surveill</source>
<volume>15</volume>
.</mixed-citation>
</ref>
<ref id="pone.0072866-Tokyo1"><label>8</label>
<mixed-citation publication-type="other">Tokyo metropolitan government (2008) Tokyo Statistical Yearbook, chapter 2-2 Population by 5 year age group and sex (1920–2005). Available: <ext-link ext-link-type="uri" xlink:href="http://www.toukei.metro.tokyo.jp/tnenkan/2008/tn08qyti0510b.htm">http://www.toukei.metro.tokyo.jp/tnenkan/2008/tn08qyti0510b.htm</ext-link>
. Accessed: 15 Aug 2013.</mixed-citation>
</ref>
<ref id="pone.0072866-Kermack1"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Kermack</surname>
<given-names>WO</given-names>
</name>
, <name><surname>McKendrick</surname>
<given-names>AG</given-names>
</name>
(<year>1927</year>
) <article-title>Contributions to the mathematical theory of epidemics</article-title>
. <source>Proc Royal Soc Ser A</source>
<volume>115</volume>
: <fpage>700</fpage>
–<lpage>721</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0072866-Dawood1"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Dawood</surname>
<given-names>F</given-names>
</name>
, <name><surname>Jain</surname>
<given-names>S</given-names>
</name>
, <name><surname>Finelli</surname>
<given-names>L</given-names>
</name>
, <name><surname>Shaw</surname>
<given-names>M</given-names>
</name>
, <name><surname>Lindstrom</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Emergence of a novel swine-origin influenza A (H1N1) virus in humans</article-title>
. <source>N Engl J Med</source>
<volume>360</volume>
: <fpage>2605</fpage>
–<lpage>2615</lpage>
.<pub-id pub-id-type="pmid">19423869</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-Germann1"><label>11</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>Ira</surname>
<given-names>M</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>J</given-names>
</name>
, <name><surname>Macken</surname>
<given-names>CA</given-names>
</name>
(<year>2006</year>
) <article-title>Mitigation strategies for pandemic influenza in the united states</article-title>
. <source>PNAS</source>
<volume>105</volume>
: <fpage>5935</fpage>
–<lpage>5940</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0072866-Cooley1"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Cooley</surname>
<given-names>P</given-names>
</name>
, <name><surname>Brown</surname>
<given-names>S</given-names>
</name>
, <name><surname>Cajka</surname>
<given-names>J</given-names>
</name>
, <name><surname>Chasteen</surname>
<given-names>B</given-names>
</name>
, <name><surname>Ganapathi</surname>
<given-names>L</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>The role of subway travel in an influenza epidemic: A new york city simulation</article-title>
. <source>Journal of Urban Health</source>
<volume>88</volume>
: <fpage>982</fpage>
–<lpage>995</lpage>
.<pub-id pub-id-type="pmid">21826584</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-Feigin1"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Feigin</surname>
<given-names>R</given-names>
</name>
, <name><surname>Baker</surname>
<given-names>C</given-names>
</name>
, <name><surname>Herwaldt</surname>
<given-names>L</given-names>
</name>
, <name><surname>Lampe</surname>
<given-names>R</given-names>
</name>
, <name><surname>Mason</surname>
<given-names>E</given-names>
</name>
, <etal>et al</etal>
(<year>1982</year>
) <article-title>Epidemic meningococcal disease in an elementary-school classroom</article-title>
. <source>N Engl J Med</source>
<volume>307</volume>
: <fpage>1255</fpage>
–<lpage>1257</lpage>
.<pub-id pub-id-type="pmid">6813733</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-Dick1"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Dick</surname>
<given-names>E</given-names>
</name>
, <name><surname>Jennings</surname>
<given-names>L</given-names>
</name>
, <name><surname>Mink</surname>
<given-names>K</given-names>
</name>
, <name><surname>Wartgow</surname>
<given-names>C</given-names>
</name>
, <name><surname>Inhorn</surname>
<given-names>S</given-names>
</name>
(<year>1987</year>
) <article-title>Aerosol transmission of rhinovirus colds</article-title>
. <source>J Infect Dis</source>
<volume>156</volume>
: <fpage>442</fpage>
–<lpage>448</lpage>
.<pub-id pub-id-type="pmid">3039011</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-Duguid1"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Duguid</surname>
<given-names>J</given-names>
</name>
(<year>1946</year>
) <article-title>The size and duration of air-carriage of respiratory droplets and droplet nucleii</article-title>
. <source>J Hyg (Lond)</source>
<volume>44</volume>
: <fpage>471</fpage>
–<lpage>479</lpage>
.<pub-id pub-id-type="pmid">20475760</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-Siegel1"><label>16</label>
<mixed-citation publication-type="other">Siegel J, Rhinehart E, Jackson M, Chiarello L (2007) 2007 Guideline for Isolation Precautions: Preventing Transmission of Infectious Agents in Healthcare Settings. The Healthcare Infection Control Practices Advisory Committee.</mixed-citation>
</ref>
<ref id="pone.0072866-Jefferson1"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Jefferson</surname>
<given-names>T</given-names>
</name>
, <name><surname>Rivetti</surname>
<given-names>A</given-names>
</name>
, <name><surname>Harnden</surname>
<given-names>A</given-names>
</name>
, <name><surname>Pietrantonj</surname>
<given-names>CD</given-names>
</name>
, <name><surname>Demicheli</surname>
<given-names>V</given-names>
</name>
(<year>2008</year>
) <article-title>Vaccines for preventing influenza in healthy children (review)</article-title>
. <source>Cochrane Database of Systematic Reviews</source>
<volume>2008</volume>
: <fpage>Art. No.: CD004879</fpage>
.</mixed-citation>
</ref>
<ref id="pone.0072866-Jefferson2"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Jefferson</surname>
<given-names>T</given-names>
</name>
, <name><surname>Pietrantonj</surname>
<given-names>CD</given-names>
</name>
, <name><surname>Rivetti</surname>
<given-names>A</given-names>
</name>
, <name><surname>Bawazeer</surname>
<given-names>GA</given-names>
</name>
, <name><surname>Al-Ansary</surname>
<given-names>LA</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Vaccines for preventing influenza in healthy adults (review)</article-title>
. <source>Cochrane Database of Systematic Reviews</source>
<volume>2010</volume>
: <fpage>Art. No.: CD001269</fpage>
.</mixed-citation>
</ref>
<ref id="pone.0072866-Jefferson3"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Jefferson</surname>
<given-names>T</given-names>
</name>
, <name><surname>Pietrantonj</surname>
<given-names>CD</given-names>
</name>
, <name><surname>Al-Ansary</surname>
<given-names>L</given-names>
</name>
, <name><surname>Ferroni</surname>
<given-names>E</given-names>
</name>
, <name><surname>Thorning</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Vaccines for preventing influenza in the elderly (review)</article-title>
. <source>Cochrane Database of Systematic Reviews</source>
<volume>2010</volume>
: <fpage>Art. No.: CD004876</fpage>
.</mixed-citation>
</ref>
<ref id="pone.0072866-Geoffery1"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Geoffery</surname>
<given-names>AW</given-names>
</name>
, <name><surname>Perter</surname>
<given-names>GS</given-names>
</name>
(<year>2010</year>
) <article-title>Vaccine epidemiology: Efficacy, effectiveness, and the translational research roadmap</article-title>
. <source>The Journal of Infectious Diseases</source>
<volume>201</volume>
: <fpage>1607</fpage>
–<lpage>1610</lpage>
.<pub-id pub-id-type="pmid">20402594</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0072866-Ministry1"><label>21</label>
<mixed-citation publication-type="other">Ministry of Health, Labour and Welfare Japan. (2009) Current standard vaccination schedule (Japanese). Available: <ext-link ext-link-type="uri" xlink:href="http://www.mhlw.go.jp/kinkyu/kenkou/influenza/dl/infu091217-01.pdf">http://www.mhlw.go.jp/kinkyu/kenkou/influenza/dl/infu091217-01.pdf</ext-link>
. Accessed: 15 Aug 2013.</mixed-citation>
</ref>
<ref id="pone.0072866-Ministry2"><label>22</label>
<mixed-citation publication-type="other">Ministry of Health, Labour and Welfare Japan. (2010) Study Committee adverse reactions after vaccination flu Fifth 2009 (Japanese). Available: <ext-link ext-link-type="uri" xlink:href="http://www.mhlw.go.jp/shingi/2010/02/s0212-6.html">http://www.mhlw.go.jp/shingi/2010/02/s0212-6.html</ext-link>
. Accessed: 15 Aug 2013.</mixed-citation>
</ref>
<ref id="pone.0072866-Ministry3"><label>23</label>
<mixed-citation publication-type="other">Ministry of Health, Labour and Welfare Japan. (2012) Study committee report concerning distribution improvement of novel influenza vaccine (Japanese). Available: <ext-link ext-link-type="uri" xlink:href="http://www.mhlw.go.jp/stf/shingi/2r9852000002oeqs-att/2r9852000002oext.pdf">http://www.mhlw.go.jp/stf/shingi/2r9852000002oeqs-att/2r9852000002oext.pdf</ext-link>
. Accessed: 15 Aug 2013.</mixed-citation>
</ref>
<ref id="pone.0072866-Inaida1"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>Inaida</surname>
<given-names>S</given-names>
</name>
, <name><surname>Yasui</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Tada</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Taniguchi</surname>
<given-names>K</given-names>
</name>
, <name><surname>Okabe</surname>
<given-names>N</given-names>
</name>
(<year>2011</year>
) <article-title>Geographic trends and spread of the pandemic (h1n1) 2009 in the metropolitan areas of japan studied from the national sentinel data</article-title>
. <source>Japanese Journal of Infectious Disease</source>
<volume>64</volume>
: <fpage>473</fpage>
–<lpage>481</lpage>
.</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 000D720 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000D720 | 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é= |texte= }}
![]() | This area was generated with Dilib version V0.6.34. | ![]() |