Article Text

Download PDFPDF

Original research
Relationship between life satisfaction and preventable hospitalisations: a population-based cohort study in Ontario, Canada
  1. Eric De Prophetis1,
  2. Vivek Goel1,2,
  3. Tristan Watson3,
  4. Laura C Rosella1
  1. 1Dalla Lana School of Public Health, University of Toronto, Toronto, Ontario, Canada
  2. 2Institute of Health Policy, Management and Evaluation, University of Toronto, Toronto, Ontario, Canada
  3. 3Populations and Public Health, ICES, Toronto, Ontario, Canada
  1. Correspondence to Dr Laura C Rosella; laura.rosella{at}utoronto.ca

Abstract

Objective To examine if low life satisfaction is associated with an increased risk of being hospitalised for an ambulatory care sensitive condition (ACSC), in comparison to high life satisfaction

Design and setting Population-based cohort study of adults from Ontario, Canada. Baseline data were captured through the Canadian Community Health Survey (CCHS) and linked to health administrative data for follow-up information.

Participants 129 467 men and women between the ages 18 and 74.

Main outcome measures Time to avoidable hospitalisations defined by ACSCs.

Results Life satisfaction was measured at baseline through the CCHS and follow-up information on ACSC hospitalisations were captured by linking participant respondents to hospitalisation records covered under a single payer health system. Within the study time frame (maximum of 14 years), 3037 individuals were hospitalised. Older men in the lowest household income quintile were more likely to be hospitalised with an ACSC. After controlling for age, sex, socioeconomic status (SES) and other behavioural factors, low life satisfaction at baseline had a strong relationship with future hospitalisations for ACSCs (HR 2.71; 95% CI 1.87 to 3.93). The hazards were highest for those who jointly had the lowest levels of life satisfaction and low household income (HR 3.80; 95% CI 2.13 to 6.73). Results did not meaningful change after running a competing risk survival analysis.

Conclusions This study demonstrates that poor life satisfaction is associated with hospitalisations for ACSCs after adjustment for several confounders. Furthermore, the magnitude of this relationship was greater for those who were more socioeconomically disadvantaged. This study adds to the existing literature on the impact of life satisfaction on health system outcomes by documenting its impact on avoidable hospitalisations in a universal health system.

  • life satisfaction
  • avoidable hospitalization
  • epidemiology
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • This is the largest study to combine detailed survey variables linked to complete hospitalisation outcomes form from health administrative databases to study the impact of life satisfaction on ambulatory care sensitive conditions (ACSCs).

  • The linkage of various databases allowed for control of a wide range of confounders in addition to objective measures of comorbidity.

  • The linkage and study design allowed measurement of baseline life satisfaction on a sample of individuals who did not have an ACSCs in the 2 years prior, which overcomes a limitation of previous studies that may have been influenced by the bidirectional relationship between poor health and life satisfaction.

  • Life satisfaction was measured at one point in time, and thus, we were unable to capture changes in life satisfaction over the study period.

  • This study did not directly address mechanisms by which life satisfaction could influence ACSCs.

Introduction

There has been broad recognition that upstream determinants have an influence on a range of health outcomes, including social determinants and risky health behaviours. In addition, subjective well-being, specifically life satisfaction is increasingly being recognised as playing an important role.1–3 Shifting from health outcomes to a health system perspective, many ambulatory care sensitive conditions (ACSCs) are preventable and result in high amounts of healthcare resources if left unmanaged.4 5 Hospitalisations for these conditions (eg, diabetes) are regarded as avoidable due to the relationship of timely access to primary care preventing complications and acute episodes.5–7 In many countries, hospitalisations for these conditions are used as an indicator to measure the effectiveness of the healthcare system.8

Positive affect (ie, positive emotions, happiness), is the degree to which an individual experiences positive emotion,9 has been shown to be independently associated with reduced risk of coronary heart disease, even after adjustment for depressive symptoms.10 Similarly, life satisfaction is associated with a reduced risk of various chronic conditions.2 11 Evidence from a recent prospective cohort study concluded lower life satisfaction was associated with a higher incidence of cancer, stroke and type 2 diabetes11 and others have shown its relationship with health behaviours.12 Positive affect is distinct from life satisfaction in that it is an emotional response describing positive emotions or happiness, which can be transitory. This is conceptually distinct from the more stable and complex measure of being satisfied with life, which is based on an individual’s self-judgements of several factors that they feel are valuable to their life.

