Abstract

Background

Recent reports show substantial geographic variation in postacute health care spending. Little is known about variation in functional outcomes after postacute rehabilitation for patients with hip fracture. We examined variation in mobility and self-care after hip fracture rehabilitation across inpatient rehabilitation facilities (IRFs), hospital referral regions (HRRs) and states.

Methods

Retrospective cohort study using data from the Centers for Medicare and Medicaid Services (CMS) from 2006 to 2009. Study sample included 149,258 records from patients 66 years and older at 1,166 IRFs located within 292 HRRs and across 50 states. Hip fracture cases were defined by CMS impairment group codes (08.11, 08.12). Hierarchical generalized linear models were used to assess discharge mobility and self-care functional status, adjusting for individual patient characteristics and the random effect of IRFs, HRRs, and states.

Results

Variation in discharge mobility status as assessed by the intraclass correlation percentage (ICC%) was 8.8% across IRFs, 4.0% across HRRs, and 1.8% across states. For self-care, the ICCs were 10.2% across IRFs, 4.8% across HRRs, and 2.4% across states. The range of discharge mobility scores (maximum functional status rating to minimum functional status rating) showed a 9.6-point difference for IRFs, 6.5 for regions, and 2.6 for states. Range of discharge self-care scores were 13.1 for IRFs, 6.8 for HRRs, and 3.4 for states.

Conclusion

Variation in functional status following postacute hip fracture rehabilitation appears to occur primarily at the level of facilities rather than geographic location.

Hip fracture is associated with increased mortality, debility, and institutionalization in older adlults (1). Postacute rehabilitation after hip fracture surgery has been shown to improve functional independence (2,3) and lower mortality (3).

An Institute of Medicine (IOM) report indicates that around 73% of the variation in health care costs across hospital referral regions (HRRs) is comprised of spending on postacute care (4,5). Geographic variation in rates of hip fracture (6–8) and other orthopedic procedures commonly associated with rehabilitation such as knee and shoulder replacement have been reported (9,10). It is not clear from these studies if regional differences in hip fracture prevalence rates, translate to quality of care parameters such as functional status. Few studies have focused on geographic variation in functional status in postacute care settings (11–13).

Patient level demographic and clinical characteristics do not fully explain observed rates of variation (14,15) in health care (16,17). While understanding that higher spending does not necessarily translate to better quality of care and outcomes (18), the Centers for Medicare and Medicaid Services (CMS) recently targeted high spending regions in an attempt to limit Medicare coverage and lower payments (14–16). These measures do not appear to have been successful to-date in reducing variation (19). The Federal Register’s proposed rules for the Medicare Inpatient Rehabilitation Facility Prospective Payment System for 2015, include the use of metrics such as mobility and self-care function at discharge as indicators of the quality of patient care (20). The purpose of our study was to examine variation in hip fracture mobility and self-care function at discharge across inpatient rehabilitation facilities (IRFs), HRRs, and states after adjusting for patient level demographic and clinical factors.

Methods

Study Participants

We used the CMS claims files and examined a national sample of Medicare fee-for-service beneficiaries receiving inpatient rehabilitation from 2006 through 2009. We included records from the following files: (a) Medicare Provider Analysis and Review (MedPAR), (b) Inpatient Rehabilitation Facility-Patient Assessment Instrument (IRF-PAI), (c) Beneficiary Summary file, and the (d) Provider of Service file (POS). A data use agreement was completed in accordance with CMS requirements. Approval was obtained from the University’s Institutional Review Board.

Our eligible sample included 185,644 records from patients, 66 years and older, with hip fracture (IRF-PAI impairment codes 08.11, 08.12) who were discharged from acute care to IRFs. These codes correspond to ICD-9 codes 820.00–820.9. The analysis was restricted to those patients with unilateral and bilateral hip fracture, and excluded individuals with pelvis fracture, multiple fractures, and femur shaft fractures (n = 31,792). We further excluded records from 125 rehabilitation facilities that treated less than 10 patients (n = 535) over the 3-year period. Patients with atypical length of stay (ie, ≤3 days, n = 3,256) and patients with missing information on covariates used in our multivariable models were also excluded (n = 803). The final sample included 149,258 patient records from 1,166 rehabilitation facilities located across 292 HRRs in 50 states.

Outcome of Interest: Variation in Discharge Mobility and Self-Care

Functional status at discharge was assessed using items from the Functional Independence Measure (FIM) instrument (21). The 18-items from the FIM instrument are included in the IRF-PAI (22). The items from the IRF-PAI are administered by clinicians within 3 days of admission and discharge, and evaluate the patient’s motor and cognitive abilities across six subscales including self-care, sphincter control, locomotion, transfers, communication, and social recognition. Scores for each functional item range from 1 (total assistance) to 7 (complete independence). We examined variation in self-care and mobility subscales at discharge as they were likely to be influenced in patients with hip fracture (23,24). The mobility subscale includes five items: (a) transfers from bed, chair, and wheelchair, (b) toilet transfers, (c) tub/shower transfers, (d) walking/wheelchair locomotion, and (e) stairs resulting in a range of scores from 5 to 35 with higher scores indicating better performance. The self-care scale is comprised of six items: (a) eating, (b) grooming, (c) bathing, (d) upper body dressing, (e) lower body dressing, and (f) toileting with scores ranging from 6 to 42. The IRF-PAI functional items have been used extensively in the literature and their reliability and validity is well established (21,25).

