Article Text

Original research
Does gradually returning to work improve time to sustainable work after a work-acquired musculoskeletal disorder in British Columbia, Canada? A matched cohort effectiveness study
  1. Esther T Maas,
  2. Mieke Koehoorn,
  3. Christopher B McLeod
  1. School of Population and Public Health, University of British Columbia, Vancouver, British Columbia, Canada
  1. Correspondence to Dr Esther T Maas, School of Population and Public Health, University of British Columbia, Vancouver, V6T 1Z3, Canada; esther.maas{at}ubc.ca

Abstract

Objective This study investigates if gradual return to work (GRTW) is associated with full sustainable return to work (RTW) for seriously injured workers with a musculoskeletal disorder (MSD), in British Columbia, Canada.

Methods This is an effectiveness study using a retrospective cohort study design. Accepted workers’ compensation lost-time claims were extracted for workers with an MSD who were on full work disability for at least 30 days, between 2010 and 2015 (n=37 356). Coarsened exact matching yielded a final matched cohort of 12 494 workers who experienced GRTW at any point 30 days post-injury and 12 494 workers without any GRTW. The association between GRTW and sustainable RTW through to end of 12 months was estimated with multivariable quantile regression.

Results Workers who were provided with GRTW experienced more time-loss days until sustainable RTW between the 2nd and 5th months after the first time-loss day (<50th quantile of time loss), but less time-loss days until sustainable RTW between the 6th and 12th months of work disability (70th quantile of time loss), with the largest effect for women, workers with soft-tissue injuries and workers in the manufacturing or trades sector (all in the 60th and 70th percentile, after 6–7 months of time loss).

Conclusions For seriously injured workers with at least 30 days of disability due to a work-acquired MSD, the effect of GRTW becomes apparent at longer disability durations (more than 6 months), with larger beneficial effects for women, workers with soft-tissue injuries and for trade and manufacturing sectors.

  • epidemiology
  • sick leave
  • musculoskeletal system
  • occupational health

Data availability statement

Data may be obtained from a third party and are not publicly available. The BC Freedom of Informant and Protection of Privacy Act [RSBC 1996] CHAPTER 165 (the legislation that authorises us to use the data for our study) does not authorise this data be access or view outside of Canada. Under our obligations under this legislation and under an information sharing agreement with WorksafeBC to use the data, we are prohibited from sharing the data outside of Canada and outside our research team. If other researchers wanted access to the data, they would need to negotiate access directly with the data provider, WorkSafeBC.The data for this project can only be accessed in Canada and making an application through https://www.popdata.bc.ca/dataaccess to WorkSafeBC and by entering into an Research Agreement with WorkSafeBC.

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.

Key messages

What is already known about this subject?

  • Gradual return to work opportunities after a period of work disability due to a musculoskeletal disorder are considered a key aspect of work disability management. There is limited evidence on the effectiveness of gradual return to work for longer term claims as a population-based intervention at a jurisdictional level, where there may be specific policies or practices that determine the likelihood and effectiveness of gradual return to work.

What are the new findings?

  • Gradual return to work was associated with sustainable return to work after 6 months of work disability for workers with a work-acquired musculoskeletal disorder, compared with similar workers without gradual return to work, in British Columbia, Canada. The positive long-term effects of gradual return to work were larger for women, workers with torso sprains and workers in the manufacturing and trades sector.

How might this impact on policy or clinical practice in the foreseeable future?

  • For seriously injured workers (off work for at least 30 days) with a musculoskeletal disorder, the effect of gradual return to work becomes more apparent in the longer term, at least 6 months after the first time-loss day. A focus on increasing the use of gradual return to work for longer term claims could improve overall work disability duration.

Introduction

Work injury rates are declining in most high-income countries, while work disability durations remained static.1–4 There is an increased interest in workplace-based return to work (RTW) interventions, like gradual RTW (GRTW),5–7 as part of work disability management.8 GRTW provides workers with the opportunity to increase working hours and workload, and to limit or modify work tasks while recovering from an injury, with the goal to return to full hours and duties.9 GRTW is intended to reduce work disability or work disability duration.7 10–13 Krause et al found that modified work facilitates RTW for temporarily and permanently disabled workers,10 and a recent systematic review from Cullen et al supported a strong association between gradual work activities and a reduction in lost time associated with work disability.7

