Skip to main content

Income inequalities in multimorbidity prevalence in Ontario, Canada: a decomposition analysis of linked survey and health administrative data



The burden of multimorbidity is a growing clinical and health system problem that is known to be associated with socioeconomic status, yet our understanding of the underlying determinants of inequalities in multimorbidity and longitudinal trends in measured disparities remains limited.


We included all adult respondents from four cycles of the Canadian Community Health Survey (CCHS) (between 2005 to 2011/12), linked at the individual-level to health administrative data in Ontario, Canada (pooled n = 113,627). Multimorbidity was defined at each survey response as having ≥2 (of 17) high impact chronic conditions, based on claims data. Using a decomposition method of the Erreygers-corrected concentration index (CErreygers), we measured household income inequality and the contribution of the key determinants of multimorbidity (including socio-demographic, socio-economic, lifestyle and health system factors) to these disparities. Differences over time are described. We tested for statistically significant changes to measured inequality using the slope index (SII) and relative index of inequality (RII) with a 2-way interaction on pooled data.


Multimorbidity prevalence in 2011/12 was 33.5% and the CErreygers was − 0.085 (CI: -0.108 to − 0.062), indicating a greater prevalence among lower income groups. In decomposition analyses, income itself accounted more than two-thirds (69%) of this inequality. Age (21.7%), marital status (15.2%) and physical inactivity (10.9%) followed, and the contribution of these factors increased from baseline (2005 CCHS survey) with the exception of age. Other lifestyle factors, including heavy smoking and obesity, had minimal contribution to measured inequality (1.8 and 0.4% respectively). Tests for trends (SII/RII) across pooled survey data were not statistically significant (p = 0.443 and 0.405, respectively), indicating no change in inequalities in multimorbidity prevalence over the study period.


A pro-rich income gap in multimorbidity has persisted in Ontario from 2005 to 2011/12. These empirical findings suggest that to advance equality in multimorbidity prevalence, policymakers should target chronic disease prevention and control strategies focused on older adults, non-married persons and those that are physically inactive, in addition to addressing income disparities directly.


An increasing number of adults in high-income settings have been diagnosed with multiple coexisting chronic conditions [1], also known as multimorbidity. This increased burden has created significant challenges in the effective provision of clinical care and has added pressures on health systems that traditionally provide highly specialized care for a single condition. Extensive research has shown that persons with multimorbidity have a greater risk of functional decline, cognitive decline and early mortality [2,3,4], are more frequent users of healthcare services, experience greater care fragmentation and longer hospital stays, and also incur higher healthcare costs [5,6,7,8,9,10,11,12].

Multiple studies from high-income settings have shown that multimorbidity is concentrated among persons from lower socio-economic ranks [13,14,15,16,17]. Multimorbidity may develop 10–15 years earlier among young and middle-aged persons in the most (vs. least) deprived areas [18]. While quantifying multimorbidity inequalities may help inform where policy action is needed, less is known about the relative contributions of determinants to socioeconomic disparities in multimorbidity [19]. These data are particularly useful to policymakers for informing future resource allocation and designing targeted interventions. Moreover, few studies have quantified longitudinal trends in measured inequalities in multimorbidity prevalence by income, education, or occupation, despite calls for action by federal agencies and professional organizations in many jurisdictions to reduce health disparities [20, 21].

To address these gaps, we undertook a comprehensive examination of inequalities in multimorbidity prevalence in Ontario, Canada. Our key research objectives were to: 1) quantify household income inequalities in multimorbidity prevalence among adults in Ontario, 2) identify the relative contribution of key determinants of multimorbidity to this inequality including socio-demographic, socio-economic, lifestyle and health system factors; and 3) assess whether these disparities were widening, and whether the drivers contributing to this inequality were changing over time.


We used (cross-sectional) survey data and linked health administrative information that is routinely collected in Ontario. The use of this data was authorized under section 45 of Ontario’s Personal Health Information Protection Act, which does not require review by a Research Ethics Board. The study is reported according to the RECORD guidelines [22].

Data and setting

Residents of Ontario, Canada aged 18 years and older that participated in any of four cycles of the Canadian Community Health Survey (CCHS) – 2005, 2007/08, 2009/10 and 2011/12 – and whose responses were linked to population-based health administrative databases were included in the study. Each survey was analyzed separately, unless otherwise stated (pooled cross-sections for longitudinal analysis).

