Skip to main content

Post-acute care referral in United States of America: a multiregional study of factors associated with referral destination in a cohort of patients with coronary artery bypass graft or valve replacement



The use of post-acute care (PAC) for cardiovascular conditions is highly variable across geographical regions. Although PAC benefits include lower readmission rates, better clinical outcomes, and lower mortality, referral patterns vary widely, raising concerns about substandard care and inflated costs. The objective of this study is to identify factors associated with PAC referral decisions at acute care discharge.


This study is a retrospective Electronic Health Records (EHR) based review of a cohort of patients with coronary artery bypass graft (CABG) and valve replacement (VR). EHR records were extracted from the Cerner Health-Facts Data warehouse and covered 49 hospitals in the United States of America (U.S.) from January 2010 to December 2015. Multinomial logistic regression was used to identify associations of 29 variables comprising patient characteristics, hospital profiles, and patient conditions at discharge.


The cohort had 14,224 patients with mean age 63.5 years, with 10,234 (71.9%) male and 11,946 (84%) Caucasian, with 5827 (40.96%) being discharged to home without additional care (Home), 5226 (36.74%) to home health care (HHC), 1721 (12.10%) to skilled nursing facilities (SNF), 1168 (8.22%) to inpatient rehabilitation facilities (IRF), 164 (1.15%) to long term care hospitals (LTCH), and 118 (0.83%) to other locations. Census division, hospital size, teaching hospital status, gender, age, marital status, length of stay, and Charlson comorbidity index were identified as highly significant variables (p- values < 0.001) that influence the PAC referral decision. Overall model accuracy was 62.6%, and multiclass Area Under the Curve (AUC) values were for Home: 0.72; HHC: 0.72; SNF: 0.58; IRF: 0.53; LTCH: 0.52, and others: 0.46.


Census location of the acute care hospital was highly associated with PAC referral practices, as was hospital capacity, with larger hospitals referring patients to PAC at a greater rate than smaller hospitals. Race and gender were also statistically significant, with Asians, Hispanics, and Native Americans being less likely to be referred to PAC compared to Caucasians, and female patients being more likely to be referred than males. Additional analysis indicated that PAC referral practices are also influenced by the mix of PAC services offered in each region.

Peer Review reports


Post-acute care (PAC) facilities provide treatment for acute-care patients following hospital discharge and are known to improve patient outcomes, readmission rates, mortality, and functional disability [1, 2]. Their usage has grown over 80% since 1996 [3], causing the U.S. Medicare’s annual PAC spending to double since 2001 [4]. Medicare spending on PAC for heart attack, congestive heart failure, and hip fracture grew 4.5–8.5% annually from 1994 to 2009, as compared to a growth of 1.5–2% per year for total spending in the U.S. [4]. Recent reports demonstrated PAC to be the largest contributor (40%) to Medicare spending variation among geographic regions [5]. For example, in 2013, Medicare spent one of every six dollars (about $60 billion) on PAC reimbursement [6].

Conditions frequently requiring PAC referral upon acute care discharge include respiratory failure, stroke, joint replacement, cardiac surgery, heart failure, and pneumonia. Services are provided to these patients through various settings, which include long-term acute care hospitals (LTCH, 428 facilities in the U.S.), inpatient rehabilitation facilities (IRF, 1165 in the U.S.), skilled nursing facilities (SNF, 16,000 in U.S.), and Certified Home Health Agencies (HHA, 33,000 in U.S.) [7]. Although these facilities play an essential role in improving acute-care patient outcomes, they are somewhat heterogeneous, poorly coordinated with acute-care hospitals, and exhibit high regional variations in usage and availability [8]. Overall, PAC is perhaps the least understood portion of the U.S. healthcare continuum, and limited research was completed on PACs’ effectiveness for the conditions and settings mentioned above.

In 2007, Heinemann [9] called for research to establish an evidence-based practice for PAC rehabilitation outcomes. Since then, many studies focused on the associations between PAC, hospital length of stay (LOS), and readmission, especially for stroke patients. Sacks et al. [2] observed positive associations between the increased use of PAC, shorter average hospital LOS, and lower risk-adjusted readmission rates. Burke et al. [10] worked on readmission from PAC facilities, identifying risk factors (e.g., impaired functional status, increased acuity) and timing (30 and 100 days) for readmission risk. Other researchers focused on variables associated with the Functional Independence Measure (FIM) score and PAC referral [11,12,13,14,15]; important determinants of discharge location included memory and comprehension, living status, and social support. The relationship between readmission and discharge location was also investigated [16, 17], which indicated that strong hospital-SNF linkages [16] and high nursing care quality [17] could mitigate readmission rates. Some studies on stroke patients also highlighted the influence of regional and facility-level variation in post-acute settings and hypothesized facility-level variation influenced rehabilitation outcome more than geographic location variation [18].

In contrast to stroke, little works are present that addressed PAC referral for cardiovascular disease (CVD) patients. CVD is the leading cause of death in the U.S. [19], and patients with CVD often require rehabilitation after cardiac surgery [20], especially those undergoing coronary artery bypass graft (CABG) or heart valve replacement (VR) accounting for over 500,000 operations annually [21, 22]. Since these patients exhibit increased risks of additional cardiac events, studies focusing on PAC rehabilitation of these patients are essential for restoring the quality of life and mitigating mortality risk [23, 24]. However, most of the existing studies on CABG patients explored variables associated with readmission and reported postoperative infection, heart failure, and cardio-pulmonary complications as the most common risk factors [25,26,27]. A Massachusetts study reported that readmission rates and patient mortality were held constant when significant reductions in acute-care LOS were accompanied by increased PAC usage [28]. Few others predicted hospital readmission following heart failure [29,30,31]; significant predictor variables included the type of valve surgery, hospital LOS, discharge location, age, and the degree of patient follow up.

