Introduction
Since the legalisation of cannabis in Canada in 2018, many studies have been undertaken to study the effect of legalisation on cannabis-related emergency department (ED) visits and hospitalisations in the general population (Auger et al., Reference Auger, Luu, Ayoub, Bilodeau-Bertrand, Lo and Low2021; Kim et al., Reference Kim, Chum, Nielsen, Allin, Penney, Rittenbach and O'Campo2022, Reference Kim, Chum, Nielsen, MacMaster, Rittenbach, Allin and Hamilton2023; Myran et al., Reference Myran, Pugliese, Tanuseputro, Cantor, Rhodes and Taljaard2022b; Yeung, Weaver, Janz, Haines-Saah, & Lang, Reference Yeung, Weaver, Janz, Haines-Saah and Lang2020). Given that past-month cannabis use is approximately twice as prevalent among adults with depression compared to those without depression (Gorfinkel, Stohl, & Hasin, Reference Gorfinkel, Stohl and Hasin2020), and heavy-use can exacerbate depressive symptoms (Degenhardt, Hall, & Lynskey, Reference Degenhardt, Hall and Lynskey2003), Canadian cannabis legalisation may lead to more adverse events in individuals with depression. This study aims to (1) ascertain the rates of cannabis-related ED-visits in depressed and non-depressed individuals in Ontario, Canada, and (2) fill the gap in the literature regarding the effect of cannabis legalisation on ED-visits in individuals who were diagnosed with depression.
Prior research shows that those with depression in the past year had more severe symptoms of marijuana use disorder based on the Diagnostic and Statistical Manual-5 (compared to those without depression), namely, they used marijuana in larger amounts or over a longer period than intended (OR 1.8, 95% CI 1.56–2.19, p < 0.001), repeatedly failed to cut down on the amount of marijuana use (OR 1.9, 95% CI 1.50–2.30, p < 0.001), spent an inordinate amount of time on the acquisition, use and recovery of marijuana (OR 1.6, 95% CI 1.41–1.93), and continued to use marijuana despite adverse social consequences (OR 3.2, 95% CI 2.41–4.30) (Dierker, Selya, Lanza, Li, & Rose, Reference Dierker, Selya, Lanza, Li and Rose2018). In an analysis of the US National Survey on Drug Use and Health annual surveys (2005–17), a higher prevalence of cannabis use in the past month among people with depression (19%) v. those without (9%) was reported (Pacek, Weinberger, Zhu, & Goodwin, Reference Pacek, Weinberger, Zhu and Goodwin2020). Another study found that depressed mood scores (using the Centre for Epidemiological Studies-Depression scale) were higher for adolescents who reported cannabis use, compared to those who did not (Hernandez et al., Reference Hernandez, Cancilliere, Graves, Chun, Lewander and Spirito2016). These studies suggest that depression may modify cannabis use such that it elevates its risk of potential adverse events relative to the general population. Increased use of recreational cannabis has been found to be associated with higher risk of cardiovascular adverse events (Singh et al., Reference Singh, Saluja, Kumar, Agrawal, Thind, Nanda and Shirani2018). There is also evidence that individuals with depression are at higher risk for cannabis-related ED-visits. Using state-wide data from Colorado from 2012 to 14, a study found that individuals with mood disorders (which includes depression and bipolar disorders) have 7.4 times higher prevalence of ED-visits for cannabis use than individuals without (Hall et al., Reference Hall, Monte, Chang, Fox, Brevik, Vigil and James2018). There is also limited evidence that more liberal cannabis laws may promote cannabis use among depressed individuals: in a study of US adults with mood or anxiety disorders, those living in states with greater access to cannabis through medical licensing laws had 21% higher prevalence of cannabis self-medication compared to those in states without these laws (Sarvet et al., Reference Sarvet, Wall, Keyes, Olfson, Cerdá and Hasin2018).
To establish the unique contribution of our study to the literature, a comprehensive literature search using MEDLINE and Scopus was conducted (1) to identify prior studies that quantified cannabis-related ED-visits in individuals with depression, and (2) the impact of cannabis legalisation on cannabis-related ED-visits among those with depression. Cannabis-related ED visits is defined as any ED visit for mental and behavioural disorder and/or poisoning due to cannabis use. See Fig. 1 for PRISMA flowchart and search terms. Only two studies met our search criteria. In a Colorado study of a single hospital system that included 4202 cannabis-related ED and urgent care visits among adolescents, a comorbid psychiatric diagnosis was made in 71% of the sample, with the highest diagnoses for depression (39%), mood disorder (22%), anxiety disorder (13%) and bipolar disorder (6%) (Wang, Davies, Halmo, Sass, & Mistry, Reference Wang, Davies, Halmo, Sass and Mistry2018). This study provides evidence that depression is likely the most common psychiatric condition among adolescents with cannabis-related ED-visits. A single-centre study involving retrospective chart review of cannabis-related ED-visits found that patients with depression accounted for 4.3% of visits 2020 (Shelton et al., Reference Shelton, Mills, Saben, Devivo, Williamson, Abbott and Monte2020). Both studies were limited by (1) convenience sampling, (2) samples drawn from a single-centre/hospital and (3) the lack of population-wide state-level data, which may increase the risk of selection bias and reduce generalisability. Furthermore, no studies investigated the potential differential impact of cannabis legalisation on depressed v. non-depressed individuals.
PRISMA flowchart and search terms.
Scopus
TITLE-ABS-KEY (“depress*” OR “dysthymia” OR “mood disorder”) AND TITLE-ABS-KEY (“cannabis” OR “marijuana”) AND ALL (“emergency service” OR “emergency department” OR “emergency room”)
Pubmed.
((depress*[Title/Abstract]) OR (dysthymia[Title/Abstract]) OR (mood?disorder[Title/Abstract])) AND ((cannabis[Title/Abstract]) OR (marijuana[Title/Abstract])) AND ((emergency?service) OR (emergency?department) OR (emergency?room))
Based on the results of our systematic search, this is the first study to characterise the risk of cannabis-related ED-visits in depressed and non-depressed populations using population-wide data, and is also the first study to compare the impact of cannabis legalisation on the cannabis-related ED-visits in depressed and non-depressed individuals. In Ontario, legalisation was multi-phased, with the legalisation of flower and herbs (phase 1) in October 2018, and the legalisation of cannabis edibles one year later (phase 2). Given the prior literature which indicates the increased vulnerability of depressed individuals regarding cannabis use, we hypothesise that (1) cannabis-related ED-visits are greater in depressed v. non-depressed individuals, and (2) the effects of recreational cannabis legalisation on cannabis-related ED-visits is greater among those with depression. Since prior studies have shown gender differences in the rates and severity of cannabis use disorders (Dierker et al., Reference Dierker, Selya, Lanza, Li and Rose2018), we tested these hypotheses for men and women separately accounting for potential differences in the effects of legalisation across men and women.
Methods
We performed a comparative interrupted time-series (CITS) analysis to examine whether the effect of Ontario's cannabis legalisation on cannabis-related ED-visits differed between depressed and non-depressed individuals.
Data source and study population
Using health administrative data accessed through ICES (formerly the Institute for Clinical Evaluative Sciences), we created a cohort that included the entire adult population (18+) eligible for Ontario Health Insurance Plan, OHIP (n = 11 156 100) on 17 October 2015, and they were followed through their healthcare records until 20 May 2021. OHIP is Ontario's universal healthcare that covers over 95% of Ontario residents, those who are not eligible includes individuals in their 3-month waiting periods and migrants with temporary status (e.g. international students, temporary workers) (Ontario Ministry of Health, 2022), which includes approximately 500 000 people (OHIP For All, 2022). Participants must be 18 or over as of 17 October 2018 (the date of cannabis legalisation), and have continuous OHIP coverage and residency in Ontario for the entire study period to be included in the study (October 2015 to May 2021).
This study complied with privacy regulations of ICES. To protect privacy, all cell sizes (i.e. the number of individuals in a categorical variable) with fewer than six individuals in all descriptive tables were suppressed and reported as n < 6. Consent was not obtained for participants for the use of their data in this study. ICES is an independent, non-profit research institute whose legal status under Ontario's health information privacy law allows it to collect and analyse healthcare and demographic data used for the study, without consent, for the purposes of health system evaluation and improvement. All patient information was anonymised and de-identified prior to analysis. Ethics approval for this study was obtained through York University (REB# 2022-254).
Population at risk
We identified individuals diagnosed with depression by a physician [OHIP diagnostic code 311 (Government of Ontario, 2015)] in a look back period up to 5 years prior to legalisation between 17 October 2013 to 16 October 2018 (n = 929 844). This look back period was established based on consultations with two clinical psychiatrists, given that depression diagnoses that are older than 5 years are likely to indicate an illness that has been resolved, while diagnoses within 5 years may indicate unresolved illness. A prior study has shown that the code has high specificity and positive predictive value for the measure of depressive disorder (i.e. 99.69% specificity and 89.62% for positive predictive value) (Alaghehbandan, MacDonald, Barrett, Collins, & Chen, Reference Alaghehbandan, MacDonald, Barrett, Collins and Chen2012).
Exposure: cannabis-related policy in Ontario
To examine the impact of recreational cannabis policy on cannabis-related ED-visits among individuals diagnosed with depression (hypothesis #2), we used three distinct time periods related to the Ontario legal cannabis availability: (1) the pre-legalisation period (October 2015–October 2018), (2) phase 1 of legalisation (October 2018–February 2020), when flower and herb were available online through a government website and at limited private physical retail locations, and (3) phase 2 (March 2020 onwards), which marked the rapid private retail expansion (removal of provincial retailer cap) alongside increased edible cannabis availability. These time periods have been used in prior studies on cannabis-related ED-visits following cannabis legalisation and commercialisation in the Ontario population (Kim et al., Reference Kim, Chum, Nielsen, Allin, Penney, Rittenbach and O'Campo2022). Cannabis-related ED-visits between March and April 2020 were censored to account for significant decreases in healthcare utilisation at the onset of the COVID-19 pandemic, as suggested by prior literature (Myran et al., Reference Myran, Cantor, Finkelstein, Pugliese, Guttmann, Jesseman and Tanuseputro2022a).
Outcome
The outcome of the study was monthly cannabis-related ED-visits from October 2015 to May 2021 in Ontario (the frequency of the event by gender for each month of the study period). We define cannabis-related ED visits as any ED visit for mental and behavioural disorders and poisonings due to cannabis use. Monthly cannabis-related ED-visits were collected from the National Ambulatory Care Reporting System (NACRS) using at least one of the following ICD-10 codes for either the primary or supplemental diagnosis: F12 (including 11 subcodes for mental and behavioural disorders related to cannabis) and T40.7 (cannabis poisoning) (Canadian Institute for Health Information, 2021a). Positive cases were identified using the primary and non-primary diagnosis codes in NACRS. This method of extracting cannabis-related acute care events conforms with the definition set out for the Canadian Institute for Health Information (Canadian Institute for Health Information, 2021b), which has been used in prior studies (Callaghan et al., Reference Callaghan, Sanches, Vander Heiden, Asbridge, Stockwell, Macdonald and Kish2021; Counsil of State and Territorial Epidemiologists, n.d.; Myran et al., Reference Myran, Pugliese, Tanuseputro, Cantor, Rhodes and Taljaard2022b; Yeung et al., Reference Yeung, Weaver, Janz, Haines-Saah and Lang2020; Yeung, Weaver, Hartmann, Haines-Saah, & Lang, Reference Yeung, Weaver, Hartmann, Haines-Saah and Lang2021), and shows a high positive predictive value (>95%) (DeYoung et al., Reference DeYoung, Chen, Beum, Askenazi, Zimmerman and Davidson2017).
Statistical analysis
To compare the depressed population with non-depressed individuals (hypothesis #1), we conducted a matching method to identify an appropriate counterfactual, i.e. a sample of the population that has not been diagnosed with depression that is matched on key variables to individuals diagnosed with depression. We used 1-to-1 nearest-neighbour propensity score matching on (1) age group, (2) sex, (3) neighbourhood socioeconomic status [using the Ontario Marginalisation Index (Public Health Ontario, 2022)], (4) residency in Northern or Southern Ontario and (5) seven Aggregated Diagnosis Groups (ADGs) clusters, produced by the Johns-Hopkins Adjusted Clinical Group® System (version 10), which are used as our comorbidity indicators (John Hopkins University, 2022). Northern v. southern indicator was included because prior research has found that individuals living in Northern Ontario have higher rates of cannabis-related ED-visits (Kim et al., Reference Kim, Chum, Nielsen, Allin, Penney, Rittenbach and O'Campo2022).
After matching, we produced age- and sex-specific rates for cannabis-related ED-visits separated for the depressed and non-depressed groups in matched and unmatched samples. For comparison of the outcome between depressed and non-depressed individuals, the exact Poisson method was used to calculate incidence rate ratios (IRR) with 95% confidence intervals (CI), which presents effect sizes. These comparisons were performed separately for men and women, and on the matched and unmatched samples.
To investigate the potential differential impact of legalisation across the depressed and non-depressed groups (hypothesis #2), we conducted a CITS analysis for cannabis-related ED-visits between depressed and non-depressed individuals using the user-written package for ITSA (Linden, Reference Linden2015):
where Yt is the monthly count of cannabis-related ED events at time-point t. Z is the dummy variable denoting cohort assignment status [depressed (=1) v. not depressed (=0)], while X is the dummy variable denoting cannabis legalisation phases, with X 1t for phase 1 of legalisation and X 2t for phase 2. Tt is the time since the start of the study at time point t and XtTt is an interaction term. Detailed descriptions of the dummy variables X 1t and X 12, and continuous variables Tt, T 1t and T 2T are presented in online Supplementary Tables S2 and S3. Betas 1, 2, 3, 8, 9 refer to parameters for the control group (i.e. non-depressed), and betas 5, 6, 7, 10, 11 are used to modify the former coefficients to estimate the parameters of interest for the depressed group. For example, change in slope in the pre-legalisation period for the non-depressed group is represented by β 1, while β 1 and β 5 are summed to represent pre-legalisation slope for the depressed group. β 4 in the above represents the difference in the level (intercept) of the outcome variable between target and comparison prior to the intervention; however, it is not used in the parameter estimation because it is not of interest in our study. See details in the online Supplementary Table S1 for the full explanations for every parameter that are used in our analyses. Four gender-stratified models were specified, where model 1 is the males-only model using the unmatched sample, model 2 is the female unmatched model, model 3 is the matched male model and model 4 is the matched female model. STATA16 (College Station, TX) was used for our data analysis. Each test (per every postestimation analysis) was two-tailed with a 5% as alpha to test statistical significance. To correct for multiple testing, we estimated the probability of making at least one type 1 error to be 53%. The overall level of significance was calculated as 0.0033 using the Bonferroni correction. We present both the matched and unmatched results since matched results compare depressed with a non-depressed group similar to the depressed group, while the unmatched results compare the depressed with the non-depressed general population. While the matched results highlight the effects of independent effects of depression, the latter highlights average (unadjusted) differences that are relevant for healthcare planning.
Sensitivity analyses
There may be a difference in healthcare utilisation patterns between individuals whose ED-visit lists cannabis as either the primary diagnosis, or as one of the secondary diagnoses. To ensure the consistency of our results, we repeated our analysis where the outcome is restricted to ED-visits with cannabis listed as the primary diagnosis (i.e. removing cases with cannabis as secondary diagnosis only).
Results
The most common ICD-10 code used for cannabis-related ED-visits was F12.1 (mental/behavioural disorder from harmful use), which was found in 52% of those with depression and 48% of those without depression. The second and third most common codes, F12.0 (mental/behavioural disorder from acute intoxication) and T40.7 (poisoning) respectively, were more prevalent among those without depression than those with depression (see Table 1 for code distribution in the matched and unmatched samples) (Table 1).
Percent of total cannabis-related ED-visits by ICD-10-CA codes, depressed v. non-depressed, in matched and unmatched sample
Note: percentage add up to more than 100% because some hospitalisation records had more than one cannabis-related diagnostic code. 2.07% of hospitalisations included 2+ cannabis-related diagnoses.
Table 2 presents the sample characteristics for unmatched and matched samples at baseline (Table 2). Association between each exposure and depression status is obtained through logistic regression. The strength and direction of the association is summarised by Somer's D, where a higher absolute value represents a stronger association. It should be noted that the matching procedure reduced the association between every covariate and depression status (as indicated by a lower Somer's D), which implies an improved covariate balance between the target and comparison groups in the matched sample. For example, in the unmatched sample, the proportion of women in the target (depressed) and comparison (non-depressed) group was 62% and 50% respectively, as opposed to 62% and 62% in the matched sample. For further information on the differences between the matched and unmatched groups in terms of the proportion of those with cannabis-related ED visits across depression and substance use disorder statuses, please see the online Supplementary Tables S4 and S5.
Descriptive characteristics of unmatched and matched study samples at baseline
a Association between each target and comparison status (i.e. depressed v. non-depressed) is obtained through logistic regression. The strength and direction of logistic regression is summarised by Somer's D, where a higher absolute value represents a stronger association.
Rates of cannabis-related ED-visits by gender and depression status in each phase of cannabis-related policy are presented in Table 3 (matched and unmatched). Depressed individuals have consistently higher rates of cannabis-related ED-visits for both matched and unmatched samples, but the difference in the rates between depressed and non-depressed individuals is greater in the unmatched sample (7–8 times higher), compared to the matched sample (4 times greater) as indicated by the IRR value found in Table 3.
Rates of cannabis-related ED-visits per 100 000 person-years in depressed and non-depressed individuals (both matched and unmatched sample)
IRR, incidence rate ratio.
While the initial disparities between depressed and non-depressed individuals are large, they appear to reduce over the course of the study period. In the matched sample in Table 3, depressed women were 4.81 (95% CI 3.02–8.00) times more likely to have cannabis-related ED-visits compared to non-depressed women during the pre-legalisation period. This rate fell to 3.80 during phase 1, and fell yet again in phase 2 to 3.13 (95% CI 2.22–4.51). This pattern of falling IRRs is repeated among men, with depressed men having cannabis-related ED-visits at a rate 4.63 times higher than non-depressed men during the pre-legalisation period, 4.20 during phase 1 and 4.03 in phase 2. A similar pattern of decreasing disparities is also seen in the unmatched sample (Table 3).
With regards to the CITS analyses, in the matched sample (Table 4), trends in ED-visits among depressed individuals are compared to those who were as similar as possible to them without depression. While there was modest month-to-month increase in the number of ED-visits in the pre-legalisation period in both depressed and non-depressed men and women (ranging from 1% to 2% per month), neither legalisation nor the introduction of edibles led to further sudden or gradual changes (as shown in Table 4). There was also no difference in the size or direction of change observed across the target and comparison groups.
Comparative interrupted time-series analyses, estimations based on a multiple group, multiple intervention design for both matched and unmatched samples
In the unmatched sample, trends in ED-visits among depressed individuals are compared to the Ontario general population (see the online Supplementary Table S1 for the complete description of the following coefficients). A similar upward pre-legalisation trend was observed in the target (β 1 + β 5) and comparison (β 1) group; however, in both men and women, the upward trajectory (β 5) was stronger in the non-depressed group, i.e. an additional 0.5% increase was seen in non-depressed men, and an additional 0.7% in non-depressed women, relative to their depressed counterparts. The immediate change (β 4) associated with phase 1 (Table 4) appeared to have divergent impact in the target and comparison groups. For men, there was an immediate 20.6% reduction (95% CI 5.4–35.8% reduction) in the depressed group, while non-depressed men saw a 23.1% increase. For women, phase 1 was not associated with any change in ED-visits in the depressed group, while the non-depressed group saw a 26.1% immediate increase (95% CI 11.1–41.1%). The difference tests (β 4) in the immediate change found that there were statistically significant differences between non-depressed and depressed individuals for both men and women.
The results of the supplementary analysis where the outcome was ED-visits with cannabis listed as the primary outcome in online Supplementary Table S6. The results show that removing cases where cannabis was secondary diagnosis only, the sensitivity test showed substantially similar findings compared to our main analyses.
Discussion
We found that individuals diagnosed with depression have consistently higher rates of cannabis-related ED-visits compared to those without for both matched and unmatched samples: approximately four times higher risk based on the matched sample (11 v. 50 events/100 000 person-years in men, 5 v. 20 events/100 000 person-years in women), which provides support for hypothesis #1. Phase 1 and 2 legalisation was not associated with significant changes in ED-visits in the matched sample; however, immediate increase was associated with phase 1 in the unmatched sample in the non-depressed group (i.e. the general population). This immediate increase (23% in men, 26% in women) in the general population is consistent with prior literature on the subject. For example, Myran et al., also based in Ontario, found an immediate increase of 12% in rates of cannabis-attributable ED-visits in the general population (Myran et al., Reference Myran, Pugliese, Tanuseputro, Cantor, Rhodes and Taljaard2022b). The difference across the studies may be attributable to different age cut-offs used, and the fact that individuals diagnosed with depression were not included in our comparison group.
On the other hand, contrary to our hypothesis #2, depressed individuals did not experience any increase in ED-visits associated with legalisation. This is the first study to examine cannabis-related acute care rates and disparities in depressed v. non-depressed individuals. Despite the prior evidence that depressed individuals are more vulnerable to cannabis-related disorders and problematic use, our study shows that this group did not experience increased risk of cannabis-related ED-visits due to the recreational cannabis legalisation. Prior research shows that a proportion of people with depression may be using cannabis for self-medication purposes (i.e. to help cope with symptoms of depression) (Wallis et al., Reference Wallis, Coatsworth, Mennis, Riggs, Zaharakis, Russell and Hale2022), and many are regular users of cannabis even prior to legalisation, and thus the legalisation for recreational use may not have a significant effect on their cannabis-related ED outcome. Further research should investigate whether self-medicating with cannabis among depressed individuals is a major risk factor for cannabis-related ED visits.
The discrepancy between the matched and the unmatched sample may be interpreted as adjusted v. unadjusted results, this is because the matching procedure balances the distribution of sociodemographic and health characteristics across the target and comparison groups, and subsequently we can estimate, ceteris paribus, the independent effects of legalisation by depression status. On the other hand, comparison to the general population (shown in our unmatched sample results), which gives the average risk of cannabis-related ED-visits in the depressed v. the average risk in the general population, can also be informative since it can better represent the disparities between the two groups regardless of their sociodemographic and health information.
The limitations of this study are as follows: (1) the depressed group that we identified would not include undiagnosed patients with depression, who may be hidden in the general population. Since they have not been clinically assessed, their symptoms may be worse or better than those diagnosed, which may bias the results towards an unknown direction. (2) Although our use of CITS may have reduced the effects of unobservable concurrent interventions, there is still possibility that the unobservable confounders have differential effects for the target and comparison groups. (3) Although we had the censoring window between March and April 2020 to account for the effects of the COVID-19 pandemic on healthcare utilisation, this could bias the results towards an unknown direction if there were differential effects between depressed and non-depressed individuals. (4) There is a possibility that our matching method may introduce dependency among the data which is a known issue with propensity score-based matching methods. However, given that unmatched and matched samples produced similar results, our use of matching is likely robust.
Conclusion
Individuals diagnosed with depression have a higher rate of cannabis-related ED-visits; however, legalisation did not appear to be associated with further elevated risk in ED-visits for such individuals. The general elevated risk of cannabis misuse in this group is still a concern, and healthcare providers may consider screening for cannabis use in individuals diagnosed with depression, and provide adequate resources and education to reduce the risk of cannabis-related acute care. Furthermore, while the legalisation of herbs and flowers (phase 1) and edibles (phase 2) was not associated with greater increases of cannabis-related acute care in the depressed v. general population, there is a need to continue monitoring this clinical group given the expansion of the cannabis market in Canada, which have led to the introduction of new edible products such as cannabis-infused coffee, beers, spirits, chocolates and desserts and new methods of delivery (e.g. cannabis patches, vapes, suppositories) that are now being marketed with greater choice than ever before. In addition, future studies may take into account the impacts of the market regulation for cannabis on depressed and non-depressed individuals. For example, while this study only included Ontario, where the online and wholesale distribution of cannabis products is centralised through a crown corporation monopoly (i.e. Ontario Cannabis Store), future studies could compare with other jurisdictions that have a decentralised approach (e.g. Alberta).
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291723000569.
Data
The data used in this study are third-party data that were provided by ICES (formerly known as the Institute for Clinical Evaluative Sciences), and were accessed remotely through their secure Data & Analytics Services. While data sharing agreements prohibit ICES from making the dataset publicly available, access may be granted to those who meet pre-specified criteria for confidential access, please contact a representative at das@ices.on.ca to request access and verification. The full dataset creation plan and underlying analytics code are available from the authors upon request, understanding that the computer programmes may rely upon coding templates or macros that are unique to ICES and are therefore either inaccessible or may require modification.
Acknowledgements
This study contracted ICES Data & Analytic Services (DAS) and used de-identified data from the ICES Data Repository, which is managed by ICES with support from its funders and partners: Canada's Strategy for Patient-Oriented Research (SPOR), the Ontario SPOR Support Unit, the Canadian Institutes of Health Research and the Government of Ontario. The opinions, results and conclusions reported are those of the authors. No endorsement by ICES or any of its funders or partners is intended or should be inferred.
Author contributions
A. C., C. K. and G. J. D. contributed to the conception and design of the study. C. K. and G. J. D. performed the statistical analysis. C. K., G. J. D., A. C. and A. N. wrote the first draft of the manuscript. All authors contributed to manuscript revision, read and approved the submitted version.
Financial support
Funding for the study was provided by Canadian Institutes of Health Research (Project Grant FRN# 173447, NPI: Antony Chum). The funding agencies had no role in the design and conduct of the study; the collection, management, analysis or interpretation of the data; the preparation, review or approval of the manuscript; or the decision to submit the manuscript for publication. The project principal investigator (Antony Chum) is supported by the Canada Research Chair Program (CRC-2021-00269).
Conflict of interest
None.
Ethical standards
Ethics approval for this study was obtained through York University (REB# 2022-254).
Patient consent for publication
Not applicable.
Permission to reproduce material from other sources
Not applicable as no materials were reproduced from other sources.
Open access
