Skip to main content
  • Research Article
  • Open access
  • Published:

Latent class analysis of multimorbidity patterns and associated outcomes in Spanish older adults: a prospective cohort study



This study sought to identify multimorbidity patterns and determine the association between these latent classes with several outcomes, including health, functioning, disability, quality of life and use of services, at baseline and after 3 years of follow-up.


We analyzed data from a representative Spanish cohort of 3541 non-institutionalized people aged 50 years old and over. Measures were taken at baseline and after 3 years of follow-up. Latent Class Analysis (LCA) was conducted using eleven common chronic conditions. Generalized linear models were conducted to determine the adjusted association of multimorbidity latent classes with several outcomes.


63.8% of participants were assigned to the “healthy” class, with minimum disease, 30% were classified under the “metabolic/stroke” class and 6% were assigned to the “cardiorespiratory/mental/arthritis” class. Significant cross-sectional associations were found between membership of both multimorbidity classes and poorer memory, quality of life, greater burden and more use of services. After 3 years of follow-up, the “metabolic/stroke” class was a significant predictor of lower levels of verbal fluency while the two multimorbidity classes predicted poor quality of life, problems in independent living, higher risk of hospitalization and greater use of health services.


Common chronic conditions in older people cluster together in broad categories. These broad clusters are qualitatively distinct and are important predictors of several health and functioning outcomes. Future studies are needed to understand underlying mechanisms and common risk factors for patterns of multimorbidity and to propose more effective treatments.

Peer Review reports


Chronic conditions and non-communicable diseases are the leading cause of morbidity and disability worldwide [1]. The co-occurrence of two or more conditions, or multimorbidity, is especially common in elderly people. Approximately two out of three persons at retirement age suffer from at least two chronic diseases [2]. The therapeutic management of multimorbidity is often complex because people suffering from multimorbidity present high treatment burden, poorer health outcomes, problems in independent living [3, 4] and higher rates of mortality [5].

The classical single-disease paradigm might not be adequate for patients with multiple chronic conditions [6]. Studies have focused on different approaches to assess multimorbidity. Counting the number of diseases has been a commonly used method [7, 8] but results could be affected by the number of chronic conditions considered. One variant is the conditional count [9, 10]. For example, people with arthritis have, on average, 3.8 co-morbid conditions, compared with 1.8 conditions among persons without arthritis [9]. Other approaches include the most prevalent pair combinations. Diabetes and cataract, or asthma and chronic obstructive pulmonary disease (COPD), are especially prevalent and associated with poor levels of quality of life and high disability [4]. However, these approaches might not adequately capture larger clusters of conditions linked to greater burden. In recent years, the co-occurrence of chronic conditions has been addressed with cluster analysis methods as they can help to identify broad comorbidity patterns.

Islam et al. [6] compared distinct methods to assess multimorbidity patterns in a community sample of elderly people, including the most frequent pairs and triplets of comorbid diseases, cluster analysis, principal component analysis and latent class analysis (LCA). They found considerable consistency across the analytic methods although some differences arose, mainly due to the underlying methodology. Cluster analysis is based on distance measures, whereas principal components and LCA are based on correlations, with similarities between the last two approaches [6].

Few studies have analyzed multimorbidity patterns using LCA. LCA is based on structural equation modeling which allows identification of latent groups based on a set of observed variables. It can be useful in describing how people are clustered according to patterns of chronic conditions and identifying the main differences between these groups in terms of sociodemographics, functioning and clinical characteristics. Whitson et al. [11] reported six latent groups in a community dwelling sample of people aged 65+. These multimorbidity groups were significantly associated with emergency department use and hospitalization over 1 year. However, the authors reported also high misclassification errors as one limitation of the study. There is still a need to understand how chronic conditions are clustered together, whether these multimorbidity patterns are reliable and valid, and determine the impact on health and functioning outcomes at baseline and longitudinally. Moreover, there is also a need for general-population studies to ensure the external validity of multimorbidity patterns [12].

The aim of this study was to investigate whether eleven common chronic conditions cluster together in a Spanish representative sample of people aged 50 and older according to their pattern of co-occurrence. LCA was conducted to identify latent groups of people based on a set of observed variables [13]. These latent groups were then described in terms of sociodemographic characteristics, co-morbid conditions, quality of life, disability, functioning and cognitive function. The second part of the analysis involved the study of the effect of these latent co-morbid groups on several outcomes measured at baseline and at three-year follow-up.


Design and setting

The present study used data from “Edad con Salud”, a longitudinal, nationally representative survey of adult, non-institutionalized people in Spain. The first wave was part of the Collaborative Research on Ageing in Europe (COURAGE in Europe) study [14]. A stratified, multistage clustered design was used and people over 50 and 80 years old were oversampled. Strata included all Autonomous Communities in Spain (except Ceuta and Melilla). Data on households were provided by the Spanish Statistical Office. A total of 4753 persons participated in face-to-face structured interviews conducted at their homes between 2011 and 2012. The final response rate at baseline was 69.9%. The second wave took place after approximately 3 years (2014–2015). 2528 participants (53.2%) completed the interview at follow-up, 259 (5.5%) had died, 862 (18.1%) declined to participate in the second assessment, and 23.2% did not participate for other reasons (e.g., unable to locate, institutionalization).

