Study Design and Setting
We conducted a retrospective cohort study using large, population-based data from the Healthcare Cost and Utilization Project (HCUP) State Inpatient Database (SID) of seven geographically-dispersed US states (Arkansas, California, Florida, Iowa, Nebraska, New York, and Utah) from 2010 through 2013. The HCUP is a family of healthcare databases that are developed through a federal-state-industry partnership and sponsored by the Agency for Healthcare Research and Quality (AHRQ). The HCUP is the largest collection of longitudinal hospital care data in the US, with all-payer, encounter-level information. The HCUP SIDs capture all hospitalizations, regardless of source of disposition, from acute care, non-federal, general and other specialty hospitals within the participating states [8]. These seven states were selected for their high data quality, geographic distribution, and mainly because their data included unique encrypted patient identifiers that enable longitudinal follow-up of specific individuals across years. The institutional review board of Massachusetts General Hospital approved this study.
Study Population
We identified all hospitalized adult patients (aged >40 years) with a principal discharge diagnosis of COPD, as defined by the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) diagnosis codes of 491.21, 491.22, 491.8, 491.9, 492.8, 493.20, 493.21, 493.22, and 496, or those with a primary diagnosis of respiratory failure (codes 518.81, 518.82, 518.84, and 799.1) and a secondary diagnosis of COPD [9, 10]. In the current analysis, we used only the first hospitalizations of the eligible patients during the study period. We also excluded patients who left the hospital against medical advice, those who died in-hospital at their index hospitalization, those who were transferred to another acute-care facility, and out-of-state residents.
Measurements
The SID contain information on the patient characteristics, including demographics (age, sex, and race/ethnicity), primary insurance type (payer), quartiles for estimated household income, patient residence, ICD-9-CM diagnosis and procedure codes, patient comorbidities (29 Elixhauser comorbidity measures and arrhythmia), hospital course (e.g., hospital length-of-stay, in-hospital death), and disposition.
Primary Exposure
The primary exposure was the development of in-hospital AKI during the index hospitalization for AECOPD, as defined by the ICD-9-CM diagnostic codes of 584.5, 584.6, 584.7, 584.8, and 584.9 in any diagnostic fields [11, 12, 13], with excluding AKI as an admission diagnosis. Additionally, we also identified AKI with the use of dialysis, defined as having both of AKI (diagnostic codes, 584.5-584.9) and hemodialysis (procedure code of 39.95 or diagnostic code of V45.1, V56.0 or V56.1) [11, 12].
Outcome Measures
The outcome measures were readmission attributable to any cause within 30 and 90 days of discharge from the index hospitalization for AECOPD. In the COPD literature, 30-day readmission rates have been investigated [9, 13, 14] in the context of the Centers for Medicare and Medicaid Services’ Hospital Readmissions Reduction Program (HRRP) [15]; 90-day readmission rates have also been recognized as an important clinical indicator [16, 17]. The secondary outcome measure was the primary discharge diagnosis of the readmission. To make data presentation and interpretation more meaningful, we consolidated the principal discharge diagnoses (>14,000 ICD-9-CM diagnosis codes) into 285 mutually exclusive diagnostic categories by using the AHRQ-defined Clinical Classifications Software (CCS) [18].
Statistical Analysis
First, we compared the patient characteristics between patients with and without AKI by using Wilcoxon rank sum test or chi-squared test, as appropriate. We also compared Kaplan-Meier curves between the two groups with the use of the log-rank test. Next, we modeled the time-to-readmission by fitting Cox proportional hazards models with generalized estimating equations accounting for patient clustering within hospitals (e.g., severity of patients, physicians’ preference in disease management) [19, 20]. The time-to-readmission for each patient was defined as the period from the discharge to when the first readmission occurred within the 30-day and 90-day follow-up periods. Patients who did not have an outcome were censored at 30 days (or 90 days) from discharge or in-hospital death during the corresponding follow-up period, whichever occurred first. We fitted Cox proportional hazards model with adjustment for potential confounders, such as age, sex, race/ethnicity, primary insurance, quartiles for median household income, residential status, length-of-stay at the index hospitalization, hospital state, and 28 Elixhauser comorbidities as well as arrhythmia [21, 22]. Furthermore, as sensitivity analyses, we repeated the analysis with stratifications by age category (40-64 years and ≥65 years) and sex as previous studies have reported age- and sex-related differences in the readmission rate after hospitalization for AECOPD [9, 13, 14]. Lastly, we compared the 30-day and 90-day readmission rates and calculated the unadjusted and adjusted hazard ratios among patients without AKI, those with AKI without dialysis use, and healthcare use and with AKI and dialysis use. We primarily conducted an available case analysis, and examined consistency with the results of complete case analysis. All analyses used STATA 14.0 (STATA Corp, College Station, TX). All P values were two-tailed, with P<0.05 considered statistically significant.