Variation in discharge scores was assessed instead of change from admission, as the later does not completely capture patient level of function (26). Change scores among patients may differ in the direction of change. Furthermore, comparisons in the magnitude of change for someone starting at the lower end of the scale to others starting higher on the scale are not appropriate and are difficult to interpret from a clinical point of view. Finally, discharge function ratings represent a metric understandable to clinicians as well as discharge planners.

Covariates

Patient level demographic covariates included age, gender and race/ethnicity, marital status, and pre-hospital living status (27–31). Age was examined as a continuous variable. Race/ethnicity was categorized as non-Hispanic white, non-Hispanic black, Hispanic, and others. We used responses for marital status (married vs unmarried [never married, separated, divorced, widowed]) and pre-hospital living status (not alone [family/relatives, friends, attendant, other] vs alone) to construct a support variable. Patients were considered to have support if they were married or were living with others prior to hospitalization. Clinical covariates included mobility or self-care ratings at admission, length of stay, and 30 comorbid conditions from the Elixhauser Comorbidity index (27,32). Additionally, we included clinical comorbidities of osteoarthritis, osteoporosis, avascular necrosis, deep vein thrombosis, and pulmonary embolism. In-patient rehabilitation facilities included freestanding hospitals/centers and units within a hospital, identified from ID codes in the MedPAR file. The region variable in our analysis was based on the Dartmouth HRRs (17). We used the Dartmouth Atlas crosswalk linked to facility ZIP codes in the Provider of Service file to identify the HRRs and states.

Statistical Analysis

Summary statistics (percentages [%], mean and standard deviation [SD]) were calculated for the sample. Multilevel mixed models with a restricted maximum likelihood estimation approach were used to assess discharge mobility and self-care function following hip fracture rehabilitation. The percentage of variance (ICC: intraclass correlation) in mobility and self-care ratings was calculated independently at the level of the rehabilitation facilities, HRRs, and states, after adjusting for their clustering (random) effect. Patient level (fixed effects) adjustments included age, gender, race/ethnicity, support, length of stay, admission self-care and mobility ratings, clinical comorbidities, and Elixhauser comorbidity index to examine discharge mobility and self-care FIM scores. Variation profiles for discharge function along with 95% Empirical Bayes confidence intervals (33,34) were constructed using adjusted two and three level models (patients–IRFs, patients–IRFs–HRR, patients–IRFs–states). Likewise, maps for multivariable adjusted predicted mean discharge mobility and self-care by state were created at two and three levels (patients–states, patients–IRFs–states). All analyses were performed with SAS version 9.4 (SAS, Inc., Cary, NC).

Results

Table 1 presents the unadjusted demographic and clinical characteristics of the study sample. Approximately 73% of the sample were women and the mean age was 81.9 (SD = 7.2) years. The sample was predominantly white (90.8%) and had support (60.9%). The most common comorbidities were hypertension (61.2%), deficiency anemia (27.3%), osteoarthritis (21.6%), osteoporosis (21.1%), chronic pulmonary disease (18.6%), arrhythmias (17.7%), and diabetes without chronic complications (16.0%). The average functional status rating at admission was 9.2 (SD = 3.2) for mobility and 19.5 (SD = 5.3) for self-care. Mean discharge functional status ratings were 19.0 (SD = 6.2) for mobility and 29.5 (SD = 6.8) for self-care. The average length of stay (LOS) was 13.6 days (SD = 4.7).

Table 1.

Characteristics of the Study Sample*

CharacteristicsStudy Sample (N = 149,256)
Age, years81.9 ± 7.2
Female (%)108,898 (73.0)
Race/ethnicity (%)
 Non-Hispanic White135,446 (90.8)
 Non-Hispanic Black5,142 (3.5)
 Hispanic/Latino6,041 (4.1)
 Others2,627 (1.8)
Married (%)57,164 (38.3)
Living with someone (%)86,694 (58.1)
Support (%)
 Married/living with someone90,948 (60.9)
 Not-married/living alone58,308 (39.1)
Elixhauser comorbidities (%)
 Hypertension91,337 (61.2)
 Deficiency anemia40,673 (27.3)
 Chronic pulmonary disease27,828 (18.6)
 Arrhythmias26,396 (17.7)
 Diabetes without chronic complications23,926 (16.0)
Clinical comorbidities (%)
 Osteoarthritis32,255 (21.6)
 Osteoporosis31,508 (21.1)
 Avascular necrosis246 (0.2)
 Deep vein thrombosis2,093 (1.4)
 Pulmonary embolus731 (0.5)
Length of stay, days13.6 ± 4.7
Admission mobility score, mean ± SD (range)9.2 ± 3.2 (5–29)
Discharge mobility score, mean ± SD (range)19.0 ± 6.2 (5–35)
Admission self care score, mean ± SD (range)19.5 ± 5.3 (6–42)
Discharge self care score, mean ± SD (range)29.5 ± 6.8 (6–42)
CharacteristicsStudy Sample (N = 149,256)
Age, years81.9 ± 7.2
Female (%)108,898 (73.0)
Race/ethnicity (%)
 Non-Hispanic White135,446 (90.8)
 Non-Hispanic Black5,142 (3.5)
 Hispanic/Latino6,041 (4.1)
 Others2,627 (1.8)