Research on life satisfaction has also recently shown to impact future high-health care use.13 Individuals with lower life satisfaction had three times the odds of being in the top 5% of healthcare resource utilisation.13 However, another important system indicator that is absent from the current life satisfaction literature is avoidable hospitalisations, such as those caused by ACSCs. Given the evidence around the interrelationships between emotional factors and chronic disease self -management14 and in particular how emotional and psychological distress can influence the ability of a patient to manage their chronic condition,15 16 more research is needed that quantitatively measures the impact of subjective well-being on important healthcare system outcomes.

On a recent search, only eight studies were identified that fit the criteria of a similar exposure to life satisfaction and hospitalisations for one or more ACSCs as an outcome. Previous studies mostly focused on singular conditions, with small sample sizes, limited follow-up time and clinical or convenience-based samples.17 18 Regarding the length of follow-up, these ranged from 6 months19 to 4 years.20 Studies with shorter follow-up times can be limited because many of the conditions regarded as ACSCs are also chronic diseases (eg, diabetes) require a longer follow-up time necessary to allow for outcomes to be observed. Furthermore, hospitalisations for ACSCs are a relatively rare event in Canada where people with an ACSC hospitalisation constitutes only 0.4% of the population aged 12–74,21 and therefore, the limited sample sizes in previous studies may not be sufficient to observe an effect. We address these limitations by conducting the largest population-based cohort study to date. The low rate of ACSCs is typical of similar health systems in Europe, the UK and Australia.

The primary objective of this study was to determine if poor life satisfaction increased the risk of being hospitalised for an ACSC in a relatively healthy baseline cohort among an adult population using linked survey and complete hospitalisation records. Our secondary objective was to determine if this association was stronger among those with lower socioeconomic status (SES).

Methods

Participants

The study was a longitudinal population-based cohort study of adult Ontario participants of the Canadian Community Health Survey (CCHS) pooled across five cycles: cycle 2.1 (2003–2004), cycle 3.1 (2005–2006), cycle (2007–2008), cycle (2009–2010) and cycle (2011–2012). The CCHS conducted surveys on a 2-year collection cycle (ie, 2003, 2005), therefore, no cycle existed for 2004 and 2006 and cycle naming conventions changed after 2005 to remove cycle numbers. The CCHS is a cross-sectional survey administered by Statistics Canada, representative of 98% of the Canadian population aged ≥12 years living in private dwellings with response rates>75%.22 23 The respondents in the CCHS survey consented to participate in the survey and have their data linked to administrative data for research purposes. Where the CCHS is a cross-sectional survey, the longitudinal aspect of this study is achieved by retrieving the outcome measure from hospital administrative databases. The data were linked to population-based health administrative data held at ICES. These data sources capture all hospitalisation records for every person living in the province of Ontario covered under the single-payer health system.

Eighty per cent of the CCHS survey respondents consented to have their data linked to the single-payer health insurance data, referred to as the Ontario Health Insurance Plan, which captures all related healthcare encounters. All survey respondents were linked to Ontario’s population registry, the Registered Persons Database (RPDB), which captures core demographic and clinical information as well as death, in addition to the Discharge Abstract Database. The analytical sample included adults (aged 18–74) who reported on life satisfaction and did not have an ACSC-related hospitalisation in the 2 years before their CCHS interview date (figure 1).

Figure 1

Flow chart of study participants. ACSC, ambulatory care sensitive condition; RPDB, Registered Persons Database.

Measures

Self-reported life satisfaction, the primary exposure variable, is captured from CCHS. The question that respondents answered regarding life satisfaction was, ‘How satisfied are you with your life in general?’ With response options being: very satisfied, satisfied, neither satisfied nor dissatisfied, dissatisfied or very dissatisfied. Due to small sample sizes within each category, we collapsed the categories of dissatisfied and very dissatisfied. Furthermore, we did not hypothesise significant conceptual differences between these two categories related to hospitalisations for ACSCs.