Studies have focused on GRTW interventions within the workplace. However, there is limited evidence for the effectiveness of GRTW for longer term disability due to serious injury (ie, at least 30 days of work disability), or as a population-based intervention at a jurisdictional level. In the case of the Canadian province of British Columbia (BC), GRTW is encouraged but not legislated. Accordingly, GRTW in BC, and many other jurisdictions, is not standardised by the disability compensation system and the provision of GRTW by employers varies substantially.14–16

In BC, almost 80% of all work disability days between 2009 and 2013 were due to a work-acquired musculoskeletal disorder (MSD).2 Besides being a major cause of work disability and reduced productivity, MSDs are often episodic and recurrent, causing longer duration claims. GRTW may be an appropriate work disability management strategy to maintain workplace interaction while in a recovery phase.8 The objective of this study is to investigate the effects of GRTW on sustainable RTW for workers with a disability duration 30 days or longer for a work-acquired MSD, in BC, Canada.

Methods

Jurisdictional context

Workers who experience a recognised work-acquired injury or disease in BC are provided with disability benefits, medical aid and rehabilitation services by WorkSafeBC (the workers’ compensation system in BC). WorkSafeBC is funded through employer-paid insurance premiums and provides short-term disability wage replacement to injured workers with the goal to RTW in a timely manner. For most workers, short-term disability payments, representing 90% of workers’ pre-injury wage, are provided until workers fully RTW. Workers who do not completely recover from their injury and who have a permanent partial impairment are eligible for vocational re-training and/or long-term disability benefits. Over 97% of the BC labour force is covered by this work disability insurance system.17

Study design and data

This is an effectiveness study using a retrospective cohort study design, with a 1 year follow-up period from the first recorded time-loss day. Administrative claims data collected by the province of BC’s Ministry of Health (Medical Registry data, Medical Services Plan (MSP) data, Hospital discharge abstract data and Pharmanet data) and WorkSafeBC (claims and RTW data) were linked at the individual level, and used to select workers with an accepted work-acquired MSD lost-time claim filed between 1 January 2010 and 31 December 2015.18–22 More information about the databases can be found in online supplemental appendix 1.

Supplemental material

The data were obtained and de-identified by Population Data BC, a multi-university, data and education resource that supports access to data on BC’s 4.6 million residents for research.23 The Behavioral Research Ethics Board at the University of British Columbia approved the study (certificate no. H17-02019).

MSDs were categorised into nine major categories using the Barell matrix for musculoskeletal injuries24 (sprains/strains, fractures, dislocations) and International Classification of Diseases (ICD)-9-CM codes for musculoskeletal diseases (dorsopathies and rheumatism (excluding the back)). Sprains/strains and fractures were divided into three body regions: (a) head and neck/spine/back/torso; (b) upper extremities; (c) lower extremities.

The cohort was restricted to the first MSD work disability claim per worker in the study period, with a follow-up period of 365 days. Most workers RTW within 1 year after the first time-loss day.25 Workers were excluded from the study for the following reasons (figure 1):

  1. Healthcare only claims (no time loss).

  2. Exclusions based on cohort definitions:

    • Non-MSD-related claims (defined using ICD-9-CM codes and National Work Injury Statistics Program Nature of Injury codes).

    • Age <15 or ≥65 years.

    • Multiple jobholders at the time of injury.

    • Claims from self-insured industry sectors. There is less RTW data available for these firms as WorkSafeBC interacts differently with these companies.

    • Work-related fatal injuries.

  3. Exclusions based on missing data on firm size, industry, wage or gender.

  4. Exclusions based on missing MSP registration date, or lack of registration in the year before and after injury.

  5. Pregnancy in year before claim.

Figure 1

Construction of a cohort of workers with work-acquired musculoskeletal disorders from compensation claims data in British Columbia between 2010 and 2015. GRTW, gradual return to work; MSD, musculoskeletal disorder; MSP, Medical Services Plan.

A total of 127 181 unique MSD claims were included in the study, including 53 102 workers (41.8%) with at least 1 day of GRTW.

