The Association of Shared Care Networks With 30-Day Heart Failure Excessive Hospital Readmissions: Longitudinal Observational Study

Background Higher-than-expected heart failure (HF) readmissions affect half of US hospitals every year. The Hospital Reduction Readmission Program has reduced risk-adjusted readmissions, but it has also produced unintended consequences. Shared care models have been advocated for HF care, but the association of shared care networks with HF readmissions has never been investigated. Objective This study aims to evaluate the association of shared care networks with 30-day HF excessive readmission rates using a longitudinal observational study. Methods We curated publicly available data on hospital discharges and HF excessive readmission ratios from hospitals in California between 2012 and 2017. Shared care areas were delineated as data-driven units of care coordination emerging from discharge networks. The localization index, the proportion of patients who reside in the same shared care area in which they are admitted, was calculated by year. Generalized estimating equations were used to evaluate the association between the localization index and the excessive readmission ratio of hospitals controlling for race/ethnicity and socioeconomic factors. Results A total of 300 hospitals in California in a 6-year period were included. The HF excessive readmission ratio was negatively associated with the adjusted localization index (β=–.0474, 95% CI –0.082 to –0.013). The percentage of Black residents within the shared care areas was the only statistically significant covariate (β=.4128, 95% CI 0.302 to 0.524). Conclusions Higher-than-expected HF readmissions were associated with shared care networks. Control mechanisms such as the Hospital Reduction Readmission Program may need to characterize and reward shared care to guide hospitals toward a more organized HF care system.


Introduction
Higher-than-expected heart failure (HF) readmission impacts approximately half of US hospitals every year, and almost every hospital has experienced it at least once in the period between 2012 and 2017.By 2030, HF is projected to affect at least 8 million people in the United States, with an incidence of 21 per 1000 people older than 65 years and an estimated cost of US $69.8 billion [1].The number of patients with HF receiving HF care and requiring advanced HF therapies such as left ventricular assisted devices (LVADs) will increase exponentially [2].Addressing higher-than-expected HF readmissions for patients with HF is needed as demand increases, with the aging population requiring improved care coordination mechanisms that promote a more organized HF care system [3].
HF is managed through a complex system that serves both affluent and vulnerable patient populations, and encompasses nonlinear interactions among primary care, general cardiology, specialized HF clinics, and tertiary and quaternary centers.The implementation of any control mechanism can produce unintended consequences if the complexity of the HF care system is not taken into consideration [4,5].Systemwide control programs such as the Hospital Reduction Readmission Program (HRRP) [6] may be a first step toward organizing the HF care system.Nonetheless, they will continue to create unintended consequences and penalize hospitals for factors beyond their control [7] unless these programs specifically foster care coordination mechanisms capable of promoting organization for HF care's complex system.Shared care integrates primary, secondary, and tertiary levels of care [8], and has been advocated as a necessary model to promote a more organized HF care system [9] such as the spoke-hub-and-node model [10].Shared care has been studied among chronic diseases [11], including HF [12], but only recently has it been advocated for by international working groups as a way to organize HF care [9], particularly among patients with advanced HF [10] such as patients with LVAD support [13].Shared care areas (SCAs) are data-driven units of care coordination captured from large-scale data on hospital discharges to patient residencies, and SCAs may explain variation in medical adherence to HF guideline-directed medical therapy [14].The localization index (LI) of an SCA is the proportion of patients who reside in the same SCA they are admitted and is a measure of local care coordination commonly used to evaluate SCAs [15].This study aims to evaluate the longitudinal association between higher-than-expected HF readmissions and the LI of SCAs both unadjusted and adjusted for racial/ethnic and socioeconomic factors.

Methods
This methods section was written according to the STROBE (Strengthening the Reporting of Observational Studies in Epidemiology) standard of writing.

Study Design, Study Setting, and Participants
This is an observational longitudinal study.All data used in this study are made publicly available by the HRRP and Office of Statewide Health Planning and Development (OSHPD).The study setting was hospitals in California during the period from 2012 to 2017.Participants were all in hospitals reported in the HRRP [6].The eligibility criteria were as follows: at least 2 repeated measures of higher-than-expected HF readmission in the HRRP and availability of discharge data from the OSHPD [16].These criteria enabled carrying out a longitudinal study that requires repeated measures and linking data from the HRRP with date from OSHPD.Between 233 and 237 hospitals in California were included depending on the year.Ethical approval was unnecessary because all data were at the hospital level and already made publicly available from both HRRP and OSHPD.All code, processed data, built networks, and data analysis resulting from this study are available on the Open Science Framework repository for this study [17].

