Using Administrative Data to Incorporate Age and Sex-Dependent Resource Use for COVID-19 Acute Care Resource Use Simulations in Ontario, Canada

Evolving COVID-19 Data in Ontario, Canada

CCM Plus is Ontario’s province-wide population-based dataset on all individuals who test positive for COVID-19. (3) CCM Plus includes individual-level data on demographics (e.g., age, sex, region), epidemiology (e.g., likely acquisition), patient characteristics (e.g., co-morbidities), acute care resource utilization (e.g., hospitalization, intensive care unit (ICU) admission, invasive mechanical ventilation (IMV)), health outcomes (e.g., mortality), and long-term care (LTC) residency. We obtained research ethics board approval from the University of Toronto.

We conducted a population-based cohort study using public health data to describe age- and sex-specific acute care use, LOS, and mortality.(4) This study showed that acute care use, including LOS, and mortality varies by age and sex, and has decreased between March and September in Ontario overall and within each groups. We include these changes and temporal trends in the updated model to capture unmeasured effects (e.g., changes in clinical practice).(4)

Demographics

The modelled population (i.e., incident cases starting March 1, 2020) takes into consideration individuals’ age and sex, informed by Ontario reported COVID-19 case demographics. In CORE 2.0, demographic parameter values (e.g., age and sex) of infected cases are sampled from distributions informed by 4-month averages between March-June, and July-October (Table 1). Day-to-day demographics fluctuate similarly to daily reported cases. We chose 4-months because after June, there was a noticeable increase in the proportion of reported COVID-19 cases in the younger age groups: 0-9, 10-19, 20-29 and 30-39 years. Using this data, an individual’s age and sex now determines the individual’s probability of being hospitalized, requiring ICU admission, mechanical ventilation, and death. A 4-month time frame balanced using specific data and computational resources required for simulation.

Time from Episode Date to Hospitalization

In CORE 2.0, we incorporated a lag period of 8 days between symptom onset and hospitalization for all cases. This lag period was informed by the mean number of days from episode date (defined as the earliest of case created, case reported, specimen, and symptom onset dates) to hospital admission date observed in Ontario.(5) We assumed that all COVID-19 individuals do not interact with the healthcare system in this period of 8 days, and do not consider variability across individuals.

Acute Care Resource Use

The probability of hospitalization, ICU admission, and ICU admissions requiring mechanical ventilation were all updated using age and sex-stratified parameter values (Table 2). For ICU admissions, there was a decrease in the proportion of ICU admissions in September within all age groups, and we therefore used the mean proportion of ICU admissions for September and October for anyone requiring ICU admission post-August.

Mortality

The probability of death for patients in each of the acute care levels were updated by age, sex, and time (Table 3). Similar to the case demographics, there was a notable change in mortality, overall, and within age groups, in May and June. Hence, we used the mean mortality from March to May for the first three months of the pandemic, and the mean mortality from June to October for incident cases from June onwards. To parameterize the model, which is in daily time-steps, mortality was converted over the duration of the average LOS in each respective level of acute care (ward, ICU without IMV, and ICU with IMV).

Length of Stay

We updated LOS (days) for each acute care level, which also changed over time by month (Table 4). When the individual is hospitalized, LOS is sampled from poisson and gamma distributions informed by monthly mean LOS and standard deviation. For patients who are admitted to ward, the average LOS was between 5.7 and 16.15 days. For patients requiring ICU care without ventilation, their average stay in the ICU was between 6.5 to 10.7 days. For patients who needed ICU care and ventilation, their average stay in the ICU was between 14.8 and 21.6 days. For any individual requiring ICU care (ventilator or no ventilator), the time spent in the ward prior and post ICU admission was 6.83 days; we assumed the 6.83 days was equally divided for the ward LOS pre- and post-ICU. Individuals requiring ventilators spend on average 1.2 days in the ICU without ventilation. Overall, for an individual requiring ventilation, their LOS includes the following distinct components: time spent in the ward prior and post to ICU admission, time in the ICU not ventilated, plus time on ventilation. Anytime during this LOS, patients may die. When an individual no longer requires a resource, either due to completion of LOS and discharge or death, the resource will become available to other individuals. Individuals who cannot access their required resource form a queue and may die depending on their disease severity.

Analysis and Calibration