The present analysis focused on people aged 50 years of age and older who completed a non-proxy interview at baseline (n = 3625). Of these, 1970 participated in an interview in the second wave. In a small proportion (n = 84), a proxy responded to the second assessment because the original participant presented some evidence of cognitive deterioration. Since proxy interviews were much shorter, we excluded these proxy interviews at follow-up, resulting in a final n of 3541 participants at baseline, of whom 1886 also participated in the second assessment.


Presence or absence of eleven chronic conditions at baseline was considered. Participants reported whether they had received a medical diagnosis during the previous 12-months of depression, arthritis, asthma, COPD, angina, stroke, hypertension, diabetes, edentulism and cataract. Other symptom questions were added to the interview based on the WHO-SAGE protocol [15] and the adapted version of the World Health Organization Composite International Diagnostic Interview (CIDI) for depression [16]. These additional questions allowed us detect undiagnosed cases. Algorithms were implemented [17, 18] and an individual was considered to have one of these conditions if he/she met criteria for at least the self-reported diagnosis, or the diagnosis made according to symptoms. The presence of hypertension was based on self-reported diagnosis or presence of systolic blood pressure ≥ 140 mmHg or diastolic blood pressure ≥ 90 mmHg [19, 20] measured at the time of the interview. Blood pressure was measured twice with less than a 1 minute interval using an arm blood pressure monitor and taking the average of the measurements. Interviewers recorded participants’ height and weight using a stadiometer and a routinely calibrated electronic weighing scale, respectively. Body Mass Index (BMI) was calculated as weight (in kilograms) divided by the square of height (in meters). A BMI of 30 or higher was used as cut-off point for obesity [21].

Self-reported demographic variables at baseline included age, gender, years of schooling, quintiles of household income (with the first quintile indicating lowest level and the fifth the highest), and marital status (never married, married or currently cohabiting, separated or divorced, and widowed).

Respondents were asked to recall a list of words three times immediately and once after a short delay which was filled with other cognitive tests (Consortium to Establish a Registry for Alzheimer’s Disease) [22]. The psychometric properties of the Word List Learning task have previously been established [23, 24]. A composite score was calculated as the sum of the number of correct words, ranging from 0 to 40, with higher scores indicating better memory. Participants were asked to name as many animals as possible in 1 minute, both at baseline and follow-up. Animal naming tasks are considered a measure of verbal fluency whereas word-list recall is regarded as a measure of verbal memory.

Disability was assessed with the 12-item version of the World Health Organization Disability Assessment Schedule 2.0 (WHODAS 2.0) [25] which evaluates functioning in six life domains. A total score was obtained as the sum of the items and transformed into a 0 to 100 scale (higher scores indicating greater disability). Presence or absence of difficulties in independent living was evaluated through activities of daily living (ADLs) and instrumental activities of daily living (IADLs). ADLs describe a set of daily self-care activities and assess the need for help with personal care activities such as eating, bathing and dressing; IADLs describe higher-level functioning considered necessary to live independently (using transportation, housekeeping, or preparing food). ADLs and IADLs difficulties were present if the person answered severe or extreme/cannot do it to any of the questions. Quality of life was measured at baseline and follow-up with the WHOQOL-AGE [26], a modified version of the World Health Organization Quality of Life instrument (WHOQOL) adapted for the elderly population. It has 13 items and a global score can be obtained, ranging from 0 (minimum quality of life) to 100 (maximum quality of life). The total number of visits to any health professional in the last year as well as whether the participant had been hospitalized in the previous 12 months (yes/no) was also recorded at baseline and at follow-up.

Statistical analysis

Latent Class Analysis (LCA) was conducted on the 3541 participants at baseline using the Stata plugin [27]. Eleven chronic health conditions (arthritis, asthma, COPD, angina, stroke, hypertension, diabetes, edentulism, cataract, depression, and obesity) were used as observed indicators. The optimal number of latent classes was determined using the adjusted Bayesian Information Criterion (BIC) [28] and the consistent Akaike Information Criterion (CAIC) [29], which have been shown to be more robust indicators of class enumeration with categorical outcomes [30]. The adjusted BIC and CAIC were used to compare several plausible models where the lowest values indicate the best fitting model.

Furthermore, interpretability and clinical judgment were used. After selecting the best model, each participant was assigned to one class according to his or her highest computed probability of membership. Average posterior probabilities above 70% indicate optimal fit [31]. Each latent class was labeled according to those chronic conditions whose prevalence exceeded the prevalence in the full cohort [11]. The final latent groups were compared in terms of co-morbidities and functioning and sociodemographic variables at baseline using chi-squared and Kruskal-Wallis tests for significance.