Married (%)57,164 (38.3)
Living with someone (%)86,694 (58.1)
Support (%)
 Married/living with someone90,948 (60.9)
 Not-married/living alone58,308 (39.1)
Elixhauser comorbidities (%)
 Hypertension91,337 (61.2)
 Deficiency anemia40,673 (27.3)
 Chronic pulmonary disease27,828 (18.6)
 Arrhythmias26,396 (17.7)
 Diabetes without chronic complications23,926 (16.0)
Clinical comorbidities (%)
 Osteoarthritis32,255 (21.6)
 Osteoporosis31,508 (21.1)
 Avascular necrosis246 (0.2)
 Deep vein thrombosis2,093 (1.4)
 Pulmonary embolus731 (0.5)
Length of stay, days13.6 ± 4.7
Admission mobility score, mean ± SD (range)9.2 ± 3.2 (5–29)
Discharge mobility score, mean ± SD (range)19.0 ± 6.2 (5–35)
Admission self care score, mean ± SD (range)19.5 ± 5.3 (6–42)
Discharge self care score, mean ± SD (range)29.5 ± 6.8 (6–42)
*

Data presented are unadjusted number (percentages) or mean values ± SD, as appropriate for each variable.

Frequencies for the top five Elixhauser comorbidity measures.

Table 1.

Characteristics of the Study Sample*

CharacteristicsStudy Sample (N = 149,256)
Age, years81.9 ± 7.2
Female (%)108,898 (73.0)
Race/ethnicity (%)
 Non-Hispanic White135,446 (90.8)
 Non-Hispanic Black5,142 (3.5)
 Hispanic/Latino6,041 (4.1)
 Others2,627 (1.8)
Married (%)57,164 (38.3)
Living with someone (%)86,694 (58.1)
Support (%)
 Married/living with someone90,948 (60.9)
 Not-married/living alone58,308 (39.1)
Elixhauser comorbidities (%)
 Hypertension91,337 (61.2)
 Deficiency anemia40,673 (27.3)
 Chronic pulmonary disease27,828 (18.6)
 Arrhythmias26,396 (17.7)
 Diabetes without chronic complications23,926 (16.0)
Clinical comorbidities (%)
 Osteoarthritis32,255 (21.6)
 Osteoporosis31,508 (21.1)
 Avascular necrosis246 (0.2)
 Deep vein thrombosis2,093 (1.4)
 Pulmonary embolus731 (0.5)
Length of stay, days13.6 ± 4.7
Admission mobility score, mean ± SD (range)9.2 ± 3.2 (5–29)
Discharge mobility score, mean ± SD (range)19.0 ± 6.2 (5–35)
Admission self care score, mean ± SD (range)19.5 ± 5.3 (6–42)
Discharge self care score, mean ± SD (range)29.5 ± 6.8 (6–42)
CharacteristicsStudy Sample (N = 149,256)
Age, years81.9 ± 7.2
Female (%)108,898 (73.0)
Race/ethnicity (%)
 Non-Hispanic White135,446 (90.8)
 Non-Hispanic Black5,142 (3.5)
 Hispanic/Latino6,041 (4.1)
 Others2,627 (1.8)
Married (%)57,164 (38.3)
Living with someone (%)86,694 (58.1)
Support (%)
 Married/living with someone90,948 (60.9)
 Not-married/living alone58,308 (39.1)
Elixhauser comorbidities (%)
 Hypertension91,337 (61.2)
 Deficiency anemia40,673 (27.3)
 Chronic pulmonary disease27,828 (18.6)
 Arrhythmias26,396 (17.7)
 Diabetes without chronic complications23,926 (16.0)
Clinical comorbidities (%)
 Osteoarthritis32,255 (21.6)
 Osteoporosis31,508 (21.1)
 Avascular necrosis246 (0.2)
 Deep vein thrombosis2,093 (1.4)
 Pulmonary embolus731 (0.5)
Length of stay, days13.6 ± 4.7
Admission mobility score, mean ± SD (range)9.2 ± 3.2 (5–29)
Discharge mobility score, mean ± SD (range)19.0 ± 6.2 (5–35)
Admission self care score, mean ± SD (range)19.5 ± 5.3 (6–42)
Discharge self care score, mean ± SD (range)29.5 ± 6.8 (6–42)
*

Data presented are unadjusted number (percentages) or mean values ± SD, as appropriate for each variable.

Frequencies for the top five Elixhauser comorbidity measures.

Table 2 shows the results for the variation (ICC%) in discharge mobility and self-care function at the level of IRFs, HRRs, and states. Without adjustments for random effects (ie, Model I) variation in discharge mobility was 8.8% across facilities, 4.0% across HRRs, and 1.8% across states. Likewise, discharge self-care ICCs for IRFs, HRRs, and states were 10.2%, 4.8%, and 2.4%, respectively.

Table 2.

Variation in Discharge Mobility and Self-Care Functional Status Across Inpatient Rehabilitation Facilities (IRFs), Hospital Referral Regions (HRRs), and States