The model starts simulations on March 1, 2020. The number of COVID-19 cases between March 1, 2020 and October 31, 2020 were taken from Ontario’s Ministry of Health data.(6) We calibrated some LOS components (ICU, non-ventilated) to reported ward and ICU (with and without ventilation) occupancy between June 15 and October 31, 2020.(7) The non-ventilated ICU LOS was calibrated only for the most recent month of October since the analysis is subject to some censoring (i.e., individuals who may not have an outcome at the end of analysis follow-up time). The simulated occupancy from CORE 2.0 is compared to the reported occupancy data for all levels of care in Figures 2a-2c. Overall, simulations for all three level of care resources compare well with reported occupancy data up until October 31, 2020. A brief summary of all parameters is in Table 5.

• Download figure
• Open in new tab

Figure 2a. Predicted and confirmed ward occupancy in Ontario based on administrative data

• Download figure
• Open in new tab

Figure 2b. Predicted and confirmed ICU occupancy in Ontario based on administrative data

• Download figure
• Open in new tab

Figure 2c. Predicted and confirmed ventilator occupancy in Ontario based on administrative data

• Download figure
• Open in new tab

Figure 3. Predicted resource use in Waterloo region between Dec 1 and 30, 2020 for A) Ward beds B) ICU beds C) Ventilators

Regional Application of CORE 2.0 in Ontario

Here we describe an example of the use of CORE 2.0 for a regional analysis of acute care resource use for the Region of Waterloo, Ontario, Canada that serves 617,870 people from the cities of Cambridge, Kitchener, Waterloo and townships of North Dumfries, Wellesley, Wilmot and Woolwich. In this example, we simulate the ward bed, ICU bed, and ventilator occupancies for 30 days starting on December 1, 2020 for three different case trajectory scenarios:

1. Growth 1.01: Growth of incident cases for 10 days at 1%, followed by contraction (reduction) of incident cases by 5% (i.e., [incident cases] _day = [incident cases] _day-1*0.95, where day = current date) for the remainder of the 30 days.

2. Growth 1.03: Growth of incident cases for 10 days at 3%, followed by contraction (reduction) of incident cases by 5% for the remainder of the 30 days.

3. Growth 1.05: growth of incident cases for 10 days at 5%, followed by contraction (reduction) of incident cases by 5% for the remainder of the 30 days.

These scenarios were chosen based on input from stakeholders from the region, and are designed to capture the potential decrease in cases after recent implementation of public health interventions to limit community transmission. The contraction rate of 5% starting at 10 days after peak are user-defined assumptions. The model used a run-in time of 8-days prior to the simulation start date (i.e., November 23, 2020), allowing the reported cases in the past 8 days to also start being hospitalized (if needed) on December 1^st, and onwards, to provide a more accurate number of occupancies required within the next 30 days. The reported cases for November 23 to 30, 2020 were obtained from the Region of Waterloo.(8) On November 30, 2020, there were 66 incident cases reported, which is what the case growth scenarios started with; all reported cases used as input is reported in Table 6.

The results reported are the ward bed, ICU bed, and ventilator occupancies in Figures 3A-3C. These are the number of people requiring this resource on specific dates, and not the new admissions (i.e., there are new admissions, and also discharges based on survival, and length of stay). In these scenarios, as the increase of incident cases grows faster, the need for all resources increases. In the most conservative scenario used to estimate the capacity required, the Region of Waterloo may require up to 29 ward beds, 10 ICU beds, and 8 ventilators for COVID-19 patients, and in the worst case growth estimate (5%), these numbers increase to 37, 13, and 10, respectively.

There are some limitations to this regional analysis. First, the model parameters have all been calibrated and validated to Ontario’s population-level use of COVID-19 resource use, length of stay, and mortality by age and sex. This does not necessarily reflect the population from the Region of Waterloo. However, the Region of Waterloo is not considered a rural area or geographically different from many of the hotspots in Ontario (e.g., Toronto, Peel Region) and so the demographics of COVID-19 cases should still apply. Second, the results represent the occupancy that would be seen in the entire Region of Waterloo region, and not necessarily at a specific hospital. These results do not consider the current ward, ICU and ventilator occupancies prior to December 1, 2020. Unfortunately, because these individuals’ trajectory in their LOS is unknown, it was not possible to add these current occupancies into the model and so these must be considered for the near-term future in addition to the new resources predicted. Despite these caveats, CORE 2.0 was able to simulate acute care resource needs for regional analysis and provide an estimate for local hospital planners for resource allocation purposes.