The cohort was restricted to those on full disability for at least 30 days after the first time-loss injury day, a threshold to distinguish between short and longer work disability duration.26 Further, the provision of GRTW within 30 days of the first time-loss day is less common for injuries with longer disability durations25 27 and workers with these injuries may be less likely to benefit from GRTW.

Outcomes

The primary outcome was calendar days from the first time-loss day until sustainable RTW (end of claim, no further disability days) up to 365 days. Secondary outcomes were (a) time-loss days until full sustainable RTW or sustainable GRTW (end of claim defined as GRTW, no further full disability days) until 365 days and (b) time-loss days until sustainable RTW up to 730 days.

Explanatory variable

GRTW was defined as temporarily changing a worker’s duties, hours and/or days of work within the same job/same employer while being on a compensation claim. Partial benefits compensate for the income loss due to the disability. This is distinct from vocational rehabilitation that would be provided to workers with permanent partial work disability.

GRTW was indicated in WorkSafeBC’s RTW event data (work disability, GRTW, non-RTW or RTW) daily after injury. GRTW was measured as ‘yes’ versus ‘no’ for at least 1 day of GRTW in each month of follow-up, starting at the second month after the first time-loss day. More information on GRTW can be found in Maas et al.27

Covariates

The following variables were included as potential confounders based on prior research and those with the greatest association with being provided GRTW (as explained below):

  • MSD, as specified above.

  • Age (15–24, 25–34, 35–44, 45–54, 55–64 years).

  • Gender (woman, man).

  • Annual wage (<$20 000, $20 000–$39 999, $40 000–59 999, >$59 999 Canadian dollars (CAD$)).

  • Occupation, classified according to Statistics Canada’s Standard Occupational Classification.28

  • History of prior claims (yes/no): at least one accepted claim for any type of injury or illness in the preceding 10 years to the MSD claim.

  • Industry sector, classified according to the WorkSafeBC industry classification structure.29

  • Size of the workers’ firm measured as fulltime equivalent (FTE) workers employed by the firm at the time of the MSD injury (<20, 20–99, 100–499, 500–999, >999 FTE).

  • Claim year was identified as start date of the claim (2010–2015).

  • Opioid, non-steroidal anti-inflammatory drugs, and skeletal muscle relaxants use within 30 days after the first time-loss day.

  • Comorbidities (see online supplemental appendix 1).

Statistical analyses

Descriptive statistics identified differences between workers with and without GRTW in the cohort of workers that were on full disability for 30 days.

A logistic regression model assessed the relative contribution of injury, sociodemographic, workplace, healthcare and temporal characteristics in predicting the likelihood of GRTW 30 days or more after the first time-loss day to inform the matching strategy. The model was built according to recommended methodological guidelines for multivariable logistic regression analysis.30 31 We have used the prediction model in addition to prior independent work to select the matching variables.27

Workers with and without GRTW were matched on the identified characteristics using coarsened exact matching (CEM).32 33 The aim of matching was to find balance across the multidimensional distribution of covariates. This reduces the degree of dependence on the estimation model for the outcome variable and reduces estimation bias.34

Quantile regression was used to estimate the number of additional time-loss days until sustainable full RTW for workers with GRTW compared with workers without GRTW at different points in the disability distribution. This method is ideal for analyses where the outcome variable can be skewed, and outliers may affect the mean.35 Quantile regression is not restricted by the proportional hazard’s assumption like cox proportional hazards, which require complex modifications when this assumption is violated. In contrast, quantile regression estimates the effect on the outcome directly instead of the hazard rate and this can be more intuitive, easier to interpret and relevant to decision-making. Based on visualisation of time loss days, the quantile regressions were specified a priori from the 10th to the 90th percentile in 10% increments. All covariates were included in the model.

The main model was used to estimate predicted time-loss days between workers with and without GRTW until full sustainable RTW up to 365 days after the first time-loss day at various points in the disability distribution, holding other independent variables at their mean values.36

To test study assumptions, secondary models were used to estimate differences in time-loss days until (a) sustainable RTW or sustainable GRTW (no further full disability days) up to 365 days (the assumption being that sustainable work accommodation is preferred over continued time-loss as a positive outcome), and (b) sustainable RTW until 730 days. For the latter model, the cohort was redefined and re-matched for claims occurring between 2010 and 2014 to enable 2 years of follow-up. Further, separate models were constructed for the unmatched cohort for comparison of results, and for the cohort of workers employed at firms who had offered GRTW at least once between 2010 and 2015 to control for workplace variability. This leaves out workers at firms that never offered GRTW in the study period and who are not ‘at risk’ of being offered GRTW.