The primary outcome variable was hospitalisations for an ACSC, which we used as a composite outcome. The ACSCs that we chose to report on are grand mal status and other epileptic convulsions, chronic obstructive pulmonary disease, asthma, diabetes, heart failure and pulmonary oedema, hypertension and angina. These conditions are in accordance with the Canadian Institute for Health Information’s (CIHI) methodology24 and this composite outcome is an established health system indicator in Canada. The CIHI indicator applies only to individuals under the age of 75 as the hospitalisations in those above the age of 75 are not as clearly avoidable through timely and effective primary care. See online supplementary table 1 for the list of included conditions and their corresponding International Classification of Diseases 10th Revision (ICD-10) codes.

Aggregate diagnosis groups (ADGs)25 were captured through administrative data as a summary measure of comorbidity and are based on the Johns Hopkins ACG System, which is a person-focused, diagnosis-based method of categorising subjects’ illnesses.22 ADGs have previously been used and validated as a reliable method of comorbidity adjustment in the Ontario population,26 and we used V.10.0.1 in this analysis. Additionally, we used the Ontario Marginalization Index (ON-Marg) as a measure of SES. The ON-Marg is a census-based, geographically derived index that was used to calculate area-level material deprivation.27 Specifically, the material deprivation dimension measures the proportion of the population within a geographical region that is low income, without high school diploma, lone parent families, receiving government transfer payments, unemployed and living in dwellings in need of repair. All other covariates were captured through self-report from the CCHS interview questions.

CCHS variables

Household income quintile categorises individuals based on their total household income in addition to the number of individuals living in the household. Individuals are then ranked from the lowest levels of household income (Q1) to the highest (Q5). Body mass index (BMI) was categorised into five categories ranging from ‘underweight’ (BMI <18.5 kg/m2) to ‘severely obese’ (BMI >34.9 kg/m2). Physical activity was based on an individual’s self-reported daily energy expenditure and further categorised into three levels: active, moderately active and inactive. Smoking status measured an individuals’ self-reported past and present smoking habits by considering both the total amount of cigarettes smoked and the type of smoker they are (eg, daily vs occasional). This variable was categorised into three levels: current smoker, former smoker and never smoker. Alcohol consumption was based on the participant’s sex and the quantity of alcohol consumed each day. This variable was then categorised into four levels: heavy drinker, moderate drinker, light drinker and never drinker. Mood disease was captured through the CCHS interview question ‘have you ever been diagnosed by a health professional for depression, bipolar disorder, mania or dysthymia?’ This variable was used to control for depression. Anxiety disorder, which was captured through the question: ‘have you ever been diagnosed by a health professional for an anxiety disorder such as a phobia, obsessive–compulsive disorder or panic disorder?’ Education level is a derived variable which indicates the highest level of education acquired by the participant; this variable was explored as a potential indicator for SES.

Statistical analysis

We calculated the distribution of demographic, socioeconomic, health status and behaviour characteristics according to ACSCs and life satisfaction. Excluding those with an ACSC in the year prior allowed for the investigation of the upstream determinants (ie, the factors associated with future development of an ACSC hospitalisation in a cohort who were without a recent hospitalisation for one of these conditions).

Cox proportional hazards models were used to estimate the hazards associated with baseline life satisfaction on the risk of being hospitalised for an ACSC. Time is defined as survey date until disease or censoring for study endpoint (max follow-up until 31 March 2017) or death. The models were used to quantify the association between life satisfaction and the hazard of being hospitalised for an ACSC using ‘very satisfied’ as the referent category. We calculated unadjusted, age-adjusted and sex-adjusted, minimally adjusted and fully adjusted models to transparently demonstrate the impact of adjustment. The minimally adjusted model controlled for age, sex and household income while the fully adjusted model included age, sex, household income, smoking status, alcohol consumption, physical activity and BMI. In order to show the fully adjusted results were not meaningfully affected by over comorbidity, mood disorders and anxiety, we ran three additional models, which controlled for ADG score, mood disease and anxiety separately. We ran these models separately in order to quantify their impact on the life satisfaction hazards in the fully adjusted model that controlled for sociodemographic and behavioural factors that we conceptualised as confounders (ie, age, sex, household income, smoking status, alcohol consumption, physical activity and BMI).

To evaluate the proportional hazards assumption for life satisfaction, a plot of log (-log (survival)) versus the logarithm of follow-up time in days was run. With this method, the proportional hazards assumption is met if the plot produces parallel curves.