Missing data at baseline and follow-up was handled using a Stata command for Imputation by Chained Equations (ICE) [32] assuming missing-at-random (MAR). The imputation model included all the variables used in the regression models plus other auxiliary variables (variables not included in the analysis model but which are potential predictors of missingness [33]) (see Additional file 1: Table S1). Since the imputed values of the continuous variables lay outside the observed data range, they were transformed to normality before imputing their values and then converted back to the original variables following multiple imputation [34]. LCA was performed without imputing missing values in the eleven indicators because LCA is not supported by mim command [35] for analyzing multiply imputed datasets. Instead, missing data in one of the indicators was tolerated and handled with a full-information maximum likelihood (FIML) technique, assuming MAR [27]. Additional file 1: Table S1 shows the proportion of missingness in each of the indicators for LCA.

Multivariable linear, logistic and Poisson regression models were computed to assess the association of each multimorbidity class with several outcomes at baseline and after 3 years of follow-up. Models for outcomes at baseline were adjusted for gender, age, years of schooling, household income and marital status (at baseline). Models for outcomes at follow-up were additionally adjusted for the same measure as the outcome at baseline. Regression models were conducted separately in one hundred imputed datasets and results combined using Rubin’s rules [36] in Stata SE version 13 (College Station, TX).

For sensitivity analysis, the same models described above were conducted again in the subsample of participants with valid values in all the variables (n = 1508) and compared with the results from the imputed datasets.


Latent classes of multimorbidity pattern

Table 1 shows the adjusted BIC and CAIC values for the two to six-class models. There was an important drop in the adjusted BIC and CAIC values from the 2-class to the three-class model. The three-class model yielded the lowest CAIC value (CAIC = 1542.11) and, although the adjusted BIC corresponded to the four-class model (adjusted BIC = 1382.76), further inspection showed that one of the four classes presented a posterior probability lower than 0.7. Given the negligible difference between three and four-class models in terms of the adjusted BIC and following the parsimonious principle, the three-class model was finally chosen over the rest.

Table 1 Comparison between models

Table 2 shows the distribution of sociodemographic, clinical and functioning characteristics in the overall sample at baseline and by multimorbidity classes. Some 63.8% of people were classified as being in the “healthy” class, with prevalence of all conditions below that observed in the whole sample. The “cardiorespiratory/mental/arthritis” class (6.2%) presented excess prevalence of depression (39.5%), arthritis (64.0%), asthma (73.2%), COPD (83.2%) and angina (40.0%). The “metabolic/stroke” class, which comprised 29.9% of the sample, showed excess prevalence of stroke (12.3%), obesity (55%), diabetes (35.1%) and hypertension (89.2%). The proportion of oral problems and cataract was high in both the “cardiovascular/mental/arthritis” and the “metabolic/stroke” classes.

Table 2 Baseline characteristics of the overall sample and people assigned to three multimorbidity classes

The average posterior probabilities for all three classes exceeded 0.7 (0.85 for the “healthy” group, 0.81 for the “cardiorespiratory/mental/arthritis” class, and 0.75 for the “metabolic/stroke” group). Participants in the “healthy” class were more likely to be well-classified (79.7% of them had 0.7 or greater posterior probability). 68.6% and 60.9% of participants presented a posterior probability equal to or higher than 0.7 in the “cardiorespiratory/mental/arthritis” and the “metabolic/stroke” class, respectively.

People in the “metabolic/stroke” and “cardiorespiratory/mental/arthritis” latent groups were significantly older than those in the “healthy” group (Table 2). There were more females in the “metabolic/stroke” group, whereas the highest proportion of widows was seen in the “cardiorespiratory/mental/arthritis” class. This last group presented fewer years of schooling. The number of chronic conditions was significantly higher in the “cardiorespiratory/mental/arthritis” class. The prevalence of limitations in ADLs and IADLs was significantly higher in the “cardiorespiratory/mental/arthritis” class. The lowest level of quality of life was observed in the “cardiorespiratory/mental/arthritis” class, followed by the “metabolic/stroke” group. Level of disability was significantly higher in the “cardiorespiratory/mental/arthritis” and “metabolic/stroke” classes, compared with the healthy group. The lowest scores in verbal memory and verbal fluency at baseline were also seen in the “cardiorespiratory/mental/arthritis” class.

Association between multimorbidity classes with outcomes at baseline and after 3 years

Table 3 displays the adjusted unstandardized coefficients for several outcomes assessed at baseline and follow-up. The “healthy” group was used as the reference group. At baseline, the two multimorbidity groups were associated with all the outcomes, except for verbal fluency. Being in the “cardiovascular/mental/arthritis” class was significantly related to lower levels of verbal memory, higher levels of disability, poorer quality of life, higher number of medical visits, higher risk of limitations in ADLs (OR = 7.91, 95%CI = 5.41–11.73, p < 0.001), IADLs (OR = 11.42, 95%CI = 8.04–16.22, p < 0.001) and more hospital admissions (OR = 3.77, 95%CI = 2.79–5.10, p < 0.001). Similarly, and compared with the healthy class, participants in the “metabolic/stroke” class were more likely to present poorer verbal memory and quality of life, higher levels of disability and number of medical visits in the last 12 months, greater odds for limitations in ADLs and IADLs, and being hospitalized more frequently. Inspection of ORs shows that the risk of presenting limitations in ADLs and IADLs and hospital admission was greater for people assigned to the “cardiovascular/mental/arthritis” group than for those assigned to the “metabolic/stroke” class.