All analyses were performed in Stata V.16.0 (Stata Corp).

Results

Descriptive statistics

The characteristics used for cohort matching are summarised in table 1.

Table 1

Descriptive statistics of workers on full work disability for 30 days with an accepted MSD lost-time claim between 2010 and 2015, by GRTW status, in the Canadian jurisdiction of British Columbia

The unmatched cohort of workers with a work-acquired MSD lost-time claim who were on work disability for at least 30 days between 2010 and 2015 comprised 37 356 workers, of which 47% of workers had experienced GRWT between 30 days and 1 year after injury. The sample comprised 65% men, had a median age of 46 years (IQR 35–54) and a median annual wage prior to injury of approximately CAD$ 42 000 (IQR 30 000–60 000). Back sprains and strains were the most common disorder type (34%). The most common industry sector was the services sector (39%), and the most common occupation was trades, transport and equipment operations (38%). The median firm size associated with an injured worker was 46 FTE (IQR 9–254). More information and clinical characteristics (healthcare utilisation and medication use) are shown in online supplemental appendix 2.

Proportionately more women were provided with GRTW (56%) compared with men (42%). The provision of GRTW among injured workers increased with increasing age, wage and year of the study. There were more workers with GRTW in larger firms, in the trade industry and on fixed job shifts. There was little variation between workers with and without GRTW in terms of clinical characteristics in the 30 days after the first time-loss day. Clinical characteristic, as one of the proxy measures for injury severity, was not found as being predictive of being offered GRTW.

Likelihood of workers to receive GRTW

Firm size, gender, claim year, shift type, annual wage, occupation and industry sector were statistically significant (p<0.05) predictors of the likelihood that a worker would be provided with GRTW based on the multivariate logistic regression model, as well univariate F tests (online supplemental appendix 3), and were used for the CEM. Occupation and industry sector were collinear and only industry sector was used for CEM.

Matching

After using CEM, the matched cohorts resulted in equal distributions in matched covariates between workers with (n=12 494) and without (n=12 494) GRTW (table 1). The multivariable imbalance measure L1 improved from 0.33 to 0.00 by CEM, indicating that balance in the matched cohorts could not be improved.

Difference in time-loss days between workers with and without GRTW

Among the matched cohort, the quantile regression results for the main model indicated that between the 10th and the 50th percentile of the disability distribution, workers with GRTW had on average 24 more work disability days until sustainable RTW compared with workers without GRTW (eg, at the 10th percentile=β 21.8; 95% CI 20.7 to 22.8) and at the 50th percentile=β 24.6; 95% CI 19.9 to 29.2) (table 2). By the 70th percentile of the distribution, however, workers with GRTW had fewer disability days than workers without GRTW (β −9.4; 95% CI −22.2 to 3.5). There was no difference in time-loss days until sustainable RTW between workers with and without GRTW at the 80th and 90th percentile.

Table 2

Adjusted quantile regression for differences in time-loss days until full sustainable RTW 1 year after the first time-loss day*

Results for the effect of GRTW were larger when sustainable GRTW was included in the outcome along with sustainable RTW (table 2). Between the 10th and 50th percentile of the distribution, workers with GRTW had on average 10 more disability days than workers without GRTW. However, workers with GRTW had fewer disability days at the 60th percentile (β −21.5 days; 95% CI –28.7 to –14.4 days), 70th (β −85.7; 95% CI –94.5 to –76.9) and 80th percentile (β −82.0; 95% CI –88.8 to –75.3). No differences in disability duration were found at the 90th percentile due to the censoring effect at 365 days.

A graphical representation of these results is presented in online supplemental appendix 8.

Results were also comparable for the matched cohort of workers followed for 2 years (online supplemental appendix 6), but the effect sizes were larger given the longer disability window. Among the cohort comprised of 23 318 workers, beneficial effects for GRTW on sustainable RTW were observed at the 60th (β −22.5 days; 95% CI –41.6 to –3.5) and 70th percentile (β −139.5; 95% CI –166.8 to −111.9); and on sustainable RTW or GRTW between the 50th and 80th percentiles (eg, β at 70th percentile −374.7 days; 95% CI –393.6 to –355.9).