Discharge Functional StatusUnadjusted ICC%Adjusted ICC%*
Model IModel IIModel III§Model IModel IIModel III§
At IRFs
 Discharge mobility10.18.38.98.87.98.2
 Reduction from Model I||17.8%11.9%10.2%6.8%
 Discharge self-care10.17.68.710.29.29.7
 Reduction from Model I||24.7%13.9%9.8%4.9%
At HRRs
 Discharge mobility5.82.24.01.0
 Reduction from Model I||62.1%75.0%
 Discharge self-care5.42.94.81.1
 Reduction from Model I||46.3%77.1%
At states
 Discharge mobility2.71.81.80.7
 Reduction from Model I||33.3%61.1%
 Discharge self-care2.31.82.40.8
 Reduction from Model I||21.7%66.7%
Discharge Functional StatusUnadjusted ICC%Adjusted ICC%*
Model IModel IIModel III§Model IModel IIModel III§
At IRFs
 Discharge mobility10.18.38.98.87.98.2
 Reduction from Model I||17.8%11.9%10.2%6.8%
 Discharge self-care10.17.68.710.29.29.7
 Reduction from Model I||24.7%13.9%9.8%4.9%
At HRRs
 Discharge mobility5.82.24.01.0
 Reduction from Model I||62.1%75.0%
 Discharge self-care5.42.94.81.1
 Reduction from Model I||46.3%77.1%
At states
 Discharge mobility2.71.81.80.7
 Reduction from Model I||33.3%61.1%
 Discharge self-care2.31.82.40.8
 Reduction from Model I||21.7%66.7%
*

Adjusted for age (years), gender (male, female), race (non-Hispanic White, non-Hispanic Black, Hispanics, others), social support (present, absent), admission self-care or mobility, comorbidities (Elixhauser and clinical comorbidities) and length of stay (days).

Model I: patients–IRFs or patients–HRRs or patients–states.

Model II: patients–IRFs–HRRs.

§Model III: patients–IRFs–states.

||Reduction from Model I = 1 – (Model II or Model III ICC%/Model I ICC%).

Table 2.

Variation in Discharge Mobility and Self-Care Functional Status Across Inpatient Rehabilitation Facilities (IRFs), Hospital Referral Regions (HRRs), and States

Discharge Functional StatusUnadjusted ICC%Adjusted ICC%*
Model IModel IIModel III§Model IModel IIModel III§
At IRFs
 Discharge mobility10.18.38.98.87.98.2
 Reduction from Model I||17.8%11.9%10.2%6.8%
 Discharge self-care10.17.68.710.29.29.7
 Reduction from Model I||24.7%13.9%9.8%4.9%
At HRRs
 Discharge mobility5.82.24.01.0
 Reduction from Model I||62.1%75.0%
 Discharge self-care5.42.94.81.1
 Reduction from Model I||46.3%77.1%
At states
 Discharge mobility2.71.81.80.7
 Reduction from Model I||33.3%61.1%
 Discharge self-care2.31.82.40.8
 Reduction from Model I||21.7%66.7%
Discharge Functional StatusUnadjusted ICC%Adjusted ICC%*
Model IModel IIModel III§Model IModel IIModel III§
At IRFs
 Discharge mobility10.18.38.98.87.98.2
 Reduction from Model I||17.8%11.9%10.2%6.8%
 Discharge self-care10.17.68.710.29.29.7
 Reduction from Model I||24.7%13.9%9.8%4.9%
At HRRs
 Discharge mobility5.82.24.01.0
 Reduction from Model I||62.1%75.0%
 Discharge self-care5.42.94.81.1
 Reduction from Model I||46.3%77.1%
At states
 Discharge mobility2.71.81.80.7
 Reduction from Model I||33.3%61.1%
 Discharge self-care2.31.82.40.8
 Reduction from Model I||21.7%66.7%
*

Adjusted for age (years), gender (male, female), race (non-Hispanic White, non-Hispanic Black, Hispanics, others), social support (present, absent), admission self-care or mobility, comorbidities (Elixhauser and clinical comorbidities) and length of stay (days).

Model I: patients–IRFs or patients–HRRs or patients–states.

Model II: patients–IRFs–HRRs.

§Model III: patients–IRFs–states.

||Reduction from Model I = 1 – (Model II or Model III ICC%/Model I ICC%).

Furthermore, we examined the change in explained variation due to adjustments for random effects (Model II: patients–IRFs–HRRs and Model III: patients–IRFs–states). At IRFs discharge mobility scores decreased 10.2% after accounting for HRRs and 6.8% when adjusting for states. In contrast among HRRs, discharge mobility scores decreased by 75.0% when adjusting for facility effects. Adjusting for facility effects among states resulted in a 61.1% decrease in discharge mobility. Estimates for variables included in the models are reported in Supplementary Tables 3 and 4.

Additionally, we computed variation plots for discharge mobility and self-care functional status (Figures 1 and 2) across facilities. The 1,166 IRFs were ranked based on their multivariable adjusted mean mobility or discharge self-care ratings, and accounting for random state level effects. Each vertical line in the figure represents the adjusted mean discharge mobility or self-care scores along with the Empirical Bayes confidence interval for each facility. For discharge mobility (Figure 1) the shaded vertical lines on the left show the 250 facilities significantly below the mean while shaded vertical lines on the right indicate the 252 facilities above the national average. Similarly, Figure 2 illustrates the variation in discharge self-care scores across facilities.