Table 3 Association between latent multimorbidity membership and outcomes at baseline and follow-up

In terms of cognitive function, verbal memory at follow-up was not significantly associated with multimorbidity groups. Only people in the “metabolic/stroke” class presented significantly lower levels of verbal fluency at follow-up, compared with the “healthy” group (b = −0.78, SE = 0.30; p < 0.01). Being in the “cardiorespiratory/mental/arthritis” or the “metabolic/stroke” groups was significantly associated with lower quality of life scores as well as a higher number of medical visits in the previous year. Persons classified in the “cardiorespiratory/mental/arthritis” group were four times more likely to have limitations in ADLs (OR = 4.05, 95%CI = 2.65–6.18) and twice as likely to present limitations in IADLs (OR = 2.17, 95%CI = 1.32–3.55) at follow-up. Membership of the “metabolic/stroke” class was significantly associated with greater risk of presenting limitations in ADLs (OR = 2.11, 95%CI = 1.64–2.72) and limitations in IADLs (OR = 1.63, 1.17–2.29) after 3 years. Both groups, compared with the “healthy” group, were at greater risk of hospital admission at follow-up.

Sensitivity analysis

Regression models were run again using complete case analysis (n = 1508). In general, the results were similar to those obtained using multiple imputation. However, when using complete cases, the “cardiorespiratory/mental/arthritis” group was not associated with lower memory scores at baseline, and being in the “metabolic/stroke” class was not significantly related to lower verbal fluency scores at follow-up. Conversely, membership of the “cardiorespiratory/mental/arthritis” group appeared to be significantly related to disability at follow-up whereas this association was not significant when analyzing the imputed data (see Additional file 2: Table S2).


This study sought to describe multimorbidity patterns using LCA in a representative sample of Spanish community-dwelling adults over 50 years old. The LCA procedure identified three latent classes of multimorbidity which were statistically and clinically distinct, based on the presence or absence of eleven chronic conditions. These clusters generated using LCA were significantly related to several health, functioning and use of health service outcomes at baseline and were still significant predictors of most of them after 3 years of follow-up.

Previous studies conducting LCA to describe latent classes of co-occurring conditions in community samples of older people have yielded mixed results as regards the number of clusters detected. Four latent classes were identified in a cross-sectional sample of 4574 Australian seniors using eleven chronic conditions [6], including cancer, Parkinson’s disease, or osteoporosis. Another study found six clusters of multimorbidity in a sample of 14,502 people aged 65 years old and over, using 13 conditions, which also included neurological diseases, cancer and osteoporosis [11]. Comparison with these studies is difficult since the results might be influenced by the number and type of diseases included in the analysis, characteristics of the sample, or how data on diseases were collected. In our study, information about cancer, osteoporosis and dementia was not available. Conversely, we included other highly prevalent conditions such as edentulism, cataract and obesity.

The majority of our sample (63.8%) was classified into the “healthy” class. This latent group has previously been reported in other studies which also conducted LCA [6, 11]. However, the proportion classified as “healthy” in our study is larger than that described in these studies. This difference could be explained by the age of participants. For example, Whitson et al. [11] found that 32.8% of their community sample was classified in the minimal disease category but the average age was older than in our study (76.4 vs. 65.7 years). Our findings support the existence of broad multimorbidity patterns. These clusters are very similar to those reported in a review where 14 studies on patterns of multimorbidity were considered [12]. Despite the fact that there was considerable heterogeneity between studies in terms of number and types of chronic conditions included or the statistical approach used, the authors concluded that there are at least three broad patterns; one comprising cardiovascular and metabolic diseases, a second one related to mental health problems and a third including musculoskeletal disorders. In our study, 30% of participants were classified under the “metabolic/stroke” category. This cluster is close to the “metabolic syndrome” which has been shown to increase the risk of stroke and diabetes [37]. The least frequent latent group was the “cardiorespiratory/mental/arthritis” class. However, it appeared to be the most severe category, with worst functioning and greatest burden. This category clusters a large number of pathologies together (i.e., angina, COPD, asthma, depression and arthritis). Previous studies also reported a cluster of angina and respiratory diseases [17]. The co-occurrence of health problems and musculoskeletal disorders has been consistently reported [38, 39] and has been associated with a constellation of comorbidities. However, the link between mental health conditions and COPD is still unclear [40].

Cataract and edentulism were both highly prevalent in the two multimorbidity clusters. In a similar study, visual impairment was associated with a history of stroke, diabetes and arthritis [4, 41]. Patients with arthritis are more likely to suffer from cataract, after adjusting for glucocorticoid intake [42]. Heart diseases have also been related to higher risk for cataract [43], supporting the theory of inflammatory pathways [44]. The presence of edentulism has been related to diabetes, coronary artery disease, hypertension, and rheumatoid arthritis [45].

