Abstract
Background Adults hospitalized for cardiovascular events are at high risk for post-discharge mortality. Hospital-based screening of health-related psychosocial risk factors is now prioritized by the Joint Commission and the National Quality Forum to achieve equitable, high-quality care. We tested our hypothesis that key patient-reported psychosocial and behavioral measures could predict post-hospitalization mortality in a cohort of adults hospitalized for a cardiovascular event.
Methods This was a prospective cohort of adults hospitalized at Vanderbilt University Medical Center. Validated patient-reported measures of health literacy, social support, disease self-management, and socioeconomic status were used as predictors of interest. Cox survival analyses of mortality were conducted over a median 3.5-year follow-up (range: 1.25 – 5.5 years).
Results Among 2,977 adults, 1,874 (63%) were hospitalized for acute coronary syndrome and 1,103 (37%) were hospitalized for acute decompensated heart failure; 60% were male; and the mean age was 53 years. After adjusting for demographic, clinical, and other psychosocial factors, mortality risk was greatest among patients who reported being unable to work due to disability (Hazard Ratio (HR) 2.36, 95% Confidence Interval (CI): 1.73-3.21), who were retired (HR 2.14, 95% CI 1.60-2.87), and who reported unemployment (HR 1.99, 95% CI 1.30-3.06) as compared to those who were employed. Patient-reported measures of disease self-management, perceived health competence and exercise frequency, were also associated with mortality risk after full covariate adjustment (HR 0.86, 95% CI 0.73-1.00 per four-point increase), (HR 0.86, 95% CI 0.77-0.96 per three-day change), respectively.
Conclusions Patient-reported measures of employment status independently predict post-discharge mortality after a cardiac hospitalization. Measure of disease self-management also have prognostic modest utility. Hospital-based screening of psychosocial risk is increasingly prioritized in legislative policy. Incorporating brief, valid measures of employment status and disease self-management factors may help target patients for psychosocial, financial, and rehabilitative resources during post-discharge transitions of care.
Introduction
Mortality after hospitalization for acute coronary events and acute decompensated heart failure is high, and ranges from 12-20%.1-2 Valid, feasible methods to identify patients at greatest risk for mortality after discharge are necessary to improve care transitions and optimally target follow-up support programs for these vulnerable groups. The contribution of medical factors including primary disease status and comorbidities on mortality risk have been long-established, as has the general importance of social and behavioral factors.3 However, only recently have hospitals and health systems begun to collect data on psychosocial and behavioral risk factors for health outcomes to inform patient care.4-5 New regulations from The Joint Commission now require hospitals to collect data on health-related psychosocial risk factors and will increase the availability and potential uses of these data.6
Patient-reported psychosocial and behavioral measures have the potential to identify those at greatest risk of post-discharge adverse outcomes, including mortality.7 Psychosocial measures include socioeconomic status, perceived social support, subjective and functional health literacy and numeracy; behavioral measures include alcohol use and aspects of disease self-management such as perceived health competence and medication adherence. These measures are known to add significant prognostic value to predictive models that estimate mortality risk among individuals living with a malignancy,7 but their prognostic utility in vulnerable populations after an acute hospitalization is less well-studied. This gap in knowledge limits the health system’s ability to optimally target resources to improve transitions of care after an acute hospitalization.
Informed by a thorough literature review of psychosocial and behavioral measures that predict mortality,8-20 we tested the hypothesis that patient-reported health literacy, social support, disease self-management, and socioeconomic status would independently predict mortality after hospitalization for acute coronary syndrome or acute decompensated heart failure, after adjustment for demographic and clinical factors.
Methods
Study design and sample
This analysis is part of the National Heart Lung and Blood Institute-funded Vanderbilt Inpatient Cohort Study (VICS), aimed to identify novel, modifiable factors that predict adverse outcomes during the post-hospitalization period.21 Briefly, VICS was a prospective cohort that included 2,977 adults hospitalized for acute decompensated heart failure or acute coronary syndromes at Vanderbilt University Hospital in Nashville, Tennessee between 2011 and 2015. The cohort is unique in its administration of a wide range of validated patient-reported measures at cohort enrollment. Eligibility was determined from medical record review conducted by a physician and research staff. Criteria for exclusion were severe cognitive impairment, unstable psychiatric illness, inability to communicate in English, current hospice status, and unstable contact information for follow-up. For this analysis, we excluded patients who died during the index hospitalization and therefore did not enter the post-discharge follow-up period. The study was approved by Vanderbilt’s Institutional Review Board, and written informed consent was obtained from enrollees.
Demographic factors collected included age, gender, years of maximum educational attainment, household income, self-reported race, and marital status. Number of hospitalizations during the prior 12 months and a validated score of total clinical comorbidity burden (the Elixhauser Comorbidity Index)22 were also obtained. The Elixhauser Comorbidity Index includes whether a participant has congestive heart failure, cardiac arrythmias, cardiac valvular disease, pulmonary hypertension, peripheral vascular disease, hypertension, paralysis, neurodegerative disordres, chronic obstructive pulmonary disease, diabetes, hypothyroidism, kidney disease, liver disease, peptic ulcer disease, acquired immunodeficiency symdrome, lymphoma, metastatic or solid malignancies, collagen vascular disease, coagulopathy, obesity, weight loss, fluid and electrolyte disorders, anemia, evidence of alcohol or illicit substance abuse, depression, or psychosis. Cause of cardiovascular hospitalization (acute coronary syndrome vs. acute decompensated heart failure) were ascertained from the medical record.
Exposures/predictors of interest
We identified validated, patient-reported measures of both functional and subjective health literacy, subjective numeracy, social support, and health behaviors as worthy of further investigation for prognostic utility in our sample. Given the known significant associations between self-efficacy and health outcomes, including mortality, across multiple chronic conditions23-24 we also tested whether a validated measure of health self-efficacy, perceived health competence, was associated with mortality in our sample. Given that health behaviors and emotion management are encompassed by a construct known in the literature as disease self-management,25 we grouped patients’ reports of health behaviors, coping, and perceived health competence into a ‘Disease Self-Management’ domain.
Functional health literacy was measured using the short form of the Test of Functional Health Literacy in Adults (sTOFHLA; scores 0-36 with higher scores indicating greater ability to interpret health-related text). This test categorizes functional health literacy as inadequate (Score 0-16), marginal (score 17-22), or adequate (scores 23-36).26 Subjective health literacy was assessed using the three-item Brief Health Literacy Screen (BHLS; scores 3-15 with higher scores indicating higher confidence in interpreting written and spoken health-related material).27 Numeracy was measured using a three-item version of the Subjective Numeracy Scale (SNS; scores 1-6 with higher scores indicating greater comfort with interpreting numbers).28
Social support was assessed using multiple measures, each of which identified a different aspect of social support. These included items from the Health and Retirement Survey and items from the Midlife Development in the United States (MIDUS) survey that measured frequency of contact and level of support from friends, family, and neighbors.29-30 We also added the Enhancing Recovery in Coronary Heart Disease Emotional Social Support Instrument (ESSI; scores 8-34 with higher scores indicating greater perceived emotional support).31
Disease self-management was assessed using patients’ frequency in engaging in health behaviors and ability to manage emotional stressors associated with living with a health condition.25 Perceived health competence was assessed using a validated, two item abbreviation of the original, eight-item Perceived Health Competence Scale, a measure of individuals’ self-confidence in managing health and achieving health goals (PHCS-2; scores 2-10 with higher scores indicating higher perceived health competence).20,32 Resilient coping behaviors, the ability to rebound from or positively adapt to significant stressors, was measured by the validated four-item Brief Resilient Coping Scale (BRCS; scores 4-20 with higher scores indicating a greater ability to cope with stress).33 Medication adherence before hospitalization was measured using a seven-item version of the Adherence to Refills and Medicines Scale (ARMS, scores 7-28 with higher scores indicating higher self-reported medication adherence).34 Tobacco and alcohol use, dietary habits, and weekly exercise frequency were measured using the Centers for Disease Control Behavioral Risk Factor Surveillance System questions, the Starting the Conversation scale, and the Exercise Vital sign scale, respectively.35-37
Socioeconomic status was assessed by self-reported annual household income, highest level of education attained, employment status (employed, not employed for wages, retired, or unemployed due to a disability), and degree of difficulty paying bills.
Outcome variable
Mortality was ascertained using the Social Security Administration’s Death Master File; information from patient charts recorded after clinical contact with family members; and obituaries. The Death Master File contains over 85 million records since 1936 and is linked to patient records by Social Security Number. The Death Master File provides accurate matches for death (in excess of 90% among American born individuals), but has had considerable omissions, particularly since 2011.38-39 We therefore supplemented Death Master File data with clinical records and obituaries. Mortality follow-up was for a minimum of 1.25 years for the last patient enrolled (December 2015), and up to 5.5 years for the initial enrollees (October 2011). Therefore, any variability in followup was not due to loss to followup but due to administrative censoring.
Statistical analyses
We summarized categorical variables with counts and percentages, and continuous variables with the median (25th, 75th) percentiles. For the primary analyses, we fitted multivariable Cox proportional hazards regression models to investigate the associations between clinical and patient-reported risk factors and time to death.
We defined six domains of variables: 1) demographics and reason for hospitalization (age, gender, race, and acute coronary syndrome vs acute decompensated heart failure), 2) health literacy and numeracy (BHLS, sTOFHLA, SNS), 3) social support (relationship status, number of close family and/or friends, ESSI), 4) disease self-management (PHCS-2, BRCS, ARMS, alcohol use, tobacco use, exercise, dietary habits), 5) socioeconomic status (income, education, employment, difficulty paying bills) and 6) health status (Elixhauser comorbidity score, hospitalizations in prior 12 months). We performed redundancy analysis to examine the potential for collinearity among all risk factors. The highest coefficients of determination corresponded to the income and age values, with values of 0.59 and 0.50, respectively. All other coefficients of determination were 0.4 or lower, and no factors were removed from consideration.
To examine the association of mortality rates and risk factors or domains of risk factors, we implemented a hierarchical modeling approach. We first fit a Base Model that only included variables in the demographics and diagnosis domain. We then added all variables in each of the five other domains (one domain at a time) to the Base Model to examine their associations while adjusting for the Base Model variables but without adjusting for all other domains. Finally, we fit a full model with all variables from all domains simultaneously.
We reported hazard ratios and 95% confidence intervals corresponding to interquartile range changes in continuous variables. To test the association between each domain and mortality rates we used likelihood ratio tests (LRT) with degrees of freedom equal to the number of parameters estimated for the domain. For the domain-specific models, the LRT compares the Base Model to a model that adds a single domain of variables to the Base Model. For the full model results, the LRT compares the full model to a model that excludes all variables within in a single domain. Additionally, in the full model, to quantify the contribution of each risk factor and domain, we calculated a measure of relative explained variation (REV) by dividing the likelihood ratio Chi-square statistic associated with each variable or domain by the total likelihood ratio Chi-square of the full model.
We summarized the missingness rate of each variable and conducted multiple imputation with chained equation, with five imputated datasets using a predictive mean matching algorithm. We conducted a single global test for the proportional hazards assumption across all variables. Because there was not a clear violation of the proportional hazards assumption (p=0.06), we chose to assume proportional hazards for all variables. We used a two-sided p < 0.05 significance level and conducted all analyses in R version 4.3.0 (R Foundation for Statistical Computing, Vienna, Austria).40-42
Results
Our sample consisted of 2,977 patients who survived to discharge; 1,874 (63%) were hospitalized for acute coronary syndrome and 1,103 (37%) were hospitalized for acute decompensated heart failure. The median age was 61, 60% were men, and 85% identified as white. The median years of education was 13, 22% had an income below $20,000, 35% were retired, and 25% reported being unable to work due to a disability (Table 1).
During the follow-up period, 616 patients died. In descriptive analyses (Table 1), we observed that patients who died tended to be older, admitted for a heart failure exacerbation, had poorer health literacy, consumed less alcohol, exercised on fewer days, had lower income, were less likely to be employed, were more likely to be unable to work due to disability, had more hospital admissions in the last year, and had higher Elixhauser comorbidity scores.
Table 2 identifies patient-reported factors predictive of mortality in each of our analyses. The Base Model corresponds to the simple model that includes only the Demographic and Study Diagnosis domain variables. The Domain-Specific Models columns includes results from models that add each of the domains (separately) to the Base Model. Finally, the Full Model columns includes the final models with all domain variables included simultaneously. Each domain (Health literacy and Numeracy, Social support, Disease self-management, and Socioeconomic status) was found to be statistically significantly associated with mortality (p<0.01) when individually added to the Base Model (middle columns of Table 2).
The fully adjusted Cox proportional hazard models contained all domains simultaneously, as well as variables in the Base Model (the last columns of Table 2). Its estimated C-index was 0.79 consistent with moderate discrimination performance. The global, likelihood ratio chi-square statistic was 652.63, and on 27 degrees of freedom; this is consistent with minimal model overfitting (i.e., 4% shrinkage).43
The socioeconomic status and Health status domains were strongly associated with mortality (p<0.001) (Table 2). The Disease self-management domain was marginally associated with mortality (p=0.03) while the Health literacy/Numeracy and Social support domains were not. Within the Disease self-management and Socioeconomic Status domains, a few individual variables showed associations with mortality. A four point increase in the Perceived Health Competence Score was associated with a 14 percent decrease in hazard for mortality (Hazard Ratio (HR): 0.86, 95% Confidence Interval (CI): 0.73-1.00) in the Full Model. Similarly a three-day increase in the number of exercise days was associated with a 14 percent decrease in hazard for mortality (HR: 0.86, 95% CI: 0.77-0.96) in the Full Model. Employment status was strongly associated with mortality risk in the Full Model. Hazard ratios for those who were not employed, retired, and unable to work due to disability were 1.99, (95% CI 1.30-3.06), 2.14 (95% CI 1.60-2.87) and 2.36 (95% CI 1.73-3.21), respectively, when compared to those who were employed.
Figures 1a and 1b show relative contributions of the individual variables and domains to the information (i.e., variability in mortality explained) contained in the full model. We observed that reason for cardiac hospitalization and health status were the strongest contributors to model information. However, eight percent of the total variability explained in the model was due to the socioeconomic status construct. Further, of all variables, employment status contributed the third most information to the model.
These figures display individual variable and domain contributions to the explained variability in the full model. We sorted the contribution by the likelihood ratio chi-square statistics corrected for the degrees of freedom and reported the relative explained variation (REV). REV was calculated using the Chi-square statistic associated with each variable or domain divided by the total Chi-square of the full model. It captures the proportion of information lost from the model if a single variable (Figure 1a) or if a domain of variables (Figure 1b) is removed from the full model.
Discussion
To our knowledge, ours is the first large cohort that has assessed the mortality risk associated with a wide range of validated patient-reported measures in a high-risk cardiovascular population. Our results support that patient-reported measures of socioeconomic status and disease self-management may independently provide prognostic utility for adults hospitalized for cardiovascular events.
The addition of patient-reported measures to clinical data improves mortality risk prediction, though much of the literature thus far has focused on patients with malignancy. In an analysis of 3,240 men diagnosed with prostate cancer, the ten-year area under the received operating curve improved from 0.721 to 0.812 when patient-reported health measures were added.44 Cardiovascular hospitalizations are common and associated with increased morbidity and mortality, but few studies have determined the prognostic utility of patient-reported measures among adults hospitalized for these events to support successful transitions post-discharge.
In our analysis, employment status was the measure of socioeconomic status with the greatest prognostic significance. Patients were specifically asked whether they were unemployed due to disability. The significant, independent associations between being unemployed due to disability and mortality are striking. While our questionnaires did not ask patients to specify the type of disability they were experiencing, disability can often include physical function decline. Prior evidence demonstrates that physical function trajectory independently predicts cardiovascular and all-cause mortality, and we are currently investigating the prognostic utility measures of physical function (physical frailty and physical vulnerability) in populations with cardiovascular disease using validated measures.45-46 While our study design precludes our ability to identify the precise biological and behavioral mechanisms by which being disabled and/or unemployed contribute to mortality, our work further supports the utility of collecting these measures among adults hospitalized for cardiovascular events to identify those who may need additional post-discharge vocational or rehabilitative resources.
In our cohort, functional health literacy and exercise frequency were each individually predictive of mortality. We and others have previously demonstrated that functional health literacy is independently predictive of mortality among adults hospitalized for cardiovascular events, presumably through an inability to manage the complex medication prescriptions and dietary requirements associated with cardiovascular disease.17,47 Similarly, it has long been established that patient-reported exercise frequency is predictive of mortality among adults with increased risk for cardiovascular disease, though the precise amount of leisure time physical activity needed to confer a mortality benefit remains a source of ongoing investigation.48 However, our cohort added to these data by using a wide range of validated measures of disease self-management, followed patients for an extended period of time, and adjusted for numerous demographic and clinical characteristics.
Our analysis also included perceived health competence, an under-investigated construct that measures an individual’s sense of self-confidence in managing health conditions. In our cohort, perceived health competence was still marginally associated with long-term mortality even after adjustment for demographics, clinical comorbidities, and other patient-reported factors. Validity and reliability of the Perceived Health Competence Scale, which measures perceived health competence, has been established using adults living with a chronic illness.20,32 Given that the validated, abbreviated version applied in this cohort contains just two questions, incorporating these questions as part of a routine patient-facing assessment during hospital admission may be a novel, efficient way to identify those patients who need psychosocial support during transitions of care to reduce post-cardiac hospitalization mortality.
This study has several strengths. Few large cohorts exist that have included the type and range of validated patient-reported measures as has been included in our sample. In focusing on patient-reported measures, the study identifies factors that can alter the post-hospitalization trajectory of a group already at high risk for adverse outcomes. An additional strength is that long-term mortality data were obtained using thorough and rigorous methods, leading to robust outcome data.
Our study also has key limitations. Non-English speaking participants were excluded due to the fact that many of the patient-reported surveys used have only been validated in English. While this study was conducted at one academic medical center, the center has a very large, three-state catchement area. Those individuals who report unemployment due to disability may represent a group that also has a high comorbidity burden. More than likely, providers must aggressively treat these patients’ clinical comorbidities in addition to refer them to psychosocial or rehabilitative interventions to reduce mortality. Other limitations include that most adults in our cohort identified as white, this reducing the generalizability of our results across within other racial and ethnic groups. Finally, our patient-reported measures were only collected at baseline, limiting our ability to test for changes in these measures over time and their impact on outcomes.
Conclusion
In a large cohort of adults hospitalized for a cardiovascular event, patient-reported measures of socioeconomic status, particularly employment, significantly and independently predicted post-discharge mortality at a median of 3.5 years after adjustment for demographic and clinical factors. Patient-reported measures of disease self-management; in particular exercise frequency and perceived health competence, were modestly associated with mortality after covariate adjustment. As collecting a large array of patient-reported measures takes time, focusing on these particular measures could be an efficient way to identify those patients who should be targeted for psychosocial and rehabilitative resources during the post-hospitalization period.
Funding
This study was supported by grants R01 HL109388 (Kripalani) and 1K23 DK1297742 (Nair). The content is solely the responsibility of the authors. The funding agencies were not involved in the design and conduct of the study; collection, management, analysis, and interpretation of the data; and preparation, review, or approval of the manuscript. The contents do not necessarily represent the official views of the NIH.
Acknowledgements
We thank our study research analysts and project coordinators: Courtney Cawthon, Catherine Couey, Monika Rizk, Hannah Rosenberg, Daniel Lewis, Blake Hendrickson, Olivia Dozier, Vanessa Fuentes, Cardella Leak, Mary Lou Jacobsen, Catherine Evans, Joanna Lee, Emily Lucianno, and Erin Acord.
Data Availability
All analytical data will be available upon request.