Distribution of discharge mobility scores across facilities.
Figure 1.

Distribution of discharge mobility scores across facilities.

Distribution of discharge self-care scores across facilities.
Figure 2.

Distribution of discharge self-care scores across facilities.

Figures 3 and 4 show the predicted mean discharge mobility and self-care scores across states with and without adjusting for random effects of facilities. California, Texas, Florida, and Alabama had below average discharge mobility as well as self-care scores. In contrast, Arkansas had discharge mobility scores above the national mean. Iowa, Ohio, and Wisconsin were in the upper tertile for discharge self-care.

Distribution of discharge mobility scores across states. (A) Discharge mobility scores adjusted for patient characteristics. (B) Discharge mobility scores adjusted for patient characteristics and facility effects.
Figure 3.

Distribution of discharge mobility scores across states. (A) Discharge mobility scores adjusted for patient characteristics. (B) Discharge mobility scores adjusted for patient characteristics and facility effects.

Distribution of discharge self-care scores across states. (A) Discharge self-care scores adjusted for patient characteristics. (B) Discharge self-care scores adjusted for patient characteristics and facility effects.
Figure 4.

Distribution of discharge self-care scores across states. (A) Discharge self-care scores adjusted for patient characteristics. (B) Discharge self-care scores adjusted for patient characteristics and facility effects.

Overall, the range of discharge functional status scores (maximum functional status rating to minimum functional status rating) was higher for IRFs when compared to regions and states. For discharge mobility, there was a 9.6-point difference for IRFs, 6.5 for HRRs, and 2.6 for states. Likewise, the range for discharge self-care scores was 13.1 for IRFs, 6.8 for HRRs, and 3.4 for states.

Discussion

In the current study, we examined variation in mobility and self-care status at the level of IRFs, HRRs, and across states in patients receiving postacute rehabilitation following hip fracture. We found greater variation in functional status at discharge among IRFs when compared to HRRs or states. The results were consistent when comparing unadjusted and adjusted multi-level models for both mobility and self-care at discharge.

Our findings of greater variation in discharge functional status at the facility level are consistent with the Donabedian model for health service evaluation and quality of care (35). Structural factors within a facility such as space, equipment, payment methods, availability of staff, and their level of expertise may influence functional outcomes (36–38). In addition, process factors such as physician practice styles, discharge planning, and continuity of care models also influence functional status outcomes (39–43). Gozalo et al. (44) found that among hip fracture patients a successful return to the community depends on the volume and staffing characteristics of the facility. Another study by Leland et al. (45) concluded that existing metrics such as 30-day readmission do not completely capture the needs of the PAC market and there is a need for additional quality measures. A previous study on variation in functional status among stroke patients in inpatient rehabilitation revealed that facility-level practice patterns and procedures influence outcomes (13).

The results of the current study are important for several reasons. First, geographic variation in functional status outcomes is understudied across postacute settings. This is certainly true for inpatient medical rehabilitation facilities. The IOM report (4) on variation in health care spending highlights the need to better understand sources of variation in health care beyond those associated with acute care and patient geography. Another important strength of our study is the focus on variation in measures of functional status as an outcome. Most previous studies examining regional variation in outcomes have focused on health care spending, for example, length of stay or specific drugs, and procedures or services (4,5,46). This is the first study we are aware of that examines regional and facility variation in patient level functional status using a national sample of patients receiving hip fracture rehabilitation. Little is currently known regarding variation in postacute recovery trajectories for persons with hip fracture across geographic regions or different rehabilitation facilities (47).

The need for national information related to variation in functional outcomes is particularly important in the context of the programs and services associated with the Affordable Care Act (48) including the development of Accountable Care Organizations, Medical Homes, site neutral payments and the CMS Bundled Payments for Care Improvement (BPCI) Initiative (49). Our findings are particularly relevant to the recent implementation of the Comprehensive Care for Joint Replacement Program (50) and the Improving Medicare Post-Acute Care Transformation (IMPACT) Act (51).

Information on how patient-centered outcomes such as self-care vary across different postacute facilities and geographic regions is becoming increasingly important as the programs of the Affordable Care Act are implemented. For example, in the bundled payment initiative the focus is on the services needed and outcomes attained by patients across an episode of care (47). Knowing what factors are related to variation in patient-centered outcomes across facilities and regions will be helpful in the successful implementation of these programs (52). Our results suggest substantial differences in variation in mobility and self-care across postacute rehabilitation facilities, HRRs, and states for Medicare beneficiaries receiving treatment following a hip fracture. Factors such as volume, provider supply, workforce differences, as well as other state programs and policies may be responsible for the observed findings (53,54). These differences are a part of our ongoing research examining rehabilitation service areas and require further exploration.

