Abstract
Introduction
Patients with complex health care needs (PCHCN) have numerous and costly health issues that place a heavy burden on health care resources (Hayes et al., 2016). Efforts in the United States to better classify PCHCN have shown that such patients are mainly frail elderly people, or individuals with multiple chronic conditions (multimorbidity). Targeting the needs of these patients, and particularly of those at high risk of specific health-related events, has raised interest in designing health care systems that can improve patients’ health outcomes while reducing the related costs (Hochman & Asch, 2017).
Multimorbidity in PCHCN poses new challenges to health care services (Marengoni et al., 2011). Various approaches have been used to elucidate the complexity of multimorbidity, based on the fundamental assumption that health outcomes in patients with multimorbidity are influenced not only by their single diseases, but also by the additional effects of interactions between them (Fortin, Stewart, Poitras, Almirall, & Maddocks, 2012). Identifying patterns in how chronic conditions occur in combinations in patients needing the most care and support is important to the efficient allocation of resources. Understanding how chronic conditions cluster together, and clarifying the impact of these disease patterns on health outcomes is essential to funding medical care and planning prevention and treatment services. Most research on multimorbidity has focused on measures adopting simple disease counts, but such approaches fail to capture specific combinations or patterns, and few studies have associated different chronic disease clusters in PCHCN with different demands on health care services with a view to better targeting these services to their needs.
Collerton and colleagues (2016) used cluster analysis to identify patterns in the burden of morbidity for a given patient, also exploring associations with their use of medication and health care. They found differences between comorbidity clusters that concerned the numbers of medicines, hospital admissions, general practitioner consultations, and general practice nurse consultations. Another study (Olaya et al., 2017) identified multimorbidity patterns using latent class analysis (LCA), finding associations between these patterns and several outcomes, including hospital admissions. The authors found two latent classes—cardiovascular/mental/arthritis and metabolic/stroke—and the former was significantly associated with more medical consultations and more hospital admissions.
A recently published review of studies on patterns of multimorbidity (Prados-Torres, Calderón-Larrañaga, Hancco-Saavedra, Poblador-Plou, & van den Akkerm, 2014) concluded for an urgent need to conduct further research of this type, essentially for the purpose of orienting future health policies, developing case management programs, and designing clinical practice guidelines to help professionals respond to the special needs of people with multimorbidity. The aim of the present study was to ascertain which disease clusters are associated with a greater need for hospital services, and with a higher mortality in a population of elderly patients with numerous health care needs, with a view to better organizing health care for these PCHCN.
Method
Context
The Italian National Health System (NHS) is a public service financed mainly by general taxation. It is grounded on values of universality, free access, freedom of choice, pluralism in provision, and equity. Regional authorities plan and organize health care facilities and activities through their regional health departments in accordance with a national health plan designed to assure an equitable provision of comprehensive care throughout the country. The regional authorities coordinate and control local health units (LHUs), each of which is a separate geographically based public company delivering public health promotion and community health services, primary care and hospital care, either with their own facilities and personnel or through outside contractors. The Veneto Regional Health Service had 21 such LHUs serving a population of about 5 million. The LHU involved in the present study was the “Azienda ex-ULSS4-Veneto,” which served a population of about 190,000 in the province of Vicenza, in northeast Italy.
To see which multimorbidity patterns in PCHCN are associated with particular demands for hospital services, our analysis was developed as follows:
a cohort of PCHCN was identified;
patients were characterized by their multimorbidity, identifying the number of chronic conditions, diseases dyads (or pairs), and clusters for each patient;
an association was sought between different multimorbidity patterns and different demands for hospital services (measured as at least one hospital admission, at least two hospital admissions, at least one preventable admissions, total number of admissions, total days in hospital), or death.
The ACG® (Adjusted Clinical Groups) System: Identifying the Sample
The ACG® System was implemented in the Veneto Region in 2012 as a tool for population risk stratification. The ACG® System is a method used internationally to characterize multimorbidity on the strength of routinely collected administrative data (e.g., hospital discharge records, pharmaceutical prescriptions, access to emergency departments, prescription charge exemptions) gathered using record linkage (The Johns Hopkins University, 2014). It relies on an algorithm that starts from individual-level diagnoses and is based on clinical judgments of likelihood (persistence or recurrence over time, demand for specialist services, hospitalization, disability or decline in quality of life, expected need for and use of diagnostic or therapeutic procedures), and is then adjusted for age and sex, to group a population by 93 mutually exclusive combinations of conditions. The ACG define clinically logical categories of patients expected to need similar levels of health care. Based on their health care resource usage, the ACG® System automatically collapses the different ACG into six Resource Utilization Bands (RUBs), which are defined as follows: 0 = nonuser or invalid diagnosis, 1 = healthy user, 2 = low morbidity, 3 = moderate morbidity, 4 = high morbidity, 5 = very high morbidity.
For the purposes of the present study, only people older than 65 years old in 2012, residing in the area served by the LHU “Azienda ULSS4-Veneto,” and characterized as PCHCN, classified as RUBs 4 and 5 were considered.
The ACG® System: Identifying Chronic Conditions
For each individual included in the study, we extracted the EDC (Expanded Diagnosis Clusters), which coincide with the clinical diagnosis that the ACG system assigns to single patients by combining different information flows. To improve the sensitivity of our model, patients with chronic conditions were also identified by means of the information available from the Pharmacy (RX)-based Morbidity Marker Groups (Rx-MGsTM), which is a drug-related diagnostic approach. A dichotomous variable was assigned to each chronic disease (taking into account either the EDC or the Rx-MGs diagnosis). Cases of neoplastic disease, Alzheimer’s disease, cardiac arrhythmia, and cerebrovascular disease were only available from the EDC codes, while cases of hyperlipidemia could only be obtained from the Rx-MGs diagnoses.
Unfortunately, there is no standard way to measure multimorbidity, so decisions concerning which morbidities to include, and how to define them are bound to be subjective to some degree, and depend strictly on the data available. This study focused on a subset of conditions including cancer, congestive heart failure (CHF), ischemic heart disease, high blood pressure (HBP), cardiac arrhythmia, cerebrovascular disease, Alzheimer’s disease, depression, asthma/bronchitis, diabetes, chronic obstructive pulmonary disease (COPD), osteoporosis, hypothyroidism, and chronic renal disease.
Identifying Outcomes
The subjects classified in 2012 as PCHCN (as explained above) were linked with the hospital discharge records in 2013 to identify patients who experienced any of the following:
at least one hospital admission;
at least two hospital admissions;
at least one preventable admission among the ones used to identify quality of care for “ambulatory care sensitive conditions,” defined as all discharges with any of the following International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) principal diagnostic codes: bacterial pneumonia, hypovolemia, urinary tract infection, angina, CHF, hypertension, asthma, COPD, uncontrolled diabetes, short-term complications of diabetes, and long-term complications of diabetes (see Appendix A);
total number of hospital admissions;
annual days in hospital.
The subjects were also linked to a mortality registry to identify patients who died in or before 2013 to assess the last outcome:
overall mortality.
Statistical Methods
Descriptive analyses were run by estimating the mean number of morbidities by sex and age group (“65-69,” “70-74,” “75-79,” “80-84,” and “85+” years old), and reporting the corresponding frequencies.
A previous study conducted an exploratory LCA to place patients in a number (K) of clinically meaningful classes of chronic diseases (Buja et al., 2017) and it has been reported in Appendix B. The number of classes was defined a priori using the Bayesian information criterion (BIC), a model selection approach that balances fit with parsimony (Schwarz, 1978). Ten different models were delineated, characterized by increasing numbers of chronic disease classes (from one to 10). Then, the model with five latent classes was chosen because it had the lowest BIC index. This result was confirmed by applying the Bootstrap LRT criterion (McLachlan, 1987) with 999 iterations: This latter method showed evidence in favor of the five-class model as opposed to the four-class model (p = .0099) or the six-class model (p = .0275).
Appendix B shows the class-conditional probabilities for each disease. Patients with ischemic heart disease, diabetes, hypertension, and hyperlipidemia seemed to be prevalent in Class 1 (which could be labeled as Metabolic-ischemic heart disease), whereas Class 2 was characterized by higher probabilities of Alzheimer’s disease, cerebrovascular disease, and depression (and was called Neurological). Congestive heart failure and fibrillation revealed higher class-conditional probabilities related to Class 3 (Heart impairment), and Class 4 included the majority of patients with asthma and chronic obstructive pulmonary disease (Cardio-respiratory diseases). Patients with neoplasms were only prevalent in the Class 5 (which was consequently labeled as Cancer). Hypertension clearly showed high class-conditional probabilities in all five classes, and this was because its prevalence in our dataset was high (87.8%). Logistic regression analyses were used to examine the associations between dichotomous dependent variables (e.g., at least one hospital admission, at least two hospital admissions, at least one preventable hospital admission, and death) and exposure variables (e.g., number of chronic conditions, or specific dyads, or latent classes). We defined dyads as pairs of diseases affecting the same patient. The results of the LCA were adopted here to shed light on the complex patterns of morbidity in our dataset. Each model was adjusted for age and sex.
A negative binomial model and a Tobit regression model, applied to censored dependent variables (Tobin, 1958), were used to study the associations between the continuous outcomes “number of hospital admissions” or “annual days in hospital” (total number of days spent in hospital in a given year) and the exposure variables (each disease dyad, or number of chronic conditions, or latent classes), respectively. Each model was adjusted for age and sex.
In the case of the “number of chronic conditions” exposure variable, for dichotomous dependent variables, trends in the proportions of the events were sought with the Cochran–Armitage test, whereas for continuous variables, the trend in the mean values of the outcomes was examined with a Mann-Kendall test.
As we were drawing a large number of comparisons, there being 91 possible dyads ([14 × 13] / 2), it became essential to correct the significance levels: We opted to do so by applying the Benjamini–Hochberg procedure (Benjamini & Hochberg, 1995).
Results
We analyzed data on a population of 190,000 residents served by the LHU: 2,028 of these people were older than 65 years old, classified as high-need, high-cost elderly according to the ACG system, and alive at the end of 2012. Table 1 shows the characteristics of this study population. More than 65% of these patients had four or more of the chronic diseases considered in this study.
Description of the Sample’s Characteristics.
Number of Chronic Conditions
Table 2 shows that the adjusted odds ratio (OR) for a range of adverse outcomes increases in a dose-response fashion with increasing numbers of chronic conditions. As regards outcomes in terms of at least one hospital admission, and at least two admissions, a significantly higher adjusted OR only emerged for the groups of patients with four or more diseases vis-à-vis the reference group with two. The number of chronic conditions was strongly associated with the outcome concerning one or more preventable hospital admissions: All the ORs were significantly higher than 1.0, and ranged from 3.28 (for patients with three diseases) to 10.62 (for those with seven or more). A significant trend emerged for each of these outcomes in relation to patients grouped by number of diseases.
Association Between the Number of Chronic Conditions and the Binary Outcomes as Measured With the Adjusted Odds Ratio and 95% Confidence Interval.
Note. Odds ratios highlighted in bold are significantly greater than 1.0.
Significant trend in proportions (Cochran–Armitage test).
Table 3 shows that having more chronic diseases coincided with a higher likelihood of spending more days in hospital per year, and of having a greater number of hospital admissions.
Effects of the Number of Chronic Conditions and the Outcomes “Number of Hospital Admission” and “Annual Days in Hospital.”
Note. IRR = incidence risk ratios.
Bold values signify the confidence intervals.
Significant trend in proportions (Cochran–Armitage test).
Dyads
Figure 1 shows the disease dyads associated with the outcome “at least one hospital admission,” including the adjusted OR and the relative confidence interval (CI). The strongest associations with this outcome were for the dyads chronic kidney disease and asthma (OR = 3.02, CI = [1.54, 6.09]), and coronary heart disease and chronic kidney disease (OR = 2.66, CI = [1.48, 4.86]). The two dyads showing a negative association with this outcome (adjusted OR less than 1.0) both included Alzheimer’s disease. Figures 2 and 3 show that chronic kidney disease, osteoporosis, and diabetes mellitus were involved in the dyads most strongly associated with the outcomes “at least two hospital admissions,” and “at least one preventable hospital admission.” In particular, the ORs for the dyad chronic kidney disease and depressive disorder were 4.11 (CI = [2.16, 7.58]) for the former, and 3.50 (CI = [1.78, 6.54]) for the latter, whereas for the dyad COPD and osteoporosis, they were 3.50 (CI = [1.85, 6.37]) and 4.98 (CI = [2.69, 8.92]), respectively. Figure 4 shows that chronic kidney disease was strongly associated with “death,” as it was included in all six disease dyads most strongly associated with this outcome.

