<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.0 20040830//EN" "http://dtd.nlm.nih.gov/publishing/2.0/journalpublishing.dtd">
<?covid-19-tdm?>
<article article-type="letter" dtd-version="2.0" xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-id journal-id-type="publisher-id">JMIR</journal-id>
      <journal-id journal-id-type="nlm-ta">J Med Internet Res</journal-id>
      <journal-title>Journal of Medical Internet Research</journal-title>
      <issn pub-type="epub">1438-8871</issn>
      <publisher>
        <publisher-name>JMIR Publications</publisher-name>
        <publisher-loc>Toronto, Canada</publisher-loc>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="publisher-id">v24i7e37744</article-id>
      <article-id pub-id-type="pmid">35679053</article-id>
      <article-id pub-id-type="doi">10.2196/37744</article-id>
      <article-categories>
        <subj-group subj-group-type="heading">
          <subject>Research Letter</subject>
        </subj-group>
        <subj-group subj-group-type="article-type">
          <subject>Research Letter</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>Using Machine Learning to Efficiently Vaccinate Homebound Patients Against COVID-19: A Real-time Immunization Campaign</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="editor">
          <name>
            <surname>Leung</surname>
            <given-names>Tiffany</given-names>
          </name>
        </contrib>
      </contrib-group>
      <contrib-group>
        <contrib contrib-type="reviewer">
          <name>
            <surname>Rahaman</surname>
            <given-names>Oindrila</given-names>
          </name>
        </contrib>
        <contrib contrib-type="reviewer">
          <name>
            <surname>Calisto</surname>
            <given-names>Francisco Maria</given-names>
          </name>
        </contrib>
      </contrib-group>
      <contrib-group>
        <contrib id="contrib1" contrib-type="author" corresp="yes" equal-contrib="yes">
          <name name-style="western">
            <surname>Kumar</surname>
            <given-names>Anish</given-names>
          </name>
          <degrees>BA</degrees>
          <xref rid="aff1" ref-type="aff">1</xref>
          <address>
            <institution>Icahn School of Medicine at Mount Sinai</institution>
            <addr-line>1 Gustave L. Levy Place</addr-line>
            <addr-line>New York, NY, 10029</addr-line>
            <country>United States</country>
            <phone>1 9085668399</phone>
            <email>anish.kumar@icahn.mssm.edu</email>
          </address>
          <ext-link ext-link-type="orcid">https://orcid.org/0000-0001-9183-7161</ext-link>
        </contrib>
        <contrib id="contrib2" contrib-type="author" equal-contrib="yes">
          <name name-style="western">
            <surname>Ren</surname>
            <given-names>Jennifer</given-names>
          </name>
          <degrees>BA</degrees>
          <xref rid="aff1" ref-type="aff">1</xref>
          <ext-link ext-link-type="orcid">https://orcid.org/0000-0001-8309-8928</ext-link>
        </contrib>
        <contrib id="contrib3" contrib-type="author">
          <name name-style="western">
            <surname>Ornstein</surname>
            <given-names>Katherine A</given-names>
          </name>
          <degrees>MPH, PhD</degrees>
          <xref rid="aff2" ref-type="aff">2</xref>
          <xref rid="aff3" ref-type="aff">3</xref>
          <ext-link ext-link-type="orcid">https://orcid.org/0000-0001-6270-4423</ext-link>
        </contrib>
        <contrib id="contrib4" contrib-type="author">
          <name name-style="western">
            <surname>Gliatto</surname>
            <given-names>Peter M</given-names>
          </name>
          <degrees>MD</degrees>
          <xref rid="aff4" ref-type="aff">4</xref>
          <ext-link ext-link-type="orcid">https://orcid.org/0000-0003-4405-4230</ext-link>
        </contrib>
      </contrib-group>
      <aff id="aff1">
        <label>1</label>
        <institution>Icahn School of Medicine at Mount Sinai</institution>
        <addr-line>New York, NY</addr-line>
        <country>United States</country>
      </aff>
      <aff id="aff2">
        <label>2</label>
        <institution>Department of Geriatrics and Palliative Medicine</institution>
        <institution>Icahn School of Medicine at Mount Sinai</institution>
        <addr-line>New York, NY</addr-line>
        <country>United States</country>
      </aff>
      <aff id="aff3">
        <label>3</label>
        <institution>Institute for Care Innovations at Home</institution>
        <institution>Icahn School of Medicine at Mount Sinai</institution>
        <addr-line>New York, NY</addr-line>
        <country>United States</country>
      </aff>
      <aff id="aff4">
        <label>4</label>
        <institution>Department of Medicine</institution>
        <institution>Icahn School of Medicine at Mount Sinai</institution>
        <addr-line>New York, NY</addr-line>
        <country>United States</country>
      </aff>
      <author-notes>
        <corresp>Corresponding Author: Anish Kumar <email>anish.kumar@icahn.mssm.edu</email></corresp>
      </author-notes>
      <pub-date pub-type="collection">
        <month>7</month>
        <year>2022</year>
      </pub-date>
      <pub-date pub-type="epub">
        <day>12</day>
        <month>7</month>
        <year>2022</year>
      </pub-date>
      <volume>24</volume>
      <issue>7</issue>
      <elocation-id>e37744</elocation-id>
      <history>
        <date date-type="received">
          <day>4</day>
          <month>3</month>
          <year>2022</year>
        </date>
        <date date-type="rev-request">
          <day>11</day>
          <month>5</month>
          <year>2022</year>
        </date>
        <date date-type="rev-recd">
          <day>4</day>
          <month>6</month>
          <year>2022</year>
        </date>
        <date date-type="accepted">
          <day>7</day>
          <month>6</month>
          <year>2022</year>
        </date>
      </history>
      <copyright-statement>©Anish Kumar, Jennifer Ren, Katherine A Ornstein, Peter M Gliatto. Originally published in the Journal of Medical Internet Research (https://www.jmir.org), 12.07.2022.</copyright-statement>
      <copyright-year>2022</copyright-year>
      <license license-type="open-access" xlink:href="https://creativecommons.org/licenses/by/4.0/">
        <p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in the Journal of Medical Internet Research, is properly cited. The complete bibliographic information, a link to the original publication on https://www.jmir.org/, as well as this copyright and license information must be included.</p>
      </license>
      <self-uri xlink:href="https://www.jmir.org/2022/7/e37744" xlink:type="simple"/>
      <kwd-group>
        <kwd>home care</kwd>
        <kwd>covid</kwd>
        <kwd>vaccination</kwd>
        <kwd>COVID-19</kwd>
        <kwd>machine learning</kwd>
        <kwd>vaccine</kwd>
        <kwd>geographic cluster</kwd>
        <kwd>patient data</kwd>
        <kwd>clustering algorithm</kwd>
        <kwd>geospatial</kwd>
        <kwd>digital surveillance</kwd>
        <kwd>public health</kwd>
        <kwd>logistic operation</kwd>
        <kwd>electronic health record integration</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec sec-type="introduction">
      <title>Introduction</title>
      <p>In March 2021, the New York City Department of Health and Mental Hygiene (NYCDOHMH) received a supply of single-dose Janssen COVID-19 vaccines to vaccinate the city’s patients who are homebound, but they needed assistance to identify and reach this population [<xref ref-type="bibr" rid="ref1">1</xref>,<xref ref-type="bibr" rid="ref2">2</xref>]. Mount Sinai Visiting Doctors (MSVD) is a large home-based primary care program serving more than 1200 patients who are homebound throughout Manhattan. We partnered with NYCDOHMH to vaccinate patients in our program. The administrative team generally schedules routine home visits based on zip codes, using 12 unique catchment areas covering all of Manhattan. This existing zoning system was inadequate for vaccination purposes for several reasons, mainly because these zones were too large. Furthermore, the additional task of manually scheduling patients by zone would be overwhelming for administrative staff, especially given the temporal constraints of the Janssen vaccine (ie, doses expired in 6 hours). In response, we developed a system to geographically cluster patients to efficiently vaccinate our homebound patients against COVID-19.</p>
    </sec>
    <sec sec-type="methods">
      <title>Methods</title>
      <sec>
        <title>Ethics Approval</title>
        <p>The Icahn School of Medicine at Mount Sinai’s Program for the Protection of Human Subjects approved and granted a waiver of consent for this secondary data analysis study (STUDY- 21-00157) which was conducted in accordance with the Helsinki Declaration.</p>
      </sec>
      <sec>
        <title>Overview</title>
        <p>We developed a software program that takes a cohort of unvaccinated patients in the MSVD Program as input and assigns each patient a cluster number. We used Python 3 (Python Software Foundation) [<xref ref-type="bibr" rid="ref3">3</xref>] to process patient addresses, group them into clusters, and export these clusters into a standard database format and onto a map. Specifically, we used an open-source implementation of a modified unsupervised K-means clustering machine learning algorithm to group patients [<xref ref-type="bibr" rid="ref4">4</xref>].</p>
        <p>The first step was to convert patient households to spatial representations. Using the Google Maps Geocoding application programming interface (API), we obtained latitude/longitude coordinates of patient residences, which served as the input data to our algorithm.</p>
        <p>These coordinates were then fed to the modified unsupervised K-means clustering algorithm. The standard K-means clustering algorithm clusters data points into <italic>K</italic> groups, while the modified algorithm allows one to impose constraints on cluster size, ensuring each group contains a number of data points within a specified range [<xref ref-type="bibr" rid="ref4">4</xref>]. This was important for our process to avoid assigning an excessive number of patients to any one cluster; although many patients live close to each other, provider routes could not exceed 6 hours. Targeting approximately 5 patients per provider and bearing in mind that not all identified patients would desire and be available when the vaccine was offered by the scheduling team, we enforced a minimum of 12 and a maximum of 15 patients per cluster.</p>
        <p>Next, we projected the data points onto a geospatial visualization for easy interpretability by the scheduling team and providers. Using the Google Maps API, we generated an HTML page that contained pins with distinct numbers and colors representing the cluster to which a patient belonged as well as the patient’s name.</p>
        <p>Our process was iterative: every week after patients were scheduled and immunized, we updated our census to reflect the remaining number of patients who were unvaccinated. The team called patients to confirm appointments and updated the patient list. We then re-ran our clustering process, integrating any changes to cluster size preferences, for the next round of immunizations for patients who are homebound (<xref rid="figure1" ref-type="fig">Figure 1</xref>).</p>
        <fig id="figure1" position="float">
          <label>Figure 1</label>
          <caption>
            <p>An example of a vaccination campaign, using mock patient data, for patients who are homebound in the Mount Sinai Visiting Doctors Program routing process. It was iterative week-to-week and included the following steps: (1) patients who were unvaccinated were identified and clustered into groups, (2) patient addresses were pinned on a map with identifying group labels, (3) patients were called for vaccine scheduling according to group number, and (4) the patient list was updated to reflect who remained unvaccinated.</p>
          </caption>
          <graphic xlink:href="jmir_v24i7e37744_fig1.png" alt-version="no" mimetype="image" position="float" xlink:type="simple"/>
        </fig>
      </sec>
    </sec>
    <sec sec-type="results">
      <title>Results</title>
      <p><xref ref-type="table" rid="table1">Table 1</xref> summarizes the demographic information of the patients and the results of our campaign between March and April 2021, which involved about 100 vaccination clusters and routes over 6 weeks. On average, we vaccinated 5.6 patients per provider per day, averaging 22.1 total patients per day [<xref ref-type="bibr" rid="ref5">5</xref>]. Each provider accomplished their assigned route within the time constraints of the Janssen vaccine, and no doses were wasted.</p>
      <table-wrap position="float" id="table1">
        <label>Table 1</label>
        <caption>
          <p>Demographics and statistics from our COVID-19 homebound vaccination campaign of Manhattan-based patients in the Mount Sinai Visiting Doctors Program.</p>
        </caption>
        <table width="1000" cellpadding="5" cellspacing="0" border="1" rules="groups" frame="hsides">
          <col width="30"/>
          <col width="30"/>
          <col width="30"/>
          <col width="780"/>
          <col width="130"/>
          <thead>
            <tr valign="top">
              <td colspan="4"> </td>
              <td>Value</td>
            </tr>
          </thead>
          <tbody>
            <tr valign="top">
              <td colspan="5">
                <bold>Patient demographics</bold>
              </td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td colspan="3">
                <bold>Patients who are homebound</bold>
              </td>
              <td>428</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td colspan="2">Average age (years)</td>
              <td>83.9</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td colspan="2">Average Elixhauser comorbidity score [<xref ref-type="bibr" rid="ref6">6</xref>]<sup>a</sup></td>
              <td>3.8</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td colspan="3">
                <bold>Sex, n (%)</bold>
              </td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td> </td>
              <td>Female</td>
              <td>323 (75.5)</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td> </td>
              <td>Male</td>
              <td>105 (24.5)</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td colspan="3">
                <bold>Racial/ethnic identity, n (%)</bold>
              </td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td> </td>
              <td>White</td>
              <td>148 (34.6)</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td> </td>
              <td>Black or African American</td>
              <td>63 (14.7)</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td> </td>
              <td>Asian</td>
              <td>15 (3.5)</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td> </td>
              <td>Hispanic</td>
              <td>108 (25.2)</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td> </td>
              <td>Other</td>
              <td>86 (20.1)</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td> </td>
              <td> </td>
              <td>Unknown</td>
              <td>8 (1.9)</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td colspan="3">Patient’s family members and caregivers, n</td>
              <td>92</td>
            </tr>
            <tr valign="top">
              <td colspan="5">
                <bold>Vaccination campaign statistics</bold>
              </td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td colspan="3">Providers vaccinating per day (n), range</td>
              <td>3-6</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td colspan="3">Average number of patients vaccinated per day</td>
              <td>22.1</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td colspan="3">Average duration of provider time spent vaccinating (hours)</td>
              <td>4.6</td>
            </tr>
            <tr valign="top">
              <td> </td>
              <td colspan="3">Average duration of individual vaccination (including transit time, vaccine administration, and 15-minute postvaccination observation time; minutes)</td>
              <td>52</td>
            </tr>
          </tbody>
        </table>
        <table-wrap-foot>
          <fn id="table1fn1">
            <p><sup>a</sup>Elixhauser scores were available for 372 of the patients who are homebound.</p>
          </fn>
        </table-wrap-foot>
      </table-wrap>
    </sec>
    <sec sec-type="discussion">
      <title>Discussion</title>
      <p>The new tool optimized logistic operations for an acute public health intervention while minimizing wasted resources. The model was later used to quickly deploy booster immunizations during the surge of the Omicron subvariant from December 2021 to January 2022. Future research will consider the ability to create routes with ordered stops given a provider’s choice of transportation and more flexibility in the cluster size depending on the population density of a given region. This feature is especially important in densely populated cities, where providers are unlikely to travel by car. Additional future work can include electronic health record integration for more streamlined access and allowing scheduling teams to directly recluster patients.</p>
      <p>Ultimately, the newly developed approach was instrumental in maximizing efficiency and minimizing vaccine waste, suggesting its potential for future use in home-based health care or other public health interventions. The success of this project further demonstrates the value of novel technological approaches in improving the efficiency of clinical operations.</p>
    </sec>
  </body>
  <back>
    <app-group/>
    <glossary>
      <title>Abbreviations</title>
      <def-list>
        <def-item>
          <term id="abb1">API</term>
          <def>
            <p>application programming interface</p>
          </def>
        </def-item>
        <def-item>
          <term id="abb2">MSVD</term>
          <def>
            <p>Mount Sinai Visiting Doctors</p>
          </def>
        </def-item>
        <def-item>
          <term id="abb3">NYCDOHMH</term>
          <def>
            <p>New York City Department of Health and Mental Hygiene</p>
          </def>
        </def-item>
      </def-list>
    </glossary>
    <fn-group>
      <fn fn-type="conflict">
        <p>None declared.</p>
      </fn>
    </fn-group>
    <ref-list>
      <ref id="ref1">
        <label>1</label>
        <nlm-citation citation-type="web">
          <person-group person-group-type="author">
            <name name-style="western">
              <surname>Otterman</surname>
              <given-names>S</given-names>
            </name>
          </person-group>
          <article-title>Homebound and waiting: vaccine push misses some older New Yorkers</article-title>
          <source>The New York Times</source>
          <year>2021</year>
          <month>04</month>
          <day>23</day>
          <access-date>2021-11-18</access-date>
          <comment>
            <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://www.nytimes.com/2021/04/23/nyregion/homebound-vaccine-nyc.html">https://www.nytimes.com/2021/04/23/nyregion/homebound-vaccine-nyc.html</ext-link>
          </comment>
        </nlm-citation>
      </ref>
      <ref id="ref2">
        <label>2</label>
        <nlm-citation citation-type="web">
          <person-group person-group-type="author">
            <name name-style="western">
              <surname>Baulkman</surname>
              <given-names>J</given-names>
            </name>
          </person-group>
          <article-title>Vaccinating homebound patients is an uphill battle</article-title>
          <source>MDedge Family Medicine</source>
          <year>2021</year>
          <month>04</month>
          <day>23</day>
          <access-date>2021-11-18</access-date>
          <comment>
            <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://www.mdedge.com/familymedicine/article/239068/coronavirus-updates/vaccinating-homebound-patients-uphill-battle">https://www.mdedge.com/familymedicine/article/239068/coronavirus-updates/vaccinating-homebound-patients-uphill-battle</ext-link>
          </comment>
        </nlm-citation>
      </ref>
      <ref id="ref3">
        <label>3</label>
        <nlm-citation citation-type="book">
          <person-group person-group-type="author">
            <name name-style="western">
              <surname>Van Rossum</surname>
              <given-names>G</given-names>
            </name>
            <name name-style="western">
              <surname>Drake</surname>
              <given-names>FL</given-names>
            </name>
          </person-group>
          <source>Python 3 Reference Manual</source>
          <year>2009</year>
          <publisher-loc>Scotts Valley, CA</publisher-loc>
          <publisher-name>CreateSpace</publisher-name>
        </nlm-citation>
      </ref>
      <ref id="ref4">
        <label>4</label>
        <nlm-citation citation-type="journal">
          <person-group person-group-type="author">
            <name name-style="western">
              <surname>Bradley</surname>
              <given-names>PS</given-names>
            </name>
            <name name-style="western">
              <surname>Bennet</surname>
              <given-names>KP</given-names>
            </name>
            <name name-style="western">
              <surname>Demiriz</surname>
              <given-names>A</given-names>
            </name>
          </person-group>
          <article-title>Constrained k-means clustering</article-title>
          <source>Microsoft Res</source>
          <year>2000</year>
          <fpage>1</fpage>
          <lpage>8</lpage>
        </nlm-citation>
      </ref>
      <ref id="ref5">
        <label>5</label>
        <nlm-citation citation-type="journal">
          <person-group person-group-type="author">
            <name name-style="western">
              <surname>Gliatto</surname>
              <given-names>P</given-names>
            </name>
            <name name-style="western">
              <surname>Franzosa</surname>
              <given-names>E</given-names>
            </name>
            <name name-style="western">
              <surname>Chavez</surname>
              <given-names>S</given-names>
            </name>
            <name name-style="western">
              <surname>Ng</surname>
              <given-names>A</given-names>
            </name>
            <name name-style="western">
              <surname>Kumar</surname>
              <given-names>A</given-names>
            </name>
            <name name-style="western">
              <surname>Ren</surname>
              <given-names>J</given-names>
            </name>
            <name name-style="western">
              <surname>DeCherrie</surname>
              <given-names>LV</given-names>
            </name>
            <name name-style="western">
              <surname>Ornstein</surname>
              <given-names>K</given-names>
            </name>
          </person-group>
          <article-title>Covid-19 vaccines for homebound patients and their caregivers</article-title>
          <source>New Engl J Med Catalyst</source>
          <year>2021</year>
          <fpage>14</fpage>
        </nlm-citation>
      </ref>
      <ref id="ref6">
        <label>6</label>
        <nlm-citation citation-type="journal">
          <person-group person-group-type="author">
            <name name-style="western">
              <surname>Elixhauser</surname>
              <given-names>A</given-names>
            </name>
            <name name-style="western">
              <surname>Steiner</surname>
              <given-names>C</given-names>
            </name>
            <name name-style="western">
              <surname>Harris</surname>
              <given-names>DR</given-names>
            </name>
            <name name-style="western">
              <surname>Coffey</surname>
              <given-names>RM</given-names>
            </name>
          </person-group>
          <article-title>Comorbidity measures for use with administrative data</article-title>
          <source>Med Care</source>
          <year>1998</year>
          <month>01</month>
          <volume>36</volume>
          <issue>1</issue>
          <fpage>8</fpage>
          <lpage>27</lpage>
          <pub-id pub-id-type="doi">10.1097/00005650-199801000-00004</pub-id>
          <pub-id pub-id-type="medline">9431328</pub-id>
        </nlm-citation>
      </ref>
    </ref-list>
  </back>
</article>