Study Outcome
The main study outcome was hospital excessive readmission ratio (ERR), which is a risk-standardized 30-day readmission ratio that adjusts for a set of patient-specific covariates such as congestive HF, renal failure, and chronic obstructive pulmonary disease [18].It is used by the HRRP to assess excess hospital readmissions and calculate hospital penalties [6].The ERR is calculated by dividing the predicted readmissions by the expected readmissions.Using a hierarchical generalized linear model, both predicted and expected readmissions are estimated using an adjusted average intercept over all hospitals, but predicted readmissions, in addition, are estimated using a hospital-specific intercept deviation from the adjusted average intercept over all hospitals.Such methodology, well documented in the Condition-Specific Readmission Measures Updates and Specifications Report from the Centers for Medicare & Medicaid Services [18], makes the ERR an appropriate instrument for comparing hospitals within and between years.

Study Variables
The main study variable was the LI, which represents the proportion of patient discharges from hospitals within the same SCA of which these patients live [19,20].A higher LI represents a homogenous SCA with localized care coordination (ie, patients tend to receive care where they live).Other study variables were

Data Sources
The ERR data used in this study was made publicly available from the HRRP [6].The ERR data of each year in the period from 2012 to 2017 (ie, fiscal year 2014 and 2019) was separately downloaded from HRRP and compiled into a single file.The Patient Origin/Market Share data was made publicly available from the OSHPD [16].Patient Origin/Market Share data are aggregated numbers of emergency department (ED) discharges among zip codes of hospitals and patient residencies.Zip Codes were converted to the Zip Code Tabulation Areas (ZCTAs) using the Zip Code to ZCTA Crosswalk made publicly available by the Uniform Data System [22].Demographic data was gathered for the state of California from the American Community Survey [21].

Uncovering Shared Care Areas and Localization Index From Hospital-Patient Discharge Data
Six yearly hospital-patient discharge networks were built from OSHPD hospital-patient ED discharges between 2012 to 2017.In a hospital-patient discharge network [15], a node is the ZCTA of a hospital or patient residency, and the link between two nodes (ie, ZCTAs) is the total number of ED discharges.For each yearly hospital-patient discharge network, SCAs were delineated using community detection algorithms.Each delineated SCA consists of a set of ZCTAs in which hospitals are embedded.A set of four diverse community detection algorithms were considered to decrease both variability and bias [23].The algorithms were Louvain [24] with resolution equal to 1, Stochastic Block Model [25,26] with degree corrected, Infomap [27] with two levels, and Speaker-Listener Label Propagation [28] with postprocessing threshold equal to 0.5

Statistical Analysis
The ERR hospitals and the LI of SCAs were integrated at each year by linking the ZCTAs of hospitals and SCAs (Table S1 and Figure S1 in Multimedia Appendix 1).A longitudinal regression was specified in which the dependent variable ERR of a hospital at time t as a function of the LI of its SCA at time t.We used a generalized estimating equation (GEE) using a Gaussian family and an exchangeable working correlation structure to account for multiple observations of ERR from the same hospital across years and SCAs [29].The estimated regression coefficients (beta) were used to measure unadjusted associations between the dependent and independent variables, and adjusted associations after controlling for racial/ethnic and socioeconomic confounders associated with HF readmission at the regional level [30].The GEE was estimated using the Statsmodels Python package [31].Additionally, hospitals were stratified based on quartiles of the LI and all covariates that were found statistically significant, and median values of ERRs and percentage of hospitals penalized were calculated for each quartile (Q1, Q2, Q3, Q4).We estimated 95% CIs using 10,000 bootstrap samples with replacement from each quartile, the estimation of CIs for medians using the Bootstrapped Python package [32].

Predicting Higher-Than-Expected Heart Failure Readmissions for Changes in Localization Index
The estimated GEE model was used to predict HF's ERRs assuming a range of changes in the LI in SCAs with distinct percentages of Black residents, the only statistically significant covariate.The differences in the LI between subsequent years were calculated for all hospitals.The 25th, 50th, and 75th percentiles were separately calculated for both positive (+q1, +q2, and +q3) and negative (-q1, -q2, and -q3) differences.The SCAs were stratified by quartiles of Black residents (Q1, Q2, Q3, and Q4).The ERR was predicted using the GEE model after each positive and negative percentile difference in the LI was applied to the stratified SCA data.