Odds ratios and relative CIs for regressions between disease dyads and the outcome “at least one hospital admission.”

Odds ratios and relative CIs for regressions between disease dyads and the outcome “at least two hospital admissions.”

Odds ratios and relative CIs for regressions between disease dyads and the outcome “at least one preventable admission.”

Odds ratios and relative CIs for regressions between disease dyads and the outcome “death.”
As shown in Figures 5 and 6, the analyses on the outcomes “number of hospital admissions” and “annual days in hospital” essentially confirmed the previous results: Alzheimer’s disease was part of all the disease dyads showing a negative association with these outcomes.

IRR and relative CIs for regressions between disease dyads and the outcome “number of hospital admissions.”

Tobit marginal effects and relative CIs for regressions between disease dyads and the outcome “annual days in hospital.”
Latent Classes
Table 4 shows that the analysis of the associations between the five latent classes and the outcomes “at least two hospital admissions,” “at least one preventable hospital admission,” and “death” did not differ significantly for patients in the Neurological and Cancer classes, whereas the other classes showed a stronger association with these outcomes (taking the Neurological class for reference). Patients classed as Heart impairment, Metabolic-ischemic, Cardio-respiratory, or Cancer carried a higher risk of being hospitalized more than once than patients in the Neurological class. There was strong evidence of patients in the Heart impairment, Ischemic, or Cardio-respiratory disease groups being significantly more likely to spend more days in hospital too.
Associations Between Latent Classes and Different Outcomes (Reference Class: “Neurological”).
Bold values signify the confidence intervals.
Discussion
This study confirms that multiple chronic conditions are associated with higher odds of hospitalization or death. More intriguingly, our analyses also showed that certain disease dyads are more strongly associated with the need for hospital services than others. LCA also identified some classes of patients as being characterized by statistically significant higher odds of hospitalization or death.
In particular, the high-need, high-cost population studied here showed a statistically significant rising trend in the odds of being admitted to hospital (“at least one admission,” “at least two or more admissions,” and “at least one preventable admissions”) associated with more numerous chronic diseases. Several studies have found that having more chronic conditions is associated with a significantly greater recourse to health care services. Notably, the number of comorbidities is weakly associated with the number of hospital admissions but strongly associated with the total number of days spent in hospital, probably because it takes longer to stabilize patients with several comorbidities once they have arrived in hospital. One study (Hernandez et al., 2009) showed that the number of chronic conditions correlated with hospital admissions irrespective of age, and that each chronic condition added to the risk of hospitalization. Our findings indicate that the risk of “at least one preventable admission” rose steadily and significantly with the number of diseases in a given patient. The connection between multimorbidity and preventable hospital admission is well known. For example, a study by Wolff, Starfield, and Anderson (2002) found that the numbers of avoidable hospital admissions and of inpatients’ preventable complications increased dramatically for patients with larger numbers of chronic conditions.
Our present results also point to a significant association between the number of chronic diseases and mortality in patients with high health care needs. This result confirms a previous report of chronic conditions having a particularly significant association with the likelihood of death, even if the study focused on inpatient mortality (Iezzoni et al., 1994). In contrast, Marengoni, von Strauss, Rizzuto, Winblad, and Fratiglioni (2009) found that it was not multimorbidity, but chronic disability that emerged as the strongest negative prognostic factor for functionality and survival in a population-based cohort of older adults in Sweden. They reported that the hazard ratio for death was the same whether people had only one or as many as four chronic diseases.
Our study identified different probabilities of hospitalization depending on the type(s) of chronic disease. To our mind, it is intriguing that the “Neurological” class of conditions taken for reference coincided with significantly lower odds of being hospitalized than any of the other classes considered. Analyzing disease dyads confirmed as much, as the dyads with the lowest odds of hospital admission always included Alzheimer’s disease. This could be partly because people with dementia are more frequently institutionalized—in the Veneto Region, this is true of almost 28% of dementia sufferers (DRG Regione Veneto 653, 28/4/2015)—and this makes them less likely to be hospitalized than patients with high health care needs who live at home (Smith, Sherlaw-Johnson, Ariti, & Bardsley, 2015). This is a novel finding, as previous population studies compared the hospital admission rates of institutionalized versus noninstitutionalized general populations, not within a selected high-need population (Godden & Pollock, 2001).
On the contrary, some studies identified a different usage of health care resources for different multimorbidity patterns. For instance, one (Collerton et al., 2016) found five clusters of patients that differed in terms of number of medicines, hospital admissions, and consultations with general practitioners and general practice nurses, but their clusters were not characterized by type of chronic disease, so a comparison with our results was not possible. The authors concluded that clustering patients by similar morbidity profiles can help to inform future health care provision. As in the present study, another report (Dong, Wressle, & Marcusson, 2013) on multimorbidity patterns in relation to the use of health services again found five clusters (vascular, cardiopulmonary, cardiac, somatic-mental, malignancy), and showed that clusters involving cardiac and pulmonary conditions were more strongly associated with hospitalization than any single morbidity. Our findings confirm as much, as our cardiopulmonary group had higher odds of all outcomes regarding the need for hospitalization.
These findings are relevant for the purposes of developing an up-to-date, valid approach to population health management that is capable of identifying subpopulations and predicting their health care service needs. In addition, identifying subgroups of PCHCN at higher risk of needing certain kinds of health care enables the use of “impactibility models” (Lewis, 2010). Such models can help to identify subsets of patients at risk most likely to benefit from preventive care. In fact, the most often-mentioned way to increase the impact of predictive risk models has been to give priority to patients with certain diagnoses, known to be particularly amenable to upstream care. These are diseases for which prompt, high-quality primary or outpatient care can reduce the risk of hospitalization (Billings et al., 1993). More broadly, some authors (Billings & Mijanovich, 2007) claim that the “business case” for population health management lies in that the welfare of potentially high-cost patients is offset by savings deriving from reducing their future hospitalizations. As our study shows, existing data resources can be used to predict with a reasonable degree of accuracy which comorbidities make patients most likely to need hospitalizing. These data point to the importance of a highly functional, high-value, proactive primary care model for individuals with multiple chronic conditions.
The main strength of our study lies in that it was population-based, thus minimizing selection bias and relying on independently collected data. The study also suffers from some limitations, however, because administrative health care records are unable to accurately register some highly prevalent chronic health conditions that may be associated with a heavy health care burden, such as chronic pain, or visual or hearing impairments. Administrative records also do not track economic and social factors that might independently influence the demand for health care services.
Conclusion
In terms of hospitalization, our results are consistent with those of other studies investigating the association between multimorbidity and recourse to hospital services. The advantage of our approach lies in that analyzing morbidity clusters or dyads provides more detailed information about the risk of hospitalization events. Unlike morbidity counts (in which all morbidities are scored equally, irrespective of the relationships between them), analyzing morbidity clusters and dyads reveals which combinations of morbidities are associated with the highest hospitalization rates. In our sample, cardiopulmonary diseases were the chronic conditions most strongly associated with hospitalization.
Footnotes
Appendix
Class-Conditional Probabilities for Each Disease Estimated With the Five Latent Classes Model.
| Disease | Latent classes |
||||
|---|---|---|---|---|---|
| Class 1 |
Class 2 |
Class 3 |
Class 4 |
Class 5 |
|
| % of respondents in each class |
|
|
|
|
|
| Cancer | 41 | 17 |
|
28 |
|
| Coronary heart disease |
|
8 | 22 | 21 | 8 |
| Alzheimer’s disease | 5 |
|
14 | 16 | 8 |
| Cerebrovascular disease | 35 |
|
30 | 17 | 16 |
| Diabetes mellitus |
|
22 | 33 | 30 | 26 |
| High blood pressure |
|
|
|
|
75 |
| Heart failure |
|
41 |
|
|
20 |
| Chronic obstructive pulmonary disease | 4 | 4 | 4 |
|
5 |
| Chronic kidney disease | 11 | 2 | 16 | 12 | 3 |
| Depressive disorder | 26 |
|
25 | 28 | 18 |
| Cardiac arrhythmia | 34 | 23 |
|
37 | 16 |
| Asthma | 12 | 4 | 15 |
|
12 |
| Osteoporosis | 20 | 16 | 19 | 19 | 19 |
| Hypothyroidism | 12 | 7 | 8 | 9 | 12 |
| Hyperlipidemia |
|
10 | 4 | 14 | 12 |
Note. The higher class-conditional probabilities are highlighted in bold.
Author Contributions
Alessandra Buja was involved in conceptualization, study design and methods, supervision, writing and approval of final article as submitted. Michele Rivera was involved in formal analysis, writing and revision of article, approval of final article as submitted. Maria Chiara Corti was involved in conceptualization, supervision, critical review of article, approval of final article as submitted. Francesco Avossa was involved in conceptualization, data collection, approval of final article as submitted. Elena Schievano was involved in conceptualization, data collection, approval of final article as submitted. Stefano Rigon was involved in conceptualization, data collection, approval of final article as submitted. Vincenzo Baldo was involved in conceptualization, review of manuscript, approval of final article as submitted. Giovanna Boccuzzo was involved in conceptualization, critical review of article, approval of final article as submitted. Mark Ebell was involved in revision of article, critical review of article, and approval of final article as submitted.
Ethical Considerations
The data analysis was performed on anonymized aggregate data with no chance of individuals being identifiable. The study complied with the Declaration of Helsinki and with Italian Law n. 196/2003 on the protection of personal data. The recent resolution n. 85/2012 of the Italian Guarantor for the Protection of Personal Data also confirmed the allowability of processing personal data for medical, biomedical, and epidemiological research, and that data concerning health status may be used in aggregate form in scientific studies.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