So far, only a limited number of studies focused on CVD patients’ relation with PAC referral following acute-care discharge. Brown et al. [24] noted that 56% of the coronary artery patients of their study were referred for rehabilitation at discharge; patients exhibiting non-ST-segment elevation myocardial infarction, comorbidities, and greater age being more likely to be referred. Dolansky et al. [20] reported the prominent factors associated with PAC use for aged cardiac patients are: race, gender, and LOS. However, none of the studies on PAC referral of CVD patients incorporated detailed hospital information, i.e., location, capacity, and other specialties in their analysis so far. The integration of these information along with patient demographic and clinical data would be useful to achieve a more detailed understanding of the drivers influencing the variation in PAC referral practice.

Previous studies highlighted that PAC referral practices vary widely [32, 33], and it led many to believe that standardized referral protocols will be beneficial from both outcome and cost perspectives [2]. However, no uniform guidelines have been established yet to help providers predict the appropriate PAC destination for CABG and VR patients. To achieve this target, first, it is imperative to understand prevailing PAC referral practices across the dimensions of geographic region and patient acuity. Considering this target and the gaps as mentioned above in the literature, the objective of this study is to examine the geographic variations in PAC referral pattern and identify the associated risk factors related to hospital characteristics, patient demographics, and clinical information for the decision of discharge location for CABG and VR. To accomplish the objective, patient-level detailed cohort data was obtained through an electronic health record system, and the association of the risk factors influencing PAC as a discharge destination was quantified. Details of the analysis, results, and discussions are presented in the following sections.


Data source

Data was extracted from the Cerner Health Facts data warehouse, which was shared by the Oklahoma State University Center for Health System Innovation (OSU-CHSI). Cerner Health Facts data warehouse is a Health Insurance Portability and Accountability Act (HIPAA) compliant, geocoded data warehouse containing comprehensive clinical records related to pharmacy, laboratory, admission, and billing collected from participating clinical facilities starting from 2000 [34]. The Cerner’s HealthFacts data warehouse contains electronic medical records for over 63 Million unique patients obtained from more than 400 U.S. Cerner clients (hospital systems). This, de-identified fully HIPAA compliant relational database covers 16 years of longitudinal episodes, connected at the patient level using a unique patient identifier. The data contains information on patient specifics, hospital specifics, doctor specifics, diagnostic/clinical information, lab, pharmacy, and billing data.

Study design

The study was a retrospective record review of 5 years’ EHR data spanned from January 2010 to December 2015 collected from the Cerner’s HealthFacts data warehouse. This data warehouse is donated to the Center for Health Systems Innovation at Oklahoma State University for medical and healthcare-related academic research. The data warehouse organized the electronic records in tables and relationships that readily allows for data extraction on ICD-9 procedure codes. The data is extracted by one of the co-authors on this paper (he is also the research director for the research center that provided the data). This study was reviewed by the Institutional Review Board of Texas A&M University and approved as an exempt study (IRB reference number IRB2016-0453 M).

Study population

The study population consisted of patients who had CABG or valve replacement (n = 14,224) surgery. International Classification of Diseases, 9th Clinical Modification (ICD-9-CM) procedure codes (10 codes) [36.10–36.17, 36.19–36.2] were used to identify the CABG patients and (21 codes) [35.00–35.04, 35.10–35.14, 35.20–35.38, 35.97, 35.99] were used to identify patients with valve replacement. In case a patient had multiple hospitalizations for the same condition in the study window, only the first admission was included to avoid potential effects of aging or readmission.

The study population was individuals with CABG or valve replacement and who were discharged alive after their index hospitalization. Patients who expired (n = 185), left against medical advice (LMA) or discharged for outpatient service were excluded from the study. The outpatient service was considered as inappropriate for this study because this study focuses on the referral to post-acute care facilities of the patient who got admitted and stayed in the acute care hospital for some days to get the required procedure done. Patients who had procedures performed before the admission date or after the discharge date (considered as incoherent data) were also excluded. These entries were considered incoherent data indicating data collection error because clearly it is not possible to have a procedure performed before admitting the hospital or after the discharge from the hospital. This study only included adult patients (> = 20 years) admitted through the emergency department or transferred from other clinical facilities or referred by a physician/HMO. This study excluded patients with length of stay > 75 days (n = 8) and age < 20 years (n = 8). Patients with missing predictor variables (n = 2685) were excluded. All these exclusions resulted in a final sample size of 14,224 patients from 49 acute care hospitals. Figure 1 summarizes the data cleaning and study cohort generation process.

Fig. 1

Flow diagram of the final sample cohort processing through data analysis

Definitions and variables