Several underlying mechanisms could explain the nonrandom association between chronic conditions. Insulin resistance has been proposed as one possible underlying mechanism explaining the strong association between metabolic syndrome and stroke, affecting metabolic processes and leading to abnormalities of vascular reactivity [37]. A change of lifestyles to reduce metabolic syndrome can help prevent stroke and other vascular diseases. Other therapeutic strategies can include targeting insulin resistance [37]. The association between respiratory diseases and coronary heart disease has previously been established [46] and could be explained by inflammation, hypoxia, or stress processes. Other environmental risk factors could include smoking or air pollution [17]. Unexpected associations between conditions (such as arthritis and respiratory diseases) should be studied in the future. It has been suggested that medication could be a risk factor for co-occurrence of certain groups of diseases [12]. Future research should focus on finding the underlying pathogenesis connecting these medical conditions, and the shared risk factors.

Our findings show that these three multimorbidity patterns are qualitatively distinct with important differences with respect to sociodemographic, clinical and functioning aspects. Being older, female and having a low educational level have been consistently associated with more risk of suffering from multimorbidity [47, 48]. In common with other studies [4], our results show that multimorbidity clusters were cross-sectionally associated with a greater degree of disability, poor functioning, lower quality-of-life levels, poor memory function and greater risk of health care visits and hospitalizations. Persons assigned to one of these two multimorbidity clusters still presented poor quality of life, were more likely to present problems in independent living, and used health services more frequently (including being hospitalized) after 3 years of follow-up.

One interesting result is the significant association between multimorbidity class membership and cognitive function. Worse verbal fluency at follow-up was only linked to membership of the “metabolic/stroke” class. Depression, which is part of the “cardiovascular/mental/arthritis” cluster, has been consistently linked to memory complaints [49] as well as worse performance in memory tests [50]. Previous longitudinal studies have also reported an association between musculoskeletal diseases, lung diseases or arthritis, with cognitive decline [51]. Vascular risk factors and vascular diseases, such as stroke, have been consistently associated with cognitive deterioration in older adults [52]. Similarly, Ganguli et al. [53] found cross-sectional associations between a history of stroke, diabetes and abdominal adiposity with worse memory and executive function in a population-based cohort of old people, although these associations were not observed after 4 years of follow-up. The authors suggested that the effect of these diseases on brain structures would be static rather than progressive [53]. Performance in memory tests might also be affected by the learning effect [54]. Longer follow-up periods are needed to observe possible cognitive decline over time associated with these multimorbidity clusters and to avoid possible learning effects.

Some limitations should be considered when interpreting our findings. In our study, some diseases were not evaluated (such as cancer or neurological diseases) and the inclusion of additional chronic conditions might have yielded some different patterns. The COURAGE protocol included a limited number of chronic conditions based on their high prevalence and impact on health outcomes. Conversely, depression and obesity were considered in our study, although they are commonly omitted in other research. Second, the presence of chronic diseases was partially based on self-reporting, and can thus be affected by measurement errors or lack of accuracy. However, the literature shows that self-reported measures of chronic diseases are widely used in large population-based studies and show reasonable accuracy [55,56,57]. Moreover, we used additional questions about symptoms during the interview, allowing us to detect undiagnosed cases. For hypertension and obesity, objective measures were obtained during the interview. Next, misclassification of persons assigned to each latent class is reasonable [31]. However, there is still some degree of uncertainty associated with latent class membership and results should be interpreted with caution [27]. Finally, the use of multiple imputations could add some bias. Nevertheless, the sensitivity analysis showed similar results. The few differences between results with complete cases and imputed data might be explained by greater change variation when using the former, and because under MAR assumption, multiple imputation should correct biases that may arise in complete cases analyses [58].


This study identified three qualitatively separate, broad multimorbidity clusters using LCA in a Spanish nationally representative sample of older adults with distinct clinical and sociodemographic characteristics. The latent classes identified presented relatively low misclassification errors, and demonstrated predictive and external validity. Multimorbidity has been consistently related to greater burden and increased use of health services. The single-disease paradigm does not seem to fit the majority of persons with more than one chronic condition. Future efforts should focus on the underlying mechanisms of these multimorbidity clusters (lifestyles, metabolic or inflammatory factors, stress, and environmental factors) and determine targets for prevention and intervention.



95% confidence interval


Activities of daily living


Bayesian Information Criterion


Body Mass Index


Corrected Akaike Information Criterion


Composite International Diagnostic Interview


Chronic obstructive pulmonary disease


Collaborative Research on Ageing in Europe


Diagnostic and Statistical Manual of Mental Disorders


Full-information maximum likelihood


Instrumental activities of daily living


Imputation by chained equations


Latent class analysis




Odds Ratio


Standard error


World Health Organization Disability Assessment Schedule 2.0