Difference in time-loss days between workers with and without GRTW, by gender and industry sector

Differences emerged in the pattern of disability duration by GRTW status were interacted by gender (table 3). Women with GRTW had fewer time-loss days at the 60th and 70th percentile of the disability distribution compared with women without GRTW and no difference after that, while men with GRTW had more time-loss days until the 80th percentile compared with men without GRTW.

Table 3

Adjusted quantile regression for differences in time-loss days until full sustainable RTW 1 year after the first time-loss day,* by gender

The largest differences in the effect of GRTW appear when interacted with injury type, for torso sprains, upper extremity fractures and rheumatism (table 4; other injury types presented in online supplemental appendix 4). For torso sprains, workers with GRTW had fewer time-loss days at the 60th (β −18.9; 95% CI –33.7 to –4.1) and the 70th percentile (β −82.2; 95% CI –101.9 to –62.3). Censoring prevented identifying differences at the 80th and 90th percentiles. For upper extremity fractures, workers with GRTW had more time-loss days at every percentile of the disability duration; and for rheumatism claims, workers with GRTW had smaller but beneficial effects at the 60th (β −13.8; 95% CI −49.7 to 22.2) and 70th percentile (β −17.6; 95% CI −65.8 to 30.6).

Table 4

Adjusted quantile regression for differences in time-loss days until full sustainable RTW 1 year after the first time-loss day,* by MSD

By industry sector (table 5), being provided with GRTW showed beneficial effects in the trade and manufacturing sectors in the 60th and 70th percentile. No differences were observed in the services or construction sectors, or among the remaining other smaller industry sectors (online supplemental appendix 5).

Table 5

Adjusted quantile regression for differences in time-loss days until full sustainable RTW 1 year after the first time-loss day,* by industry sector

Other comparisons

The matched cohort of workers who were followed-up for 2 years (online supplemental appendix 6) comprised 23 318 workers and showed large beneficial effects for GRTW at the 60th (β −22.5; 95% CI –41.6 to –3.5) and 70th percentile (β −139.5; 95% CI –166.8 to −111.87). Larger effects were shown between the 50th and 80th percentile in favour of GRTW when sustainable GRTW is considered a positive outcome (online supplemental appendix 6).

Online supplemental appendix 7 provides the results for the matched cohort of workers limited to employers who offered GRTW at least once between 2010 and 2015 (n=19 568 workers). The overall conclusions were similar, but the effect sizes were attenuated, as expected, since the contrast of workers with and without GRTW was decreased with the removal of workers at employers who never offered GRTW.

Discussion

Main results

The objective was to investigate if GRTW facilitates sustainable RTW for workers with a work-acquired MSD, who are off work for at least 30 days, between 2010 and 2015 in BC, Canada. First, there were differences in work disability duration by GRTW status over the longer term, even after accounting for a wide range of sociodemographic, economic and injury predictors. Workers with GRTW had fewer disability (time-loss) days at the 60th and 70th percentile of the disability duration. This is an important finding for work disability management of the benefits of GRTW beyond typical work disability windows for most MSDs. Second, workers with GRTW had fewer disability days when both sustained RTW or GRTW were considered a positive outcome. It is important for work disability management systems to consider sustained work accommodations for workers with long duration claims and to reflect the real-world application of GRTW as an intermediate state in the RTW process. However, there is a group of workers who have significant injuries and permanent disability who may not be able to RTW to their previous job, and for who GRTW is not effective. These workers require vocational rehabilitation and, in some cases, will not RTW. Third, the effect of GRTW differed by MSD type, with the greatest differences observed for soft-tissue disorders, which may be more likely to resolve after initial treatment and compared with chronic inflammatory disorders.

Another notable finding is the larger effect of GRTW for women. This effect was likely attributable to the lower likelihood of sustainable RTW for women who were not provided with GRTW over the longest disability durations, as shown in other studies.37–39 By industry, larger effects were found for the trade and manufacturing sectors, sectors that are more likely to offer GRTW to their injured workers.27 These sectors generally have better access to resources and more possibilities to offer modified hours or duties compared with other industry sectors.