The dependent variable, discharge destination, was obtained from the initial encounter table. The categories of discharge destination were: (1) Discharged to home, (2) Discharged to home health care service (HHC) (3) Discharged to skilled nursing facility (SNF) (4) Discharged to long term care hospitals (LTCH) (5) Discharged to inpatient Rehabilitation facility (IRF) (6) Discharged to others. Discharged to ‘others’ included several miscellaneous discharge locations, and the number of patients discharged to these locations was very low. The miscellaneous locations were discharged to another short-term hospital, discharged within this institution to Medicare-approved swing bed, discharged to court/ law enforcement/jail, discharged to a designated cancer center or children’s hospital, discharged to a federal health care facility, discharged to a psychiatric hospital, and unknown. All of these miscellaneous locations are binned into one category ‘others’ to bring clarity in our analytical model.

In the analytical model, in total, 29 independent variables were considered. The independent variables were categorized into five categories like hospital location (census region), provider/hospital characteristics, patient demographics, related factors of PAC referral discharge, and comorbidity and diagnosis information. Table 1 provides a list of the 29 variables considered in this study.

Table 1 Variables considered for the discharge location analytical model

This study considered census division of the hospital in the regional (hospital location) category and hospital bed size range, the teaching facility affiliation, and Hospital status (urban or rural) are categorized under provider/hospital characteristics. Demographic variables include age, marital status (married, divorced, single and widowed), race (Caucasian, African American, Asian, Hispanic, Native American and others), gender. Other predictor variables related to PAC referral discharge were the length of stay and the Charlson Index. The length of stay in the hospital was identified by the number of hospital days.

Based on the frequencies of the CABG and VR procedures in the study population, we identified four CABG procedures and two VR procedures, accounting for approximately 90.9% (12935) patients’ reasons for hospitalization. These include coronary bypass surgery for two arteries (n = 4496), coronary bypass surgery for three arteries (n = 3133), coronary bypass surgery for one artery (n = 2743), coronary bypass surgery for four or more arteries (n = 1108), open and other replacement of aortic valve with tissue graft (n = 874), and open and other replacement of aortic valve (n = 581). These six events were considered as independent binary variables to facilitate examination of the effect of these specific cardiac procedures in the discharge decision.

Comorbidity is defined as the coexistence of additional diseases or disorders in the same person with a specific index disease [35]. To assess the contribution of comorbid conditions in the discharge location (PAC referral), we examined if the patient had suffered from atrial fibrillation (ICD9–427.31), hypertension (ICD9–401.9), coronary atherosclerosis (ICD9–414.01), intermediate coronary syndrome (ICD9–411.1), hyperlipidemia (ICD9–272.4), acute posthemorrhagic anemia (ICD9–285.1), acute myocardial infarction (ICD9–410.71), tobacco use disorder (ICD9–305.1), diabetes mellitus without complication (ICD9–250), acute kidney failure (ICD9–584.9), pulmonary collapse (ICD9–518), congestive heart failure (ICD9–428) and unspecified anemia (ICD9–285.9). These 13 comorbid diagnoses were selected for assessment because they were the most frequent common comorbidities in the study population. However, the Charlson comorbidity index was used to capture the overall effect of comorbidities in each patient [36].

Descriptive analysis and model development

The primary focus of this study is the analysis of patient discharge location (PAC referral). Analyses included descriptive statistics for discharge location (Fig. 2) and exploratory analysis (univariate and bivariate analyses). Variables with large numbers of missing values and outliers were excluded. Chi-square tests were performed for categorical variables to test for differences in distribution of discharge locations among patients. Variables with p-value less than 0.1 [37] in the bivariate test were included as candidates in the multinomial logistic regression model. Percentages and medians with interquartile ranges are recorded for categorical and continuous variables in Table 3 in Appendix. The likelihood ratios for all variables are also reported in Table 4 in Appendix.

Fig. 2

Distribution of discharge location

Regression analysis has been widely used in healthcare and medical research in different predictive models specially in the field of disease prediction [38, 39], patient outcome prediction (i.e. readmission, mortality) [25, 27] and so on. Multinomial logistic regression is a popular method used for predicting a response variable with more than two categories (i.e. Home, LTCH, SNF, IRF, HHC). In this study, multinomial logistic regression was used to develop the analytical model, and ‘Home’ was used as the reference category. ‘Home’ was selected as the reference category because this category represented the highest percentage (40.96%) of the discharge destinations. To reduce bias in estimation of such analytical models, the clustering effect of patients within facilities within geographic regions is emphasized to consider performing multilevel mixed model [18, 40]. Therefore, we tested the multilevel analysis approach considering random effects from the census division and found that the difference between single and multilevel results is negligible. For example, the difference between single and multilevel model misclassification errors is 0.21% only. Also, the Akaike Information Criterion (AIC) difference between two models is 0.417%. It implies that there is none or minimal clustering effect of census divisions in our dataset. So, we only considered a single-level analysis approach in this work. The model’s accuracy was calculated based on multiclass receiver operating characteristic (ROC) value and overall misclassification error. A 10-fold cross-validation of the model was conducted to assess model overfitting. We contrasted the mean misclassification error of cross-validation with the misclassification error of the model developed with the entire cohort. A flowchart describing the methodology used to develop and validate the model is shown in Fig. 3.

Fig. 3

Flowchart of the methodology followed to develop the multinomial predictive model

The relative risk ratio (RR), the p-value, and the 95% likelihood confidence intervals of the predictor variables for each category are reported in Table 4 in Appendix section. The data analysis and all the statistical tests were carried out in R version 3.2.3, an open-source package from the R Foundation for Statistical Computing [41].