Our results must be interpreted in the context of the study’s limitations. First, we examined only records for Medicare fee-for-service beneficiaries, and our findings are not generalizable to the overall Medicare population. Second, our observations were restricted to IRFs and did not include data from other postacute care settings such as home health, long-term acute care hospitals or skilled nursing facilities. Third, we did not explore facility specific factors such as urban versus rural designations, teaching versus nonteaching, or number of beds. Finally, we used the Dartmouth based HRRs, originally developed for acute care settings as our region level geographic variable for rehabilitation facilities. This unit of analysis may not be ideal for variation across postacute settings. Quality of care among PAC settings is currently a topic of active research. Compared to patients treated at SNFs, hip fracture patients receiving postacute care at IRFs have demonstrated better functional outcomes despite shorter length of stay (55,56). In previous research, people treated at IRFs showed lower readmission rates when compared to other PAC settings (57). In contrast, Deutsch et al. (58) found similar or better outcomes in hip fracture patients receiving care at SNFs. Further research is needed to better understand the redistribution of patients across PAC settings and the potential changes in outcomes due to healthcare reform.

In conclusion, we observed greater variation in postacute hip fracture mobility and self-care at discharge within IRFs compared to geographical regions as defined by HRRs and states. Additional research is required to better understand the factors responsible for differences in functional status outcomes across both facilities and geographic regions. This is particularly true for outcomes in postacute care where substantial changes are planned in service delivery models over the next decade.

Supplementary Material

Supplementary data are available at The Journals of Gerontology, Series A: Biological Sciences and Medical Sciences online.

Funding

This research was supported in part by grants from the NIH, National Institute of Child Health and Human Development (NICHD), the National Center for Medical Rehabilitation Research (NCMRR: K01-HD068513, R24-HD065702, and R01-HD069443); the National Institute on Aging (P30-AG024832); the Agency for Healthcare Research and Quality (AHRQ: R01-HS024711 and R24-HS022134), and the National Institute on Disability Independent Living and Rehabilitation Research (NIDILRR: H133P110012).

Conflict of Interest

The authors report no conflicts of interest related to this manuscript.

Acknowledgment

Srinivas Teppala, MBBS, MPH was a postdoctoral fellow in the Division of Rehabilitation Sciences at the University of Texas Medical Branch during the completion of this study.

References

1.

Tajeu
GS
Delzell
E
Smith
W
et al. .
Death, debility, and destitution following hip fracture
.
J Gerontol A Biol Sci Med Sci
.
2014
;
69
(
3
):
346
353
.

2.

Salpakoski
A
Tormakangas
T
Edgren
J
et al. .
Effects of a multicomponent home-based physical rehabilitation program on mobility recovery after hip fracture: a randomized controlled trial
.
J Am Med Dir Assoc
.
2014
;
15
(
5
):
361
368
.

3.

Torpilliesi
T
Bellelli
G
Morghen
S
et al. .
Outcomes of nonagenarian patients after rehabilitation following hip fracture surgery
.
J Am Med Dir Assoc
.
2012
;
13
(
1
):
81
85
.

4.

Institute of Medicine
.
Interim Report of the Committee on Geographic Variation in Health Care Spending and Promotion of High-Value Care
.
The National Academies Press
;
2013
. http://www.nap.edu/openbook.php?record_id=18308. Accessed
June 30, 2014
.

5.

Newhouse
JP
Garber
AM
.
Geographic variation in Medicare services
.
N Engl J Med
.
2013
;
368
(
16
):
1465
1468
.

6.

Bacon
WE
Smith
GS
Baker
SP
.
Geographic variation in the occurrence of hip fractures among the elderly white US population
.
Am J Public Health
.
1989
;
79
(
11
):
1556
1558
.

7.

Jacobsen
SJ
Goldberg
J
Miles
TP
Brody
JA
Stiers
W
Rimm
AA
.
Regional variation in the incidence of hip fracture. US white women aged 65 years and older
.
JAMA
.
1990
;
264
(
4
):
500
502
.

8.

Hinton
RY
Lennox
DW
Ebert
FR
Jacobsen
SJ
Smith
GS
.
Relative rates of fracture of the hip in the United States. Geographic, sex, and age variations
.
J Bone Joint Surg Am
.
1995
;
77
(
5
):
695
702
.

9.

Buntin
MB
Colla
CH
Deb
P
Sood
N
Escarce
JJ
.
Medicare spending and outcomes after postacute care for stroke and hip fracture
.
Med Care
.
2010
;
48
(
9
):
776
784
.

10.

FitzGerald
JD
Boscardin
WJ
Ettner
SL
.
Changes in regional variation of Medicare home health care utilization and service mix for patients undergoing major orthopedic procedures in response to changes in reimbursement policy
.
Health Serv Res
.
2009
;
44
(
4
):
1232
1252
.

11.

Kane
RL
Lin
WC
Blewett
LA
.
Geographic variation in the use of post-acute care
.
Health Serv Res
.
2002
;
37
(
3
):
667
682
.

12.

Reistetter
TA
Karmarkar
AM
Graham
JE
et al. .
Regional variation in stroke rehabilitation outcomes
.
Arch Phys Med Rehabil
.
2014
;
95
(
1
):
29
38
.

13.

Reistetter
TA
Kuo
YF
Karmarkar
AM
et al. .
Geographic and facility variation in inpatient stroke rehabilitation: multilevel analysis of functional status
.
Arch Phys Med Rehabil
.
2015
;
96
(
7
):
1248
1254
.

14.