A joint-effects model was used to test if the relationship between life satisfaction and avoidable hospitalisations varied by SES. A joint-effects variable, which contains each combination of life satisfaction and household income, was included in the model while controlling for age, sex, smoking status, alcohol consumption, physical activity and BMI.

We ran sensitivity analyses by rerunning the fully adjusted model with a consistent survival time of 5 years. Additionally, models for individuals who did not have an event in the first 2 years of the study were run to control for undocumented comorbidity. The subdistribution hazards model, which was initially developed by Fine and Gray,28 was run to test the possibility of death behaving as a competing risk.

Survey and bootstrap sampling weights provided by Statistics Canada were applied in all descriptive and survival regression analyses to account for the complex survey design and to maintain population representativeness.29 The bootstrap sampling weights were applied using balanced repeated replication, in order to properly calculate CIs. Finally, all statistical analyses were performed in 2018 and 2019 using SAS V.9.4.

Patient involvement

Patients were not involved in the development of the research question, outcome measures, recruitment, design or the implementation of the study objectives. Furthermore, no patients were consulted on the interpretation of results, and there are no plans to disseminate the results of this study to the relevant participants or their communities.

Results

After combining the five cycles of data linked to the RPDB and excluding those less than 18, older than 74, in multiple survey cycles, missing life satisfaction, had an ACSC-related hospitalisation in the 2 years before the start of the study resulted in 129 467 individuals remaining. Those who experienced an ACSC-related hospitalisation 2 years before their survey interview date were removed to examine the impact of life satisfaction on future hospitalisations for an ACSC and reduce the possibility of reverse causality (figure 1). In the case of people in multiple survey cycles (n=1589), we used the earliest survey response. Each respondent was followed for a maximum of 14 years or until the study end date, after which we determined that 3037 individuals had experienced an ACSC-related hospitalisation.

The distributions of baseline characteristics according to life satisfaction categories are shown in table 1. Those with the lowest levels of life satisfaction (very dissatisfied) compared with the highest level of life satisfaction (very satisfied) had a lower proportion in the youngest (18–34 years) age group (14.5% vs 30.6%), more likely to have less than secondary education (12.9% vs 4.62%), had a greater proportion in the lowest income quintile (38.9% vs 9.5%), higher comorbidity levels (ADG score 10.1 vs 3.42), and higher disease-related risk factors such as smoking (43.8% vs 17.4%) and physical inactivity (70.6% vs 40.7%) (table 1).

Table 1

Weighted* distributions of baseline characteristics across five levels of life satisfaction (n=129 467)

The distributions of all the cohort characteristics according to ACSC status are show in table 2. Those that had an ACSC during the follow-up compared with those that did not were more likely to be in the older age group (65–74 years) (21.6% vs 9.55%), more likely to have less than secondary education (19.2% vs 6.01%), had a greater proportion in the lowest income quintile (20.7% vs 13.5%), higher comorbidity levels (ADG score 11.8 vs 3.86) and higher disease-related risk factors such as smoking (36.7% vs 21.5%) and physical inactivity (61.4% vs 48.2%) (table 2).

Table 2

Weighted* distributions of baseline characteristics according to hospitalisation for an ACSC (n=129 467)

Life satisfaction had a strong unadjusted relationship with hospitalisation for ACSC (unadjusted Kaplan-Meier curves are shown in online supplementary figure 1). Table 3 presents the unadjusted and adjusted HRs for the relationship between life satisfaction and ACSC-related hospitalisations, which include models that adjusted for age and sex, and then further adjusted for socioeconomic and lifestyle factors. Although full adjustment does reduce the size of the effect, the HR of an individual with the lowest levels of life satisfaction (dissatisfied and very dissatisfied combined) compared with those who were very satisfied is 2.71 (95% CI 1.87 to 3.93). The observed relationship follows a dose–response pattern, or in other words, the HRs increase in size for each decreasing level of life satisfaction. For example, in the fully adjusted model, the HR for the middle life satisfaction category (neither satisfied nor dissatisfied) was 1.71 (95% CI 1.36 to 2.14) while the satisfied category produced an HR of 1.32 (95% CI 1.15 to 1.50). Furthermore, the additional analyses (table 4) that controlled for ADG score, mood disease and anxiety separately did not substantially reduce the observed effect sizes with the added adjustment of ADG score having the largest impact (HR of 2.42, 95% CI of 1.68 to 3.51). Finally, as the plot of survival by time according to life satisfaction produced parallel curves, the proportional hazards assumption was satisfied.