Analysis and model interpretation

The final cohort of the study population had a mean age of 63.5 ± 11.81 years (mean ± sd) with 10,234 (71.9%) male and 11,946 (84%) Caucasian. Figure 2 describes the distribution of the discharge destination. The most dominant discharge location was Home (40.97%), followed by HHC (36.75%), SNF (12.10%), IRF (8.2%), and LTCH (1.16%). Table 3 in Appendix summarizes demographic characteristics, hospital information, and information related to cardiac events and comorbidities for each discharge location. Discharge location ‘others’ does not represent any specific PAC type. Therefore, in the rest of the result section, discharge location category ‘others’ is ignored while discussing the model insights.

Factors related to discharge destination selection

Table 4 in Appendix provides the significance of the factors associated with discharge destination in the multinomial logistic regression. A p-value of 0.05 was used as a threshold to distinguish significant variables. The relative risk ratio and 95% confidence interval (CI) limits are also provided in Table 4 in Appendix. The association of the factors related to discharge location is stated in the following paragraphs.

  1. (i)

    Regional (Hospital Location)

The location of the hospital, captured as census division, was found to be strongly associated with the selection of discharge location. For census division, ‘East South Central’ was chosen as the reference category. Patients in West South Central are around 13 times more likely to discharge to LTCH, 9 times more likely to HHC, and 3 times more likely to SNF or IRF compared to patients in East South Central. Patients from the Middle Atlantic are around 6 times more likely to discharge to HHC, 3 times to SNF, 2 times to IRF, and 2 times to LTCH compared to home than patients from East South Central. Figure 4 summarizes the risk ratios for the nine census divisions.

  1. (ii)


Fig. 4

The relative risk ratio of 9 census divisions for 4 types of PAC

The bed size of the index hospital was also found to be a significant predictor of discharge destination. The bed size range of 500+ was chosen as the reference category in the model. Compared to the 500+ bed size hospitals, those with 300–499 bed size are 40 and 70% less likely to discharge patients to SNF and HHC, respectively, and around 3 times more likely to discharge to IRF compared to Home (Table 4 in Appendix). Patients from 200 to 299 bed size range hospitals are 50, 70, and 50% less likely to discharge to SNF, HHC, and IRF, respectively, compared to Home. Patients admitted to hospitals with bed size range 6–99 are less likely to be discharged to HHC and LTCH. Figure 5 summarizes the variation of the RR values for different bed size range.

Fig. 5

The relative risk ratio of different bed size range for 4 types of PAC

Whether a facility is a teaching hospital or not was also a significant factor of discharge destination. Hospitals with teaching are less likely to discharge patients to PAC compared to home. No significant difference was found in referral to HHC, IRF, and LTCH between urban and rural hospitals.

  1. (iii)

    Patient Demographic

Gender was found to be significant for discharge location. Females are more likely to be discharged to PAC than males. The likelihood of a female patient being discharged to SNF and LTCH is twice that of males. Further, Asians are around two times more likely to be referred to HHC compared to Caucasians, and single, divorced, and widowed patients are 2 to 3 times more likely to be discharged to SNF, IRF, and LTCH compared to married. Age is another significant predictor in the discharge destination referral, with the likelihood of PAC referral increasing with age.

  1. (iv)

    Related factors of PAC referral discharge

Length of stay and Charlson comorbidity index were also significant predictor variables for the decision of discharge location. Patients with longer length of stay and higher comorbidity index were more likely to be released to a PAC facility compared to Home.

  1. (v)

    Comorbidity and Diagnosis information

Patients with valve replacement exhibited higher discharge rates to PAC facilities compared to non-valve replacement. CABG and VR patients diagnosed with coronary atherosclerosis and tobacco use disorder are less likely to discharge to a PAC facility compared to Home. However, CABG or VR Patients diagnosed with acute kidney failure are 2 times more likely to discharge to LTCH. The associations of other individual comorbid diagnoses were not found to be significant.

Predictive power of the model

The average 10-fold cross-validated predictive accuracy of the model is 62.6% considering overall misclassification error. The average cross-validation (CV) accuracy (62.5%) is consistent with the accuracy based on the complete cohort. The standard deviation of the accuracy based on CV is very low (0.015) indicating that model is very stable to data/sample variations. The misclassification error in each of the CV runs did not differ significantly from the findings in the complete cohort. Figure 6 shows the multiclass ROC curves for every discharge location category along with overall ROC for the model. pROC package from R was used to analyze and compare the multiclass ROC curves for 6 discharge locations [42]. The area under the overall ROC curve (AUC) is 0.685, and the AUC for Home, IRF, LTCH, SNF, HHC, and others are 0.72, 0.53, 0.52, 0.58, 0.72, and 0.46, respectively.

Fig. 6

Receiver operating characteristics (ROC) curve for multiclass prediction model with multinomial logistic regression