Fisher
ES
Wennberg
DE
Stukel
TA
Gottlieb
DJ
Lucas
FL
Pinder
EL
.
The implications of regional variations in Medicare spending. Part 2: health outcomes and satisfaction with care
.
Ann Intern Med
.
2003
;
138
(
4
):
288
298
.

15.

Fisher
ES
Wennberg
DE
Stukel
TA
Gottlieb
DJ
Lucas
FL
Pinder
EL
.
The implications of regional variations in Medicare spending. Part 1: the content, quality, and accessibility of care
.
Ann Intern Med
.
2003
;
138
(
4
):
273
287
.

16.

Newhouse
JP
Garber
AM
.
Geographic variation in health care spending in the United States: insights from an Institute of Medicine report
.
JAMA
.
2013
;
310
(
12
):
1227
1228
.

17.

The Dartmouth Group
.
The Dartmouth Atlas of Health Care
.
2013
. http://www.dartmouthatlas.org/publications/reports.aspx. Accessed
June 30, 2014
.

18.

Weinstein
MC
Skinner
JA
.
Comparative effectiveness and health care spending—implications for reform
.
N Engl J Med
.
2010
;
362
(
5
):
460
465
.

19.

Zhang
Y
Baik
SH
Fendrick
AM
Baicker
K
.
Comparing local and regional variation in health care spending
.
N Engl J Med
.
2012
;
367
(
18
):
1724
1731
.

20.

Federal Register
.
Medicare Program; Inpatient Rehabilitation Facility Prospective Payment System for Federal Fiscal Year 2015
.
2014
. https://www.federalregister.gov/articles/2014/05/07/2014–10321/medicare-program-inpatient-rehabilitation-facility-prospective-payment-system-for-federal-fiscal. Accessed
June 6, 2014
.

21.

Ottenbacher
KJ
Hsu
Y
Granger
CV
Fiedler
RC
.
The reliability of the functional independence measure: a quantitative review
.
Arch Phys Med Rehabil
.
1996
;
77
(
12
):
1226
1232
.

22.

Research Data Assistance Center
.
Inpatient Rehab Facility-Patient Assessment Instrument
. Centers for Medicaid & Medicare Services;
2013
. http://www.resdac.org/cms-data/files/irf-pai. Assessed
June 30, 2014
.

23.

Foss
NB
Kristensen
MT
Kehlet
H
.
Prediction of postoperative morbidity, mortality and rehabilitation in hip fracture patients: the cumulated ambulation score
.
Clin Rehabil
.
2006
;
20
(
8
):
701
708
.

24.

Kristensen
MT
Bandholm
T
Foss
NB
Ekdahl
C
Kehlet
H
.
High inter-tester reliability of the new mobility score in patients with hip fracture
.
J Rehabil Med
.
2008
;
40
(
7
):
589
591
.

25.

Heinemann
AW
Linacre
JM
Wright
BD
Hamilton
BB
Granger
C
.
Relationships between impairment and physical disability as measured by the functional independence measure
.
Arch Phys Med Rehabil
.
1993
;
74
(
6
):
566
573
.

26.

Streiner
DL
Norman
GR
Cairney
J.
Health Measurement Scales: A Practical Guide to Their Development and Use
. 5th ed.
New York, NY
:
Oxford University Press
;
2015
.

27.

Cornwall
R
Gilbert
MS
Koval
KJ
Strauss
E
Siu
AL
.
Functional outcomes and mortality vary among different types of hip fractures: a function of patient characteristics
.
Clin Orthop Relat Res
.
2004
;(
425
):
64
71
.

28.

Graham
JE
Chang
PF
Berges
IM
Granger
CV
Ottenbacher
KJ
.
Race/ethnicity and outcomes following inpatient rehabilitation for hip fracture
.
J Gerontol A Biol Sci Med Sci
.
2008
;
63
(
8
):
860
866
.

29.

Ortiz-Alonso
FJ
Vidan-Astiz
M
Alonso-Armesto
M
et al. .
The pattern of recovery of ambulation after hip fracture differs with age in elderly patients
.
J Gerontol A Biol Sci Med Sci
.
2012
;
67
(
6
):
690
697
.

30.

Reistetter
TA
Graham
JE
Deutsch
A
Markello
SJ
Granger
CV
Ottenbacher
KJ
.
Diabetes comorbidity and age influence rehabilitation outcomes after hip fracture
.
Diabetes Care
.
2011
;
34
(
6
):
1375
1377
.

31.

Semel
J
Gray
JM
Ahn
HJ
Nasr
H
Chen
JJ
.
Predictors of outcome following hip fracture rehabilitation
.
PM R
.
2010
;
2
(
9
):
799
805
.

32.

Elixhauser
A
Steiner
C
Harris
DR
Coffey
RM
.
Comorbidity measures for use with administrative data
.
Med Care
.
1998
;
36
(
1
):
8
27
.

33.

Carlin
BP
Louis
TA.
Bayes and Empirical Bayes Methods for Data Analysis
. 2nd ed.
Boca Raton, FL
:
Taylor & Francis
;
2010
.

35.

Hoenig
H
Lee
J
Stineman
M
.
Conceptual overview of frameworks for measuring quality in rehabilitation
.
Top Stroke Rehabil
.
2010
;
17
(
4
):
239
251
.