Descriptive Statistics of Heart Failure Hospital Readmissions in the United States and California
The ERR is calculated every year by the HRRP for the approximately 2700 to 2900 hospitals in the United States, from which 233 to 237 hospitals are from California (Table 1).Overall, approximately half of US hospitals are penalized, and this percentage has not changed during the study period between 2012 to 2017.The ERR (and the percentage of hospitals penalized) of US hospitals have remained approximately constant during the study period, from 1.0013 (49.76%) in 2012 to 1.0016 (48.94%) in 2017.The ERR (and the percentage of hospitals penalized) of hospitals in California increased from 0.9914 (49.36%) to 1.0087 (56.12%).In 2017, the percentage of hospitals penalized in California (56.12%, 95% CI 49.75%-62.29%) is slightly higher than that among all hospitals in the United States (48.91%, 95% CI 47.06%-50.76%).Although not statistically significant, the ERR SD appears to be decreasing over the years.

Association of the Excessive Readmission Ratio and Localization Index
The results of the quartile analysis indicate that the ERR of hospitals was negatively associated with the LI (Table 2) as well as with the percentage of Black residents (Table 3).In 2017, for instance, the ERR of hospitals in SCAs with the lowest quartile (Q1) of the LI was 1.03 (95% CI 1.02-1.04)with 65.7% (95% CI 59.4%-72.0%) of hospitals penalized.In SCAs with the highest quartile (Q4) of the LI; however, the median ERR was 0.98 (95% CI 0.97-0.99)with only 43.1% (95% CI 35.3%-51.0%) of hospitals penalized.From 2012 to 2017, the disparities between the ERR and percentage of hospitals penalized among SCAs belonging to the lowest (Q1) and highest LI (Q4) quartiles has increased mainly because of increases in the ERR and percentage of hospitals penalized within SCAs in the lowest LI quartile (Q1).Similarly, in 2017, the ERR of hospitals in SCAs with the lowest quartile (Q1) of Black residents was 0.99 (95% CI 0.98-1.0)with 45.2% (95% CI 38.2%-52.2%) of hospitals penalized.In SCAs with the highest percentage of Black residents quartile (Q4), however, the median ERR was 1.03 (95% CI 1.02-1.04)with 67.6% (95% CI 60.7%-74.6%) of hospitals penalized.The percentage of Black residents is slightly higher in SCAs with lower localization (Table S4 in Multimedia Appendix 1).The results of the regression analysis (Figure 1 and Table 4) indicate that the ERR of hospitals was negatively associated with the adjusted and unadjusted LI of their SCAs (eg, ERRs were lower when hospitals were located in SCAs where more patients received care close to where they resided) according to both unadjusted (β=-.0717;P<.001) and adjusted (β=-.0495;P=.049) coefficients when the regression was controlled for racial/ethnic and socioeconomic covariates.The percentage of Black residents in the SCA was the only covariate with a statistically significant association according to the regression coefficient (β=.3892;P<.001).The results can be separately analyzed for each community detection algorithm (Table S3, Multimedia Appendix 1), and the Stochastic Block Model uncovered SCAs with the LI anomalously lower and was not considered in the final analysis.The results can be separately analyzed for each community detection algorithm for ERR (Table S5 in Multimedia Appendix 1), percentage of hospitals penalized (Table S6 in Multimedia Appendix 1), and the percentage of Black residents (Table S7 in Multimedia Appendix 1).

Principal Findings
Regional variation in health care delivery is a ubiquitous phenomenon [3,19], and the HRRP may have differently impacted almost 3000 US hospitals depending on their state.The main finding in this study is that higher-than-expected HF hospital readmissions are associated with the share care networks in which hospitals are embedded.Specifically, hospitals within SCAs with a high LI are associated with lower ERRs than hospitals within SCAs with lower LIs.The LI represents the proportion of patient discharges from hospitals within the same SCA of which these patients live.The LI is widely used as a measure of care coordination and unwarranted health care variation [4,19], but to our knowledge, this is the first documentation of its association with HF higher-than-expected readmissions.In this study, the LI is ultimately derived from the shared care discharge networks.In SCAs with a high LI, discharges are localized with a lower proportion of discharges of patients from other SCAs.Not only has shared care been advocated as an appropriate model to organize HF care [9,10], but partnerships among community physicians and local hospitals have been identified as hospital strategies to reduce 30-day HF readmission [33].Characterizing shared care networks provides a road map for hospitals to work together, improving their shared care network as a whole instead of focusing on their hospital penalties.
Though the HRRP is a nationwide effort to reduce higher-than-expected hospital readmissions, it has also created unintended consequences in the complex system of HF care by penalizing hospitals for issues beyond their control, leaving them without specific guidance on how to improve and focusing on punishment instead of process improvements [7].Patients with HF should be managed as a continuum of care within the primary, secondary, and tertiary level of care, promoting timely patient referrals and delivering care within a strong working relationship [9].Integrated HF care will improve care coordination that influences patient outcomes.The features identified that result in improved shared care include liaisons between levels of care and institutions, shared professional education, and medication optimization.Comprehensive pathways across primary, secondary, and tertiary care and institutions should be developed and implemented considering patients and health care providers in the design of these pathways [34].