This study revealed that 54.5% of CABG patients and 73.3% VR patients were discharged with some PAC care. This finding seems reasonable because VR procedures are associated with more complexity than CABG. VR patients experience frequent complications after surgery that result in arrhythmias and unspecified heart failure [20]. For those receiving PAC, most were referred to HHC (relative proportion 63.1%), which is consistent with Dolansky et al. [20], who stated that surgery patients require less recovery care than non-surgical medical patients requiring lesser need for PAC. In our study population, IRF and LTCH were infrequently used as only 9.4% patients were sent to IRF and LTCH combined. This is also reasonable for the CABG or VR patients as they typically require little daily physical or occupational therapy (> 3 h) [20], which is a necessary admission criterion to discharge to IRF. Further, the average length of stay in our study population was 10 days, which does not meet the admission criteria of LTCH (more than 25 days for LTCH admission [43]).

Geographic variation of PAC use was significant, which is consistent with the existing literature [32, 33]. Picone et al. [44] hypothesized that the rate of PAC referral for cardiac patients aged 65 or more is positively correlated with the number of PAC facilities per 10,000 people, which our results partially support. Compared to other divisions (see Table 2), West South-Central exhibits higher relative capacity for both LTCH and HHC compared to the mean (LTCH: 2.27% vs 1.45%, HHC: 57.1% vs 37.83%) and higher relative referral to LTCH and HHC compared to the mean (LTCH: 1.80% vs 1.09%, HHC: 45.8% vs 31.59%). Similarly, West North Central exhibits both higher relative SNF capacity and relative SNF referral (69.73% vs 56.39, 20.4% vs 14.36%, respectively).

Table 2 PAC capacity and referral rate across census divisions

However, this capacity effect does not always hold. For example, patients in the East North Central are more likely to be referred to IRF as compared to other divisions (12.2% vs 6.72%), even though the relative capacity is lower than average (3.87% vs 4.32%). Further, among divisions, Mountain exhibits high relative HHC capacity (45.29% vs 37.83%) with low relative HHC referral (2.5% vs 31.59%). For capacity and referral profiles within divisions, Middle Atlantic exhibits high SNF capacity (67.69%) with low SNF referral (16.00%) and low HHC capacity (24.78%) with high HHC referral (68.3%). Overall, these results strongly indicate that, while PAC capacities are sometimes positively associated with PAC referral, other significant underlying factors exist that may contravene the capacity effect. Although researchers conjecture causes such as practice styles, service quality, insurance coverage, and acute / PAC business relationships [32] for these underlying factors, geographic variation in PAC referral is not yet clearly understood.

Our findings suggest that hospital characteristics also affect PAC referrals significantly. Smaller hospitals are more likely to refer patients to SNF (Table 4 in Appendix Referent 500 beds: bed size < 5, 6–99: SNF RR ratios: 3.2, 2.3, respectively), while larger hospitals are more likely to refer to HHC (Table 4 in Appendix: bed size < 5, 6–99: HHC RR ratios; 0.4, 0.5, respectively). Teaching hospitals are less likely to refer to PAC across all PAC types (Table 4 in Appendix Referent Non-Teaching: Teaching Hospital: SNF RR: 0.2; HHC RR: 0.4; IRF RR: 0.4; LTCH RR: 0.1).

Length of stay and comorbidity are both correlated with PAC referral, which is consistent with past findings [32, 45, 46]. Hospital length of stay is important because early discharge can contribute to less control over the patient’s condition and more reliance on PAC use [32]. Our study indicates that total comorbidity (Charlson Index) is a better predictor than specific comorbid conditions. This means that overall health complexity has more influence on referral than individual comorbid conditions. As comorbidity increases, the patient is more likely to be referred to SNF, IRF, or LTCH than to HHC (Table 4 in Appendix Charlson: HHC RR 1.1; SNF RR 1.3; IRF RR 1.3; LTCH RR 1.3). This is consistent with studies on PAC referral for patients with stroke and hip replacement [45, 46]. However, tobacco users or patients with a smoking history are less likely to be referred to PAC, which contradicts the results reported by Brown et al. [24]. Our analysis also indicated that CABG and VR patients with coronary atherosclerosis are less likely to be discharged to the PAC facilities (SNF, IRF, LTCH).

Female patients are more likely to be referred to PAC than are males (Table 4 in Appendix Referent Male: RR > 1 for all PAC categories, SNF RR 2.0), which is consistent with Suaya et al. [47], and older patients are more likely to be referred to PAC than younger (Table 4 in Appendix Age: RR > 1 for all PAC categories, SNF RR 1.1). These results are consistent with the cardiac study of Dolansky et al. [20], but again contradict the findings of Brown et al. [24], who finds that younger cardiac patients are more likely to be referred to PAC (Age Referent < 50: 66–80 Odds Ratio (OR) 0.9; > 80 OR 0.7). We note that, although the average ages of our dataset and Brown’s are very similar, Brown et al. considers not only CABG and VR, but also myocardial infarction, percutaneous coronary intervention (PCI), stable angina, and heart transplant, which could account for these differences.

Race also influences PAC referral, with Caucasians being more likely to use SNF (Table 4 in Appendix SNF Referent Caucasian: Asian RR 0.4; Native American RR 0.2; Hispanic RR 0.9; African American RR 0.9); Asians and African Americans being more likely to use HHC (Table 4 in Appendix HHC: Asian RR 1.6; African American RR 1.2); and African Americans and Hispanics being more likely to use LTCH (Table 4 in Appendix LTCH: African American RR 1.6; Hispanic RR 1.2). Our results are generally consistent with the review of Cortes and Arthur [48], although they do not consider multiple categories of PAC. Explanations for these racial disparities in the referral practice include cultural practices, education, and language [48]. However, similar demographics-driven disparities are also observed in cardiac rehabilitation practice in Canada. Studies on cardiac rehabilitation referral on using Canadian dataset [49] also reported women, ethnocultural minorities and people with low income to be less referred for cardiac rehabilitation despite their higher needs.