Table 3

Multivariable AHRs and 95% CIs from proportional hazards regression (n=129 467)

Table 4

Multivariable AHRs and 95% CIs from proportional hazards regression after additionally adjusting for ADG score, mood disease and anxiety (n=129 467)

Regarding the joint effects model (table 5), individuals who identified as having both low life satisfaction and low household income produced a fully adjusted HR of 3.80 (95% CI 2.13 to 6.73). Therefore, in comparison to the fully adjusted model presented in table 3 (HR of 2.71), poor SES individuals are at an increased risk of being hospitalised when reporting low life satisfaction.

Table 5

Multivariable AHRs and 95% CIs from joint effects models (n=129 467)

A range of sensitivity analyses were conducted to test the robustness of the study findings. Both the 5-year survival model (where everyone was limited to a 5-year follow-up) and the model which excluded individuals who had an event within the first 2 years did not change the effect sizes to a significant degree. The HRs for the lowest levels of life satisfaction compared with those who were very satisfied were 2.74 and 2.77, respectively (online supplementary table 2). In the competing risk analysis (online supplementary table 3), the unadjusted HRs produced in the subdistribution hazard model were similar to those produced in the final model (4.38 vs 4.51, respectively).

Discussion

This study focuses on how life satisfaction can impact health system indicators such as avoidable hospitalisations in a general population cohort without a recent ACSC. We investigated this relationship and accounted for a wide variety of sociodemographic and behavioural risk factors. We saw a robust association that poor life satisfaction had a strong independent relationship with future ACSC hospitalisations. The lowest levels of life satisfaction (dissatisfied and very dissatisfied) being associated with almost a three times higher hazard of an avoidable hospitalisation compared with those who were very satisfied after accounting for several sociodemographic and behavioural confounders.

Previous studies have linked life satisfaction and related exposures (eg, positive affect or happiness) to health outcomes such as stroke11 and heart disease.10 Life satisfaction has also been shown to be associated with a wide variety of health behaviours.30 31 For instance, one study found that those who exercised more were generally happier.32 Furthermore, life satisfaction has also been shown to be experienced differently across categories of SES.33 Due to the detailed survey, variables available from the survey data and the linkage of these data to complete hospitalisation outcomes from health administrative databases were able to adjust for these health behaviours and measures of SES in our analysis as well as examine how the influence of life satisfaction is strengthen or lessened across levels of SES.

This study addresses an important gap in the literature by providing a robust population sample size and examining how life satisfaction related to a meaningful health system outcome. There are few studies that addresses life satisfaction or other forms of subjective well-being and their relationship to hospitalisations for ACSCs. Furthermore, the existing research is limited by small sample sizes and limited follow-up times to capture the relatively rare ACSC event. This study also has more direct implications for the health system given ACSCs are defined as conditions for which hospitalisations should be prevented, given timely and effective access to primary care.7 Considering the preventable nature of these conditions, hospitalisations for ACSCs are an ineffective use of healthcare resources and insight into the risk factors for these conditions can help improve health system functioning.

Unlike previous studies, we measured baseline life satisfaction on a sample who did not have an ACSC hospitalisation in the 2 years prior while also presenting analyses that additionally controlled for comorbidity. These aspects of the study help mitigate the possibility that that poor life satisfaction could have been the result of the bidirectional relationship between poor health and life satisfaction. A possible explanation for the observed results is that individuals who experience poor life satisfaction tend to have higher rates of depression, given its observed relationship with poor health outcomes.34 35 To address this, we further adjusted for mood disease and found that this has little effect on attenuated the observed association. The survey did not include continuous measures that could capture subthreshold levels or undiagnosed anxiety or depression.

Limitations