Finally, the importance of workplace factors being most predictive for providing GRTW speaks to the need to provide supports and resources to workplaces to improve workplace practices in accommodating workers while they recover.

Strengths and limitations

Major strengths of this study are the population-based sample that represents all compensated lost-time MSDs in BC from 2010 to 2015, as well as the use of longitudinal and comprehensive health data from the pre-injury and post-injury period. We also implemented a robust analytical approach that minimised bias due to differences in the distribution of covariates between workers with and without GRTW. CEM improved the balance in the distribution of covariates to achieve robust inferences.

Regardless of the strengths, reliance on administrative data may be subject to misclassification and information bias. However, derivation of variables from the administrative data was detailed, and any bias is hypothesised to be non-differential with a conservative effect on the observed findings. Reliance on administrative data did not allow us to get a more granular understanding of the type of GRTW, and this study provides the reader with a general effect of GRTW on a jurisdictional level, compared with specific practices.

Restricting the analyses to claims with more than 30 days of wage replacement may lead to collider-stratification bias.40 However, the construction of the cohort was done to minimise selection bias, including by thorough variable construction, transparency of exclusions, selecting workers who were on work disability for at least 30 days and might benefit from GRTW, and the use of a matching strategy to create comparable cohorts with and without GRTW.

Implications for research and practice

This study provides evidence for policymakers and occupational health professionals concerned with work disability management and the effect of GRTW in facilitating both sustainable RTW or sustainable GRTW. GRTW was effective for MSD claims with a longer duration representing moderate to severe injuries and claims with sustainable GRTW as a positive outcome, with larger effects observed for women, workers with soft-tissue injuries and for major industrial sectors. A focus on increasing the use of GRTW for long-term claims among employers in BC, and similar jurisdictions, could reduce work disability duration for those off work for 30 days.

Future research would benefit from addressing the issue of more refined or detailed measures of GRTW and sustained RTW outcomes. For example, we were not able to identify the type of GRTW that was most effective.

Conclusion

We found a positive effect of GRTW on sustainable RTW at more than 6 months of disability duration, among seriously injured workers with a work-acquired MSD in BC, Canada. We recommend endorsement of GRTW for injured workers by employers and workers’ compensation systems to reduce lost time for MSDs, including additional resources and supports to help workplaces to offer accommodations and the strengthening of jurisdictional policies to support this as a standard component of disability management.

Data availability statement

Data may be obtained from a third party and are not publicly available. The BC Freedom of Informant and Protection of Privacy Act [RSBC 1996] CHAPTER 165 (the legislation that authorises us to use the data for our study) does not authorise this data be access or view outside of Canada. Under our obligations under this legislation and under an information sharing agreement with WorksafeBC to use the data, we are prohibited from sharing the data outside of Canada and outside our research team. If other researchers wanted access to the data, they would need to negotiate access directly with the data provider, WorkSafeBC.The data for this project can only be accessed in Canada and making an application through https://www.popdata.bc.ca/dataaccess to WorkSafeBC and by entering into an Research Agreement with WorkSafeBC.

Ethics statements

Patient consent for publication

Ethics approval

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

Acknowledgments

Esther Maas is supported by a Canadian Institutes of Health Research (CIHR) postdoctoral fellowship. Mieke Koehoorn was supported in part by a CIHR Chair in Gender, Work and Health. Christopher McLeod is supported by a CIHR New Investigator Award and a Michael Smith Foundation for Health Research Scholar Award. All authors are supported with research operating funds from WorkSafeBC through the Partnership for Work, Health and Safety. All inferences, opinions and conclusions drawn in this manuscript are those of the authors, and do not reflect the opinions or policies of the Data Steward(s).

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors All authors were involved in the conceptualisation of the project. ETM and CBM were involved in the analyses. ETM drafted the manuscript. All authors edited the manuscript.

  • Funding This study is funded by a Canadian Institutes of Health Research postdoctoral fellowship (MRE-152431).

  • Disclaimer All inferences, opinions and conclusions drawn in this dissertation are those of the author, and do not reflect the opinions or policies of the Data Stewards.

  • Competing interests None declared.

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

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.