This study is, however, not without limitations; major limitations include the following: First, we had limited information on patient socioeconomic status and potentially important variables such as discharge condition and discharge medication of the patients. This information can be used to validate the discharge locations referred to for patients in the dataset. Second, limited information was available regarding the PAC facilities; information on the PAC location, quality, and type of service, patient adherence, and length of stay at the PAC would have facilitated analysis of the patient metrics and outcome. Finally, since data on number of beds and health professionals working in the PACs were not available, the number of PAC facilities was used as a proxy variable for PAC capacity. A future extension to this study can consider true enrollment of patients into the PAC and analyze discrepancies between referral and enrollment patterns. Further studies should also include larger datasets, especially including more hospitals from each census region, and multilevel mixed modeling should be performed for the analysis to reduce the clustering effects of patients within facilities within geographic regions.


In this study, factors linked to PAC referral following acute care discharge was investigated using an EHR-extracted CABG and VR patient dataset. Our findings concluded that the regional location of the hospitals and hospital capacity (bed size) influenced the patient discharge practice. Disparities in PAC availability vis-a-vis referral across different U.S. census regions were also observed. The racial and gender-based disparity was also statistically significant, with Asians, Hispanics, and Native Americans being less likely to be referred to PAC compared to Caucasians, and female patients being more likely to be referred than males. Though patients diagnosed with relevant comorbid conditions were, in most cases, likely to be discharged to PAC facilities after the CABG or VR procedure, tobacco disorder and coronary atherosclerosis patients were less likely to be referred to PAC. These findings can help the clinicians to streamline the discharge planning process early in the patient’s acute care stay, and thereby, facilitate discharge processes, care coordination, and transition of care, following surgery. In future, inclusion of supportive information from the PAC facilities could allow to account for the PAC effectiveness and result in more robust and insightful findings.

Availability of data and materials

The datasets analyzed in this study were based on sensitive EHR data and is not publicly available to protect patient privacy. Researchers can contact corresponding author for details of obtaining the data use agreement to have access to the data.



Area under the Curve


Coronary artery bypass graft


Confidence interval


Cross validation


Cardiovascular disease


Electronic health record


Functional independence measure


Home health care


Health Insurance Portability and Accountability Act


International Classification of Diseases, 9th Clinical Modification


Intensive care unit


Institutional Review Board


Inpatient rehabilitation facilities


Left against medical advice


Long term care facilities


Oklahoma State University Center for Health System Innovation


Post-acute care


Receiver operating characteristic


Relative risk ratio


Skilled nursing facilities