We acknowledge that this study has some limitations and interpretive cautions. First, this study was an observational study and although we controlled for several potential confounders and excluded those with a recent history of ACSC at baseline through linkage, we cannot rule out the possibility of unmeasured or residual confounding. We note, however, that the effect sizes are large, and this study did control for many more confounders than previous studies through a combination of survey data in addition to health administrative variables to capture comorbidity. Second, we measured life satisfaction at one point in time (survey interview date). Research on life satisfaction has shown it to be consistent over time36; however, it is possible that life satisfaction could have changed during the study period. Therefore, we only capture the effect from that initial time point and cannot account for the influence of changes in life satisfaction that happen following the baseline assessment. Furthermore, there are other instruments that can be used to assess life satisfaction that were not available in this survey. Finally, life satisfaction is a subjective measure. It has been shown to be an accurate and robust measure, it is still up to the individuals to judge and reflect on their life satisfaction. This means that its meaning can differ based on the individual, which can result in variation in the exposure. Despite these limitations, this study has still provided an essential contribution to the literature by being one of the first assessments of life satisfaction on avoidable hospitalisation that uses a longitudinal population-based study design while measuring this exposure before the hospitalisation outcome.

Mechanisms

This study did not directly address mechanisms by which life satisfaction can influence avoidable hospitalisations; however, these have been studied by others. Out of the existing literature, there are three main mechanisms that could help explain the observed relationship: behavioural, health service use and biological. Subjective well-being has been shown to impact a variety of health behaviours, such as increased physical activity37–40 and reduced smoking habits.40 However, many of these studies lacked proper adjustment of confounding variables such as SES and psychological distress.40 Regarding health service use, a study from 2014 found that participants who were identified as having a greater purpose in life were more likely to receive preventative health services such as mammograms or colonoscopies.41 What remains to be seen regarding this mechanism is whether the use of these services reflects access to primary care, or the decisions made by the individuals themselves to seek these services. However, a recent population-based cohort study noted how hospitalisations for ACSCs could not be explained by a lack of access to primary care,42 and therefore, it is possible that this phenomenon is a result of individual decisions that are influenced by their satisfaction with life. Finally, biological mechanisms could also play a contributing role. A recent meta-analysis noted how individuals who were identified as having greater psychological well-being experienced favourable lipid profiles.43 However, these associations were largely mitigated once behavioural characteristics were taken into account.43 Our study did not compare the extent to which life satisfaction is associated with ACSC compared with other types of hospitalisation. This was out of scope for this current study, but in a future analysis could help provide insight into mechanisms that are specific to preventable hospitalisations. Our study did not examine whether ACSC hospitalisations were reduced for Ontario hospitals, but instead focused on individual’s risk based on their life satisfaction measure. In future studies, analyses at the system level (ie, ACSC hospitalisations in Ontario hospitals) would supplement these findings to better understand the health system impact.

Conclusions

This study demonstrates that life satisfaction is associated with hospitalisations for ACSCs, even after controlling for socioeconomic characteristics, health behaviours, comorbidities and mental health factors. Furthermore, more socioeconomically deprived individuals were shown to be at an increased risk. While governments plan to improve the sustainability and functioning of their health systems, there is a greater need to understand social supports that can improve life satisfaction costly and preventable conditions such as ACSCs. The findings of this study suggest that broader considerations, such as life satisfaction, can potentially influence avoidable hospitalisations, a burden to individuals and healthcare systems. Further research in this area may contribute to the development of wide-ranging approaches to target a potentially avoidable burden on the health system.

References

Footnotes

  • Twitter @LauraCRosella

  • Contributors LCR and VG conceptualised the study, TW created the database and EDP planned the analysis, analysed the data and drafted and revised the paper. LCR is the guarantor. All authors contributed to the manuscript and revised the draft paper. All have approved the final version for publication.

  • Funding This research was supported by a Canada Research Chair held by LCR (CRC-950-23072). This study was supported by ICES, which is funded by an annual grant from the Ontario Ministry of Health and Long-Term Care.

  • Disclaimer The opinions, results and conclusions reported in this paper are those of the authors and are independent from the funding sources.

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval This project received ethics approval from the Health Sciences Research Ethics Board at the University of Toronto (Ref # 36123, 15 August 2018).

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement Data may be obtained from a third party and are not publicly available. The data set from this study is held securely in coded form at ICES. While data sharing agreements prohibit ICES from making the data set publicly available, access may be granted to those who meet prespecified criteria for confidential access (available at www.ices.on.ca/DAS). The full data set creation plan and underlying analytical code are available from the authors on request, understanding that the programmes may rely on coding templates or macros that are unique to ICES.