Abstract
In response to unprecedented surges in the demand for hospital care during the SARS-CoV-2 pandemic, health systems have prioritized patients with COVID-19 to life-saving hospital care to the detriment of other patients. In contrast to these ad hoc policies, we develop a linear programming framework to optimally schedule elective procedures and allocate hospital beds among all planned and emergency patients to minimize years of life lost. Leveraging a large dataset of administrative patient medical records, we apply our framework to the National Health Service in England and show that an extra 50,750–5,891,608 years of life can be gained compared with prioritization policies that reflect those implemented during the pandemic. Notable health gains are observed for neoplasms, diseases of the digestive system, and injuries and poisoning. Our open-source framework provides a computationally efficient approximation of a large-scale discrete optimization problem that can be applied globally to support national-level care prioritization policies.
Similar content being viewed by others
Main
Health systems worldwide are struggling to provide hospital treatment during the surges in demand for emergency care caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, despite efforts to increase capacity1,2,3,4. To manage demand, many countries have prioritized patients with coronavirus disease 2019 (COVID-19), canceled elective (that is, planned) procedures5,6,7 and rationed access to life-saving care for all patients3,4. For example, policies in Italy, in March 2020, involved prioritizing intensive care to patients with COVID-19 above 70 years who previously had no more than one admission per year for a chronic illness (for example, exacerbated chronic obstructive pulmonary disease, advanced neoplasms and congestive heart failure)8. In England, the cancellation of non-urgent elective surgeries after 17 March 2020 was combined with the prioritization to critical care (CC) of those with high capacity to benefit as signaled by a low frailty score5,6. Even though several vaccines are now available, as the pandemic progresses, hospitals are still challenged by capacity constraints to schedule elective surgeries3,4,9.
Policies that prioritize patients with COVID-19 over other patients with greater capacity to benefit may result in immediate or delayed deaths and morbidity among both non-COVID-19 patients and patients with COVID-19 and increase the financial burden on health systems as delays in planned treatments may accelerate disease progression and the need for more costly interventions later10,11. Also, when implemented, these policies generate a backlog of non-COVID-19 patients in need of care12,13 that require prioritization rules that differ from pre-pandemic ones in order to be managed, since heterogeneity in disease progression over the postponement period might change their relative priority when compared with other patients. In England, the National Health Service (NHS) Confederation has projected waiting lists to reach 9.8 million by the end of 202114, highlighting how essential it is to identify ways to prioritize care and prevent hospitals from being overwhelmed by the various constraints posed by the pandemic15.
The optimal scheduling of hospital care has been studied albeit focusing on specific diseases16,17, types of care18,19,20 or individual hospitals21. These studies deploy methodologies that allow for a faithful modeling of the operations of an individual hospital or a specific care setting and thus enable prioritization of subgroups of patients, but they do not scale to the health system of an entire country and all patients in need of care, which is the focus of our paper.
In many health systems worldwide, hospital scheduling relies on a set of prioritization rules, which seldom look at the entirety of the health system and therefore lead to suboptimal decisions. The need to mitigate this gap in evidence has been explicitly raised during the pandemic. In England, the Nuffield Council on Bioethics issued a national statement on the need for clearer and “nationally developed and coordinated guidance on how decisions about the allocation of constrained resources should be made,”22 highlighting also the need for prioritization to consider all disease groups and their capacity to benefit. Therefore, in contrast to the approach followed in the literature, we propose a framework to optimally schedule elective care and allocate hospital capacity to elective and emergency patients in general and acute (G&A) and CC across all disease areas while accounting for national capacity and needs.
Our approach relies on a linear programming (LP) formulation of the problem, which we present in abridged form in Methods and in its entirety in the Supplementary Information.
We apply our model to the NHS in England between 2 March 2020 and 1 March 2021, with the aim of minimizing years of life lost (YLL) under alternative scenarios. These consider capacity constraints, demand for emergency care and projected COVID-19 emergencies under a set of pandemic trajectories that reflect varying stringency of control strategies, projected with a susceptible–exposed–infected–recovered (SEIR) model of SARS-CoV-2 transmission. We compare the outcomes of our optimized schedules to those of a set of simulated prioritization policies that reflect those implemented by various health systems during the pandemic, including: (1) postponement of elective procedures; (2) prioritization to CC based on frailty; and (3) re-scheduling patients to elective care using pre-pandemic prioritization rules.
Results
Model overview
Optimized schedules (OS) are formulated with a deterministic LP model that optimally schedules the admission of patients to hospital across all diseases and nationally. Specifically, the key decision variables of the model are which patients to admit and when (admission scheduling), and which patients to allocate to CC in case of capacity shortages. In the OS, emergency patients are always admitted if capacity is available, while elective admissions are scheduled weekly over the 52-week horizon with the objective of minimizing YLL. Figure 1 offers a conceptual input–output overview of the LP formulation.
When demand exceeds capacity, OS ration care. First, patients in need of elective care remain on the waiting list. Second, patients requiring CC are treated in G&A (with different transition probabilities) until CC capacity becomes available. Third, during the peaks of the pandemic and under some scenarios, emergency capacity is not sufficient to accommodate all patients. Given that we cannot rely on pre-pandemic data to estimate the implications of denied emergency care, our analysis considers two boundary cases. Our LP model conservatively assumes that all patients who are denied emergency admission die (upper-bound case). Alternatively, we use the solution of our LP model but re-evaluate the YLL under the assumption that the denied patients’ probabilities and outcomes remain as if they would be cared for in hospital (lower-bound case). By modeling both extreme cases, we fully characterize the spectrum of reasonable outcomes.
We compare our OS with standard policies (SP) that simulate prioritization policies implemented in England during the pandemic. SP manage surges in demand by postponing elective admissions. Patients with postponed procedures are kept waiting and admitted once there is capacity according to a first-in-first-out rule, that is, following an order that preserves their assigned priority when care was originally scheduled. Different disease groups are admitted in equal proportions. Under SP, CC is prioritized for non-frail patients. All other assumptions are as in OS.
We model four SP informed by elective care postponements and actual admission policies in England. SP1 and SP3 assume that elective procedures were canceled between 17 March and 29 April 2020 (weeks 3–8), as actually occurred7. SP2 and SP4 additionally trigger cancellations anytime over the intervention horizon when the projected COVID-19 incidence exceeds specific trigger points. SP1 and SP2 assume 100%, SP3 and SP4 75% elective cancellations.
The incremental cost-effectiveness ratio (ICER, a standard metric used by policymakers to decide on interventions to implement in health systems) of OS versus SP is calculated across all disease groups and ages as (where YLG is years of life gained):
Model data
To forecast non-COVID-19 elective and emergency needs (CC and G&A), we use administrative patient medical records (Hospital Episode Statistics23 (HES)) from all public acute hospitals in England between January 2015 and February 2020. To model care pathways of hospital patients with COVID-19, we use medical records from 614 patients admitted with SARS-CoV-2 infection to Imperial College Healthcare NHS hospitals between 25 February 2020 and 5 April 202024. YLL are calculated using standard life tables25. Patients are individually costed using the National Cost Collection dataset from 2015 to 201926 matched to HES data at Health Resource Group (HRG) level. Staff numbers are obtained from the NHS Electronic Staff Records for 2020. G&A beds are calculated using the March 2020 extract of the Quarterly Bed Availability and Occupancy Dataset27. CC beds are obtained from the Critical Care Monthly Situation Reports dataset for February 202028. Monthly emergency admissions are obtained from the A&E Attendances and Emergency Admissions dataset from NHS England Statistics for March to June 202029.
Scenario analysis
For both OS and SP as well as under both the upper- and the lower-bound assumptions, we account for uncertainty by considering alternative scenarios based on capacity constraints, inflows of emergency admissions and pandemic trajectories. For the capacity constraints, we consider that capacity either remains at pre-pandemic levels or is expanded to reflect hospital interventions introduced to increase total capacity by 16,500 beds (for example, field hospitals) and 38,462 staff30 (for example, recruitment of retired and student medical staff), see Supplementary Section 3.2. For the inflows of emergency admissions, we consider that emergency admissions either remain as forecasted or that there is a reduction in the forecasted number of patients in need of emergency care (to reflect behavioral changes due to the pandemic and its mitigation strategies as well as potentially increased mortality at home). We reduce our forecasted emergency needs by 34%, using A&E attendance data to estimate the proportion of the reduction in emergency admissions throughout the pandemic (Supplementary Section 3.2). For the pandemic trajectories, and to capture a range of scenarios characterizing the severity of the pandemic and the pressure on hospital capacity, we consider the following reproduction numbers (Rt) : Rt = 1.2, which yields a second-wave peak in hospital occupancy similar to the peak of the first wave; and Rt = 1.1, which generates a notably milder second wave of the pandemic, and a peak hospital occupancy akin to seasonal flu. These are combined with two alternative control interventions that impose lockdowns on 1 December 2020 (early LD) and 1 January 2021 (late LD).
The combination of these different assumptions leads to the following scenarios: baseline early-LD and late-LD scenarios assume an effective reproduction number of Rt = 1.1, pre-pandemic staff and bed capacity, and no reduction in emergencies. Best-case early- and late-LD scenarios assume Rt = 1.1 plus a 34% reduction in emergencies, and capacity expanded. Worst-case early- and late-LD scenarios assume Rt = 1.2, no reduction in emergencies and pre-pandemic capacity (Supplementary Section 8).
Scenario results
Except where flagged, YLL are reported as the minimum (which occurs in either the best-case early-LD or the best-case late-LD scenarios) and the maximum (which occurs in either the worst-case early-LD or the worst-case late-LD scenarios) observed YLL differences across all scenarios considering both the upper- and lower-bound results. Values reported refer to the comparison of OS with SP1. Comparisons with SP2–4 are in Supplementary Fig. 6.
YLL are greater under SP than OS in total (50,750–5,891,608 YLL across the scenarios) and for most disease groups and scenarios (Fig. 2). In per-capita terms and across all diseases, the largest gains accrue in the worst-case LD scenario with the OS offering up to 1.58 and 0.82 additional years of life for patients aged <25 and aged 25–64, respectively. The International Statistical Classification of Diseases and Related Health Problems (ICD) groups with the largest gains are neoplasms (C00–D48: 18,111–144,719 YLG; with the highest observed per-capita gains observed for young patients amounting to 1.1 YLG), digestive system diseases (K00–K93: 3,935–606,016 YLG; highest observed per-capita gains for young patients amounting to 1.2 YLG), and injuries and poisoning (S00–T98: 4,818–742,851 YLG; highest observed per-capita gains for young patients amounting to 1.6 YLG). Large differences are also observed for diseases of the circulatory system (I00–I99: −253–368,797 YLG; highest observed per-capita gains for young patients amounting to 2.2 YLG) and for respiratory diseases (J00–J99: −385–1,131,860 YLG; highest observed per-capita gains for young patients amounting to 1.9 YLG), see Fig. 2. OS prioritize these patients (to emergency care and to CC) over patients with COVID-19, resulting in lower YLL than SP. The highest per-capita gains of OS over SP occur in the worst-case late-LD scenario for young (<25) patients with COVID-19 (3.8 YLG) and for diseases of the circulatory system (2.2 YLG). The OS penalize the elderly patients with COVID-19 with 3.4 YLL compared with SP.
The reduction in YLL achieved by OS is greater in scenarios where capacity constraints are more restrictive (worst-case upper bound), suggesting that an optimal prioritization of electives is increasingly beneficial as resources become scarce. When capacity enables accommodating all emergencies or there is scope to invest in extra capacity for emergencies beyond existing levels (lower bound), there is less difference between OS and SP because all patients can receive care. The notable health gains of OS do not come at increased costs in baseline and best-case scenarios (Fig. 3 and Table 1). For scenarios in which OS is costlier than SP, small extra spending is required to save many lives, with the extra costs being associated with an increased number of elective admissions and shifting from low-priority to costlier high-priority patients. In the worst-case scenarios, OS are cost effective for thresholds between £57 and £1,070 per YLL (Fig. 3 and Table 1).
OS accommodate more elective admissions than SP since there is no blanket postponement of elective procedures, and admission scheduling is determined by the patient’s probability of survival and likelihood of needing emergency care while waiting for elective care. When capacity constraints are less restrictive (baseline and best-case early LD/late LD, worst-case early LD), OS exhibit fewer non-COVID-19 patients requiring emergency care than SP (Table 1). The ICDs for which the difference in elective (emergency) admissions is the highest (lowest) are digestive system diseases (K00–K93), neoplasms (C00–D48), genitourinary system disease (N00–N99), and injuries and poisoning (S00–T98) (Fig. 4 for SP1 and Supplementary Fig. 7 for SP2–4).
In the worst-case scenarios, capacity constraints are more stringent than in the baseline and best-case scenarios, and more patients require emergency care while waiting for elective care under both OS and SP. Due to resource constraints, differences in emergency admissions between OS and SP are smaller in the worst-case late-LD scenario (Table 1). Notice that there are no denials in emergency admissions in the best-case scenarios for either model. However, this is not the case in any other scenario, where emergency patients are denied admission due to capacity constraints. The number of emergency admission denials is higher for SP than OS (Table 1).
At baseline, OS deny CC to few patients (2.2% and 2.6% over the 52 weeks under early LD and late LD, respectively), most of which (83% and 85%) are patients with COVID-19 aged 65+ (Supplementary Fig. 9). In the worst-case scenarios, the share of patients denied CC increases to 4.8% and 6.4% under early LD and late LD, with 92% and 83% of patients denied CC stemming from COVID-19 aged 65+ (Fig. 5). Therefore, OS show that YLL are minimized when non-COVID-19 patients are prioritized for CC over patients with COVID-19. In the worst-case late-LD OS, patients with neoplasms and diseases of the circulatory system are prioritized over patients with COVID-19 aged 65+ (see Supplementary Fig. 10 for G&A bed utilization). In the best-case scenarios, OS admit to CC all patients requiring it, while SP deny CC to patients during the weeks following the cancellation of elective procedures (Table 1).
Discussion
In England, fears over premature mortality and morbidity associated with the postponement of elective procedures prompted policy discussions on how to accommodate patients needing elective care4. NHS England directed hospitals to sustain elective care over the winter 2020–2021. This directive is proving difficult to meet with demands from patients with COVID-19 surging31,32. Our findings are timely for the NHS as they have the potential to support re-scheduling delayed elective procedures while coping with further increases in patients with COVID-19.
While there is uncertainty about our estimates, reflected in the variability of the orders of magnitude of the reported gains of the OS over SP across the different scenarios, the scenarios and respective parameters reflect a broad spectrum of stringency of capacity constraints. As such, gains of OS over SP are characterized across plausible states of the world. Importantly, our results highlight some high-level commonalities across all scenarios that are useful to inform the formulation of national guidance akin to the guidance developed by the National Institute for Excellence Health and Care Excellence (NICE). Our findings outline life-saving key prioritization principles that can be embedded in national policies: (1) consider the relative capacity to benefit of patients with different diseases and ages; (2) postpone elective treatments for which disease progression is mild and that have lower chances of being admitted to emergency care; (3) prioritize access to emergency care and critical care based on capacity to benefit, rather than by default prioritizing one disease (for example, COVID-19) over other patients (for example, non-COVID-19)2. Where data are not available, these general principles can be used by policymakers globally to rapidly formulate national guidance in the future when national health services experience severe demand shocks such as those experienced during the pandemic. Our model is open source, and it is solved within minutes on a standard computer33,34.
While our objective is to inform the development of national guidance analogous to NICE guidance, our model does not capture variations between health authorities in provision or efficiency of implementation of the prioritization rules. However, the model can be applied to smaller geographical units that can be of use in decentralized health systems (for example, Italy or Canada). With regards to single hospitals capacity planning, while national guidance could in principle be mandatory, we would argue that—in the case of our modeling—localities should have autonomy to adapt and refine recommendations depending on the opportunities and constraints offered by local service delivery systems and population needs if, for example, they have excess capacity for certain treatments that would otherwise remain unused.
On the other hand, we do not envisage our model being used at a hospital level, as this would require detailed modeling of hospital-level capacity (by treatment) and patient-level constraints, which in turn would require discrete decision variables that would lead to much harder optimization problems.
This study has several limitations. The impact of COVID-19 on staff shortages and infection-control measures (for example, ward closures) are not modeled, which probably underestimates its impact on hospital capacity. Referrals for elective care are assumed to remain constant at pre-pandemic levels, but it is likely that primary care attendance and referrals were reduced during the pandemic35. This has two implications for the analysis: (1) we do not account for YLL due to reduced care-seeking behavior by patients; (2) the types of referral may differ because of the pandemic, thus impacting costs and YLL of policies.
Patients within each group are considered homogenous in clinical severity and disease progression, and we do not consider competing risks that could change prioritization rules. Also, patients scheduled for elective care may be heterogeneous on disease progression and probability of dying. Without an indicator of severity valid across ICDs, we account for latent severity by modeling transition probabilities as a function of waiting times and through some patients being forecasted to need care earlier than others based on historical data. Thus, this analysis may underestimate costs and YLL.
Some patients who have been scheduled for elective care may subsequently die of a hospital-acquired COVID-19 infection; this is not accounted for in our study. If granular data were available, this could be easily accounted for in the modeling of the transition probabilities. At present, the significance of hospital-acquired infections in the context of our findings is unclear. On the one hand, one could argue that due to capacity shortages, hospitals would use elective wards to admit more emergency patients with COVID-19, thus possibly increasing the risk of hospital-acquired infections. On the other hand, many patients with COVID-19 in the United Kingdom have been allocated to dedicated hospitals (for example, the Nightingale and field hospitals), and the infection-control programs and protocols developed during the pandemic have helped to reduce hospital COVID-19 transmissions36, thus mitigating the risk of nosocomial COVID-19 in hospitalized patients.
We do not consider treatment preferences of patients, the public or medical professionals37. While our model largely disregards the impact on health inequalities, it embraces the equity principle that a life year gained is of equal value, regardless of who receives it, which is akin to what is considered by NICE for priority setting. Socioeconomic deprivation and minority ethnicity is associated with increased COVID-19 mortality risk38,39. Further research is needed to ensure that OS do not inadvertently increase health inequalities, and are acceptable to clinicians, patients and the public. If data are available, the model can incorporate health inequalities.
Our objective function minimizes YLL because that is aligned with many prioritization policies routinely used by the NHS in England. For example, the underlying assumption in NICE guidance, followed also internationally, is that the limited NHS budget should be used to maximize health outcomes, albeit using quality adjusted life years (QALYs) rather than YLL. Although QALYs are used by many health systems for economic evaluations of health treatments, data are not routinely available for all diseases and health states for patients in hospital care settings. If data were available, it would in principle be possible to adapt the model to any other objective function, including the maximization of QALYs as well as any other equity considerations.
While these are important caveats that can impact the estimated YLL and costs, they are likely to affect OS and SP in a similar way, thus not impacting their comparison.
The presented model shows that changes in prioritization rules can minimize the detrimental health impact of unprecedented hospital capacity shortages during the pandemic. It operationalizes the principles of best use of limited resources underlying the management of the English health system40,41. More generally, the model is of relevance to health systems globally seeking to prioritize hospital care, substantially improving on short-sighted measures that focus on patients with COVID-19 to the detriment of the health of other patients.
Methods
Optimized schedules
OS are obtained from the solution of an LP model. In the following, we present an abridged model that contains the key decision variables, constraints and input parameters (cf. Supplementary Table 2) of the LP formulation. A more detailed description of the entire optimization model is relegated to Supplementary Section 1.1.
Our LP model can be interpreted as a multiperiod inventory model with five classes of non-negative decision variables: (1) w, tracking patients waiting for care; (2) z, recording which patients are admitted to hospital in a particular week; (3) z′, tracking the initial state of admitted patients (G&A, CC or death because of denied care due to capacity shortages); (4) y, tracking the state of hospitalized patients; (5) x, mapping the transitions of patients across states \(s,s^\prime \in {{{\mathcal{S}}}} = \left\{ {G,C,G^ \ast ,H,D} \right\}\), where G, C and G* correspond to patients in G&A, CC and G&A after having been denied CC, respectively, while H and D record the recovered and deceased patients.
The abridged LP formulation reads as follows.
The model minimizes the total YLL over a 52-week planning horizon \(t \in {{{\mathcal{T}}}} = \left\{ {1, \ldots ,52} \right\}\) across all patient groups \(p \in {{{\mathcal{P}}}}\), each of which comprises patients of a particular disease type, age group and admission type \(a \in {{{\mathcal{A}}}} = \left\{ {n,e} \right\}\), where n and e refer to elective and emergency admissions, respectively. The YLL is calculated as the specific YLL per patient group \(\lambda _p\) multiplied by the number of deaths during care (\(y_{tpa}^D\)) and due to denial of care (\(z_{tpa}^{{\prime} D}\)) (equation (1)). Other objective functions, such as the minimization of the healthcare costs, or combinations of YLL and healthcare costs, can be readily employed by modifying equation (1) accordingly.
At the beginning of week t, an exogenous inflow of patients in need of care (\(\phi _{tp}^a\)) is observed for every patient group \(p \in {{{\mathcal{P}}}}\) and admission type \(a \in {{{\mathcal{A}}}}\). Some of the elective patients are admitted to hospital during week t (\(z_{tpn}\)), while the other elective patients remain in the waiting queue or require emergency care (with probability \(\pi _{w,p}^e\)) (equation (2)). Incoming emergency patients (\(\phi _{tp}^e\)) are always admitted if capacity is available. In case of capacity shortages, admission to hospital might be denied to patients in need of emergency care; in the LP model, we conservatively assume that these patients (\(z_{tpe}^{{\prime} D}\)) die due to the lack of the required emergency treatments (equation (3)) (this assumption is relaxed in the sensitivity analysis).
When admitted to hospital, the fraction \(\pi _{z,tpa}^G\) (\(\pi _{z,tpa}^C\)) of patients requiring a G&A (CC) bed is recorded by the decision variables \(z_{tpa}^{{\prime}G}\) (\(z_{tpa}^{{\prime}C}\)) (equation (4)). The number of patients of patient group \(p \in {{{\mathcal{P}}}}\) and admission type \(a \in {{{\mathcal{A}}}}\) in a given state \(s \in {{{\mathcal{S}}}}\) at the end of week t, \(y_{t + 1,pa}^s\), is equal to the sum of the patients who remained in state s during that week (\(x_{tpa}^{ss}\)) and the transitions from other states \(s^\prime \in {{{\mathcal{S}}}}\) to s during week t (\(x_{tpa}^{s^\prime s}\)) (equation (5)). Transitions across states evolve according to the transition probabilities \(\pi _{y,pa}^{s^\prime s}\), and they are recorded by the decision variables \(x_{tpa}^{s^\prime s}\) (equation (6)). Admission to CC is an exception to this: equation (7) enforces that in case of capacity shortages, CC might be denied to a patient in need; the patient then transitions to a designated G&A state (G*) where s/he evolves according to a new set of transition probabilities until CC capacity becomes available.
The abridged model is a simplified model of the OS and does not account for the possibility to send patients requiring CC to G&A (the designated state G*) directly upon admission, it assumes that all admitted patients are hospitalized for at least one week, and it disregards the resource constraints of the health system. Those considerations are accounted for by our OS, and they are presented in Supplementary Section 1.1, which also contains more detailed explanations. The solution time of the LP model scales polynomially in the size of the input data, that is, the number of time periods, the number of patient groups as well as the number of resources. Note in particular that the solution time is unaffected by the number of patients, which makes our model particularly suited to policy advice at a national level.
Model inputs
Our model relies on weekly projections of key estimates over the 52-week horizon: patients requiring elective care (Supplementary Section 3.1), non-COVID-19 patients and patients with COVID-19 requiring emergency care (Supplementary Sections 3.1 and 4); probabilities of transitioning to various states once admitted (for example, discharged, to CC, to G&A, died) (Supplementary Section 5.2), and probabilities of patients waiting for elective care needing emergency care (Supplementary Section 5.1). Additional inputs are: forecasts of the number of patients waiting for elective care; hospitalized patients in week 0; costs and YLL of all patients (Supplementary Sections 6.1 and 6.2). Capacity constraints on the supply side are given by the maximum number of G&A/CC beds and staff (senior and junior doctors, and nurses) and recommended staff-to-bed ratios42,43. To reflect historical bed utilization rates, we assume all available capacity can be used.
Patients in need of care are categorized into groups based on primary diagnosis and age resulting in 42 and 45 disease-age groups for non-COVID-19 elective and emergency admissions, respectively, and 3 groups for COVID-19 emergency admissions (Supplementary Section 2). A binary frailty score is calculated for each patient (Supplementary Section 3)44.
Non-COVID-19 patients are forecasted using time-series methods assuming historical hospital bed utilization rates (Supplementary Section 3.1). Patients with COVID-19 are projected with a SEIR model using compartments for three age groups and degrees of severity, hospitalizations and deaths. The basic reproduction number, seed time of the epidemic, start of lockdown and reduction in transmission due to non-pharmaceutical interventions are calibrated to hospital occupancy data (Supplementary Section 4)45.
The disease progression of patients waiting for elective care may be such that they transition to the cohort requiring emergency care with a probability estimated as a function of waiting time (days) using a Kaplan–Meier estimator (Supplementary Section 5.1). Once patients are admitted for either elective or emergency care to either G&A or CC, they can transition to any of the following states: (1) remain in their current state (G&A or CC); (2) move to CC or G&A; (3) be discharged alive; or (4) die. The transition probabilities are estimated using multinomial logistic regressions, conditional on waiting time for electives (Supplementary Section 5.2).
We compute mean treatment costs for each patient group by averaging across the individual average cost of every non-COVID-19 patient in English hospitals. As the cost of patients with COVID-19 is yet unknown, we determine their HRGs in three London hospitals using the 2015–2019 national cost schedule (Supplementary Section 6.1). We calculate aggregate YLL by adding the YLL for all deaths within each age group, averaging age specific remaining life expectancies across all ages within the group (Supplementary Section 6)25.
Ethical approval
The study was approved by the NHS Digital Independent Group Advising on the Release of Data (IGARD) committee (application reference: DARS-NIC-276970). No other ethical approval was required.
Data availability
Due to the data sharing agreement, data cannot be made publicly available with this study. Individual-level data used for the analysis of COVID-19 outcomes are not publicly available and only accessible to selected researchers. Data from Hospital Episode Statistics and Electronic Staff Records are available to all researchers upon making an application to NHS Digital. Applications can be made through the NHS’s Data Access Request Service (https://digital.nhs.uk/services/data-access-request-service-dars). A&E Attendances and Emergency Admissions data are freely available at https://www.england.nhs.uk/statistics/statistical-work-areas/ae-waiting-times-and-activity/. Bed Availability and Occupancy data are freely available at https://www.england.nhs.uk/statistics/statistical-work-areas/bed-availability-and-occupancy/. Critical Care Bed Capacity data are freely available at https://www.england.nhs.uk/statistics/statistical-work-areas/critical-care-capacity/. NHS Cost Collection data are freely available at
https://www.england.nhs.uk/national-cost-collection/. All model input data are documented in the Supplementary Information and made available open source at https://github.com/HFAnalyticsLab/overflow_analysis. Source Data for Figs. 2–5 and for Supplementary Figs. 3–10 are provided with this paper.
References
Birnbaum, M. & Morris, L. With coronavirus exploding in Europe, hospitals calculate how long until they hit capacity. The Washington Post (2 November 2020); https://www.washingtonpost.com/world/europe/covid-coronavirus-europe-hospitals-capacity/2020/10/31/fe074b56-188a-11eb-8bda-814ca56e138b_story.html
Emanuel, E. J. et al. Fair allocation of scarce medical resources in the time of COVID-19. N. Engl. J. Med. 382, 2049–2055 (2020).
Coronavirus: Weston General Hospital halts admissions. BBC News (25 May 2020); https://www.bbc.co.uk/news/uk-england-somerset-52796589
Cullen, L. Covid-19: Elective surgery cancelled at Craigavon Hospital. BBC News NI (4 November 2020); https://www.bbc.co.uk/news/uk-northern-ireland-54801326
COVID-19 Rapid Guideline: Critical Care in Adults NICE Guideline No. 159 (National Institute for Health and Care Excellence, 12 February 2021); https://www.ncbi.nlm.nih.gov/books/NBK566886/
Important and Urgent—Next Steps on NHS Response to COVID-19 (NHS Improvement and NHS England, 2020); https://www.england.nhs.uk/coronavirus/wp-content/uploads/sites/52/2020/03/urgent-next-steps-on-nhs-response-to-covid-19-letter-simon-stevens.pdf
Operating Framework for Urgent and Planned Services in Hospital Settings During COVID-19 (NHS England, 2020); https://covidlawlab.org/wp-content/uploads/2020/06/Operating-framework-for-urgent-and-planned-services-within-hospitals.pdf
Raccomandazioni di Etica Clinica per l’Ammissione a Trattamenti Intensivi e per la Loro Sospensione 207-2011 (Società italiana di anestesia analgesia rianimazione e terapia intensiva, 2020); https://www.quotidianosanita.it/allegati/allegato2675063.pdf
Cullen, L. Coronavirus: three of NI’s five health trusts cancel planned surgery. BBC News NI (12 November 2020); https://www.bbc.co.uk/news/uk-northern-ireland-54921936
Griffin, S. COVID-19: waiting times in England reach record highs. Br. Med. J. 370, m3557 (2020).
Rosenbaum, L. The untold toll—the pandemic’s effects on patients without COVID-19. N. Engl. J. Med. 382, 2368–2371 (2020).
Campbell, D. NHS hospital waiting lists could hit 10 million in England this year. The Guardian (2020). at <https://www.theguardian.com/society/2020/jun/10/nhs-hospital-waiting-lists-could-hit-10-million-in-england-this-year>
Davey, M. Covid-19 Victoria: hospitals advised to reduce elective surgery as doctors warn of overcrowding. The Guardian (9 July 2020); https://www.theguardian.com/australia-news/2020/jul/09/covid-19-victoria-hospitals-advised-to-reduce-elective-surgery-as-doctors-warn-of-overcrowding
Road to recovery. NHS Confederation https://www.nhsconfed.org/news/2020/06/road-to-recovery (2020).
Hunter, D. Trying to “Protect the NHS” in the United Kingdom. N. Engl. J. Med. 383, e136 (2020).
Déry, J. et al. A systematic review of patient prioritization tools in non-emergency healthcare services. Syst. Rev. 9, 227 (2020).
Valente, R. et al. How to prioritise and optimise access to elective surgery in the COVID-19 pandemic: an interhospital alliance pathway based on the SWALIS 2020 model. A single-cohort feasibility and pilot study. SSRN https://doi.org/10.2139/ssrn.3622349 (2020).
Cayirli, T. & Veral, E. Outpatient scheduling in health care: a review of literature. Prod. Oper. Manag. 12, 519–549 (2003).
Ahmadi-Javid, A., Jalali, Z. & Klassen, K. J. Outpatient appointment systems in healthcare: a review of optimization studies. Eur. J. Oper. Res. 258, 3–34 (2017).
Cardoen, B., Demeulemeester, E. & Beliën, J. Operating room planning and scheduling: a literature review. Eur. J. Oper. Res. 201, 921–932 (2010).
Gupta, D. & Denton, B. Appointment scheduling in health care: challenges and opportunities. IIE Trans. 40, 800–819 (2008).
Statement: The Need for National Guidance on Resource Allocation Decisions in the COVID-19 Pandemic (The Nuffield Council on Bioethics, 2021); https://www.nuffieldbioethics.org/news/statement-the-need-for-national-guidance-on-resource-allocation-decisions-in-the-covid-19-pandemic
Hospital Episode Statistics (HES). NHS Digital https://digital.nhs.uk/data-and-information/data-tools-and-services/data-services/hospital-episode-statistics (2020).
Perez-Guzman, P. N. et al. Clinical characteristics and predictors of outcomes of hospitalized patients with coronavirus disease 2019 in a multiethnic London National Health Service Trust: a retrospective cohort study. Clin. Infect. Dis. https://doi.org/10.1093/cid/ciaa1091 (2020).
Past and Projected Period and Cohort Life Tables, 2018-Based, UK (Office for National Statistics, 2020); https://www.ons.gov.uk/peoplepopulationandcommunity/birthsdeathsandmarriages/lifeexpectancies/bulletins/pastandprojecteddatafromtheperiodandcohortlifetables/1981to2068
NHS England: national cost collection for the NHS. NHS Digital https://www.england.nhs.uk/national-cost-collection/ (2020).
Statistics: bed availability and occupancy. NHS Digital https://www.england.nhs.uk/statistics/statistical-work-areas/bed-availability-and-occupancy/ (2020).
Statistics: critical care bed capacity and urgent operations cancelled. NHS Digital https://www.england.nhs.uk/statistics/statistical-work-areas/critical-care-capacity/ (2020).
Statistics: A&E attendances and emergency admissions. NHS Digital https://www.england.nhs.uk/statistics/statistical-work-areas/ae-waiting-times-and-activity/ (2020).
McCabe, R. et al. Adapting hospital capacity to meet changing demands during the COVID-19 pandemic. BMC Med. 18, 329 (2020).
Elective Letter (NHS England, 20 August 2020); https://www.england.nhs.uk/publication/elective-letter/
Campbell, D. & Quinn, B. Mass cancellations of NHS operations inevitable this winter, say doctors. The Guardian (23 October 2020); https://www.theguardian.com/society/2020/oct/23/mass-cancellations-of-nhs-operations-inevitable-this-winter-say-doctors
jdaeth274 jdaeth274/overflow_analysis: publication release. Zenodo https://doi.org/10.5281/zenodo.5079551 (2021).
Moret, S. ImperialCollegeLondon/OptimalScheduling4COVID: v1.0. Zenodo https://doi.org/10.5281/ZENODO.5071220 (2021).
Watt, T., Firth, Z., Fisher, R., Thorlby, R. & Kelly, E. Use of Primary Care During the COVID-19 Pandemic (The Health Foundation, 2020); https://www.health.org.uk/news-and-comment/charts-and-infographics/use-of-primary-care-during-the-covid-19-pandemic
Rhee, C. et al. Incidence of nosocomial COVID-19 in patients hospitalized at a large US academic medical center. JAMA Netw. Open 3, e2020498 (2020).
COVID-19 Documents (The Federation of Surgical Specialty Associations, 2020); https://fssa.org.uk/covid-19_documents.aspx
Deaths Involving COVID-19 by Local Area and Socioeconomic Deprivation: Deaths Occurring Between 1 March and 31 July 2020 (Office for National Statistics, 2020); https://www.ons.gov.uk/peoplepopulationandcommunity/birthsdeathsandmarriages/deaths/bulletins/deathsinvolvingcovid19bylocalareasanddeprivation/deathsoccurringbetween1marchand31july2020
Updating Ethnic Contrasts in Deaths Involving the Coronavirus (COVID-19), England and Wales: Deaths Occurring 2 March to 28 July 2020 (Office for National Statistics, 2020); https://www.ons.gov.uk/peoplepopulationandcommunity/birthsdeathsandmarriages/deaths/articles/updatingethniccontrastsindeathsinvolvingthecoronaviruscovid19englandandwales/deathsoccurring2marchto28july2020
Next steps on the NHS five year forward view. Br. Med. J. https://doi.org/10.1136/bmj.j1678 (2017).
Our Principles (NICE, 2020); https://www.nice.org.uk/about/who-we-are/our-principles
Safe Medical Staffing. (Royal College of Physicians, 2020); https://www.rcplondon.ac.uk/projects/outputs/safe-medical-staffing
Staffing Levels: Advice Guides (Royal College of Nursing, 2020); https://www.rcn.org.uk/get-help/rcn-advice/staffing-levels
Soong, J. et al. Quantifying the prevalence of frailty in English hospitals. BMJ Open 5, 8456 (2015).
Statistics: COVID-19 hospital activity. NHS Digital https://www.england.nhs.uk/statistics/statistical-work-areas/covid-19-hospital-activity/ (2020)
D’Aeth, J. et al. Optimal scheduling for COVID [source code]. Code Ocean https://doi.org/10.24433/CO.0290795.v1 (2021).
Acknowledgements
We thank the Data Management Team at the Health Foundation for their work to prepare the data extract and manage information governance. This work uses data provided by patients and collected by the NHS as part of their care and support. We also thank C. Donnelly, S. van Esland and S. Riley for comments and suggestions. G.F. acknowledges funding from Jan Wallanders and Tom Hedelius Foundation and the Tore Browaldh Foundation, funding from MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1), jointly funded by the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth and Development Office (FCDO), under the MRC/FCDO Concordat agreement, part of the EDCTP2 program supported by the European Union; and acknowledges funding by Community Jameel. K.H. and N.F. were partially funded by the National Institute for Health Research (NIHR) Health Protection Research Unit in Modelling and Health Economics, a partnership between Public Health England, Imperial College London and London School of Hygiene and Tropical Medicine (grant code NIHR200908); and also acknowledge funding from the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1) and by the Community Jameel. P.N.P.-G. acknowledges funding from the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1; and by the Community Jameel. D.H. acknowledges part funding from the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1). J.C.D. acknowledges funding from the Wellcome Trust (reference 102169/Z/13/Z). S.M. acknowledges partial support from the Swiss National Science Foundation (SNSF) under grant number P2ELP2_188028. S.G. was funded by the Imperial College President’s PhD Scholarship. F.G. and S.R.D. were funded by the Health Foundation as part of core staff member activity. This research was funded in whole, or in part, by the Wellcome Trust (grant number 102169/Z/13/Z). For the purpose of open access, the author (J.C.D.) has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission.
Author information
Authors and Affiliations
Contributions
All authors contributed to all stages of the work and to the original draft, review and editing. M.M. coordinated the different components of the work. J.C.D., S.G., F.G., D.H., E.K., K.L., S.M. and D.R. (listed in alphabetical order) contributed equally as first authors of this work and carried out the analyses. S.G., E.K. and S.M. were responsible for the modeling, optimization and simulation developments. F.G., D.H., D.R., K.L. and P.N.P.-G. were the data curators. M.M., G.F. and W.W. are senior authors and contributed equally. M.M. is the guarantor of the paper and accepts full responsibility for the overall content. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.
Corresponding author
Ethics declarations
Competing interests
All authors have completed the Unified Competing Interest form (available on request from the corresponding author). K.L. is employed by Flatiron Health but employment started after the research finished. The other authors declare no competing interests.
Additional information
Peer review information Nature Computational Science thanks the anonymous reviewers for their contribution to the peer review of this work. Handling editor: Ananya Rastogi, in collaboration with the Nature Computational Science team.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Supplementary Information
Supplementary Sections 1–10, Figs. 1–10 and Tables 1–11.
Supplementary Data 1
Source data used for Supplementary Fig. 3.
Supplementary Data 2
Source data used for Supplementary Fig. 4.
Supplementary Data 3
Source data used for Supplementary Fig. 5.
Supplementary Data 4
Source data used for Supplementary Fig. 6.
Supplementary Data 5
Source data used for Supplementary Fig. 7.
Supplementary Data 6
Source data used for Supplementary Fig. 8.
Supplementary Data 7
Source data used for Supplementary Fig. 9.
Supplementary Data 8
Source data used for Supplementary Fig. 10.
Source data
Source Data Fig. 2
Statistical source data for Fig. 2.
Source Data Fig. 3
Statistical source data for Fig. 3.
Source Data Fig. 4
Statistical source data for Fig. 4.
Source Data Fig. 5
Modeling results, reported in eight different sheets corresponding to the eight different subplots in Fig. 5.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
D’Aeth, J.C., Ghosal, S., Grimm, F. et al. Optimal national prioritization policies for hospital care during the SARS-CoV-2 pandemic. Nat Comput Sci 1, 521–531 (2021). https://doi.org/10.1038/s43588-021-00111-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s43588-021-00111-1
This article is cited by
-
Disruption of hospital care during the first year of the COVID-19 pandemic impacted socioeconomic groups differently: population based study using routine registration data
BMC Health Services Research (2024)
-
Were the unfinished nursing care occurrence, reasons, and consequences different between COVID-19 and non-COVID-19 patients? A systematic review
BMC Nursing (2023)