Valve replacement


  1. 1.

    Jubelt LE, Goldfeld KS, Chung W, Blecker SB, Horwitz LI. Changes in discharge location and readmission rates under bundled payment. JAMA Intern Med. 2016;176:115–7.

    Article  Google Scholar 

  2. 2.

    Sacks GD, Lawson EH, Dawes AJ, Weiss RE, Russell MM, Brook RH, et al. Variation in hospital use of postacute care after surgery and the association with care quality. Med Care. 2016;54:172–9.

    Article  Google Scholar 

  3. 3.

    Medicare’s post-acute care: Trends and ways to rationalize payments. Rep. to Congr. Medicare Paym. Policy. 2015. 163. Available from: Accessed 21 Apr 2018.

  4. 4.

    Mechanic R. Post-acute care — the next frontier for controlling medicare. N Engl J Med. 2014;370:692–4.

    CAS  Article  Google Scholar 

  5. 5.

    Newhouse J, Garber A, Graham R. Contributors to Geographic Variation in Medicare Spending. Interim Rep. Comm. Geogr. Var. Heal. Care Spend. Promot. High-Value Heal. Care Prelim. Comm. Obs. . 2013. Available from: Accessed 21 Apr 2018.

  6. 6.

    Rau J. Medicare Seeks To Curb Spending On Post-Hospital Care. Kaiser Heal. News. 2013; Available from: Accessed 21 Apr 2018.

  7. 7.

    Association AH. TrendWatch:maximizing the value of post-acute care T. 2010. Available from:

    Google Scholar 

  8. 8. Health IT in Long-term and Post Acute Care. 2013. Available from: Accessed 21 Apr 2018.

  9. 9.

    Heinemann AW. State-of-the-science symposium on postacute rehabilitation: setting a research agenda and developing an evidence base for practice and public policy - executive summary. Am J Speech-Language Pathol. 2007;16:290–4.

    Article  Google Scholar 

  10. 10.

    Burke RE, Whitfield EA, Hittle D, Joon MS, Levy C, Prochazka AV, et al. Hospital Readmission From Post-Acute Care Facilities: Risk Factors, Timing, and Outcomes. J Am Med Dir Assoc. 2016;17:249–55.

    Article  Google Scholar 

  11. 11.

    Mokler J. Predicting discharge destination for patients with severe motor stroke : important functional tasks. Neurorehabil Neural Repair. 2000;14:181–5.

    CAS  Article  Google Scholar 

  12. 12.

    Agarwal V, Mcrae MP, Bhardwaj A, Teasell RW. A model to aid in the prediction of discharge location for stroke rehabilitation patients. Arch Phys Med Rehabil. 2003;84:1703–9.

  13. 13.

    Barbara J. Determinants of discharge destination for stroke patients. Rehabil Nurs. 2004;29:154–63.

    Article  Google Scholar 

  14. 14.

    Brauer SG, Bew PG, Kuys SS, Lynch MR, Morrison G. Prediction of discharge destination after stroke using the motor assessment scale on admission : a prospective , multisite study. Arch Phys Med Rehabil Vol. 2007;89:1061–5.

    Article  Google Scholar 

  15. 15.

    Frank M, Engelter S. Prediction of discharge destination after neurological rehabilitation in stroke patients. Eur Neurol. 2010;63(4):4227–33.

    Article  Google Scholar 

  16. 16.

    Rahman M, Foster AD, Grabowski DC, Zinn JS, Mor V. Effect of hospital – SNF referral linkages on rehospitalization. Health Serv Res. 2013;48(6pt1):1898–919.

    Article  Google Scholar 

  17. 17.

    Fehnel CR, Lee Y, Wendell LC, Thompson BB, Potter NS, Mor V. Post – acute care data for predicting readmission after ischemic. J Am Hear Assoc. 2015;4:e002145.

    Google Scholar 

  18. 18.

    Reistetter TA, Kuo YF, Karmarkar AM, Eschbach K, Teppala S, Freeman JL, et al. Geographic and facility variation in inpatient stroke rehabilitation: multilevel analysis of functional status. Arch Phys Med Rehabil. 2015;96:1248–54.

    Article  Google Scholar 

  19. 19.

    Prevention. C for DC and. Hospital stays grow shorter, heart disease leading cause of hospitalization. 2001. Available from: Accessed 21 Apr 2018.

  20. 20.

    Dolansky MA, Xu F, Zullo M, Shishehbor M, Moore SM, Rimm AA. Post-acute care services received by older adults following a cardiac event: a population-based analysis Mary. J Cardiovasc Nurs. 2011;25:342–9.

    Article  Google Scholar 

  21. 21.

    Li Y, Cai X, Mukamel DB, Cram P. Impact of length of stay after coronary bypass surgery on short- term readmission rate: an instrumental variable analysis. Med Care. 2014;51:45–51.

    CAS  Article  Google Scholar 

  22. 22. Heart Valve Surgery, Heart Disease and Stroke Statistics. Available from: [cited 2017 May 4]

  23. 23.

    Aroesty JM. Patient education: Recovery after coronary artery bypass graft surgery (CABG) (Beyond the Basics). UpToDate. 2017. Available from: [cited 2017 May 4]

    Google Scholar 

  24. 24.

    Brown TM, Hernandez AF, Bittner V, Cannon CP, Ellrodt G, Liang L, et al. Predictors of cardiac rehabilitation referral in coronary artery disease patients findings from the American Heart Association ’ s get with the guidelines program. Am Coll Cardiol Foundation. 2009;54:515–21.

    Article  Google Scholar 

  25. 25.

    Hannan EL, Zhong Y, Lahey SJ, Culliford AT, Gold JP, Smith CR, et al. 30-Day Readmissions After Coronary Artery Bypass Graft Surgery in New York State. JCIN. Elsevier Inc. 2011;4:569–76.

    Google Scholar 

  26. 26.

    Price JD, Romeiser JL, Rosengart TK, Gnerre JM. Risk analysis for readmission after coronary artery bypass surgery : developing a strategy to reduce readmissions. J Am Coll Surg. 2013;216:412–9 American College of Surgeons.

    Article  Google Scholar 

  27. 27.

    Stewart RD, Campos CT, Jennings B, Lollis SS, Levitsky S, Lahey SJ. Predictors of 30-day hospital readmission after coronary artery bypass. Ann Thorac Surg. 2000;70:169–74.

    CAS  Article  Google Scholar 

  28. 28.

    Bohmer RMJ, Newell J, Torchiana DF. The effect of decreasing length of stay on discharge destination and readmission after coronary bypass operation. Surgery. 2002;132:10–5.

    Article  Google Scholar 

  29. 29.

    Pack QR, Priya A, Lagu T, Pekow PS, Engelman R, Kent DM, et al. Development and validation of a predictive model for short- and medium-term hospital readmission following heart valve surgery. J Am Heart Assoc. 2016;5:e003544.

    Article  Google Scholar 

  30. 30.

    AbdelRahman SE, Zhang M, Bray BE, Kawamoto K. A three-step approach for the derivation and validation of high-performing predictive models using an operational dataset: congestive heart failure readmission case study. BMC Med Inform Decis Mak. 2014;14:41.

    Article  Google Scholar 

  31. 31.

    Inouye S, Bouras V, Shouldis E, Johnstone A, Silverzweig Z, Kosuri P. Predicting readmission of heart failure patients using automated follow-up calls. BMC Med Inform Decis Mak. 2015;15:22.

    Article  Google Scholar 

  32. 32.

    Kane RL, Lin W, Blewett LA. Geographic variation in the use of post-acute care. Health Serv Res. 2002;37:667–82.

    Article  Google Scholar 

  33. 33.

    Fisher E, Skinner J. Making sense of geographic variations in health care: the new iom report. Heal Aff Blog. 2013; Available from: Accessed 15 Apr 2018.

  34. 34.

    Center for Health Systems Innovation. Available from: Accessed 21 Apr 2018.

  35. 35.

    Valderas JM, Starfield B, Sibbald B, Salisbury C, Rloand M. Defining comorbidity: implications for understanding health and health services. Ann Fam Med. 2009;7:357–63.

    Article  Google Scholar 

  36. 36.

    Charlson M, Sax F, MacKenzie C, Fields S, Braham R, Douglas RJ. Assessing illness severity: does clinical judgment work? J Chronic Dis. 1986;39:439–52.

    CAS  Article  Google Scholar 

  37. 37.

    Torabipour A, Arab M, Zeraati H, Rashidian A. Multivariate analysis of factors influencing length of hospital stay after coronary artery bypass surgery in Tehran , Iran. Acta Med Iran. 2016;54:124–33.

    PubMed  Google Scholar 

  38. 38.

    Freeze J, Erraguntla M, Verma A. Data Integration and Predictive Analysis System for Disease Prophylaxis : Incorporating Dengue Fever Forecasts. In: Proc. 51st Hawaii Int. Conf. Syst. Sci; 2018. p. 913–22.

    Google Scholar 

  39. 39.

    Erraguntla M, Zapletal J, Lawley M. Framework for infectious disease analysis : a comprehensive and integrative multi-modeling approach to disease prediction and management. Health Informatics J. 2019;25:1170–87.

    Article  Google Scholar 

  40. 40.

    Song J, Zhao H, Pan C, Li C, Liu J, Pan Y. Risk factors of chronic periodontitis on healing response: a multilevel modelling analysis. BMC Med Inform Decis Mak. 2017;17:1–9.

    CAS  Article  Google Scholar 

  41. 41.

    RCT. R: A language and environment for statistical computing. 2015. Available from:

    Google Scholar 

  42. 42.

    Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, Müller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.

    Article  Google Scholar 

  43. 43.

    Stefanacci RG. Admission criteria for facility-based post-acute services. Ann Long-Term CareClinical Care Aging. 2015;23:18–20.

    Google Scholar 

  44. 44.

    Picone G, Wilson RM, Chou S. Analysis of hospital length of stay and discharge destination using hazard functions with unmeasured heterogeneity. Health Econ. 2003;12:1021–34.

    Article  Google Scholar 

  45. 45.

    Bond J, Gregson B, Lecouturier MSJ, Rousseau N, Rodgers H. Predicting place of discharge from hospital for patients with a stroke or hip fracture on admission. J Health Serv Res Policy. 2000;5:133–9.

    CAS  Article  Google Scholar 

  46. 46.

    Pablo PD, Losina E, Phillips CB, Fossel AH, Mahomed N, Lingard EA, et al. Determinants of discharge destination following elective total hip replacement. Arthritis Rheum. 2004;51:1009–17.

    Article  Google Scholar 

  47. 47.

    Suaya JA, Shepard DS, Normand ST, Ades PA, Prottas J, Stason WB. Health services and outcomes research use of cardiac rehabilitation by medicare beneficiaries after myocardial infarction or coronary bypass surgery. Circulation. 2007;116:1653–62.

    Article  Google Scholar 

  48. 48.

    Cortés O, Arthur HM. Determinants of referral to cardiac rehabilitation programs in patients with coronary artery disease: a systematic review. Am Heart J. 2006;151:249–56.

    Article  Google Scholar 

  49. 49.

    Grace SL, Turk-Adawi K, Santiago de Araújo Pio C, Alter DA. Ensuring cardiac rehabilitation access for the majority of those in need: a call to action for Canada. Can J Cardiol. 2016;32:S358–64.

    Article  Google Scholar 