World Health Organization Quality of Life instrument


  1. Wang H, Naghavi M, Allen C, Barber RM, Bhutta ZA, Carter A, et al. Global, regional, and national life expectancy, all-cause mortality, and cause-specific mortality for 249 causes of death, 1980–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet. 2016;388(10053):1459–544. 

  2. Kirchberger I, Meisinger C, Heier M, Zimmermann A-K, Thorand B, Autenrieth CS, et al. Patterns of multimorbidity in the aged population. Results from the KORA-age study. PLoS One. 2012;7(1):e30556.

  3. Wallace E, Salisbury C, Guthrie B, Lewis C, Fahey T, Smith SM. Managing patients with multimorbidity in primary care. BMJ. 2015;350:h176.

  4. Garin N, Olaya B, Moneta MV, Miret M, Lobo A, Ayuso-Mateos JL, et al. Impact of multimorbidity on disability and quality of life in the Spanish older population. PLoS One. 2014;9(11):e111498.

  5. Lee TA, Shields AE, Vogeli C, Gibson TB, Woong-Sohn M, Marder WD, et al. Mortality rate in veterans with multiple chronic conditions. J Gen Intern Med Springer. 2007;22(Suppl 3):403–7.

  6. Islam MM, Valderas JM, Yen L, Dawda P, Jowsey T, McRae IS, et al. Multimorbidity and Comorbidity of chronic diseases among the senior Australians: prevalence and patterns. PLoS One Pub Lib Sci. 2014;9(1):e83783.

  7. Fortin M, Bravo G, Hudon C, Vanasse A, Lapointe L. Prevalence of multimorbidity among adults seen in family practice. Ann Fam Med. American Academy of Family Physicians. 2005;3:223–8.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Fried LP, Bandeen-Roche K, Kasper JD, Guralnik JM, Fried LP, Guralnik JM, et al. Association of Comorbidity with disability in older women. J Clin Epidemiol Elsevier. 1997;52:27–37.

    Article  Google Scholar 

  9. Verbrugge LM, Lepkowski JM, Konkol LL. Levels of disability among U.S. adults with arthritis. J. Gerontol. Oxford University Press. 1991;46:S71–83.

    Article  CAS  PubMed  Google Scholar 

  10. John R, Kerby DS, Hennessy CH. Patterns and impact of Comorbidity and multimorbidity among community-resident American Indian elders. Gerontologist Oxford University Press. 2003;43:649–60.

    Article  PubMed  Google Scholar 

  11. Whitson HE, Johnson KS, Sloane R, Cigolle CT, Pieper CF, Landerman L, et al. Identifying patterns of multimorbidity in older Americans: application of latent class analysis. J Am Geriatr Soc. 2016;64:1668–73.

  12. Prados-Torres A, Calderón-Larrañaga A, Hancco-Saavedra J, Poblador-Plou B, van den Akker M, Barnett K, et al. Multimorbidity patterns: a systematic review. J Clin Epidemiol. 2014;67:254–66. 

  13. Formann AK, Kohlmann T. Latent class analysis in medical research. Stat Methods Med Res. 1996;5:179–211. 

  14. Leonardi M, Chatterji S, Koskinen S, Ayuso-Mateos JL, Haro JM, Frisoni G, et al. Determinants of health and disability in ageing population: the COURAGE in Europe project (collaborative research on ageing in Europe). Clin Psychol Psychother. 21:193–8. 

  15. Kowal P, Chatterji S, Naidoo N, Biritwum R, Fan W, Lopez Ridaura R, et al.Data resource profile: the World Health Organization study on global AGEing and adult health (SAGE). Int J Epidemiol. 2012;41:1639–49.

  16. Kessler RC, Ustün TB. The world mental health (WMH) survey initiative version of the World Health Organization (WHO) composite international diagnostic interview (CIDI). Int J Methods Psychiatr Res. 2004;13:93–121. 

  17. Garin N, Koyanagi A, Chatterji S, Tyrovolas S, Olaya B, Leonardi M, et al. Global multimorbidity patterns: a cross-sectional, population-based, multi-country study. J Gerontol A Biol Sci Med Sci. 2016;71:205–14.

  18. American Psychiatric Association. Diagnostic and statistical manual of mental disorders. 5th ed. Washington DC: American Psychiatric Association; 2001.

    Google Scholar 

  19. Basu S, Millett C. Social Epidemiology of Hypertension in Middle-Income Countries: determinants of prevalence, diagnosis, treatment, and control in the WHO SAGE study. Hypertension. 2013;62(1):18–26.

  20. Mancia G, Fagard R, Narkiewicz K, Redón J, Zanchetti A, Böhm M, et al. 2013 ESH/ESC Guidelines for the management of arterial hypertension: the Task Force for themanagement of arterial hypertension of the European Society of Hypertension (ESH) and of the European Society of Cardiology (ESC). J Hypertens. 2013;31:1281–357.

  21. World Health Organization. WHO: Global Database on Body Mass Index [Internet]. Available from:

  22. Moms JC, Heyman A, Mohs RC, Hughes JP, van Belle G, Fillenbaum G, et al. The consortium to establish a registry for Alzheimer’s disease (CERAD). Part I. Clinical and neuropsychological assesment of Alzheimer’s disease. Neurology. 1989;39:1159–65.

  23. Andel R, CA MC, Murdock GA, Fiske A, Wilcox RR, Gatz M. Performance on the CERAD word list memory task: a comparison of university-based and community-based groups. Int J Geriatr Psychiatry. 2003;18:733–9. 

  24. Chen P, Ratcliff G, Belle SH, Cauley JA, DeKosky ST, Ganguli M. Cognitive tests that best discriminate between presymptomatic AD and those who remain nondemented. Neurology. 2000;55:1847–53. 

  25. Rehm J, Üstün TB, Saxena S, Nelson CB, Chatterji S, Ivis F, et al. On the development and psychometric testing of the WHO screening instrument to assess disablement in the general population. Int J Methods Psychiatr Res. 1999;8:110–22.

  26. Caballero FF, Miret M, Power M, Chatterji S, Tobiasz-Adamczyk B, Koskinen S, et al. Validation of an instrument to evaluate quality of life in the aging population: WHOQOL-AGE. Health Qual Life Outcomes. 2013;11:177.

  27. Lanza ST, Dziak JJ, Huang L, Wagner AT, Collins LM. LCA Stata plugin users’ guide (Version 1.2). University Park: The Methodology Center, Penn State; 2015.

  28. Schwarz G. Estimating the dimension of a model. Ann Stat. 1978;6:461–4.

  29. Hurvich CM, Tsai C-L. Regression and time series model selection in small samples. Biometrika. 1989;76:297–307.

  30. Nylund KL, Asparouhov T, Muthén BO. Deciding on the number of classes in latent class analysis and growth mixture modeling: a Monte Carlo simulation study. Struct Equ Model A Multidiscip J. 2007;14:535–69. 

  31. Nagin DS. Group-based modeling of development [internet]. Cambridge, MA: Harvard University Press; 2005. 

  32. Royston P, White I. Multiple imputation by chained equations (MICE): implementation in Stata. J Stat Softw. 2011;45:1–20. 

  33. Hayati Rezvan P, Lee KJ, Simpson JA. The rise of multiple imputation: a review of the reporting and implementation of the method in medical research. BMC Med Res Methodol. 2015;15:30. 

  34. Lunt M. A Guide to Imputing Missing Data with Stata (Revision: 1.4) [Internet]. 2011. Available from:

  35. Royston P, Carlin JB, White IR. Multiple imputation of missing values: new features for mim. Stata J. 2009;9:252–64.

    Google Scholar 

  36. Rubin DB. Multiple imputation for nonresponse in surveys. New York: J. Wiley & Sons; 1987.

    Book  Google Scholar 

  37. Arenillas JF, Moro MA, Dávalos A. The Metabolic Syndrome and Stroke. Stroke. 2007;38(7):2196–203.

  38. Prados-Torres A, Poblador-Plou B, Calderón-Larrañaga A, Gimeno-Feliu LA, González-Rubio F, Poncel-Falcó A, et al. Multimorbidity patterns in primary care: interactions among chronic diseases using factor analysis. PLoS One. 2012;7:e32190. 

  39. Schäfer I, von Leitner E-C, Schön G, Koller D, Hansen H, Kolonko T, et al. Multimorbidity patterns in the elderly: a new approach of disease clustering identifies complex interrelations between chronic conditions. PLoS One. 2010;5:e15941. 

  40. Nussbaumer-Ochsner Y, Rabe KF. Systemic manifestations of COPD. Chest. 2011;139:165–73.

    Article  PubMed  Google Scholar 

  41. Garin N, Olaya B, Lara E, Moneta MV, Miret M, Ayuso-Mateos JL, et al. Visual impairment and multimorbidity in a representative sample of the Spanish population. BMC Public Health. 2014;14:815. 

  42. Falsarella GR, Coimbra IB, Barcelos CC, Costallat LT, Carvalho OM, Coimbra AM. Prevalence and factors associated with rheumatic diseases and chronic joint symptoms in the elderly. Geriatr Gerontol Int. 2013;13:1043–50. 

  43. Nemet AY, Vinker S, Levartovsky S, Kaiserman I. Is cataract associated with cardiovascular morbidity? Eye. 2010;24:1352–8.

  44. Lindblad BE, Håkansson N, Philipson B, Wolk A, Ceriello A, Motz E, et al. Metabolic syndrome components in relation to risk of cataract extraction: a prospective cohort study of women. Ophthalmology. 2008;115:1687–92.

  45. Felton DA. Edentulism and Comorbid factors. J Prosthodont. 2009;18:88–96.

  46. Müllerova H, Agusti A, Erqou S, Mapel DW. Cardiovascular Comorbidity in COPD: systematic literature review. Chest. 2013;144:1163–78.

    Article  PubMed  Google Scholar 

  47. Marengoni A, Angleman S, Melis R, Mangialasche F, Karp A, Garmen A, et al. Aging with multimorbidity: a systematic review of the literature. Ageing Res Rev. 2011;10:430–9.

    Article  PubMed  Google Scholar 

  48. Agur K, McLean G, Hunt K, Guthrie B, Mercer S. How does sex influence multimorbidity? Secondary analysis of a large nationally representative dataset. Int J Environ Res Public Health. 2016;13:391. 

  49. Sousa M, Pereira A, Costa R, Sousa M, Pereira A, Costa R, et al. Subjective memory complaint and depressive symptoms among older adults in Portugal. Curr Gerontol Geriatr Res. 2015;2015:296581.

  50. Ganguli M. Depression, cognitive impairment and dementia: why should clinicians care about the web of causation? Indian J Psychiatry. 2009;51(Suppl 1):S29–34. 

  51. Caracciolo B, Gatz M, Xu W, Marengoni A, Pedersen NL, Fratiglioni L. Relationship of subjective cognitive impairment and cognitive impairment no dementia to chronic disease and multimorbidity in a nation-wide twin study. J Alzheimers Dis. 2013;36:275–84.

  52. Beeri MS, Ravona-Springer R, Silverman JM, Haroutunian V. The effects ofcardiovascular risk factors on cognitive compromise. Dialogues Clin Neurosci. 2009;11:201–12.

  53. Ganguli M, Fu B, Snitz BE, Unverzagt FW, Loewenstein DA, Hughes TF, et al. Vascular risk factors and cognitive decline in a population sample. Alzheimer Dis Assoc Disord. 2014;28:9–15.

  54. Salthouse TA. When does age-related cognitive decline begin? Neurobiol Aging. 2009;30:507–14.

  55. Dal Grande E, Fullerton S, Taylor AW. Reliability of self-reported health risk factors and chronic conditions questions collected using the telephone in South Australia, Australia. BMC Med Res Methodol. 2012;12:108.

  56. Huerta JM, José Tormo M, Egea-Caparrós JM, Ortolá-Devesa JB, Navarro C. Accuracy of self-reported diabetes, hypertension, and Hyperlipidemia in the adult Spanish population. DINO study findings. Rev Española Cardiol. 2009;62:143–52. 

  57. Martin LM, Leff M, Calonge N, Garrett C, Nelson DE. Validation of Self-Reported Chronic Conditions and Health Services in a Managed Care Population. Am J Prev Med. 2000;18(3):215–8.

  58. Sterne JAC, White IR, Carlin JB, Spratt M, Royston P, Kenward MG, et al. Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ. 2009;338:b2393. 

  59. Bull FC, Maslin TS, Armstrong T. Global physical activity questionnaire (GPAQ): nine country reliability and validity study. J Phys Act Health. 2009;6:790–804. 

  60. Cruz-Jentoft AJ, Baeyens JP, Bauer JM, Boirie Y, Cederholm T, Landi F, et al. Sarcopenia: European consensus on definition and diagnosis: report of the European working group on sarcopenia in older people. Age Ageing. 2010;39:412–23. 

  61. Miret M, Caballero FF, Chatterji S, Olaya B, Tobiasz-Adamczyk B, Koskinen S, et al. Health and happiness: cross-sectional household surveys in Finland, Poland and Spain. Bull World Health Organ. 2014;92:716–25.