RenderX
The association of ERRs with shared care networks, however, seems to vary depending on the ethnic/racial and socioeconomic composition of SCAs.In this study, ERR is positively associated with the percentage of Black residents in the SCA.Ethnic/racial disparities may contribute to HF hospital readmissions [20,30,33,35], and HF readmission rates are consistently higher for Black patients [35][36][37].In a previous case-control study [30], after matching maximum penalty hospitals as cases to their nearest nonpenalty hospitals as controls, the authors found that maximum penalty hospitals were more likely than controls to be in counties with low socioeconomic status.
The regional variation on the impact of the HRRP raises the following question: how much HF higher-than-expected readmissions are related to hospital-specific performance, and how much it is related to issues beyond the control of a hospital?Additionally, the increased association of the ERR with the LI in SCAs with increasingly higher percentages of Black residents raises the following question: how can improved shared care networks reduce HF disparities among underserved and marginalized groups?Our findings will hopefully motivate cluster randomized clinical trials [38] to evaluate how improved shared care models will reduce hospital readmissions and overall costs, increase adherence to guideline-directed medical therapy, and improve clinical outcomes such as survival and development of chronic conditions.

Limitations
The HRRP is a nationwide program, but our study only considered hospitals in California because large-scale hospital-specific discharge data at the ZCTA level is not publicly available to examine all US hospitals.Our finding only applies to higher-than-expected HF readmissions, and the generalization to conditions other than HF (eg, acute myocardial infarction, pneumonia, and chronic obstructive pulmonary disease) will require further investigation.The primary outcome used in our study, the ERR, is a ratio between two hospital-level regressions that can be used across heterogeneous hospitals but has little inherent variability.In its current version, our study neglects to model the interactions between SCAs, which deserves further investigation.Although our study assumes that the ERR can be used to compare different hospitals as it accounts for a plethora of factors associated with the hospital-level HF readmissions at the individual level, our findings should be interpreted at the hospital level.

Conclusions
Shared care models have been advocated for in HF care but have not been explicitly characterized and rewarded by nationwide control programs such as the HRRP or health systems.In this study, we evaluated the association of higher-than-expected HF readmissions with shared care networks by curating publicly available large-scale hospital-level data on HF ERRs from Medicare HRRP as well as hospital-patient discharges from OSHPD.HF ERRs of hospitals were associated with the LI of the SCAs in which they were embedded, even after controlling for socioeconomic disparities.The HRRP, health systems, and hospitals should characterize and reward models of shared care practices for promoting the necessary integration capable of producing a sustainable and equitable HF care system.The higher-than-expected HF readmission of hospitals was associated with the shared care networks in which hospitals were embedded and the ethnic/racial composition of their SCAs.Hospitals should collectively work to systematically improve their shared care networks for improved HF care.
Improved shared care networks of HF care could mitigate higher-than-expected HF readmissions, especially among underserved and marginalized groups, and translate into economic benefits.Implementation of this model will require collaboration between providers and hospital administrations.Future clinical trials will be needed to evaluate the impact of systematic implementation of improved shared care models of HF to improve higher-than-expected HF readmissions.

Figure 1 .
Figure 1.Central illustration: association of heart failure excessive readmission with shared care networks.Hospitals are embedded in shared care areas (SCAs), which are data-driven units of care coordination emerging from the discharge networks among hospitals.The localization index is the proportion of patient discharges from hospitals within the same SCA in which these patients live.The heart failure ERRs of hospitals are associated with the SCA localization index in which they are embedded.

Table 1 .
Descriptive statistics of excessive readmission ratio (ERR) and percentage of hospitals penalized in the United States and California.

Table 2 .
Excessive readmission ratios (ERRs) for hospitals in California by the localization index (LI) quartile.

Table 3 .
Excessive readmission ratios (ERRs) for hospitals in California by percentage of Black residents in the shared care area.CIs estimated by 10,000 bootstrap samples with replacement.

Table 4 .
Results of the generalized estimating equations regression for excessive readmission ratios.

Table 5 .
Predictions of excessive readmission ratios (ERRs) and percentage of hospitals penalized based on changes in localization index (LI).