Key PointsQuestion
Is initiation of pulmonary rehabilitation after hospitalization for chronic obstructive pulmonary disease (COPD) associated with better survival?
Findings
In this retrospective observational study that included 197 376 Medicare beneficiaries discharged after hospitalization for COPD, initiation of pulmonary rehabilitation within 3 months of discharge, compared with later or no initiation of pulmonary rehabilitation, was significantly associated with lower risk of mortality at 1 year (hazard ratio, 0.63).
Meaning
These findings support current guideline recommendations for pulmonary rehabilitation after hospitalization for COPD, although the potential for residual confounding exists and further research is needed.
Importance
Meta-analyses have suggested that initiating pulmonary rehabilitation after an exacerbation of chronic obstructive pulmonary disease (COPD) was associated with improved survival, although the number of patients studied was small and heterogeneity was high. Current guidelines recommend that patients enroll in pulmonary rehabilitation after hospital discharge.
Objective
To determine the association between the initiation of pulmonary rehabilitation within 90 days of hospital discharge and 1-year survival.
Design, Setting, and Patients
This retrospective, inception cohort study used claims data from fee-for-service Medicare beneficiaries hospitalized for COPD in 2014, at 4446 acute care hospitals in the US. The final date of follow-up was December 31, 2015.
Exposures
Initiation of pulmonary rehabilitation within 90 days of hospital discharge.
Main Outcomes and Measures
The primary outcome was all-cause mortality at 1 year. Time from discharge to death was modeled using Cox regression with time-varying exposure to pulmonary rehabilitation, adjusting for mortality and for unbalanced characteristics and propensity to initiate pulmonary rehabilitation. Additional analyses evaluated the association between timing of pulmonary rehabilitation and mortality and between number of sessions completed and mortality.
Results
Of 197 376 patients (mean age, 76.9 years; 115 690 [58.6%] women), 2721 (1.5%) initiated pulmonary rehabilitation within 90 days of discharge. A total of 38 302 (19.4%) died within 1 year of discharge, including 7.3% of patients who initiated pulmonary rehabilitation within 90 days and 19.6% of patients who initiated pulmonary rehabilitation after 90 days or not at all. Initiation within 90 days was significantly associated with lower risk of death over 1 year (absolute risk difference [ARD], –6.7% [95% CI, –7.9% to –5.6%]; hazard ratio [HR], 0.63 [95% CI, 0.57 to 0.69]; P < .001). Initiation of pulmonary rehabilitation was significantly associated with lower mortality across start dates ranging from 30 days or less (ARD, –4.6% [95% CI, –5.9% to –3.2%]; HR, 0.74 [95% CI, 0.67 to 0.82]; P < .001) to 61 to 90 days after discharge (ARD, –11.1% [95% CI, –13.2% to –8.4%]; HR, 0.40 [95% CI, 0.30 to 0.54]; P < .001). Every 3 additional sessions was significantly associated with lower risk of death (HR, 0.91 [95% CI, 0.85 to 0.98]; P = .01).
Conclusions and Relevance
Among fee-for-service Medicare beneficiaries hospitalized for COPD, initiation of pulmonary rehabilitation within 3 months of discharge was significantly associated with lower risk of mortality at 1 year. These findings support current guideline recommendations for pulmonary rehabilitation after hospitalization for COPD, although the potential for residual confounding exists and further research is needed.
Approximately 15.5 million individuals in the US reported having chronic obstructive pulmonary disease (COPD) in 2015, and exacerbations resulted in more than 1.5 million emergency department visits and nearly 700 000 hospitalizations annually between 2001 and 2012.1,2 The period after hospital discharge is marked by high levels of health care utilization and elevated risk of mortality. In 2008-2014, among Medicare beneficiaries, mortality within 1 year of discharge was estimated to be 26%.3
Quiz Ref IDPulmonary rehabilitation involves exercise training and self-management education to improve physical and psychological well-being.4 It is effective at relieving dyspnea, increasing exercise tolerance, and improving health-related quality of life for individuals with all stages of COPD, and may also reduce health care utilization.4-6 Prior research has shown that few patients with COPD complete a course of pulmonary rehabilitation because of lack of physician referral, lack of access to facilities, and a variety of patient-related barriers.7 While pulmonary rehabilitation has traditionally been viewed as beneficial for symptomatic, medically stable patients, recent attention has turned to its role after an exacerbation, when patients are at high risk of experiencing a self-reinforcing spiral of dyspnea-related deconditioning.8,9 Although the mechanisms are uncertain, 2 recent meta-analyses of randomized trials found that initiation of pulmonary rehabilitation soon after COPD exacerbation was associated with reduced risk of readmission and death.10,11 However, the total number of patients studied remains small. Current guidelines recommend that patients begin pulmonary rehabilitation within 3 to 4 weeks of a COPD exacerbation.4,12,13
Because results of randomized trials are not always generalizable to routine settings, this study investigated outcomes associated with pulmonary rehabilitation among a national sample of older adults hospitalized for COPD.14-16 The primary objective was to determine the association between initiation of pulmonary rehabilitation within 90 days of discharge and 1-year survival.
Design, Setting, Patients
This study was approved by the institutional review board of Baystate Health, which granted a waiver of informed consent. We conducted a retrospective cohort study of all fee-for-service Medicare beneficiaries older than 65 years who were hospitalized for COPD in 2014, with a final date of follow-up of December 31, 2015, the most recent year for which data were available. We defined the cohort using International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes in accordance with methods used by the Centers for Medicare & Medicaid Services (eTable in the Supplement). We identified the first hospitalization in the year as the index hospitalization. Quiz Ref IDTo create an inception, or new-user, cohort, we excluded beneficiaries who had received pulmonary rehabilitation during the previous year.17,18 To emulate criteria used in previous trials we excluded patients with dementia, metastatic cancer, and acute myocardial infarction/acute coronary syndrome. We also excluded patients we thought would be unlikely to initiate pulmonary rehabilitation within 90 days of discharge and to have a high risk of mortality. These included patients with an index hospital stay longer than 31 days, those transferred to another acute care hospital, individuals discharged to skilled nursing facilities who remained in the facility for more than 30 days, and individuals who died within 30 days of discharge (Figure 1).
Participation in Pulmonary Rehabilitation
We used the Medicare outpatient and carrier files to identify patients who received pulmonary rehabilitation, using common procedure coding system codes G0424 (COPD-specific pulmonary rehabilitation) and G0237, G0238, and G0239 (nonspecific pulmonary rehabilitation). We counted the number of days from the index hospitalization discharge to the first session and the number of sessions completed within the year after discharge.
The primary outcome was death due to any cause within 1 year of discharge from the index hospitalization. Vital status was obtained from the Medicare Beneficiary Summary File, which includes a date of death for all beneficiaries but not cause of death.
Patient and Hospital Factors
Sociodemographic characteristics included age, sex, race/ethnicity, Medicaid dual eligibility (a proxy for lower socioeconomic status), current tobacco use, and distance to the nearest pulmonary rehabilitation facility.19 Race/ethnicity was included as we have noted variation in access to and use of pulmonary rehabilitation associated with race in prior research.20 Race/ethnicity was determined using an enhancement to the Social Security Administration categorical race variable, which uses an algorithm to assign beneficiaries to Hispanic or Asian based on last name when race was unknown. Asian/Pacific Islander and American Indian/Alaska Native were grouped into the “other” category because of small numbers. We assessed individual comorbidities using Agency for Healthcare Research and Quality (AHRQ) comorbidity software and computed a longitudinal Charlson Comorbidity Index score.21,22 We calculated each patient’s risk of frailty using methods described by Segal et al.23 Markers of COPD severity included a count of hospitalizations in the prior 12 months, receipt of home oxygen within 90 days before admission, and receipt of mechanical ventilation during the index hospitalization. Additionally, we assessed emergency department visits and readmissions after hospital discharge but before the initiation of pulmonary rehabilitation. Hospital characteristics included number of beds, Census region, teaching status, and rurality.
Our primary analysis focused on the association between initiation of pulmonary rehabilitation within 90 days of an index hospital discharge and 1-year survival; additional analyses explored the relationships of timing of initiation and the number of sessions completed with survival. Given the large size of our sample, we calculated absolute standardized differences (ASDs) to compare patients who initiated pulmonary rehabilitation within 90 days of discharge with those who never initiated or who initiated between days 91 and 365 after discharge.24 We developed a nonparsimonious propensity model to predict initiation of pulmonary rehabilitation within 90 days of discharge using generalized estimating equations (GEE) logistic regression, accounting for clustering of patients within hospitals. The model included all patient sociodemographic characteristics, Medicaid dual eligibility, tobacco use, comorbidities, frailty, markers of disease severity, features of the index admission, and selected interaction terms. In propensity model development, we excluded patients discharged from hospitals where pulmonary rehabilitation was not provided because such patients had no possibility of receiving it. We then applied model coefficients to estimate a propensity for treatment for patients at all hospitals.
Association Between Initiation of Pulmonary Rehabilitation and Survival
To assess the association between initiation of pulmonary rehabilitation and survival, we modeled time from hospital discharge to death, using Cox regression. We evaluated the proportionality assumption by adding an interaction term for pulmonary rehabilitation with time; the test for proportionality was not met. We then modeled survival with time-varying pulmonary rehabilitation exposure. This approach attributed days alive after discharge, but before pulmonary rehabilitation initiation, as “unexposed” or “untreated” time; it also addressed nonproportionality. Patients who initiated pulmonary rehabilitation after 90 days, but before 365 days, contributed time until initiation toward survival among those who never initiated and then were censored when they started pulmonary rehabilitation. We computed an E-value to quantify the minimum strength of association necessary for an unmeasured confounder to negate any significant association observed between pulmonary rehabilitation and survival.25,26
We performed additional analyses to evaluate whether our primary analytic strategy was robust to different methodologic approaches and to explore potential interactions associated with selected patient characteristics. First, rather than treating pulmonary rehabilitation as a time-varying exposure, we restricted the cohort to individuals who survived at least 90 days from discharge to eliminate overlap of exposure and outcome assessment periods. We fit GEE logistic regression models for mortality within the year after discharge, adjusting for unbalanced covariates and propensity to start pulmonary rehabilitation. Second, within this same cohort we fit models using stabilized inverse probability of treatment weighting (SIPTW), as well as standardized mortality rate weighting (SMR).27,28
Third, we repeated this analysis adding interaction terms with pulmonary rehabilitation initiation for age group, frailty, comorbid heart failure, tertiles of overall comorbidity burden, and receipt of home oxygen before the index admission. Fourth, we conducted a propensity-matched analysis, matching each individual who initiated pulmonary rehabilitation within 90 days of discharge to an individual with a similar propensity who had not initiated within 90 days using a greedy matching algorithm. To avoid immortal time bias, we required that the matched control be alive on the day of the pulmonary rehabilitation participant’s first session.29 Kaplan-Meier curves were then used to compare mortality of those who initiated pulmonary rehabilitation within 90 days of discharge with those matched who never initiated or did so after 90 days. Differences in mortality within 1 year of discharge among the matched pairs were assessed using the McNemar test. In addition, in the propensity-matched cohort, we used Cox regression to model time from pulmonary rehabilitation initiation to death, adding adjustment for hospital readmissions and emergency department visits that occurred after the index hospitalization but before initiation of pulmonary rehabilitation.
Timing of Pulmonary Rehabilitation Initiation
We carried out exploratory analyses to investigate whether the timing of pulmonary rehabilitation initiation in relation to hospital discharge moderates its association with survival. Recognizing that mortality risk is highest in the early days after discharge, and because it is common for patients to experience readmission after an index discharge, we computed the number of days between a patient’s first pulmonary rehabilitation visit and the closest COPD discharge. Then, using our full study cohort, we developed a Cox model with time-varying exposure to pulmonary rehabilitation, adjusting for all patient factors including propensity to start within 90 days and adding interaction terms allowing for a change in hazard for initiation after 30 days and after 60 days.
Number of Pulmonary Rehabilitation Sessions Completed
In an additional exploratory analysis, to evaluate a potential dose-response relationship, we created a new cohort restricted to individuals who initiated pulmonary rehabilitation and survived at least 90 days after discharge. We defined the exposure period as the first 90 days after the most recent COPD discharge and assessed outcomes in the subsequent 9 months. We counted the number of sessions completed in the exposure period and then modeled survival from 91 days to death (or censoring at 1 year), evaluating the number of sessions completed as a continuous variable.
Files received from the CMS have minimal missing data. A few fields include an “unknown” category; however, use was rare, and cases categorized as unknown were excluded from analysis.
We present absolute risk differences from the adjusted models for all analyses. No adjustment for multiple comparisons was made; thus, the findings of the timing of initiation and number of sessions completed analyses should be interpreted as exploratory.
All statistical testing was 2-sided, using a .05 level of significance. All analyses were performed using SAS version 9.4 (SAS Institute Inc) and figures were created using Stata release 16 (StataCorp).
Quiz Ref IDOf 197 376 patients assessed from 4446 hospitals, 2721 (1.5%) initiated pulmonary rehabilitation within 90 days of discharge and 3161 (1.6%) initiated within days 91 and 365. Compared with those who never participated or who initiated pulmonary rehabilitation after 90 days of discharge, those who initiated within 90 days were younger (mean, 74.5 vs 77 years; ASD, 0.36), more often men (47.6% vs 41.3%; ASD, 0.13), more often non-Hispanic white (92.6% vs 85.1%; ASD, 0.24), and lived closer to a pulmonary rehabilitation facility (mean, 5.8 vs 9.8 miles; ASD, 0.35). Those who initiated pulmonary rehabilitation within 90 days of discharge had less comorbidity (mean, 3.5 vs 4.2; ASD, 0.25), a lower risk of frailty (mean, 13% vs 20%; ASD, 0.51), and were more likely to have no prior-year admissions (61.9% vs 52.4%; ASD, 0.24) but were more likely to receive home oxygen before hospitalization (39.4% vs 31.7%; ASD, 0.16) (Table 1). A total of 38 302 beneficiaries (19.4%) died within a year of hospital discharge, including 198 (7.3%) among those who initiated pulmonary rehabilitation within 90 days and 38 104 (19.6%) among those who initiated after 90 days or not at all. Eighteen hospitals (0.4%) with unknown region (with 25 patients total) were excluded from analyses (Table 2).
Association Between Initiation of Pulmonary Rehabilitation and Survival
Quiz Ref IDIn the full study cohort, attributing time before initiation of pulmonary rehabilitation as untreated time, and adjusting for unbalanced covariates and propensity to start pulmonary rehabilitation, participation was significantly associated with a lower risk of death within 1 year of discharge (absolute risk difference [ARD], –6.7% [95% CI, –7.9% to -5.6%]; hazard ratio [HR], 0.63 [95% CI, 0.57 to 0.69]; P < .001). The E-value for this model was 2.1 (eFigure in the Supplement). Multiple alternative analytic approaches restricted to survivors of the first 90 days yielded similar effect estimates (Figure 2). There was no statistically significant heterogeneity of association based on age (P = .85 for interaction), frailty (P = .13 for interaction), or comorbid heart failure (P = .50 for interaction), but there was a significant interaction of pulmonary rehabilitation with prior home oxygen use (among users: crude mortality rates, 8.8% vs 15.9%; ARD, –5.7% [95% CI, –7.4% to –3.5%]; odds ratio [OR], 0.60 [95% CI, 0.49 to 0.75]; P < .001; among nonusers: crude mortality rates, 4.5% vs 12.7%; ARD, –6.8% [95% CI, –8.0% to –5.4%]; OR, 0.43 [95% CI, 0.34 to 0.54]; P < .001) and comorbidity burden (low: crude mortality rates, 2.6% vs 10.7%; ARD, –7.6% [95% CI, –8.6% to –6.2%]; OR, 0.27 [95% CI, 0.19 to 0.39]; P < .001; medium: crude mortality rates, 6.1% vs 12.5%; ARD, –5.0% [95% CI, –6.7% to –2.8%]; OR, 0.57 [95% CI, 0.43 to 0.75]; P < .001; high: crude mortality rates, 13.8% vs 18.6%; ARD, –3.8% [95% CI, –6.7% to –0.5%]; OR, 0.76 [95% CI, 0.59 to 0.97]; P = .03).
In the sample matched on propensity and survival to pulmonary rehabilitation initiation, a high level of covariate balance was achieved, with the exceptions of distance to pulmonary rehabilitation facility, hospital size, and urban location (Table 1 and Table 2). Kaplan-Meier curves showed pulmonary rehabilitation initiation within 90 days was significantly associated with lower mortality compared with initiation after 90 days or not at all (P < .001), as did the McNemar test (P < .001) (Figure 3). Mortality within 1 year of pulmonary rehabilitation initiation was 7.3% among those who initiated pulmonary rehabilitation within 90 days and 14.1% among those matched who initiated after 90 days or not at all. Modeling the time from initiation of pulmonary rehabilitation in the propensity-matched cohort, pulmonary rehabilitation was significantly associated with lower risk of death (ARD, –6.8% [95% CI, –8.4% to –5.2%]; HR, 0.50 [95% CI, 0.42 to 0.59]; P < .001).
Among individuals who initiated pulmonary rehabilitation within 90 days of discharge from their most recent COPD hospitalization, 1020 (37.5%) initiated rehabilitation within 30 days and 115 (11.3%) died within 1 year; 1075 (39.5%) between 31 and 60 days, with 59 (5.5%) deaths within 1 year, and 626 (23.0%) between 61 and 90 days, with 24 (3.8%) deaths within 1 year. In exploratory analysis, allowing a change in hazard at different start dates relative to discharge, initiation of pulmonary rehabilitation within 30 days (ARD, –4.6% [95% CI, –5.9% to –3.2%]; HR, 0.74 [95% CI, 0.67 to 0.82]; P < .001), 31 to 60 days (ARD, –10.6% [95% CI, –12.4% to –8.4%]; HR, 0.43 [95% CI, 0.34 to 0.54]; P < .001), and 61 to 90 days (ARD, –11.1% [95% CI, –13.2% to –8.4%]; HR, 0.40 [95% CI, 0.30 to 0.54]; P < .001) were each significantly associated with lower mortality compared with initiation after 90 days or no initiation (Figure 2).
Number of Sessions Completed
Quiz Ref IDAmong the group of patients who initiated pulmonary rehabilitation within 90 days of the index hospitalization, 2689 (98.8%) survived at least 90 days from discharge and 166 (6.2%) died within 1 year. These patients completed a median of 9 sessions (interquartile range, 4-14) during the 90-day period from their most recent COPD discharge. In exploratory analysis evaluating the number of sessions as a continuous factor, after adjusting for age, comorbidity, prior home oxygen use, and frailty, every 3 additional sessions (a suggested weekly dose) in the first 90 days was significantly associated with lower mortality (HR, 0.91 [95% CI, 0.85 to 0.98]; P = .01).
In this large, population-based cohort of more than 190 000 Medicare fee-for-service patients hospitalized for COPD in the US, initiation of pulmonary rehabilitation within 90 days of discharge, while rare, was significantly associated with better 1-year survival compared with initiation after 90 days or not at all. The strength of association was similar to estimates reported in recent meta-analyses of randomized trials. In exploratory analyses, initiation of pulmonary rehabilitation was significantly associated with lower mortality across a range of potential start times relative to discharge, and there was a statistically significant association between the number of sessions completed and survival. These observations corroborate the results of multiple small randomized trials within routine clinical settings and thus lend additional support to current guideline recommendations produced by professional societies in the US and Europe that recommend pulmonary rehabilitation after hospitalization for COPD.
A 2016 Cochrane meta-analysis pooled the results of 6 trials, involving 670 patients, in which pulmonary rehabilitation after an exacerbation of COPD was compared with conventional care and in which mortality was reported. While early pulmonary rehabilitation had a pooled OR for mortality of 0.68 (95% CI, 0.28 to 1.67), the nature of the intervention varied between trials, and results demonstrated high levels of heterogeneity and study bias.10 A more recent meta-analysis of 13 clinical trials, involving 801 participants, reported lower mortality associated with early pulmonary rehabilitation that reached statistical significance (relative risk, 0.58 [95% CI, 0.35 to 0.98]).11 By contrast, this study included more than 2700 patients treated with pulmonary rehabilitation, allowing the production of more precise estimates, and the investigation of heterogeneity related to patient characteristics. Additionally, because of uncertainty surrounding the optimal timing of initiation,31 and with current guidelines recommending initiation within 3 to 4 weeks of an exacerbation, this study explored the association between pulmonary rehabilitation start date relative to discharge date and mortality. In these analyses, estimates were fairly similar across a set of start dates ranging from less than 30 days to more than 60 days, allaying concern that delayed initiation may be ineffective. A unique finding from this analysis is the significant association between the number sessions completed within the first 90 days of discharge and 1-year survival.
Given the significant association of postexacerbation pulmonary rehabilitation with lower mortality suggested by this and prior studies, it is notable that less than than 2% of patients initiated treatment within 90 days of hospital discharge. To put this into broader context, noninvasive ventilation and long-term oxygen therapy are the only treatments shown to improve survival for patients requiring hospitalization for COPD. As such, the results of this study reinforce the importance of developing more effective strategies for increasing participation in rehabilitation.32 Furthermore, earlier research has demonstrated that this is especially true for women, members of racial and ethnic minority groups, and individuals with lower socioeconomic status.22
This study has several limitations. First, treatment assignment was not random. Although this study adjusted for numerous potential confounders, and while several alternative analytic approaches yielded similar results, there is likely residual bias due to unmeasured confounding. Individuals who initiated pulmonary rehabilitation were younger and less frail, and may have differed in other ways that contributed toward their better outcomes, a phenomenon known as healthy user bias.33 Second, this study accounted for multiple factors that serve as proxies for disease severity (eg, use of home oxygen, receipt of mechanical ventilation); however, pulmonary function tests results or physiologic measures such as those used to compute the BODE (Body Mass Index, Airflow Obstruction, Dyspnea and Exercise Capacity) index, a tool for predicting the risk of mortality among patients with COPD, were not available.34 Nevertheless, the E-value of 2.1 in the primary analysis, a measure of the minimum strength of association with both treatment and outcome for an unmeasured confounder to render the result nonsignificant, suggests that the observed association between pulmonary rehabilitation and 1-year survival was robust. To provide additional perspective, for an unmeasured confounder to negate the association observed between pulmonary rehabilitation and survival, it would need to have a stronger relationship to mortality than receipt of home oxygen and also would need to be much more unevenly distributed across the treatment and control groups. Third, because this study relied on claims data, information about the individual components of pulmonary rehabilitation delivered during sessions, such as education, exercise training, nutrition counseling, or smoking cessation were not available. Therefore, it was not possible to determine whether and how these components contributed to the observed mortality association. Potential mechanisms through which pulmonary rehabilitation may improve survival include increased physical activity, medication adherence, and improved care coordination.35 Because pulmonary rehabilitation programs are heterogeneous, these results should be interpreted as representing the association of treatment with survival across all US programs.36
Fourth, more patient-centered outcomes, such as exercise capacity or quality of life, were not available. Fifth, receipt of physical therapy or cardiac rehabilitation as an alternative to pulmonary rehabilitation was not assessed. To the extent that such interventions are beneficial and were disproportionately received by patients who did not initiate pulmonary rehabilitation, this would have biased the findings toward the null. Sixth, the study was restricted to the 70% of Medicare beneficiaries in the traditional fee-for-service program who were older than 65 years. Caution is thus warranted before extrapolating the findings to younger patients or those in Medicare Advantage Plans. Seventh, because of lack of availability of cause of death, it cannot be established whether pulmonary rehabilitation is associated with a lower risk of death because of COPD or because of other causes.
Among fee-for-service Medicare beneficiaries hospitalized for COPD, initiation of pulmonary rehabilitation within 3 months of discharge was significantly associated with lower risk of mortality at 1 year. These findings support current guideline recommendations for pulmonary rehabilitation after hospitalization for COPD, although the potential for residual confounding exists and further research is needed.
Corresponding Author: Peter K. Lindenauer, MD, MSc, Institute for Healthcare Delivery and Population Science, University of Massachusetts Medical School—Baystate, 3601 Main St, Springfield, MA 01199 (peter.lindenauer@baystatehealth.org).
Accepted for Publication: March 16, 2020.
Author Contributions: Drs Lindenauer and Pekow had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Concept and design: Lindenauer, Stefan, Pekow, Mazor, Lagu, ZuWallack.
Acquisition, analysis, or interpretation of data: Lindenauer, Stefan, Pekow, Priya, Spitzer, Pack, Pinto-Plata, ZuWallack.
Drafting of the manuscript: Lindenauer, Pekow, Priya.
Critical revision of the manuscript for important intellectual content: Stefan, Pekow, Mazor, Priya, Spitzer, Lagu, Pack, Pinto-Plata, ZuWallack.
Statistical analysis: Pekow, Priya.
Obtained funding: Lindenauer.
Administrative, technical, or material support: Lindenauer, Spitzer, Lagu.
Supervision: Lindenauer, Pekow, Lagu, Pinto-Plata.
Conflict of Interest Disclosures: Dr Lagu reported serving as a consultant for the Yale Center for Outcomes Research and Evaluation, under contract to the Centers for Medicare & Medicaid Services. Dr Pack reported receiving a grant from the National Heart, Lung, and Blood Institute (K23HL135440-01A1). No other disclosures were reported.
Funding/Support: The research was supported by the National Heart, Lung, and Blood Institute of the National Institutes of Health under grants K24 HL132008-01 and 5R01HL133046-02.
Role of the Funder/Sponsor: The National Heart, Lung, and Blood Institute had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
2.Ford
ES. Hospital discharges, readmissions, and ED visits for COPD or bronchiectasis among US adults: findings from the Nationwide Inpatient Sample 2001-2012 and Nationwide Emergency Department Sample 2006-2011.
Chest. 2015;147(4):989-998. doi:
10.1378/chest.14-2146PubMedGoogle ScholarCrossref 3.Lindenauer
PK, Dharmarajan
K, Qin
L, Lin
Z, Gershon
AS, Krumholz
HM. Risk trajectories of readmission and death in the first year after hospitalization for chronic obstructive pulmonary disease.
Am J Respir Crit Care Med. 2018;197(8):1009-1017. doi:
10.1164/rccm.201709-1852OCPubMedGoogle ScholarCrossref 4.Spruit
MA, Singh
SJ, Garvey
C,
et al; ATS/ERS Task Force on Pulmonary Rehabilitation. An official American Thoracic Society/European Respiratory Society statement: key concepts and advances in pulmonary rehabilitation.
Am J Respir Crit Care Med. 2013;188(8):e13-e64. doi:
10.1164/rccm.201309-1634STPubMedGoogle ScholarCrossref 7.Cox
NS, Oliveira
CC, Lahham
A, Holland
AE. Pulmonary rehabilitation referral and participation are commonly influenced by environment, knowledge, and beliefs about consequences: a systematic review using the Theoretical Domains Framework.
J Physiother. 2017;63(2):84-93. doi:
10.1016/j.jphys.2017.02.002PubMedGoogle ScholarCrossref 18.Cox
E, Martin
BC, Van Staa
T, Garbe
E, Siebert
U, Johnson
ML. Good research practices for comparative effectiveness research: approaches to mitigate bias and confounding in the design of nonrandomized studies of treatment effects using secondary data sources: the International Society for Pharmacoeconomics and Outcomes Research Good Research Practices for Retrospective Database Analysis Task Force Report—Part II.
Value Health. 2009;12(8):1053-1061. doi:
10.1111/j.1524-4733.2009.00601.xPubMedGoogle ScholarCrossref 19.National Academies of Sciences, Engineering, and Medicine. Accounting for Social Risk Factors in Medicare Payment: Identifying Social Risk Factors. National Academies Press; 2016.
28.Kurth
T, Walker
AM, Glynn
RJ,
et al. Results of multivariable logistic regression, propensity matching, propensity adjustment, and propensity-based weighting under conditions of nonuniform effect.
Am J Epidemiol. 2006;163(3):262-270. doi:
10.1093/aje/kwj047PubMedGoogle ScholarCrossref 31.Greening
NJ, Williams
JEA, Hussain
SF,
et al. An early rehabilitation intervention to enhance recovery during hospital admission for an exacerbation of chronic respiratory disease: randomised controlled trial.
BMJ. 2014;349:g4315. doi:
10.1136/bmj.g4315PubMedGoogle ScholarCrossref 32.Rochester
CL, Vogiatzis
I, Holland
AE,
et al; ATS/ERS Task Force on Policy in Pulmonary Rehabilitation. An official American Thoracic Society/European Respiratory Society policy statement: enhancing implementation, use, and delivery of pulmonary rehabilitation.
Am J Respir Crit Care Med. 2015;192(11):1373-1386. doi:
10.1164/rccm.201510-1966STPubMedGoogle ScholarCrossref 33.Brookhart
MA, Patrick
AR, Dormuth
C,
et al. Adherence to lipid-lowering therapy and the use of preventive health services: an investigation of the healthy user effect.
Am J Epidemiol. 2007;166(3):348-354. doi:
10.1093/aje/kwm070PubMedGoogle ScholarCrossref 34.Celli
BR, Cote
CG, Marin
JM,
et al. The body-mass index, airflow obstruction, dyspnea, and exercise capacity index in chronic obstructive pulmonary disease.
N Engl J Med. 2004;350(10):1005-1012. doi:
10.1056/NEJMoa021322PubMedGoogle ScholarCrossref 36.Spruit
MA, Pitta
F, Garvey
C,
et al; ERS Rehabilitation and Chronic Care, and Physiotherapists Scientific Groups; American Association of Cardiovascular and Pulmonary Rehabilitation; ATS Pulmonary Rehabilitation Assembly and the ERS COPD Audit team. Differences in content and organisational aspects of pulmonary rehabilitation programmes.
Eur Respir J. 2014;43(5):1326-1337. doi:
10.1183/09031936.00145613PubMedGoogle ScholarCrossref