- Open Access
Outcomes after early and delayed rehabilitation for exacerbation of chronic obstructive pulmonary disease: a nationwide retrospective cohort study in Japan
Respiratory Research volume 18, Article number: 68 (2017)
The effectiveness of early pulmonary rehabilitation (PR) for exacerbation of chronic obstructive pulmonary disease (COPD) remains controversial. The present study aimed to compare the outcomes between early and delayed PR for exacerbation of COPD, using a national inpatient database.
Using the Japanese Diagnosis Procedure Combination database, we examined patients who were transported to hospital for exacerbation of COPD, received PR during hospitalisation, and were discharged to their home. The patients were divided into those who received early PR (defined as PR starting within 48 h of admission) and those who received delayed PR. The outcomes included 90-day readmission, length of stay (LOS), and activities of daily living (Barthel index ≥15) at discharge. Multiple imputation was used for missing data. To assess the associations between early PR and the outcomes, we used risk-adjusted treatment effects and instrumental variable methods.
We identified 12,572 eligible patients, including 8459 patients with delayed PR and 4113 with early PR. In the risk-adjusted treatment effect models, the early PR group had lower proportion of 90-day readmission (risk difference, −3.4%; 95% CI, −5.7% to −1.5%) and shorter LOS (−9.8 days; 95% CI, −10.8 days to −8.7 days) than the delayed PR group. There was no significant difference in activities of daily living at discharge between the two groups. The instrumental variable analyses showed similar results.
In this national database study, early PR was associated with reduced 90-day readmission and shortened LOS in patients with exacerbation of COPD.
The effects of early pulmonary rehabilitation (PR) for unstable COPD patients remain controversial. Several previous small-size randomised controlled trials (RCT) showed that early PR reduced readmission and all-cause mortality, and improved patient exercise capability and quality of life . However, recent RCTs produced conflicting results [5, 6].
Although clinical RCTs remain the gold standard for assessing the efficacy of healthcare services, they can only measure the “efficacy” of an intervention under ideal and controlled circumstances. Therefore, the “effectiveness” of early PR for unstable COPD patients remains unclear in routine clinical settings with more heterogeneous populations and less-standardized treatment protocols.
Using a national inpatient database in Japan, the present study aimed to estimate the real-world effectiveness of early PR for exacerbation of COPD compared with delayed PR.
For the present retrospective cohort study, we used the Diagnosis Procedure Combination database, a Japanese national inpatient database. The database contains administrative claims data and discharge abstracts, and has information on dates and doses of drugs used and daily records of examinations and procedures (including rehabilitation). The database includes the following data: type of admission (planned or unplanned); ambulance service use; and patients’ main diagnoses, comorbidities at admission, and complications after admission. All diagnoses are coded with International Classification of Diseases and Related Health Problems 10th revision (ICD-10) codes. A previous validation study showed good sensitivity and excellent specificity of diagnoses in the database . The database also contains the following detailed patient information: age; sex; body height and weight; smoking index; severity of dyspnoea at admission, based on the Hugh-Jones dyspnoea scale (grades I to V) ; level of consciousness at admission, based on the Japan Coma Scale (JCS) ; and activities of daily living (ADL) at admission and discharge, based on the Barthel index (0–20) .
We retrospectively collected patients who were admitted to hospitals with a diagnosis of COPD (ICD-10 codes: J41–J44) as the main diagnosis or diagnosis at admission and discharged between 1 July 2010 and 31 December 2013. The readmission records of the identified patients were followed from 1 July 2010 to 31 March 2014. We included patients who met all of the following criteria: (i) primary diagnosis of COPD at admission; (ii) unplanned admission by ambulance service use; (iii) receipt of PR during hospitalisation; (iv) and discharge to their home. For patients with two or more hospitalisations, we only included data for the initial hospitalisation. We excluded patients who were hospitalised for more than 180 days. We compared the basic characteristics of the included and excluded patients.
We defined early PR as any type and intensity of PR started within 48 h of admission. Delayed PR was defined as any type of PR started after 48 h of admission. In Japan, rehabilitation is reimbursed by public health insurance. In general, geriatric PR includes physical therapy that mainly focuses on improving physical function, i.e., PR that provides early ambulation and adaptive or assistive exercises to assist patients in standing, balancing, and walking better .
The primary outcome was readmission within 90 days of discharge. We used length of stay (LOS) as a secondary outcome variable. Because LOS had a skewed distribution, we used the log-transformed LOS. A previous study showed that a Barthel Index score of <15 was a strong predictor of in-hospital mortality in patients with COPD . Therefore, the present study used Barthel Index ≥ 15 at discharge as another secondary outcome.
We compared the following covariates between the early and delayed PR groups: age; sex; Hugh-Jones dyspnoea scale; JCS on admission; Charlson comorbidity index calculated by recorded ICD-10 codes ; ADL at admission; smoking index (defined as number of cigarettes smoked per day multiplied by number of years smoked); use of corticosteroid on day of admission (converted to equivalent dose of prednisolone ); use of mechanical ventilation on day of admission; use of oxygen on day of admission; intensive care unit admission; and all-cause pre-admission history within 180 days prior to admission. According to a previous study , we divided body mass index into the following five categories: <18.5, 18.5–22.9, 23.0–24.9, 25.0–29.9, and ≥30 kg/m2. Because a previous study showed that annual case-volume of COPD patients was associated with outcomes , we included annual hospital volume of COPD admissions in the analysis. We also included population density of patient neighbourhood area and distance between patient residence and the hospital as proxies for the local organization of health care provision. To understand the details of the rehabilitation programme, we summarized data on total PR sessions (hours) and PR duration (days) for each group.
Risk-adjusted treatment effect
We estimated the risk-adjusted treatment effect of early PR as the difference between the risk-adjusted outcomes of each treatment arm [19, 20] (known as predictive margins, model-adjusted means, or g-formula). For this estimation, we conducted the following steps: (i) construction of a prediction model for each outcome including all covariates (except total PR sessions and duration of PR) and early PR as the predictors; (ii) setting of the PR status for each arm; (iii) calculation of predicted probabilities or values with the distribution of the covariates in our study population; and (iv) calculation of differences in the predicted probabilities (or values) between the arms. We estimated standard errors with the cluster bootstrap method .
Instrumental variable analysis
Because the risk-adjusted treatment effects could be biased by unmeasured confounders, we conducted instrumental variable analyses. In general, instrumental variables meet the following criteria: (i) not associated with patient background characteristics; (ii) associated with treatment selection; and (iii) not directly associated with outcomes [22, 23]. For this study, we used the differential distance (DD) [24, 25] as an instrumental variable. DD was defined as the difference between the distance from patient home to nearest hospital (d1) and the distance from patient home to nearest hospital conducting early PR for half of COPD patients (d2); that is, DD equals d1–d2. We divided DD into two categories: 0 and >0. The adjusted treatment effect was estimated by the ordinary least square (OLS) model with all covariates and the two-stage least square (2SLS) model with all covariates and DD. The validity of the instrumental variable was tested by F-statistics and the Hausman specification test. The null hypothesis for the F-statistics was that the instrumental variable was not associated with treatment selection. We investigated whether the instrumental variable met the above criterion (ii). In the present study, the null hypothesis for the Hausman test was that early PR was not endogenous. When the null hypothesis was rejected, early PR was regarded as endogenous, and we adopted the results of the 2SLS model. When the null hypothesis was not rejected, early PR was not regarded as endogenous, and we adopted the results of the OLS model because the variance of the instrumental variable estimator is generally larger than that of the OLS estimator .
We performed multiple imputation for missing data on several variables, because complete-case analyses (excluding all patients with missing data) can lead to biased results. We created 50 multiple imputed datasets by chained equations with the mice package of the R Software . All the covariates (except total PR sessions and duration of PR) and outcome variables were included in the data imputation process. We combined all the results and variances based on Rubin’s Rule . The P-value in each imputed dataset was combined by z-transformation . We also performed complete-case analyses for comparison.
Among all of the patients hospitalized owing to COPD exacerbation during the study period (n = 45,899), we excluded those who did not receive rehabilitation, those with LOS >180 days, and those who were discharged to a place other than home. Finally, we identified 12,572 eligible patients during the study period. After exclusion of patients with missing data, there were 6955 patients (Fig. 1).
Appendix 1 Table 5 shows the difference between the included and excluded patients. The excluded patients were more likely to have severe consciousness disorders, a lower BMI, severe physical dependence, a lower Barthel index and a lower Hugh-Jones dyspnoea scale score.
Table 1 presents the characteristics of the eligible patients in the early PR group (n = 8459) and delayed PR group (n = 4113). The patients in the early PR group were significantly older, had lower Hugh-Jones dyspnoea scale scores, and were more likely to be treated in hospitals with larger annual hospital volumes of COPD patients than those in the delayed PR group. The average number of total PR sessions was slightly higher and the duration of PR was slightly shorter in the early PR group, but the standardized differences were <10%, indicating a well-balanced distribution between the groups.
Appendix 2 Table 6 shows the patient characteristics in the groups with DD of 0 km or >0 km. The patient characteristics were well-balanced with standardized differences of <10.
Table 2 shows the crude outcomes in the delayed and early PR groups.
Table 3 shows the risk-adjusted treatment effects of early PR on the outcomes. Compared with the delayed PR group, the early PR group had lower 90-day readmission (risk difference, −2.1%; 95% confidence interval (CI), −3.7% to −0.5%) and shorter LOS (difference in LOS, −9.8 days; 95% CI, −10.7 days to −8.8 days) with adjustment for patient characteristics. There was no significant difference in the proportions of Barthel index ≥15 at discharge between the early and delayed PR groups (risk difference, −0.5%; 95% CI, −2.2% to 1.1%).
Table 4 shows the adjusted coefficient estimators of each outcome for early PR in the OLS and 2SLS models. The partial F-statistics from the first-stage regression for each model indicated that DD had sufficient strength for predicting early PR. Early PR was significantly associated with 90-day readmission (coefficient = −0.021; 95% CI, −0.036 to −0.005; P = 0.009) in the OLS model, but not in the 2SLS model. The Hausman specification test did not reject the null hypothesis that early PR was exogenous, and we adopted the results of the OLS model. Early PR was significantly associated with log-transformed LOS in both the OLS model (coefficient = −0.424; 95% CI, −0.463 to −0.386; P < 0.001) and 2SLS model (coefficient = −0.934; 95% CI, −1.156 to −0.712; P < 0.001). The Hausman specification test rejected the null hypothesis that early PR was exogenous, and we adopted the results of the 2SLS model. Early PR was not significantly associated with ADL at discharge in both the OLS model (coefficient = −0.009; 95% CI, −0.025 to 0.007; P = 0.264) and 2SLS model (coefficient = 0.059; 95% CI, −0.057 to 0.174; P = 0.3188). The Hausman specification test did not reject the null hypothesis that early PR was exogenous, and we adopted the results of the OLS model.
The results of the complete-case analyses are shown in Appendix 3 Tables 7 and 8. These results were similar to those in the imputed analyses.
In this national database study, we showed that early PR was associated with lower proportion of 90-day readmission and shorter LOS in patients with exacerbation of COPD. There was no significant difference in the proportions of Barthel index ≥15 between the early PR and delayed PR groups. Our instrumental variable analyses showed that early PR was significantly associated with reduced 90-day readmission and shortened LOS, but not significantly associated with ADL at discharge.
Because our dataset has several missing values, we conducted multiple imputation analyses. Although our analyses hypothesised the missing at random assumption, our complete-case analyses had similar results to the imputed analyses. This indicates the robustness of our analyses.
Previous small-size RCTs showed that early PR improved patient exercise capacity measured by the 6-min walking test [4, 29] and readmission rate [4, 30]. However, recent large RCTs did not show significant improvement of readmission rate [5, 6],, but did find escalation of the 1-year mortality rate . These disparities between the studies can be explained by the following two aspects: (i) early PR did not have any effect because of short LOS and training duration and (ii) some other factors (e.g. outpatient PR) could have confounded or mediated between the early and delayed PR groups . Our findings showed that early PR was associated with shorter LOS and lower proportion of 90-day readmission after discharge, possibly because the average LOS in Japan is generally longer than those in other developed nations. However, the effect size on 90-day readmission (about 2%) was relatively smaller than that in a previous RCT  and expected in the sample size calculation (15%). These findings suggest that the sample sizes in the previous RCTs were too small to have sufficient statistical power for estimating the efficacy of early PR.
As expected, early PR was associated with shorter LOS and lower 90-day readmission rate without worsened ADL status at discharge. Prolonged hospitalisation can reduce patient exercise capacity and may lead to a high probability of readmission . Our results showed that early PR could prevent such adverse effects of prolonged hospitalisation.
Our results indicate that early PR can reduce LOS and readmission rate. Short LOS can improve patient quality of life and reduce hospitalisation costs. A previous study showed that readmission was one of the prognostic factors for COPD . Although the American Thoracic Society recommended early PR for unstable COPD patients, this recommendation was based on an old small-size RCT . Our findings provide new evidence for early PR in unstable COPD patients. Meanwhile, although early PR reduced 90-day readmission, the effect size was smaller than that in a previous RCT and expected in the sample size calculation . Other factors may affect the readmission rate, such as post-discharge outpatient rehabilitation .
The population excluded from our study had more severe consciousness disorders, lower BMIs, and more severe physical dependence. The excluded population included those who did not receive rehabilitation and those who were discharged to a place other than home, and many of these people may have been bedridden.
The present study has some strengths. First, we used a large nationwide inpatient database. Patients with exacerbation of COPD who want to participate in randomised rehabilitation trials are not common in daily clinical settings . A previous RCT could not recruit a sufficient population to detect the effects of early PR . Second, we used instrumental variable and missing value imputation analyses to adjust unmeasured confounders for treatment selection and bias from missing values. Third, we used real-world data from a nationwide inpatient database. The present study verified the overall effectiveness of early rehabilitation compared with non-early rehabilitation for patients with acute exacerbation of COPD in a nationwide, real-world clinical setting.
There are also some limitations to the present study. First, our data did not contain post-discharge long-term outcomes. Therefore, we cannot completely detect post-discharge readmission and deaths. Second, we excluded patients who did not receive PR. This may be limit the generalisability of our study. Third, because of data limitations, we could not analyse the details of the PR programmes. The standardized differences in total PR sessions and duration of PR between the early and delayed groups were <10%. Fourth, the Barthel index, one of our secondary outcomes, may not be the most appropriate measure to evaluate the effect of PR. However, the database did not include other short-term outcomes .
In summary, early PR was associated with reduced 90-day readmission and shortened LOS without worsened ADL in patients with exacerbation of COPD. These findings suggest that early PR should be conducted in patients with exacerbation of COPD.
Two-stage least square
Activities of daily living
Chronic obstructive pulmonary disease
International Classification of Diseases and Related Health Problems 10th revision (ICD-10)
Japan Coma Scale
Length of stay
Ordinary least square
Randomised controlled trials
Pitta F, Troosters T, Probst VS, Spruit MA, Decramer M, Gosselink R. Physical activity and hospitalization for exacerbation of COPD. Chest. 2006;129:536–44.
Spruit MA, Singh SJ, Garvey C, ZuWallack R, Nici L, Rochester C, et al. An official American Thoracic Society/European Respiratory Society statement: key concepts and advances in pulmonary rehabilitation. Am J Respir Crit Care Med. 2013;188:e13–64.
Garcia-Aymerich J, Farrero E, Felez MA, Izquierdo J, Marrades RM, Anto JM. Risk factors of readmission to hospital for a COPD exacerbation: a prospective study. Thorax. 2003;58:100–5.
Puhan MA, Gimeno-Santos E, Scharplatz M, Troosters T, Walters EH, Steurer J. Pulmonary rehabilitation following exacerbations of chronic obstructive pulmonary disease. Cochrane Database Syst Rev. 2011;10:CD005305.
Ko FW, Dai DL, Ngai J, Tung A, Ng S, Lai K, et al. Effect of early pulmonary rehabilitation on health care utilization and health status in patients hospitalized with acute exacerbations of COPD. Respirology. 2011;16:617–24.
Greening NJ, Williams JE, Hussain SF, Harvey-Dunstan TC, Bankart MJ, Chaplin EJ, 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.
Yamana H, Horiguchi H, Fushimi K, Yasunaga H. Comparison of procedure-based and diagnosis-based identifications of severe sepsis and disseminated intravascular coagulation in administrative data. J Epidemiol. 2016;26:530–7.
Hugh-Jones P, Lambert AV. A simple standard exercise test and its use for measuring exertion dyspnoea. Br Med J. 1952;1:65–71.
Shigematsu K, Nakano H, Watanabe Y. The eye response test alone is sufficient to predict stroke outcome--reintroduction of Japan Coma Scale: a cohort study. BMJ Open. 2013. doi:10.1136/bmjopen-2013-002736.
Mahoney FI, Barthel DW. Functional evaluation: the Barthel index. Md State Med J. 1965;14:61–5.
Momosaki R, Yasunaga H, Matsui H, Horiguchi H, Fushimi K, Abo M. Effect of early rehabilitation by physical therapists on in-hospital mortality after aspiration pneumonia in the elderly. Arch Phys Med Rehabil. 2015;96:205–9.
Hasegawa W, Yamauchi Y, Yasunaga H, Sunohara M, Jo T, Matsui H, et al. Factors affecting mortality following emergency admission for chronic obstructive pulmonary disease. BMC Pulm Med. 2014;14:151.
Quan H, Li B, Couris CM, Fushimi K, Graham P, Hider P, Januel JM, et al. Updating and validating the Charlson comorbidity index and score for risk adjustment in hospital discharge abstracts using data from 6 countries. Am J Epidemiol. 2011;173:676–82.
Yamauchi Y, Yasunaga H, Matsui H, Hasegawa W, Jo T, Takami K, et al. Comparison of in-hospital mortality in patients with COPD, asthma and asthma-COPD overlap exacerbations. Respirology. 2015;20:940–6.
Yamauchi Y, Hasegawa W, Yasunaga H, Sunohara M, Jo T, Takami K, et al. Paradoxical association between body mass index and in-hospital mortality in elderly patients with chronic obstructive pulmonary disease in Japan. Int J Chron Obstruct Pulmon Dis. 2014;9:1337–46.
Dres M, Tran TC, Aegerter P, Rabbat A, Guidet B, Huchon G, et al. Influence of ICU case-volume on the management and hospital outcomes of acute exacerbations of chronic obstructive pulmonary disease. Crit Care Med. 2013;41:1884–92.
Austin PC. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivariate Behav Res. 2011;46:399–424.
Austin PC. Balance diagnostics for comparing the distribution of baseline covariates between treatment groups in propensity-score matched samples. Stat Med. 2009;28:3083–107.
Williams R. Using the margins command to estimate and interpret adjusted predictions and marginal effects. Stata J. 2012;12:308–31.
Wooldridge JM. Econometric analysis of cross section and panel data. 2nd ed. Cambridge: MIT Press; 2010.
Sanchez G. How do I obtain bootstrapped standard errors with panel data? Stata. 2011. http://www.stata.com/support/faqs/statistics/bootstrap-with-panel-data/. Accessed 4 Oct 2016.
Angrist JD, Pischke J-S. Mostly harmless econometrics : an empiricists companion. Princeton, NJ: Princeton University Press; 2013.
Brookhart MA, Rassen JA, Schneeweiss S. Instrumental variable methods in comparative safety and effectiveness research. Pharmacoepidemiol Drug Saf. 2010;19:537–54.
McClellan M, McNeil BJ, Newhouse JP. Does more intensive treatment of acute myocardial infarction in the elderly reduce mortality? Analysis using instrumental variables. JAMA. 1994;272:859–66.
Newhouse JP, McClellan M. Econometrics in outcomes research: the use of instrumental variables. Annu Rev Public Health. 1998;19:17–34.
Buuren S. Groothuis-Oudshoorn K: mice: multivariate imputation by chained equations in R. J Stat Software. 2011;45.
Rubin DB. Multiple imputation after 18+ years. J Am Stat Assoc. 1996;91:473–89.
Licht C. New methods for generating significance levels from multiply-imputed data. 2010. Bamberg, https://opus4.kobv.de/opus4-bamberg/frontdoor/index/index/docId/270. Accessed 1 Nov 2016.
Kirsten DK, Taube C, Lehnigk B, Jorres RA, Magnussen H. Exercise training improves recovery in patients with COPD after an acute exacerbation. Respir Med. 1998;92:1191–8.
Seymour JM, Moore L, Jolley CJ, Ward K, Creasey J, et al. Outpatient pulmonary rehabilitation following acute exacerbations of COPD. Thorax. 2010;65:423–8.
Maddocks M, Kon SS, Singh SJ, Man WD. Rehabilitation following hospitalization in patients with COPD: can it reduce readmissions? Respirology. 2015;20:395–404.
Man WD, Kon SS, Maddocks M. Rehabilitation after an exacerbation of chronic respiratory disease. BMJ. 2014;349:g4370.
Vitacca M, Paneroni M, Baiardi P, De Carolis V, Zampogna E, Belli S, et al. Development of a Barthel Index based on dyspnea for patients with respiratory diseases. Int J Chron Obstruct Pulmon Dis. 2016;11:1199–206.
This study was funded by Grants-in-Aid for Research on Policy Planning and Evaluation (Grant Numbers: H27-Policy-Designated-009 and H27-Policy-Strategy-011) from the Japanese Ministry of Health, Labour and Welfare.
Availability of data and materials
Because individual privacy could be compromised, the datasets analysed during the current study are not publicly available. But they are available from the corresponding author on reasonable request.
HM and HY contributed to the conception and design of the study. KF and HY contributed to the data collection. HM and HY contributed to the data analysis. All authors contributed to the data interpretation and drafting of the manuscript for important intellectual content. All authors read and approved the final manuscript.
All authors have completed declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous 3 years; and no other relationships or activities that could appear to have influenced the submitted work.
Consent for publication
Ethics approval and consent to participate
This study was approved by the Institutional Review Board of The University of Tokyo (Approval Number: 3501-(1)). Because all data were de-identified, the requirement for patient informed consent was waived.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Matsui, H., Jo, T., Fushimi, K. et al. Outcomes after early and delayed rehabilitation for exacerbation of chronic obstructive pulmonary disease: a nationwide retrospective cohort study in Japan. Respir Res 18, 68 (2017) doi:10.1186/s12931-017-0552-7
- COPD exacerbation
- Pulmonary rehabilitation
- Clinical epidemiology