Download references


The authors would like to acknowledge the support of the Texas A&M High-performance supercomputing facility to conduct data analysis and statistical analysis in R studio. This work was conducted using the data extracted from the Cerner Corporation’s HealthFacts data warehouse of electronic medical records. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the Cerner Corporation.


No specific funding was received for this study. The first author was supported by funding from the Texas A&M Engineering Experiment Station (TEES) during this research. However, the funders were not involved in the study design, data collection, analysis, decision to publish, or production of this manuscript.

Author information




IS participated in data analysis, model development, and interpretation, manuscript preparation. ML, ME, HK, and DD supervised the research, developed ideas, assisted in result interpretation, drafting and revising the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ineen Sultana.

Ethics declarations

Ethics approval and consent to participate

Approval for this study was taken from the Institutional Review Board (IRB) of Texas A&M University (IRB number 2016-0453 M) where need for Informed consent was waived by IRB. Besides, data sharing agreement with Oklahoma State University-Center for Health Systems Innovation (OSU-CHSI) was executed.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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



Table 3 Summary statistic of the variables included in the multinomial logistic regression model across the discharge destinations
Table 4 Risk ratios, p-values and 95% CI of the predictor variables in the model

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

Verify currency and authenticity via CrossMark

Cite this article

Sultana, I., Erraguntla, M., Kum, HC. et al. Post-acute care referral in United States of America: a multiregional study of factors associated with referral destination in a cohort of patients with coronary artery bypass graft or valve replacement. BMC Med Inform Decis Mak 19, 223 (2019).

Download citation


  • Post-acute care
  • Patient discharge
  • Referral
  • Length of stay
  • Coronary artery bypass