36.

Gage
B
Morley
M
Constantine
R
et al. .
Examining Relationships in an Integrated Hospital System: Final Report
. Research Triangle Institute;
2008
. http://aspe.hhs.gov/health/reports/08/examine/report.pdf. Accessed
November 16, 2016
.

37.

Aiken
LH
Smith
HL
Lake
ET
.
Lower Medicare mortality among a set of hospitals known for good nursing care
.
Med Care
.
1994
;
32
(
8
):
771
787
.

38.

Hartz
AJ
Krakauer
H
Kuhn
EM
et al. .
Hospital characteristics and mortality rates
.
N Engl J Med
.
1989
;
321
(
25
):
1720
1725
.

39.

Magaziner
J
Simonsick
EM
Kashner
TM
Hebel
JR
Kenzora
JE
.
Predictors of functional recovery one year following hospital discharge for hip fracture: a prospective study
.
J Gerontol
.
1990
;
45
(
3
):
M101
M107
.

40.

Shortell
SM
Evashwick
C
.
The structural configuration of U.S. hospital medical staffs
.
Med Care
.
1981
;
19
(
4
):
419
430
.

41.

Shortell
SM
LoGerfo
JP
.
Hospital medical staff organization and quality of care: results for myocardial infarction and appendectomy
.
Med Care
.
1981
;
19
(
10
):
1041
1055
.

42.

Davis
KM
Koch
KE
Harvey
JK
Wilson
R
Englert
J
Gerard
PD
.
Effects of hospitalists on cost, outcomes, and patient satisfaction in a rural health system
.
Am J Med
.
2000
;
108
(
8
):
621
626
.

43.

Roy
A
Heckman
MG
Roy
V
.
Associations between the hospitalist model of care and quality-of-care-related outcomes in patients undergoing hip fracture surgery
.
Mayo Clin Proc
.
2006
;
81
(
1
):
28
31
.

44.

Gozalo
P
Leland
NE
Christian
TJ
Mor
V
Teno
JM
.
Volume matters: returning home after hip fracture
.
J Am Geriatr Soc
.
2015
;
63
(
10
):
2043
2051
.

45.

Leland
NE
Gozalo
P
Christian
TJ
et al. .
An examination of the first 30 days after patients are discharged to the community from hip fracture postacute care
.
Med Care
.
2015
;
53
(
10
):
879
887
.

46.

Detsky
AS
.
Regional variation in medical care
.
N Engl J Med
.
1995
;
333
(
9
):
589
590
.

47.

Brock
J
Mitchell
J
Irby
K
et al. .
Association between quality improvement for care transitions in communities and rehospitalizations among Medicare beneficiaries
.
JAMA
.
2013
;
309
(
4
):
381
391
.

48.

111th Congress
.
Patient Protection and Affordable Care Act PL 111-148
.
2010
. http://www.healthcare.gov/law/index.html. Accessed
November 16, 2016
.

49.

Bundled Payments for Care Improvement (BPCI) Initiative: General information
. Centers for Medicare and Medicaid Services;
2016
. https://innovation.cms.gov/initiatives/bundled-payments/. Accessed
June 7, 2016
.

50.

Comprehensive Care for Joint Replacement Model
. Centers for Medicare and Medicaid Services;
2016
. https://innovation.cms.gov/initiatives/cjr/. Accessed
June 7, 2016
.

51.

Improving Medicare Post-Acute Care Transformation Act of 2014 S.2553
.
2014
. https://www.congress.gov/bill/113th-congress/senate-bill/2553/committees Accessed
November 15, 2016
.

52.

Mechanic
R
Tompkins
C
.
Lessons learned preparing for Medicare bundled payments
.
N Engl J Med
.
2012
;
367
(
20
):
1873
1875
.

53.

Baicker
K
Chandra
A
.
Medicare spending, the physician workforce, and beneficiaries’ quality of care
.
Health Aff (Millwood)
.
2004
;Suppl Web Exclusives:
W4
W97
.

54.

Gage
B
Moon
M
Chi
S
.
State-level variation in Medicare spending
.
Health Care Financ Rev
.
1999
;
21
(
2
):
85
98
.

55.

Herbold
JA
Bonistall
K
Walsh
MB
.
Rehabilitation following total knee replacement, total hip replacement, and hip fracture: a case-controlled comparison
.
J Geriatr Phys Ther
.
2011
;
34
(
4
):
155
160
.

56.

Munin
MC
Seligman
K
Dew
MA
et al. .
Effect of rehabilitation site on functional recovery after hip fracture
.
Arch Phys Med Rehabil
.
2005
;
86
(
3
):
367
372
.

57.

Riggs
RV
Roberts
PS
Aronow
H
Younan
T
.
Joint replacement and hip fracture readmission rates: impact of discharge destination
.
PM R
.
2010
;
2
(
9
):
806
810
.

58.

Deutsch
A
Granger
CV
Fiedler
RC
et al. .
Outcomes and reimbursement of inpatient rehabilitation facilities and subacute rehabilitation programs for Medicare beneficiaries with hip fracture
.
Med Care
.
2005
;
43
(
9
):
892
901
.

Author notes

Decision Editor: Stephen Kritchevsky, PhD

Supplementary data