The CCHS is administered by Statistics Canada to Canadians aged ≥12 years living in private dwellings, and is representative of 98% of the Canadian population. The CCHS does not include into its sampling frame persons living on Reserves or Crown Lands, institutionalized residents or members of the Canadian Forces. Detailed methodology of the survey and sampling strategy are described elsewhere [23, 24]. Permanent residents of Ontario (2011 population 12.8 million) are covered by a universal health insurance program that covers the costs for most physician and hospital services. Immigrants receive services after a three-month waiting period. Patient encounters with the healthcare system are recorded in administrative databases. Survey and administrative data were linked deterministically at the individual level using unique encoded identifiers and analyzed at the Institute for Clinical Evaluative Sciences (ICES) in Toronto, Ontario. The databases utilized in this study included the Ontario Health Insurance Program claims database (OHIP), Discharge Abstract Database (DAD), Ontario Drug Benefits database (ODB), Registered Persons Database (RPDB), and the Client Agency Program Enrolment database (CAPE) (see:


We used health administrative data (OHIP, DAD, and ODB data) to determine the prevalence of 17 high-impact chronic conditions for each individual at time of their survey response. Consistent with previous studies of multimorbidity in Ontario [1, 12, 14, 25,26,27], these conditions were selected based on their economic impact and population burden in the general population. The 17 conditions included: acute myocardial infarction, asthma, cancers, cardiac arrhythmia, chronic coronary syndrome, chronic obstructive pulmonary disorder, congestive heart failure, diabetes, hypertension, mood and anxiety disorders, other mental illnesses (schizophrenia, delusions and other psychoses, personality disorders and substance abuse), osteoarthritis, osteoporosis, renal failure, rheumatoid arthritis and stroke (excluding transient ischemic attack). Where available, validated algorithms were applied to the health administrative data to ascertain cases, including AMI, asthma, COPD, CHD, dementia, diabetes and hypertension [28,29,30,31,32,33,34]. As in most other studies [18, 19, 35] we defined multimorbidity from this data as the co-occurrence of 2 or more (of these 17) conditions within the same individual, prevalent at the time of survey response.

We used household income from CCHS responses to measure socio-economic status (all respondents are asked: “what is your best estimate of the total income received by all household members, from all sources, before taxes and deductions, in the past 12 months?”). We focused specifically on (household) income as a marker of socioeconomic status to measure inequality because the rich-poor income gap is widening in most countries [36, 37], the strong association relating income to population health is well documented [38], and these gradients are potentially amenable to policy. Because of the sensitive nature of income reporting (resulting in missing or misclassified data), we used the imputed income responses provided by Statistics Canada. This method has been described extensively elsewhere [39]. These data were aggregated to quintiles for analysis and reporting.

Independent variables were selected a priori based on previous literature reporting associations with multimorbidity prevalence from similar high-income jurisdictions [13, 18, 40, 41]. Socio-demographic variables included age group (18–34 years, 35–49, 50–64, 65–74 and 75-plus), sex (men and women), rurality (urban, suburban or rural residence - based on the Rurality Index of Ontario [42], from RPDB data), marital status (married vs. divorced, separated, widowed or single), and immigrant status (born in Canada vs. not born in Canada). In addition to income, socio-economic variables (from CCHS data) included the respondent’s education level (no post-secondary vs. some post-secondary or higher). Health-related lifestyle characteristics included level of physical activity (active, moderately active, inactive), smoking status (heavy smoker, light smoker, former smoker, non-smoker) and body-mass index (BMI, underweight, normal weight, overweight, obese, following the classifications defined by the World Health Organization [43]). Health system variables included health region of the respondent’s residence (i.e., Local Health Integration Network [LHIN]) and enrolment in a primary model of care, which was used a proxy for access to primary care services (Family Health Group [FHG], Family Health Network [FHN], Family Health Organization [FHO], other, and not enrolled [CAPE data]) [44, 45]. Values of survey variables were selected based on the distribution of responses and to reflect meaningful categories. We excluded individuals from analysis if they were missing information (non-response) for any of these variables because there are currently no methods for combining the multiple imputation estimates specific to decomposition analyses.

Statistical analyses

For each cycle, we transformed income data into cumulative rank probabilities (ridit scores) ranging from 0 (highest income) to 1 (lowest income), with values reflecting the midpoint of the cumulative proportion of the (weighted) population in each income group. Using this variable we quantified inequalities in multimorbidity for each CCHS cycle, first, using the Concentration Index (C) [46]. This relative measure of inequality considers the entire distribution of income, and can be written as:

$$ C=\frac{2}{n\mu}\sum \limits_{i=1}^n{y}_i{R}_i-1 $$

where C is the concentration index, μ is the (weighted) population mean of the outcome (here, multimorbidity prevalence), y is the outcome mean of the ith individual, and R i is the individuals rank in the income distribution. Values range from −1 to + 1, corresponding to a concentration of the outcome (multimorbidity prevalence) among the poorer or wealthier population, respectively, where larger values (approaching −1 or + 1) indicate greater inequality, and a value of zero indicates equality. Because multimorbidity is dichotomous, we applied the Erreygers correction (CErreygers), which multiplies C by 4 times the weighted mean of the health outcome [47]. Standard errors and 95% confidence intervals (CI) for CErreygers were derived using the methods described by O’Donnell et al. [48].

Wagstaff et al. [49] showed that C of a continuous health outcome can be decomposed into a set of determinants, a methodology that has since been extended to dichotomous variables [48, 50, 51]. In comparison to traditional analyses, this decomposition method allows for the explanation of the measured health inequality across the entire distribution of socio-economic status. To perform the decomposition, we specified a probit model with marginal effects [48], including all independent variables previously listed in the regression. Income itself is included as a determinant in the decomposition regression to prevent overestimation of the contribution of all other factors to measured inequality [52]. For any determinant to contribute to measured inequality, that determinant must be associated with the health outcome (i.e., have a non-zero regression coefficient) and also be unequally distributed by socio-economic status (i.e., have a non-zero concentration index). A positive (negative) contribution means that that determinant is associated with greater (lower) income inequality in multimorbidity. In addition to the proportional contribution of each determinant is a residual, representing the amount inequality in multimorbidity that cannot be explained by variation of each determinant across the income distribution. Negative residual values are possible, suggesting that determinants entered into the decomposition explain all measured inequality. We present decomposition results from the most recent survey, and describe differences in relative contribution of determinants to the baseline cycle.

Multiple sensitivity analyses were performed to assess the robustness of these decomposition analyses. Analyses were repeated using different reference categories in the multivariable probit model [50], performing the decomposition using a logit model approach (i.e., decomposing the natural logarithm of the odds of multimorbidity) [51], decomposing the Wagstaff-corrected C [53], and using a continuous log of household income term in the model (rather than quintiles). For each sensitivity analysis, we assessed any changes in the rank-order of the contribution of determinants measured inequality. Additionally, to ensure the measured inequality was robust to missing data, we re-calculated the C-Erreygers for each survey, including those with missing determinant data.

Lastly, to confirm changes in inequality over the study period, we quantified the relative index of inequality (RII) and slope index of inequality (SII). These measures are recommended for assessing trends in health disparities [54, 55]. First, we estimated the RII in multimorbidity prevalence within each cycle using a Poisson regression model [56] including covariates for ridit score, age group and sex. The exponential of the ridit score estimate is equivalent to the RII, which can be interpreted as the proportion of multimorbidity prevalence that differs between the highest and lowest incomes. An RII > 1 is indicative of pro-rich inequality on the relative scale (with larger values suggesting a greater concentration of multimorbidity among the poor). To asses trends, we used the methods outlined by Ernstsen et al. [57] using pooled survey data (with survey-specific ridit scores) and included covariates for ridit score, age group, sex, CCHS cycle and the 2-way ridit*cycle interaction. A positive (negative) and statistically significant interaction term is indicative of increasing (decreasing) inequality, on the relative scale. These methods were repeated using a linear probability model to determine the SII, reflecting the absolute difference in multimorbidity prevalence between the highest and lowest incomes. An SII > 0 is indicative of pro-rich inequality on the absolute scale. We also assessed trends in absolute inequality trends using pooled data. As a sensitivity test, we confirmed the statistical significance of these analyses adjusting for all measured determinants. We applied bootstrap sampling weights using balanced repeated replication (n = 500) to analyses. All data management was conducted using SAS Enterprise Guide version 6.1 (SAS Institute, Inc.) and all analyses were performed using Stata/MP version 13.1 (StataCorp).


Of 134,395 Ontario respondents across 4 cycles of linked CCHS and administrative data, 122,002 (91% of sample) were adults, and 113,627 (93% of adults) had complete information across all variables for analysis. A higher proportion of adult respondents excluded (vs. included) from analyses were age 75+ years, female, from the lowest income quintile and had prevalent multimorbidity. Table 1 describes the characteristics of respondents included in analyses for each CCHS cycle. The numbers outside and inside the brackets correspond to sample counts and population weighted proportions, respectively. Crude multimorbidity prevalence increased steadily with each cycle, from 26.4% in the 2005 survey to 33.5% in the 2011/12 survey. Crude prevalence estimates by household income quintile are shown in Fig. 1, and by determinant in the Additional file 1: Table S2. Of note, prevalence in the lowest vs. highest household income quintiles, respectively, were 34.7% vs. 20.8% in 2005 and 37.8% vs. 27.6% in 2011/12.

Table 1 Characteristics of adult (age ≥ 18 years) CCHS survey respondents linked to health administrative data: 2005 to 2011/12
Fig. 1
figure 1

Longitudinal trends in multimorbidity by household income over four cycles of the CCHS

The measured CErreygers are shown in the last row of Table 1: values ranged from − 0.114 (CI: -0.132, − 0.096) in 2005 to − 0.085 (CI: -0.108, − 0.062) in 2011/12. The negative and statistically significant CErreygers indicates a greater concentration of multimorbidity among persons with low household income, although the smaller CErreygers value in the 2011/12 survey suggests that inequality was less pronounced.

Table 2 summarizes the results of the CErreygers decomposition of multimorbidity prevalence in the most recent (2011/12) and baseline (2005) surveys. Results from the 2007/08 and 2009/10 surveys are available in the Additional file 1: Table S3. Data include marginal effects from probit regressions, elasticities of each determinant (βk*x̅k/μ), the CErreygers of each determinant, and the relative contribution of each determinant to measured inequality. For example, for persons aged 75+ in the 2011/12 CCHS cycle: the difference in conditional probability of multimorbidity (vs. 18–34 age group) was 77.3%; this group represented 7.4% of the total population (from Table 1) but persons were concentrated in lower income levels (Ck*(4*μ) = − 0.061); and their contribution to measured income inequality was 55.6%.

Table 2 Decomposition of the Erreygers-corrected Concentration Index of multimorbidity prevalence among adults in Ontario: 2005 vs. 2011/12

Individual determinants with the largest contribution to the observed pro-rich inequality in the 2011/12 CCHS cycle included low household income (57.2 and 18.8% for quintiles 1 and 2, respectively), older age (55.6 and 29.9% for ages 75+ and 65–74 years, respectively), being divorced, separated, widowed or single (15.2%), and physical inactivity (10.8%). The relative contribution of each of these determinants to measured inequality increased from the decomposition analysis conducted using the 2005 CCHS survey data. For example, older ages (75+ years) were less concentrated among the low income (vs. 2005) but the elasticity (i.e., the % change in multimorbidity associated with a % change in the determinant) increased between cycles; for non-married persons, the elasticity increased, as did the concentration of non-married individuals among lower income groups. In contrast, variables that contributed negatively to inequality in the latter survey included being middle-aged (− 51.6% and − 12.2% for 50–64 and 35–49 years, respectively), having higher household income (− 10.5% for quintile 4), enrolment in a FHO (− 11.5%) and being born outside of Canada (− 7.8%). Of note, the growth in FHO enrolment between surveys (from 3.3 to 42.1%) corresponded with a greater concentration among higher household income (CErreygers from − 0.012 to 0.093); the opposite pattern was found for enrolment in a FHG (increased in enrolment from 23.6 to 27.7% and shift in CErreygers from 0.029 to − 0.042).

Aggregating these percentages by determinant (Fig. 2), income explained 69.0% of the measured inequality in 2011/12 (an increase of + 30.2% from 2005 CCHS survey), followed by age at 21% (− 16.5%), marital status at 15.2% (+ 8.6%), and physical activity at 10.9% (+ 8.6%). In contrast, determinants that contributed the largest to decreasing inequality were immigration at − 7.8 (− 7.2%) and primary care enrolment at − 7.4 (− 5.0%). In 2011/12, the modeled determinants explained − 0.093 of the total inequality (− 0.085), whereas in 2005 the determinants explained − 0.107 of the total inequality (− 0.114).

Fig. 2
figure 2

Summary of the relative contribution of key determinants to income inequality in multimorbidity in Ontario, 2005 to 2011/12

Table 3 shows the age-sex adjusted relative (RII) and absolute (SII) inequality measures for each survey, as well as the longitudinal trend analysis. The prevalence rate ratio (RII) in 2011/12 comparing high vs. low income quintile was 1.32 (CI: 1.19, 1.46) in 2011/12. The corresponding prevalence rate difference (SII) was 9.1% (CI: 5.7, 12.4%). In pooled regression analyses, longitudinal trends in inequality were not statistically significant over the study period on either the absolute or relative scale, suggesting no change in measured income inequality (p for interaction = 0.405 for RII and 0.443 for SII).

Table 3 Trends in income inequalities in multimorbidity prevalence among adults in Ontario, 2005 to 2011/12

These findings were robust to the defined sensitivity analyses (Additional file 1: Table S3).


This study confirms a moderate, statistically significant pro-rich household income gradient in multimorbidity prevalence in Ontario, Canada. This gap was observed on both relative and absolute scales, and has persisted from 2005 to 2011/12. Our results are consistent with measured wealth and socioeconomic inequalities in chronic conditions among adults observed in low and middle-income countries [19, 56, 58]. These studies considered only a small number of conditions in defining chronic disease burden. Our findings are also consistent with a study of multimorbidity (among 52 conditions) in deprived areas in the Basque Country of Spain [59]. Although socioeconomic differences in common chronic conditions among adult populations have been previously reported [60, 61], few studies have considered an aggregated count of multiple high impact chronic conditions or measured inequality with robust methods that account for the entire socio-economic distribution.

These findings highlight the unequal distribution of multimorbidity across socio-economic ranks. An important contribution of the current study to existing literature is the measurement of the relative contributions of multimorbidity’s determinants to measured income inequalities from a high-income setting. The decomposition method is useful for policy makers who seek to reduce health inequalities by targeting determinants that contribute the most to observed inequalities. Recently, Kunna et al. [19] decomposed wealth-related inequalities in multimorbidity prevalence (defined as 2 or more of 7 chronic conditions) among adults aged 50 years and older in China (middle-income setting) and in Ghana (low-income setting). Considering a similar set of measured determinants, the authors found that largest contributors to inequality were wealth, age, and education in China and BMI, wealth and rurality in Ghana. We found that the largest contributor to household income inequality in multimorbidity prevalence was income itself. This finding has been observed for similar decomposition studies of income-related inequalities in general health status from North America [62] and Europe [63]. Our data also show that the proportion of inequality in multimorbidity due to income increased over time. The growth in the rich-poor income gap in Canada, like may comparable OECD nations, has been well documented [37].

The strong negative association between increasing age and multimorbidity prevalence we observed is well established in existing literature [16]. The positive contribution of age to measured inequality in our study decreased steadily between cycles, mainly due to changes among the 50–64 year age group, who became more concentrated among higher household incomes. Among other socio-demographic variables, the contribution of marital status to measured inequality increased. Similar to age, where the change to the elasticity was minor, the increase in contribution was driven by a greater concentration of divorced, separated, widowed or single individuals in lower household income levels.

Physical activity was the strongest lifestyle factor contributing to observed inequalities in multimorbidity. We found that prevalence of multimorbidity increased over the study period (from 29.4 to 39.5%) among inactive persons in Ontario, and that inactivity was increasingly concentrated among lower incomes. Socioeconomic gradients in inactivity have been documented in other high-income jurisdictions [64], although less is known regarding changes across the socioeconomic distribution over time. Body-mass index had only a small negative contribution to inequality overall, and obesity in particular was relatively equally distributed across income. This contrasts findings by Hajizahen et al. [65], who report a small but significant concentration of obesity among the poor in Canada. Their data, however, exclude older adults (aged 65+ years). Our findings show that heavy and former smoking were unequally distributed across income but had small (although statistically significant) regression effects and small corresponding elasticities, also resulting in a minor contribution to inequality. The variability in the contribution of lifestyle determinants to inequalities in multimorbidity from previous work from low and middle-income jurisdictions [19] suggest that interventions meant to improve equality need to be targeted to local contexts. Other unhealthy lifestyle factors that are potentially amenable to public health intervention, such as healthy diet and alcohol consumption, were not included in our multivariable decomposition analyses since their relationships with multimorbidity prevalence are less certain [40].

The contribution of health system variables to overall inequality, including health region of residence and enrolment in a primary care model, were small. All participants of the study had universal access to medically necessary care; however, enrolment into a primary care model is voluntary. The negative contributions for this determinant in particular suggest that holding all else constant, persons enrolled in these models (particularly FHNs) are less likely to be multimorbid and of higher household income, which has been previously documented [44, 45]. Lastly, across all decomposition analyses, being born outside of Canada was negatively associated with multimorbidity prevalence (to varying effect sizes), and negatively contributed to multimorbidity inequality. A healthy immigrant effect for chronic disease incidence has been previously documented in Canada, although with diminishing associations as immigrants’ time in Canada lengthened [66, 67]. The effect of time since immigrating on socioeconomic inequalities in multimorbidity (and other health outcomes) is an avenue for future research.

The results from the decomposition were robust to multiple sensitivity analyses. Other strengths of the research include the use of a large sample size, survey weights that are reflective of Ontario’s population, and high survey response rates. Our results are likely to be generalizable to other high-income jurisdictions with comparable universal health care systems. Validated algorithms were also used to ascertain disease status; linked survey and administrative data enabled this inclusion. There are limitations of this research that are worth noting, however. Variables were defined using routinely-collected health administrative and survey data which are subject to known limitations [22, 23]. For health administrative data in particular, ascertainment of chronic disease is dependent having a place of residence and on health system contact (physician visits or hospital admissions), which are less likely among vulnerable, low-SES individuals. This would underreport multimorbidity prevalence and the overall measured inequality. Each survey is cross-sectional by design, and therefore, causality cannot be inferred from the data or by the decompositions. Inequality trends were assessed over a 7-year period from biennial cross-sectional surveys; the inclusion of additional surveys could be an avenue for future research, once linked to administrative data. Importantly, interpretation of decomposition results rely on a correct model specification; Pseudo R-squared values from the multivariable probit models were 0.22 (2005 CCHS) and 0.23 (2011/12 CCHS) suggesting the variables included in the regression explain only a portion of the variance observed in the outcome. Lastly, we were only able to include program enrolment as a marker of primary care access. However, a high-performing primary care sector can support equitable health outcomes by promoting seamless access to services across the care continuum, including preventive care and ongoing chronic disease management [68,69,70]. As such, other access-related determinants not measurable with our data may contribute to income inequalities in multimorbidity prevalence. Additional research that includes primary care capacities could help inform policymakers, and play a pivotal role for reducing overall disparities.


This study confirms a pro-rich income gap in multimorbidity prevalence that has persisted over time. This gap is driven primarily by income, increasing age, marital status (not being married or common-law) and physical inactivity. Additional research to uncover the root causes of co-existing chronic conditions and to develop and evaluate the effectiveness of prevention strategies are warranted to Acknowledgements.

We thank IMS Brogan Inc. for use of their Drug Information Database.



Body mass index


Concentration index


Client Agency Program Enrolment database


Canadian Community Health Survey

CErreygers :

Erreygers-corrected concentration index


Discharge Abstract Database


Institute for Clinical Evaluative Sciences


Local Health Integration Network


Ontario Drug Benefit database


Organisation for Economic Co-operation and Development


Ontario Health Insurance Program claims database


Relative index of inequality


Slope index of inequality


  1. Koné Pefoyo AJ, Bronskill SE, Gruneir A, Calzavara A, Thavorn K, Petrosyan Y, et al. The increasing burden and complexity of multimorbidity. BMC Public Health. 2015;15:415. BioMed Central Ltd

    Article  PubMed Central  Google Scholar 

  2. Rijken M, van Kerkhof M, Dekker J, Schellevis FG. Comorbidity of chronic diseases: effects of disease pairs on physical and mental functioning. Qual Life Res. 2005;14:45–55.

    Article  PubMed  Google Scholar 

  3. Tinetti ME, McAvay GJ, Chang SS, Newman AB, Fitzpatrick AL, Fried TR, et al. Contribution of multiple chronic conditions to universal health outcomes. J Am Geriatr Soc. 2011;59:1686–91.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Nunes BP, Flores TR, Mielke GI, Thumé E, Facchini LA. Multimorbidity and mortality in older adults: a systematic review and meta-analysis. Arch Gerontol Geriatr. 2016;67:130–8. Elsevier Ireland Ltd

    Article  PubMed  Google Scholar 

  5. Bähler C, Huber CA, Brüngger B, Reich O. Multimorbidity, health care utilization and costs in an elderly community-dwelling population: a claims data based observational study. BMC Health Serv Res. 2015;15:23. BioMed Central Ltd

    Article  PubMed  PubMed Central  Google Scholar 

  6. Glynn LG, Valderas JM, Healy P, Burke E, Newell J, Gillespie P, et al. The prevalence of multimorbidity in primary care and its effect on health care utilization and cost. Fam Pract. 2011;28:516–23.

    Article  PubMed  Google Scholar 

  7. Wolff JL, Starfield B, Anderson G. Prevalence, expenditures, and complications of multiple chronic conditions in the elderly. Arch Intern Med. 2002;162:2269–76.

    Article  PubMed  Google Scholar 

  8. Smith SM, Soubhi H, Fortin M, Hudon C, O'Dowd T. Managing patients with multimorbidity: systematic review of interventions in primary care and community settings. BMJ. 2012;345:e5205.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Kuo RN, Lai M-S. The influence of socio-economic status and multimorbidity patterns on healthcare costs: a six-year follow-up under a universal healthcare system. Int J Equity Health. 2013;12:69.

    Article  PubMed  PubMed Central  Google Scholar 

  10. van Oostrom SH, Picavet HSJ, de Bruin SR, Stirbu I, Korevaar JC, Schellevis FG, et al. Multimorbidity of chronic diseases and health care utilization in general practice. BMC Fam Pract. 2014;15:61. BioMed Central

    Article  PubMed  PubMed Central  Google Scholar 

  11. Zulman DM, Pal Chee C, Wagner TH, Yoon J, Cohen DM, Holmes TH, et al. Multimorbidity and healthcare utilisation among high-cost patients in the US veterans affairs health care system. BMJ Open. 2015;5:e007771. British Medical Journal Publishing Group

    Article  PubMed  PubMed Central  Google Scholar 

  12. Gruneir A, Bronskill SE, Maxwell CJ, Bai YQ, Kone AJ, Thavorn K, et al. The association between multimorbidity and hospitalization is modified by individual demographics and physician continuity of care: a retrospective cohort study. BMC Health Serv Res. 2016;16:1–9. 2nd ed. BMC Health Serv Res

    Article  Google Scholar 

  13. Roberts KC, Rao DP, Bennett TL, Loukine L, Jayaraman GC. Prevalence and patterns of chronic disease multimorbidity and associated determinants in Canada. Health Promot Chronic Dis Prev Can. 2015;35:87–94.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Lane NE, Maxwell CJ, Gruneir A, Bronskill SE, Wodchis WP. Absence of a socioeconomic gradient in older Adults’ survival with multiple chronic conditions. EBioMedicine. 2015;2:2094–100.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Salisbury C, Johnson L, Purdy S, Valderas JM, Montgomery AA. Epidemiology and impact of multimorbidity in primary care: a retrospective cohort study. Br J Gen Pract. 2011;61:12–21.

    Article  Google Scholar 

  16. Violán C, Foguet-Boreu Q, Roso-Llorach A, Rodriguez-Blanco T, Pons-Vigués M, Pujol-Ribera E, et al. Burden of multimorbidity, socioeconomic status and use of health services across stages of life in urban areas: a cross-sectional study. BMC Public Health. 2014;14:530–13. Available from:

    Article  PubMed  PubMed Central  Google Scholar 

  17. Mercer SW, Watt GCM. The inverse care law: clinical primary care encounters in deprived and affluent areas of Scotland. Ann Fam Med. 2007;5:503–10.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Barnett K, Mercer SW, Norbury M, Watt G, Wyke S, Guthrie B. Epidemiology of multimorbidity and implications for health care, research, and medical education: a cross-sectional study. Lancet. 2012;380:37–43.

    Article  PubMed  Google Scholar 

  19. Kunna R, Sebastian MS, Williams JS. Measurement and decomposition of socioeconomic inequality in single and multimorbidity in older adults in China and Ghana: results from the WHO study on global AGEing and adult health (SAGE). Int J Equity Health. 2017;16(1):79.

  20. Public Health Agency of Canada. Reducing Health Inequalities: A Challenge For Our Times. 2011. Available from:

    Google Scholar 

  21. Department of Health. Tackling health inequalities: a programme for action. 2003. Available from:

    Google Scholar 

  22. Benchimol EI, Smeeth L, Guttmann A, Harron K, Moher D, Petersen I, et al. The REporting of studies conducted using observational routinely-collected health data (RECORD) statement. PLoS Med. 2015;12:e1001885–22. Public Library of Science.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Béland Y. Canadian community health survey--methodological overview. Health Rep. 2002;13:9–14.

    PubMed  Google Scholar 

  24. Statistics Canada. Canadian Community Dent Health Survey - Annual Component (CCHS). 2018 [cited June 19 2018]. Available from:

  25. Mondor L, Maxwell CJ, Bronskill SE, Gruneir A, Wodchis WP. The relative impact of chronic conditions and multimorbidity on health-related quality of life in Ontario long-stay home care clients. Qual Life Res. 2016;25:2619–32.

    Article  PubMed  Google Scholar 

  26. Mondor L, Maxwell CJ, Hogan DB, Bronskill SE, Gruneir A, Lane NE, et al. Multimorbidity and healthcare utilization among home care clients with dementia in Ontario, Canada: a retrospective analysis of a population-based cohort. PLoS Med. 2017;14:e1002249. Brayne C, editor Public Library of Science

    Article  PubMed  PubMed Central  Google Scholar 

  27. Thavorn K, Maxwell CJ, Gruneir A, Bronskill SE, Bai Y, Koné Pefoyo AJ, et al. Effect of socio-demographic factors on the association between multimorbidity and healthcare costs: a population-based, retrospective cohort study. BMJ Open. 2017;7:e017264–15. British Medical Journal Publishing Group

    Article  PubMed  PubMed Central  Google Scholar 

  28. Austin PC, Daly PA, Tu JV. A multicenter study of the coding accuracy of hospital discharge administrative data for patients admitted to cardiac care units in Ontario. Am Heart J. 2002;144:290–6. Elsevier

    Article  PubMed  Google Scholar 

  29. Gershon AS, Wang C, Guan J, Vasilevska-Ristovska J, Cicutto L, To T. Identifying patients with physician-diagnosed asthma in health administrative databases. Can Respir J. 2009;16:183–8.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Gershon AS, Wang C, Guan J, Vasilevska-Ristovska J, Cicutto L, To T. Identifying individuals with Physcian diagnosed COPD in health administrative databases. COPD. 2009;6:388–94.

    Article  PubMed  CAS  Google Scholar 

  31. Schultz SE, Rothwell DM, Chen Z, Tu K. Identifying cases of congestive heart failure from administrative data: a validation study using primary care patient records. Chronic Dis Inj Can. 2013;33:160–6.

    PubMed  CAS  Google Scholar 

  32. Jaakkimainen RL, Bronskill SE, Tierney MC, Herrmann N, Green D, Young J, et al. Identification of physician-diagnosed Alzheimer“s disease and related dementias in population-based administrative data: a validation study using family physicians” electronic medical records. J Alzheimers Dis. 2016;54:337–49. IOS Press

    Article  PubMed  Google Scholar 

  33. Hux JE, Ivis F, Flintoft V, Bica A. Diabetes in Ontario: determination of prevalence and incidence using a validated administrative data algorithm. Diabetes Care. 2002;25:512–6.

    Article  PubMed  Google Scholar 

  34. Tu K, Campbell NR, Chen Z-L, Cauch-Dudek KJ, McAlister FA. Accuracy of administrative databases in identifying patients with hypertension. Open Med. 2007;1:e18–26.

    PubMed  PubMed Central  Google Scholar 

  35. Diederichs C, Berger K, Bartels DB. The measurement of multiple chronic diseases--a systematic review on existing multimorbidity indices. J Gerontol. 2011;66:301–11.

    Article  Google Scholar 

  36. Organization for Economic Co-operation and Development. Growing Income Inequality in OECD Countries: What drives it and how can policy tackle it? OECD forum on tackling inequality: Paris, 2011.

    Google Scholar 

  37. Organization for Economic Co-operation and Development. Divided We Stand: Why Inequality Keeps Rising: OECD Publishing; 2012. Available from:

  38. Pickett KE, Wilkinson RG. Income inequality and health: a causal review. Soc Sci Med. 2015;128:316–26. Elsevier Ltd

    Article  PubMed  Google Scholar 

  39. Yeung CW, Thomas S. Income Imputation for the Canadian Community Health Survey. Canada: Statistics Canada, Methodology Branch; 2012.

  40. Fortin M, Haggerty J, Almirall J, Bouhali T, Sasseville M, Lemieux M. Lifestyle factors and multimorbidity: a cross sectional study. BMC Public Health. 2014;14:553–8.

    Article  Google Scholar 

  41. Wikström K, Lindström J, Harald K, Peltonen M, Laatikainen T. Clinical and lifestyle-related risk factors for incident multimorbidity: 10-year follow-up of Finnish population-based cohorts 1982-2012. Eur J Intern Med. 2015;26:211–6.

    Article  PubMed  Google Scholar 

  42. Kralj B. Measuring Rurality - RIO2008 BASIC: Methodology and Results. Toronto; 2009. [cited 2016 Sep]. Accessed June 2018.

  43. World Health Organization. Obesity: Preventing and managing the global epidemic - Report of a WHO consultation. Geneva; World Health Organization; 2000. p. 1–268.

  44. Glazier RH, Kopp A, Schultz SE, Kiran T, Henry DA. All the right intentions but few of the desired results: lessons on access to primary care from Ontario’s patient enrolment models. Healthcare Q. 2012;15:17–21.

    Article  Google Scholar 

  45. Glazier RH, Zagorski BM, Rayner J. Comparison of Primary Care Models in Ontario. ICES Investigative Report. Toronto; 2012. p. 1–41. Available from:

  46. Wagstaff A, Paci P, Van Doorslaer E. On the measurement of inequalities in health. Soc Sci Med. 1991;33:545–57.

    Article  PubMed  CAS  Google Scholar 

  47. Erreygers G. Correcting the concentration index. J Health Econ. 2009;28:504–15.

    Article  PubMed  Google Scholar 

  48. van Doorslaer E, Wagstaff A, Lindelow M. Analyzing health equity using household survey data: A Guide to Techniques and Their Implementation. Washington, D.C; 2008. p. 1–234. Available from:

  49. Wagstaff A, van Doorslaer E, Watanabe N. On decomposing the causes of health sector inequalities with an application to malnutrition inequalities in Vietnam. J Econ. 2003;112:207–23.

    Article  Google Scholar 

  50. Yiengprugsawan V, Lim LL, Carmichael GA, Dear KB, Sleigh AC. Decomposing socioeconomic inequality for binary health outcomes: an improved estimation that does not vary by choice of reference group. BMC Res Notes. 2010;3:57–5.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Hosseinpoor AR, van Doorslaer E, Speybroeck N, Naghavi M, Mohammad K, Majdzadeh R, et al. Decomposing socioeconomic inequality in infant mortality in Iran. Int J Epidemiol. 2006;35:1211–9.

    Article  PubMed  Google Scholar 

  52. Wagstaff A, Bilger M, Sajaia Z, Lakshin M. Health equity and financial protection. Washington, D.C; 2011. Available from:

  53. Wagstaff A. The bounds of the concentration index when the variable of interest is binary, with an application to immunization inequality. Health Econ. 2005;14:429–32.

    Article  PubMed  Google Scholar 

  54. Mackenbach JP, Kunst AE. Measuring the magnitude of socio-economic inequalities in health: an overview of available measures illustrated with two examples from Europe. Soc Sci Med. 2003;44:757–71.

    Article  Google Scholar 

  55. Keppel K, Pamuk E, Lynch J, Carter-Pokras O, Kim Insun, Mays V, et al. Methodological issues in measuring health disparities. Vital Health Stat 2. 2005;141:1–16.

    Google Scholar 

  56. Hosseinpoor AR, Bergen N, Mendis S, Harper S, Verdes E, Kunst A, et al. Socioeconomic inequality in the prevalence of noncommunicable diseases in low- and middle-income countries: results from the world health survey. BMC Public Health. 2012;12:474. BioMed Central

    Article  PubMed  PubMed Central  Google Scholar 

  57. Ernstsen L, Strand BH, Nilsen SM, Espnes GA, Krokstad S. Trends in absolute and relative educational inequalities in four modifiable ischaemic heart disease risk factors: repeated cross-sectional surveys from the Nord-Trøndelag health study (HUNT) 1984-2008. BMC Public Health. 2012;12:266. BioMed Central

    Article  PubMed  PubMed Central  Google Scholar 

  58. JE-O A. Inequalities in multimorbidity in South Africa. Int J Equity Health. 2013;12:64.

    Article  Google Scholar 

  59. Orueta JF, García-Álvarez A, Alonso-Morán E, Vallejo-Torres L, Nuño-Solinis R. Socioeconomic variation in the burden of chronic conditions and health care provision--analyzing administrative individual level data from the Basque Country, Spain. BMC Public Health. 2013;13:870. BioMed Central

    Article  PubMed  PubMed Central  Google Scholar 

  60. Dalstra J. Socioeconomic differences in the prevalence of common chronic diseases: an overview of eight European countries. Int J Epidemiol. 2005;34:316–26.

    Article  PubMed  CAS  Google Scholar 

  61. Korda RJ, Paige E, Yiengprugsawan V, Latz I, Friel S. Income-related inequalities in chronic conditions, physical functioning and psychological distress among older people in Australia: cross-sectional findings from the 45 and up study. BMC Public Health. 2014;14:741. BioMed Central

    Article  PubMed  PubMed Central  Google Scholar 

  62. McGrail KM, van Doorslaer E, Ross NA, Sanmartin C. Income-related health inequalities in Canada and the United States: a decomposition analysis. Am J Public Health. 2009;99:1856–63.

    Article  PubMed  PubMed Central  Google Scholar 

  63. van Doorslaer E, Koolman X, Jones AM. Explaining income-related inequalities in doctor utilisation in Europe. Health Econ. 2004;13:629–47.

    Article  PubMed  Google Scholar 

  64. Farrell L, Hollingsworth B, Propper C, Shields MA. The socioeconomic gradient in physical inactivity: evidence from one million adults in England. Soc Sci Med. 2014;123:55–63. Elsevier Ltd

    Article  PubMed  Google Scholar 

  65. Hajizadeh M, Campbell MK, Sarma S. Socioeconomic inequalities in adult obesity risk in Canada: trends and decomposition analyses. Eur J Health Econ. 2014;15:203–21.

    Article  PubMed  Google Scholar 

  66. Beiser M. The health of immigrants and refugees in Canada. Can J Public Health. 2005;96(Suppl 2):S30–44.

    PubMed  Google Scholar 

  67. McDonald JT, Kennedy S. Insights into the “healthy immigrant effect”: health status and health service use of immigrants to Canada. Soc Sci Med. 2004;59:1613–27.

    Article  PubMed  Google Scholar 

  68. Starfield B. Primary care: an increasingly important contributor to effectiveness, equity, and efficiency of health services. SESPAS report 2012. Gac Sanit. 2012;26(Suppl 1):20–6.

  69. Aggarwal M, Hutchison B. Toward a primary care strategy for Canada [internet]. 2012. Available from:

    Google Scholar 

  70. Starfield B, Shi L, Macinko J. Contribution of primary care to health systems and health. Milbank Q. 2005;83:457–502. Blackwell Publishing, Inc.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We thank IMS brogan Inc. for use of their drug information Database.


This study was funded by the Health System Performance Research Network (HSPRN) and supported by the Institute for Clinical Evaluative Sciences (ICES), which is funded by an annual grant from the Ontario Ministry of Health and Long-Term Care (MOHLTC). The opinions, results and conclusions reported in this paper are those of the authors and are independent from the funding sources. No endorsement by ICES or the MOHLTC is intended or should be inferred. Parts of this material are also based on data and information compiled and provided by the Canadian Institute for Health Information (CIHI). However, the analyses, conclusions, opinions and statements expressed herein are those of the authors, and not necessarily those of CIHI.

Availability of data and materials

The data from this study is held securely in coded form at the Institute for Clinical Evaluative Sciences (ICES). While data sharing agreements prohibit ICES from making the data set publicly available, access may be granted to those who meet pre-specified criteria for confidential access, available at (contact:

Author information

Authors and Affiliations



All authors contributed to the conception and design of the research study. LM generated the analytical cohort, performed all statistical analyses, and prepared the original draft of the manuscript. All co-authors contributed intellectual content to subsequent revisions. All co-authors have approved the final draft for submission and take public responsibility for it. All co-authors have read, and confirm that they meet ICMJE criteria for authorship.

Corresponding author

Correspondence to Walter P. Wodchis.

Ethics declarations

Ethics approval and consent to participate

We used survey data and linked health information that is routinely collected in Ontario that did not require informed consent from study participants. The use of this data was authorized under section 45 of Ontario’s Personal Health Information Protection Act, which does not require review by a Research Ethics Board.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Table S1. Multimorbidity prevalence (standard deviation) by determinant. Table S2. Decomposition results from 2007/08 and 2009/10 CCHS surveys. Table S3. Sensitivity Results. (DOCX 49 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mondor, L., Cohen, D., Khan, A.I. et al. Income inequalities in multimorbidity prevalence in Ontario, Canada: a decomposition analysis of linked survey and health administrative data. Int J Equity Health 17, 90 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: