Long-term healthcare resource use and cost associated with COVID-19 disease from a health system perspective. An equity-focused population-based cohort study
1Institute of Health Policy, Management, and Evaluation, University of Toronto, Toronto, ON, Canada
2Health Systems and Policy Research Collaborative Centre, University Health Network, Toronto, ON, Canada
3ICES, Toronto, ON, Canada
4Dalla Lana School of Public Health, University of Toronto, Toronto, ON, Canada
5MAP Centre for Urban Health Solutions, Li Ka Shing Knowledge Institute, Unity Health Toronto, Toronto, ON, Canada
6Department of Medicine and Institute of Medical Science, Temerty Faculty of Medicine, University of Toronto, Toronto, ON, Canada
7Public Health Ontario, Toronto, ON, Canada
*Corresponding Author: Adeteju Ogunbameru, Health Systems and Policy Research Collaborative Centre, University Health Network, Toronto General Hospital, Eaton Building, 10th Floor, 200 Elizabeth Street, Toronto, ON M5G 2C4, Email: adeteju.ogunbameru@mail.utoronto.caAbstract
SARS-CoV-2 strained Ontario’s health system, with social determinants of health (SDH) underexplored in cost analyses. We examined COVID-19 attributable healthcare resource use and costs from the Ontario health system perspective using health administrative data. We conducted a cohort study, matching 162,633 SARS-CoV-2-exposed individuals 1:1 to unexposed individuals. We calculated 10-day per-person mean attributable costs (2023 CAD) across care phases (pre-diagnosis, acute, post-acute, terminal), stratified by individual and area-level SDH. Among exposed individuals (mean age 40.4 years, 50.7% female), 6% were hospitalized, 1.3% admitted to critical care, and 2% died within 360 days. Mean (SD) person acute phase cost was $244 ($235-$253) and higher among males, recent immigrants, individuals living in low-income neighbourhoods and neighbourhoods with a higher proportion of crowded households. Extrapolating to the population level of 166,801 exposed individuals, the mean total survival-adjusted 360-day cost was $436 million. COVID-19 increased healthcare costs, disproportionately burdening marginalized communities.
Article notes
Competing Interest Statement
The authors have declared that no competing interests exist
Funding Statement
Yes
Introduction
From January 2020 to March 2024, 1,720,202 Ontario residents were reported to have been infected with SARS-CoV-2 (1). The COVID-19 pandemic significantly impacted the health of Ontarians, prompting the provincial government to invest over $10 billion in healthcare initiatives between 2020 and 2023, including funding for population and public health programs (e.g., testing and vaccination), personal protective equipment, additional hospital and long-term care beds, and support for frontline workers (2,3)
Preexisting social determinants of health (SDH) were amplified during the COVID-19 pandemic. In Ontario, residents of low-income and materially deprived neighbourhoods experienced higher rates of SARS-CoV-2 infection and mortality, while individuals who were recent immigrants or employed in essential jobs faced greater exposure and barriers to care (4). SDH have been shown to contribute to disparities in healthcare cost, e.g., an Ontario study found that individuals with stroke living in the most materially deprived neighbourhoods had significantly higher healthcare costs than those in the least deprived areas (11). underscoring how SDH contribute to disparities in health outcomes and resource use.
Stratifying healthcare costs by SDH provides insight into how healthcare resources are utilized across populations, which is critical for equitable health system planning.This approach is particularly important in the context of COVID-19, where marginalized populations have experienced higher infection rates, severe disease, and differential access to care, potentially leading to variations in healthcare utilization and associated costs (5).
Previous studies from England and Ontario have demonstrated that COVID-19-related healthcare costs vary by individual-level characteristics such as age, income level and co-morbidities (6–8). However, these studies lacked comprehensive stratification by SDH, limiting the understanding of cost disparities among marginalized groups. For example, Yang et al. (2023) reported substantial costs among people with severe COVID-19, particularly for older adults and those with comorbidities, but did not examine SDH (6). Tsui et al. (2022) compared initial healthcare costs among marginalized populations in British Columbia and Ontario, highlighting variation driven by differences in provincial healthcare systems, with acute care hospitalizations as the primary cost driver, yet long-term costs were not captured (7). Similarly, Sander et al. (2024) quantified short- and long-term healthcare costs in Ontario, identifying acute care hospitalizations as the main contributor, but provided limited data on SDH (9). These studies underscore the need for a deeper understanding of how SDH influence COVID-19 costs to inform healthcare delivery.
Our objective is to provide a comprehensive population-based analysis of COVID-19-attributable healthcare resource use and cost, stratified by SDH, from the Ontario health system perspective.
Materials and methods
Our study methods followed the Reporting of Studies Conducted using Observational Routinely Collected Data guidelines (10).
Study Design and Setting
We conducted a population-based, retrospective cohort study to estimate short- and long-term costs attributable to COVID-19 from the health system perspective (Ontario Ministry of Health) using an incidence-based costing approach. All data were sourced from administrative datasets housed at ICES, an independent, non-profit research institute. ICES links datasets using encoded identifiers, enabling health system evaluations without requiring patient consent under Ontario privacy laws (11). The ICES repository includes health service records of approximately 14 million Ontario residents (12).
Population
The study included Ontario residents of all ages from January 1, 2016, to December 31, 2020, eligible for the Ontario Health Insurance Plan and with valid ICES identifiers. We excluded individuals with missing demographic data (e.g., sex, postal code, birth date), those aged >110 years, individuals with nosocomial COVID-19 infection, and long-term care residents, due to their distinct healthcare utilization patterns (13).
The exposed cohort comprised individuals with a positive polymerase chain reaction (PCR) test for SARS-CoV-2 in the COVID-19 Integrated Testing Data (C19INTGR) dataset between January 1 and December 31, 2020. The index date was defined as the date of the first positive PCR test.
The unexposed cohort consisted of Ontario residents identified from January 1, 2016, to December 31, 2018, with no positive COVID-19 test reported in C19INTGR dataset in 2020. We used a historical cohort to mitigate potential biases due to contamination (people who may have had COVID-19 but were not tested) and changes in health service utilization during 2020 (e.g., surgery cancellations) (14). Unexposed individuals were selected by randomly sampling 50% of the Ontario population from the Registered Persons Database (RPD) and were randomly assigned an index date based on the exposed cohort’s index date distribution.
Both cohorts were followed for one year post-index or until death.
Cohort selection is shown in Figure 1.
Covariates and Data Sources
We linked demographic, clinical, and social determinants of health (SDH) data from C19INTGR and other ICES datasets using unique encoded identifiers (15). Baseline characteristics included demographic and clinical factors (age, sex, comorbidity, rurality) and social determinants of health (SDH), defined at the individual level (e.g., immigration status) and area level (e.g., neighbourhood income, household crowding, proportion of essential workers).
Comorbidity was measured as the Resource Utilization Band based on the Johns Hopkins Aggregated Diagnosis Groups two years prior to index date (16). Age, sex, and rurality were sourced from RPD. Immigration status was obtained from the Immigration, Refugees, and Citizenship Canada Permanent Resident Database in ICES.
For area-level SDH, we considered variables previously shown to be associated with increased odds of SARS-CoV-2 infection (17) or higher healthcare costs in Ontario (18,19). We included neighbourhood income quintile, % visible minority living in an area, household crowding index (number of persons per room), and % working in essential services in an area. Area-level SDH data were sourced from the 2016 Canadian Census at the Dissemination Area (DA) level, a small, relatively stable geographic unit with an average population of 400–700 persons. DA-level data were assigned to individuals based on residential postal codes (20).
All baseline characteristics were assessed at the index date. Details and definitions for each variable are included in Appendix Table A1.
Healthcare utilization and cost data were sourced from multiple datasets (15) including the Canadian Institute for Health Information’s Discharge Abstract Database (acute inpatient hospitalization), Ontario Mental Health Reporting System (mental health), Continuing Care Reporting System (complex continuing care), National Rehabilitation System (inpatient rehabilitation), National Ambulatory Care Reporting System (same day surgeries), Ontario Health Insurance Plan (physician services), Ontario Drug Benefit (prescription drugs), and Office of the Registrar General Death registry.
Data was assessed between March 1, 2022 and December 31, 2024. Authors had access to information that could identify individual participants during or after data collection.
Outcomes
Our primary outcome were COVID-19-attributable costs by phase of care, standardized to mean 10-day per-person costs in 2023 Canadian dollars. Addtional outcomes included all-cause mortality and resource use (e.g., hospitalizations, intensive care unit (ICU) admissions) among the exposed cohort.
Matched Cohort
We matched each exposed individual to one unexposed individual using hard- and propensity score matching.
We calculated propensity scores using greedy matching with a caliper width of 0.2 standard deviations (SD) (21). Covariates included in the propensity score model were individual-level (rurality, immigration status) and area-level variables (marginalization indices, housing, occupation, and visible minority status) measured at index date. Hard matching was performed on sex, same age, index date (±60 days), and logit of the propensity score. Balance between groups was assessed using weighted standardized differences with <0.1 indicating negligible imbalance (21). Variables included in propensity score model are described in Appendix Table A1.
For stratified analyses, the variable of interest was excluded from the propensity score model and included in the hard match alongside sex, age, index date (± 60 days) and logit of the propensity score.
To estimate terminal phase costs, we matched deceased exposed individuals to unexposed individuals of the same age or older to better account for predictors of cost prior to death. Other variables included in this match were sex, logit of the propensity score with a caliper width of 0.2 SD (calculated at index), and death date ±30 days. Terminal phase costs were stratified by early (≤60 days post-index) and late deaths (>60 days post-index).
Analysis
Healthcare resource use (hospitalizations, ICU admissions, mortality) is reported for exposed individuals prior to matching. Continuous variables are reported as mean (SD) and median (IQR), categorical variables as frequencies and percentages.
Post-matching, costs, standardized to mean 10-day per-person cost, were calculated and plotted over 360 days, pre- and post-index date. Four care phases were defined using joinpoint analysis, guided by clinical expertise (SM and BS) and the WHO’s post-COVID-19 condition framework (22): pre-diagnosis (20 days pre-index), acute (80 days post-index), post-acute (days 81–360), and terminal (60 days prior to death date) (Appendix Figure A1).
Phase-specific attributable costs were calculated as the mean cost differences (with 95% confidence intervals (95% CI)) between matched exposed and unexposed individuals using a generalized estimating equation model.
Costs were stratified by individual- and area-level SDH (age, sex, immigration status, household crowding index, %visible minority living in an area, and % essential workers living in an area).
Phase-specific costs were combined with crude survival data to estimate survival-adjusted mean (minimum and maximum) per-person costs at 30, 90, 180, and 360 days (23), and stratified by sex.
All analyses were conducted using SAS Enterprise Guide 8.3 (SAS Institute, Inc.).
Ethics approval
This study was approved by the University Health Network Research Ethics Board. ICES is a prescribed entity under Ontario’s Personal Health Information Protection Act (PHIPA). Section 45 of the PHIPA authorizes the ICES to collect personal health information, without consent, for the purpose of analysis or compiling statistical information with respect to the management of, evaluation or monitoring of the allocation of resources to, or planning for all or part of the health system. Projects that use data collected by ICES under section 45 of PHIPA and use no other data are exempt from the REB review. The use of the data in this project is authorized under Section 45 and approved by the ICES Privacy and Legal Office.
Results
A cohort of 166,801 individuals with confirmed COVID-19 met the study inclusion criteria (Figure 1).
Their mean (SD) age was 40.38 (19.70) years, 50.6% were female. Most (96%) lived in urban areas, 21.4% were high healthcare resource users, 46.2% resided in low-income neighbourhoods (quintiles 1& 2), 8.3% were recent immigrants (migrated to Canada within five years), and 42.1% lived in crowded dwellings (2 or more persons per room). The unexposed cohort of 6,659,844 individuals had a mean (SD) age of 40.5 (22.7) years, 50.6% were female. Of these, 89.3% lived in urban areas, 17.1% were high healthcare resource users, 39.1% resided in low-income neighbourhoods (quintiles 1&2), 3.4% were recent immigrants, and 22.6% lived in crowded dwellings (2 or more persons per room).
Baseline characteristics are presented in Table 1.
Healthcare Utilization
Within the 360-day follow-up, 5.9% of individuals with confirmed COVID-19 were hospitalized, 1.3% required ICU admission, and 2% died (all-cause mortality).
Among individuals admitted to ICU, 49.9% were aged >65 years, 65.4% were male, 89.6% were within the moderate to very high healthcare resource utilization band, and 57.2% resided in low-income neighbourhoods (quintiles 1& 2).
Among individuals who died, 48.8% were aged >80 years, 55.1% were male, 51.1% resided in low-income neighbourhoods (quintiles 1&2), and 63.3% died within 60 days of confirmed COVID-19 (early deaths).
The median (IQR) hospital length of stay for late deaths was 239 days (55, 627).
Of the 166,801 exposed individuals who met the study inclusion criteria, 97.5% were matched 1:1 to unexposed individuals with standardized differences < 0.1, indicating good balance (Table 1).
Healthcare Costs
Full cohort Analysis
The mean (95% CI) total COVID-19 attributable 10-day per-person phase-specific cost for the pre-diagnosis, acute, post-acute and terminal care phase (combined early and late deaths) were $11 ($6-$15), $244 ($235-$253), $20 ($17-$24), and $4,753 ($4,419-$5,088), respectively (Table 2)
Cost components per phase are shown in Figure 2.
In the pre-diagnosis care phase, hospitalizations, home care services and emergency department visits were the main contributors to the mean phase cost, each accounting for approximately 27% of the mean total phase cost at $3 ($0-$6), $3 ($2-$4), and $3 ($2-$3), respectively.
In the acute care phase, hospitalizations were the largest cost component at $185 ($177-$192), accounting for 76% of the mean total cost. Emergency department visits contributed 6% to the mean total phase cost, and rehabilitation contributed 2%.
In the post-acute phase, hospitalizations contributed 40% to the mean phase costs at $8 ($6-$11), followed by outpatient visits (15%) at $3 ($2-$4), and home care services (10%) at $2 ($2-$3).
In the terminal phase, mean (95% CI) attributable costs were higher for early deaths ($5,518 [$5,099–$5,937]) than for late deaths ($3,549 [$3,001–$4,097]). For early deaths, hospitalizations accounted for 100% of the mean total cost ($5,575 [$5,189–$5,962]). For late deaths, hospitalizations accounted for 88.9% of the mean total cost at $3,152 ($2,642–$3,662), with home care costs ($6 [−$41 to $54]) and complex care costs ($14 [−$87 to $117]) being other major cost categories. (Table 2).
Stratified Analysis
Males had higher costs in the acute ($293 [278-$308]) and terminal phases ($5,690 [$5,206–$6,174]) compared to females (acute: $197 [$186- $207], terminal: ($3,569 [$3,129–$4,008]) (Table 3).
Individuals aged ≥70 years had higher costs compared to other age groups in the pre-diagnosis, acute care, post-acute care and terminal phases (Table 4).
Recent immigrants had lower costs in the pre-diagnosis, acute care, and post-acute care phases compared to long-term residents: pre-diagnosis $0 (-$9- $11) vs $11 ($6-$16), acute care $118 ($98- $137) vs $254 ($244-$263), and post-acute care phase $18 ($12-$24) vs $22 ($18-$26) (Table 5).
Individuals living in the lowest-income neighbourhoods (quintile 1) had the highest mean cost in the pre-diagnosis ($135 [$126, $144]), acute ($405 [$384, $425]), and post-acute ($132 [$125, $138]) phases compared to other quintiles, while those living in the highest-income neighbourhoods (quintile 5) had the lowest costs (pre-diagnosis: $96 ($88-$105), acute: $285 ($263-$306), post-acute: $98 ($92-$105) (Table 6).
Individuals living in neighbourhoods with crowded housing (2 or more persons per room) had higher costs in the acute phase ($248 [$233-$263]) compared to those living in non-crowded housing neighbourhoods ($239 [$227- $251]) (Table 7)
No consistent directional trend was observed in phase-specific mean attributable costs by proportion of visible minorities or essential workers in a neighbourhood (Appendix Table A2 & A3)
The mean (minimum, maximum) total survival-adjusted 360-day cost per person was $3,297 ($3,096, $3,501), with males having higher costs ($3,618 [$3,316, $3,918]) than females ($2,727 [$2,478, $2,974]). Extrapolating to the population of 166,801 exposed individuals, the total survival-adjusted 360-day cost was estimated at $436 million ($408 –$462 million) (Table 8).
Discussion
Our study estimated healthcare costs attributable to confirmed COVID-19 in Ontario, Canada, over a 360-day period post-index, stratified by SDH.
Costs followed a typical U-shaped pattern (24): high during the acute phase, declining during the post-acute phase, and rising again during terminal care. Acute phase costs were the highest, consistent with the intensive resource demands of moderate to severe disease requiring hospitalization, ICU admission, and in some cases, mechanical ventilation (25).
Terminal phase costs for individuals with confirmed COVID-19 who died within 60 days of index (early deaths) were substantially higher than for those who died after 60 days (late deaths). Higher costs for early deaths likely reflect the use of resource-intensive interventions (e.g., mechanical ventilation and intensive care) during the acute phase of moderate to severe COVID-19 (26). In contrast, individuals with late deaths (>60 days) had longer hospital stays (median: 239 days [IQR: 55–627]), with hospitalization costs approximately 56.5% of those for early deaths, but incurred higher attributable home care and complex care costs, reflecting a shift to less costly care settings.
Compared to previous Ontario studies, our attributable mean cost estimate for the acute phase $244 was substantially lower than the $ 4,491 reported by Tsui et al (7) but similar to the data reported by Sander et al. $240 (9). The discrepancies between our findings and Tsui et al. could be due to differences in the costing window and the study time period. Tsui et al. estimated short-term COVID-19 healthcare costs using a 0–120-day costing window around the index date, based on cases from the first pandemic wave (January–June 2020). In contrast, both our study and Sander et al. used a much longer 360-day follow-up period with detailed 10-day costing intervals, capturing costs across multiple pandemic waves and enabling the inclusion of long-COVID and longer-term health-system impacts (7,9). Similar to Sander et al.(9), we observed sustained healthcare costs up to 360 days post-diagnosis, indicating a significant long-term burden of COVID-19. This aligns with evidence of persistent health issues associated with COVID-19, commonly referred to as post-COVID-19 condition, which lasts for more than 12 weeks after the initial infection (27). However, unlike Sander et al., which focused on broader cost implications without stratifying by social determinants (9), our stratified analysis provides novel insights into how SDH were associated with variations in healthcare costs attributable to COVID-19. We found that sex, age, immigration status, income, and household crowding were associated with variation in phase-specific attributable costs. These findings align with Canadian epidemiological studies reporting higher rates of severe disease and mortality among these groups during the COVID-19 pandemic (13).
Multiple factors may explain the higher per-person healthcare costs observed among residents of low-income neighbourhoods. Individuals living in these areas often experience structural barriers to timely and effective healthcare, including longer travel times to healthcare facilities, fewer available healthcare providers, and limited access to preventive services (28). These barriers can potentially lead to delayed care with disease progressing to moderate or severe phase at the time of presentation, increasing reliance on acute-phase services such as emergency departments, hospitalizations, and ICU services (7,8). From our findings, hospitalization was the largest contributor to acute-phase costs among marginalized populations, including individuals living in low-income neighbourhoods. Consistent with our findings, a recent meta-analysis including over two million individuals reported higher odds of hospitalization among socioeconomically disadvantaged populations compared with advantaged groups during the COVID-19 pandemic (29).
Our study has some limitations. First, the C19INTGR dataset likely underestimates COVID-19 infections in Ontario due to incomplete testing, particularly during the early stages of the pandemic when testing was restricted to symptomatic individuals meeting specific criteria (30). As a result, the dataset potentially over-represents individuals with moderate to severe COVID-19 infections, who were more likely to be tested and incur higher healthcare costs, while those with asymptomatic or mild infections are underrepresented. Supporting this, a study conducted in Toronto, Ontario estimated that for each laboratory-confirmed symptomatic case between March and May 2020, there were approximately 18 undetected community infections (31), most of which were likely asymptomatic and required minimal to no healthcare services. Thus, our estimated attributable COVID-19 costs may be an overestimate. It is important to note that testing patterns varied by socioeconomic status: early in the pandemic, individuals in high-income neighbourhoods were more likely to be tested (4); however, testing later increased in low-income neighbourhoods due to higher infection rates and targeted public health efforts (4). This shift may have influenced our findings, as the initially higher testing rates in high-income neighbourhoods could have underestimated cost differences, while the subsequent increase in testing in low-income neighbourhoods could have contributed to greater healthcare use and costs observed in these populations. Second, our findings reflect the first pandemic year (2020) with limited preventive or treatment options and higher hospitalization rates, potentially overestimating costs compared to later waves characterized by vaccine and treatment availability, and milder variants. Furthermore, the area-level SDH variables (e.g., income, household crowding index, percentage of essential workers) were measured at the neighbourhood level using data at the DA level. While DA-level data serve as useful proxies in the absence of individual-level SDH data, they may not capture within-area and individual-level variability and could introduce bias. Therefore, our findings reflect area-level associations and should not be interpreted as individual-level effects. Lastly, while positive pre-diagnosis costs may reflect residual confounding bias, healthcare visits often preceded testing due to constrained testing capacity in 2020, aligning with early-pandemic healthcare-seeking patterns.
Our study has several strengths. We provide comprehensive estimates of the health and healthcare burden of confirmed COVID-19, stratified by key area-level SDH. Additionally, we matched exposed individuals to a historical unexposed cohort to account for healthcare delivery changes (e.g., surgery cancellations) in Ontario from May 2020 to February 2022, minimizing overestimation of COVID-19-attributable costs (14). We used a phase-of-care approach to mitigate the effect of censoring during the study follow-up period (32). Lastly, our phases of care periods and their associated costs can be mapped as health states in decision-analytic models to inform economic evaluations of COVID-19 interventions.
Our study highlights the substantial healthcare costs associated with COVID-19, particularly during acute and terminal phases, as well as the ongoing costs up to 360 days post-infection, indicating long-term sequelae. Our stratified analysis further supports the need for structural policy interventions that address the social and environmental conditions contributing to higher healthcare costs in marginalized communities. For example, enhancing healthcare infrastructure in low-income and crowded housing neighborhoods may reduce delays in care and prevent disease progression. Incorporating social risk indicators (e.g., income, housing density) into resource allocation frameworks can also help ensure that future health system planning and emergency responses are more equitable.
Conclusion
Our study highlights significant healthcare costs attributable to confirmed COVID-19 in Ontario, concentrated in the acute and terminal phases and sustained up to 360 days post-infection. Costs varied by social determinants, with higher costs observed among marginalized populations, including those residing in low-income neighbourhoods, crowded housing, recent immigrants, and older adults. Our findings highlight systemic inequities in disease burden and resource utilization, emphasizing the need for early and tailored intervention strategies (e.g., vaccination, access to testing). Our findings can inform policies addressing existing disparities, optimizing resource distribution, and strengthening health system preparedness for future public health crises.
Acknowledgements
This document used data adapted from the Statistics Canada 2016 Census Area Profiles and the Statistics Canada Postal CodeOM Conversion File, which is based on data licensed from Canada Post Corporation, and/or data adapted from the Ontario Ministry of Health Postal Code Conversion File, which contains data copied under license from ©Canada Post Corporation and Statistics Canada.
Parts of this material are based on data and/or information compiled and provided by Statistics Canada, the Canadian Institute for Health Information (CIHI), Ontario Health (OH), and the Ontario Ministry of Health, as well as Immigration, Refugees and Citizenship Canada (IRCC), current to September 2020. The analyses, conclusions, opinions, and statements expressed herein are solely those of the authors and do not reflect those of the funding or data sources, no endorsement is intended or should be inferred.
We thank the Toronto Community Health Profiles Partnership for providing access to the Ontario Marginalization Index. The authors thank Hong Lu for her coding contributions during the analysis and Patricia Aluko for proofreading the study manuscript.
Data availability
The dataset from this study is held securely in a coded form at ICES. While legal data-sharing agreements between ICES and data providers (e.g., healthcare organizations and government) prohibit ICES from making the dataset publicly available, access may be granted to those who meet pre-specified criteria for confidential access, available at www.ices.on.ca/DAS (email: das@ices.on.ca). The full dataset creation plan and underlying analytic code are available from the authors upon request, with the understanding that computer programs may rely upon coding templates or macros that are unique to ICES and are therefore either inaccessible or may require modification.
Conflict of Interest
None