Download references


Not applicable.


This work was supported by the European Community’s Seventh Framework Programme [grant number 223071-COURAGE Study]; the Instituto de Salud Carlos III-FIS [grant numbers PS09/00295, PS09/01845, PI12/01490, and PI13/00059. Projects PI12/01490 and PI13/00059 have been co-funded by the European Union European Regional Development Fund (ERDF) “A Way to Build Europe”]; the European Union’s Horizon 2020 Research and Innovation Programme [grant agreement 635,316]; and the Centro de Investigación Biomédica en Red de Salud Mental (CIBERSAM), Instituto de Salud Carlos III.

Availability of data and materials

The data analyzed for the current study are not public available but are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



BO, MVM, FFC, ST, JLAM, JM: study design; BO, IB: acquisition of data; BO, MVM, FFC: data analysis; BO, ST, IB, JMH: interpretation of the data; BO, MVM, FFC, ST, IB: preparation of the manuscript; ST, IB, JLAM, JM: critical revision of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Beatriz Olaya.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all participants and ethical approvals were granted by the Ethics Review Committees of Hospital Universitario de la Princesa (Madrid) and Parc Sanitari Sant Joan de Déu (Barcelona).

Consent for publication

Not applicable.

Competing interest

The authors declare that they have no competing interest.

Publisher’s Note

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

Additional files

Additional file 1: Table S1.

Proportion of missingness for indicators of latent classes and variables included in the multiple imputation model [59,60,61]. (DOC 76 kb)

Additional file 2: Table S2.

Association between latent multimorbidity membership with outcomes at baseline and follow-up in the completers (n = 1508). (DOCX 17 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

Olaya, B., Moneta, M.V., Caballero, F. et al. Latent class analysis of multimorbidity patterns and associated outcomes in Spanish older adults: a prospective cohort study. BMC Geriatr 17, 186 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: