Abstract
Contemporary comparisons of mortality in matched hemodialysis and peritoneal dialysis patients are lacking. We aimed to compare survival of incident hemodialysis and peritoneal dialysis patients by intention-to-treat analysis in a matched-pair cohort and in subsets defined by age, cardiovascular disease, and diabetes. We matched 6337 patient pairs from a retrospective cohort of 98,875 adults who initiated dialysis in 2003 in the United States. In the primary intention-to-treat analysis of survival from day 0, cumulative survival was higher for peritoneal dialysis patients than for hemodialysis patients (hazard ratio 0.92; 95% CI 0.86 to 1.00, P = 0.04). Cumulative survival probabilities for peritoneal dialysis versus hemodialysis were 85.8% versus 80.7% (P < 0.01), 71.1% versus 68.0% (P < 0.01), 58.1% versus 56.7% (P = 0.25), and 48.4% versus 47.3% (P = 0.50) at 12, 24, 36, and 48 months, respectively. Peritoneal dialysis was associated with improved survival compared with hemodialysis among subgroups with age <65 years, no cardiovascular disease, and no diabetes. In a sensitivity analysis of survival from 90 days after initiation, we did not detect a difference in survival between modalities overall (hazard ratio 1.05; 95% CI 0.96 to 1.16), but hemodialysis was associated with improved survival among subgroups with cardiovascular disease and diabetes. In conclusion, despite hazard ratio heterogeneity across patient subgroups and nonconstant hazard ratios during the follow-up period, the overall intention-to-treat mortality risk after dialysis initiation was 8% lower for peritoneal dialysis than for matched hemodialysis patients. These data suggest that increased use of peritoneal dialysis may benefit incident ESRD patients.
Hemodialysis and peritoneal dialysis differ profoundly, but randomized comparisons have so far proven impossible.1 Despite obvious limitations, observational studies represent the next best design for survival comparisons. Although numerous studies have made peritoneal-dialysis-to-hemodialysis survival comparisons,2–16 findings have not been entirely consistent.17 These survival comparisons are particularly salient in the United States, where mortality rates for hemodialysis patients are much higher than in Europe and in Japan,18 and peritoneal dialysis is used relatively little.19,20 Because both modalities continue to evolve greatly from year to year, up-to-date survival comparisons may help inform the modality choices patients and physicians make when considering dialysis therapy. Comparative mortality studies of contemporary hemodialysis and peritoneal dialysis patients with similar comorbidity burdens are lacking.
We used a matched-pair, retrospective cohort design to compare survival in adult patients initiating dialysis in the United States in 2003 with matching based on propensity of initial peritoneal dialysis use. Matching may reduce the influence of patients using one modality, who share few (if any) measured characteristics with patients using the other modality, because of contraindication or nonclinical forces.21 Our primary objective was to compare survival of incident dialysis patients, treated initially with hemodialysis or peritoneal dialysis, by intention-to-treat analysis in the matched-pair cohort. Secondary objectives were to compare survival across modalities in subsets defined by age, cardiovascular disease, and diabetes and to assess the sensitivity of inference to follow-up commencement (at dialysis initiation, at 90 days thereafter) and modality exposure definition (intention to treat, as-treated).
Results
The unmatched cohort included 98,875 adult patients, of whom 87,277 (88.3%) survived for 90 days after dialysis initiation. The proportions of patients receiving peritoneal dialysis at dialysis initiation (day 0) and at day 90 were 6.8% and 8.1%, respectively.
Table 1 compares patient characteristics by dialysis modality in the unmatched cohort at day 0. With multiple logistic regression, patient characteristics associated (P < 0.05) with higher likelihood of peritoneal dialysis included Native American race, non-Mexican Hispanic and non-Hispanic ethnicity, GN as primary ESRD cause, serum albumin and hemoglobin, presence of hypertension, and tobacco use. Factors associated (P < 0.05) with lower likelihood of peritoneal dialysis included age; male sex; African-American and other race; hypertension, other urologic disease, and other disease as primary ESRD cause; dual eligibility for Medicare and Medicaid; body mass index and blood urea; and presence of congestive heart failure, cerebrovascular disease, peripheral vascular disease, diabetes, chronic obstructive pulmonary disease, malignant neoplasm, alcohol or drug dependence, and inability to ambulate or transfer.
Patient characteristics by modality at dialysis initiation before and after matching
Table 1 also shows patient characteristics by dialysis modality for the matched-pair cohort. Standardized differences between hemodialysis and peritoneal dialysis were <2% of 1 SD for all covariates (data not shown). For 351 peritoneal dialysis patients (5.2%), no matching hemodialysis patient was identified (Table 1). With multiple logistic regression, discriminating characteristics (P < 0.05) of matched versus unmatched peritoneal dialysis patients included age; African-American race; hypertension as primary ESRD cause; dual eligibility for Medicare and Medicaid; serum albumin, blood urea, and hemoglobin; and presence of hypertension.
Cumulative mortality in the matched-pair cohort was 41.9% over a mean follow-up of 2.3 years. In the primary, intention-to-treat analysis of survival from day 0, cumulative survival was higher for peritoneal dialysis patients than for hemodialysis patients, with a corresponding hazard ratio (HR) for peritoneal dialysis versus hemodialysis of 0.92 [95% confidence interval (CI) 0.86 to 1.00, P = 0.04; Figure 1]. Cumulative survival probabilities for peritoneal dialysis and hemodialysis were 93.2% versus 88.6% (P < 0.01) at 6 months, 85.8% versus 80.7% (P < 0.01) at 12 months, 71.1% versus 68.0% (P < 0.01) at 24 months, 58.1% versus 56.7% (P = 0.25) at 36 months, and 48.4% versus 47.3% (P = 0.50) at 48 months. The HR was <1 (P < 0.01) in the first year (Table 2). In an as-treated analysis, the HR for peritonealdialysis versus hemodialysis was 0.80 (95% CI 0.75 to 0.85, P < 0.01), and the HR was <1 (P < 0.01) in the first year (Table 2).
Intention-to-treat in the matched cohort showed lower death risk in PD when follow up began at initiation of dialysis. Risks were similar when follow-up began at day 90. HD, hemodialysis; PD, peritoneal dialysis.
Peritoneal-dialysis-to-hemodialysis HRs of death in the matched cohort
In a secondary analysis of survival from day 90 (Figure 1), overall survival was similar for peritoneal dialysis and hemodialysis patients (HR 1.05, 95% CI 0.96 to 1.16, P = 0.23); survival probabilities for peritoneal dialysis and hemodialysis patients were 96.5% versus 94.8% (P < 0.01) at 6 months, 88.6% versus 86.9% (P = 0.05) at 12 months, 73.6% versus 73.9% (P = 0.77) at 24 months, 60.1% versus 62.3% (P = 0.10) at 36 months, and 49.5% versus 51.8% (P = 0.22) at 48 months. The HR was >1 (P = 0.03) in the second year (Table 2). In an as-treated analysis, the HR for peritoneal dialysis versus hemodialysis was 0.98 (95% CI 0.90 to 1.07, P = 0.68), and the HR was <1 (P = 0.03) in the first year and >1 (P = 0.03) in the third year (Table 2).
Regarding prespecified subgroups, interaction analysis of survival from day 0 and day 90 showed that the association of dialysis modality with survival was modified by baseline age ≥65 years, cardiovascular disease, and diabetes (P < 0.01 for each). In intention-to-treat analysis of survival from day 0, overall HRs were <1 (favoring peritoneal dialysis, P < 0.01) for patients aged <65 years and for patients without cardiovascular disease and patients without diabetes (Figure 2). In intention-to-treat analysis of survival from day 90, HRs were >1 (favoring hemodialysis, P < 0.01) for patients with cardiovascular disease and patients with diabetes (Figure 3).
Matched cohort, peritoneal-dialysis-to-hemodialysis HRs from dialysis day 0 in subgroups. CVD, cardiovascular disease.
Matched cohort, peritoneal-dialysis-to-hemodialysis HRs from dialysis day 90 in subgroups.
Figure 2 shows mortality HRs for intention-to-treat and as-treated analyses of survival from day 0 within subgroups and individual years of follow-up. HRs in the first year were <1 for all subgroups (P ≤ 0.01) for both intention-to-treat and as-treated analyses. In intention-to-treat analysis, the HR in the second year exceeded 1 for patients aged ≥65 years (P < 0.01). In as-treated analysis, HRs were <1 in the second year for patients aged <65 years and for patients without cardiovascular disease and patients without diabetes (P ≤ 0.01). HRs in the third year exceeded 1 for age ≥65 years and diabetes (P ≤ 0.02).
Figure 3 shows mortality HRs for analyses of survival from day 90. In intention-to-treat analysis, HRs in the first year were <1 for patients without cardiovascular disease and patients without diabetes (P < 0.05). HRs in the second and third years exceeded 1 with age ≥65 years, cardiovascular disease, and diabetes (P < 0.05). In as-treated analysis, HRs in the first year were <1 for age <65 years, no cardiovascular disease, and no diabetes (P < 0.01). HRs in the second and third years exceeded 1 for age ≥65 years and cardiovascular disease (P ≤ 0.01); in addition, the HR in the third year exceeded 1 for diabetes (P < 0.01).
Discussion
In this study, we compared survival by initial modality for adult patients beginning dialysis in the United States in 2003. Compared with hemodialysis counterparts, peritoneal dialysis patients in the United States were fewer and tended to be younger with less comorbidity. To balance prognostically important baseline characteristics, we matched hemodialysis and peritoneal dialysis patients on the basis of propensity of initial peritoneal dialysis use. In the matched cohort, peritoneal-dialysis-to-hemodialysis survival comparisons were complicated and partly dependent on subgroup identity and time frame of follow-up.
On average, peritoneal dialysis patients in the United States are younger and have less comorbidity than hemodialysis patients. Regression adjustment, the usual technique for controlling such differences in risk-factor distributions across the dialytic modalities, can be expected to perform well when all relationships between those factors and risk of outcome are correctly specified. This is unlikely in practice because of complex interactions and unmeasured risk factors. Propensity score analysis is an alternative to regression adjustment. The propensity score is defined here as the probability of initial peritoneal dialysis use, given all covariates, and is alone sufficient for control of measured confounding.22 Propensity scores can be used for stratification, adjustment, or matching, but stratification has been shown to admit residual confounding.23 Hypothetically, if some hemodialysis patients at very high risk for death and other peritoneal dialysis patients at very low risk for death shared few, if any, characteristics, simple adjustment might introduce a bias favoring peritoneal dialysis; matching by propensity scores reduces the chance of introducing this bias. For 351 peritoneal dialysis patients (5.2%) in our cohort, no matching hemodialysis patient was identified, and a substantial number of hemodialysis patients were not paired with a peritoneal dialysis patient. Plausibly, unmatched patients may possess unmeasured indications (or contraindications) for one modality.
To our knowledge, ours is only the third study to use propensity score analysis in a survival comparison of peritoneal dialysis versus hemodialysis. In an analysis of 1041 participants in the Choices for Healthy Outcomes in Caring for ESRD Patients study, with modality defined at approximately 10 weeks after dialysis initiation, Jaar et al. reported a propensity-score-adjusted HR of death for peritoneal dialysis versus hemodialysis of 1.47 (P = 0.32) in the first year and 2.05 (P = 0.03) in the second year, with little evidence of effect modification by age, cardiovascular disease, or diabetes.10 The cohort was notable for relatively low mortality among hemodialysis patients (24.4% over mean follow-up time of 2.4 years). In a propensity-score-matched cohort of 16,791 patients initiating dialysis in Australia and New Zealand between 1991 and 2005, McDonald et al. reported HRs of death for peritoneal dialysis versus hemodialysis of 0.99 (P = 0.80) from day 90 to 365 and 1.35 (P < 0.01) from day 366 on.2 In tandem with our secondary analysis of survival from day 90, these studies provide little evidence of an overall association between modality and survival in the first year, but stronger evidence of increased relative risk among peritoneal dialysis patients in subsequent years. However, these studies exclude patients who died within the first months of renal replacement therapy.
We found that survival was modestly longer for patients initiating renal replacement therapy with peritoneal dialysis than for matched hemodialysis patients. However, relative mortality hazards were functions of age, cardiovascular disease, diabetes, time frame of follow-up, and exposure definition (intention to treat, as-treated). In general, for patients expected to survive longer (e.g., aged <65 years, without cardiovascular disease, without diabetes), peritoneal dialysis was associated with longer survival than hemodialysis. For patients whose survival was expected to be shorter, the opposite tendency was apparent. Previous studies have reported modality-survival effect modification by age and diabetes.2,4,8,11–16 Attenuation and reversal of longer survival for peritoneal dialysis patients during the first year of dialysis have also been described elsewhere.2,7,8,10,11,15,16 Sensitivity of conclusions to study design choices likely partially reflects nonrandom switching from one modality to another.24 Our findings are also compatible with the hypothesis that longer survival may be associated with a dual-modality approach of initiation with peritoneal dialysis, followed by timely transfer to hemodialysis.25 Although first-year mortality rates for hemodialysis patients have declined very little in the past 15 years, mortality rates after the first year have declined steadily.19 Notably, we found that modality-related survival estimates depended on whether follow-up began at day 0 or day 90. This observation is compatible with the hypothesis that imbalances of unmeasured comorbid conditions, especially those leading to early demise, may be responsible. One could also reasonably argue that this difference reflects more extensive predialysis education for peritoneal dialysis patients, because observational evidence suggests that predialysis education is associated with lower risk of hospitalization.26 On the other side, catheter dependency clearly remains high among U.S. hemodialysis patients during the first 3 months, and catheter use is associated with high morbidity and mortality rates.27 Differential rates of residual kidney function loss may also be involved28; although we adjusted for baseline eGFR, serial eGFR levels were unfortunately not available.
Peritoneal dialysis use in the United States was lower than in many Asian and European countries in 2003.19 In this study, only 6.8% of adults used peritoneal dialysis at dialysis initiation, continuing a 2-decade decline from 15.8% in 1985 to 14.4% in 1990, 13.8% in 1995, and 8.1% in 2000.19 Several causes of this decline have been postulated. First, the numbers of peritoneal dialysis patients treated at individual U.S. dialysis facilities are often modest. For example, in 2003 there were 27,345 prevalent peritoneal dialysis patients and 4591 dialysis facilities, a ratio of 6.0 patients per facility, considerably lower than the ratio of 61.6 hemodialysis patients per facility.29 A low number of peritoneal dialysis patients per facility has been associated with increased rates of technique failure in the United States30,31 and in Canada.32 Plausibly, high rates of technique failure could lead to erosion of peritoneal dialysis expertise in a facility, diminished enthusiasm for considering peritoneal dialysis in incipient dialysis patients, and limited exposure to peritoneal dialysis among trainee health care professionals.33 Reimbursement considerations may also discourage peritoneal dialysis use in the United States20 because revenue from injectable medications is typically higher for hemodialysis than for peritoneal dialysis patients.34
Choice of dialysis modality may become a more pressing issue in years to come. The U.S. dialysis population is projected to reach nearly 405,000 in 2010 and to exceed 534,000 in 2020.33 Globally, the dialysis population is predicted to reach almost 2.1 million in 2010,35 and with expectations of a worldwide shift in age distribution36 and rapid growth in diabetes prevalence,37 further population growth is likely. In the setting of rising disease burden, the need to consider cost-effectiveness of different dialysis modalities is self-evident. In 2006, the respective Medicare expenditures for hemodialysis and peritoneal dialysis were $71,889 and $53,327 per patient-year.19 For segments of the dialysis population in which peritoneal dialysis outcomes compare favorably with hemodialysis outcomes, these cost differentials are nontrivial. In the future, the evolution of costs attributable to facility acquisition, maintenance, and staffing, along with the bundling of dialysis services in capitated payment systems (e.g., with the addition of injectable medications in the Medicare bundled payment for dialysis beginning in 2011), may factor into decisions concerning dialysis modality.
This study has limitations. Because dialysis modalities were not randomly assigned, causality could not be assessed; propensity matching is not a substitute for randomization. The study was retrospective and based on registry data. Prognostic covariates and biomarkers were measured only at dialysis initiation. The Medical Evidence Report [Centers for Medicare & Medicaid (CMS) form 2728] served as the predominant source of baseline data. This instrument likely has lower sensitivity for identifying comorbid conditions, compared with sensitivities in prospective studies in which precise disease definitions are rigorously applied,38 and does not ascertain the severity of conditions. Moreover, biomarker data were incomplete and serum albumin assay type was infrequently recorded; the effect of incomplete data cannot be fully assessed with multiple imputation, which rests on the unverifiable assumption of data missing at random.39 The net effect of misclassified comorbidity and biomarker data may be residual confounding, possibly resulting in bias favoring peritoneal dialysis. Before 2005, the Medical Evidence Report did not record early nephrologist referral or initial vascular access type, factors that may account for the early survival advantage of peritoneal dialysis. Prognostically important, nonrandom processes resulting in modality changes were not considered in survival analyses. In this regard, because peritoneal dialysis patients are more likely to change dialytic modality than hemodialysis patients,40 as-treated analyses may plausibly be biased in favor of peritoneal dialysis.41 Peritoneal dialysis patients are also more likely to undergo kidney transplant than hemodialysis patients,42 so informative censoring is a possibility; our study conclusions were insensitive to censoring (or not) at the time of transplant. Finally, generalizability to other countries or other eras is uncertain because we studied U.S. patients initiating dialysis in a single calendar year.
Despite its limitations, this study may have clinically relevant features and may help inform dialytic modality choices in the United States. The study represents the experience of the entire incident dialysis population in the United States, and the available sample size allowed us to use an analytical technique (propensity matching) rarely encountered to date in peritoneal-dialysis-to-hemodialysis survival comparisons. In short, the technique is an intuitive approach to a long-standing conundrum: If a substantial number of hemodialysis patients in the United States share no (or very few) similarities with peritoneal dialysis patients, what is the most equitable way to compare survival (or any other outcome) for the two modalities? Increased use of peritoneal dialysis among U.S. adults may be justified, particularly in subgroups with favorable outcomes, but in lieu of randomized comparisons, individual study results should be interpreted cautiously and subjected to rigorous replication.
Concise Methods
Participants and Matching
We analyzed data from the CMS ESRD database. The cohort consisted of adult patients who initiated dialysis in the United States in 2003. We included only those patients who immediately began hemodialysis or peritoneal dialysis according to data from the CMS Medical Evidence Report (form CMS 2728), which physicians and dialysis centers are required to submit for all patients initiating maintenance dialysis in the United States; from the CMS Standard Information Management System, to which dialysis centers are required to report modality changes monthly; and from Medicare claims.
We considered in our analyses the covariates in Table 1, ESRD Network, and payer status (Medicare as primary payer, Medicare as secondary payer, Medicare Advantage, private payer, unknown); these data are all recorded on the Medical Evidence Report. Estimated GFR was calculated using the four-variable Modification of Diet in Renal Disease equation.43 We initially excluded patients with missing information on age, sex, race, or ethnicity (223 of 99,925), and then pediatric (age <18 years) patients (827 of 99,925). Among included patients, missing values remained for primary ESRD cause (3.8%), serum albumin (26.9%), body mass index (0.7%), blood urea (6.0%), estimated GFR (1.8%), and hemoglobin (4.2%). We imputed values for these covariates using multiple imputation,28 with five draws for each missing value from a multivariate normal model of all listed covariates, plus modality.
We estimated the propensity score for peritoneal dialysis prescription at the start of follow-up with a logistic regression model that included all listed covariates as predictors.22 We constructed matched pairs with a greedy matching algorithm23 by randomly selecting for each peritoneal dialysis patient with score p a corresponding hemodialysis patient with score between p − 0.1 and p + 0.1. We also considered a matching caliper of 0.05; analytical results were practically identical and are not shown. After matching, standardized differences were used to assess similarity in baseline characteristics across the dialytic modalities.44
Statistical Analysis
Follow-up continued until the earliest occurrence of death, kidney transplant, loss to follow-up, or December 31, 2006. Death was ascertained from the CMS Death Notification Form (form CMS 2746) and the Social Security Death Index, and transplant was ascertained from CMS Standard Information Management System and Medicare claims. The primary analysis was intention-to-treat survival by dialytic modality from initiation of dialysis (day 0). Secondary analyses included (1) survival analysis from day 90, by modality at day 90; (2) as-treated analysis, with initial modality redefined at first modality switch ≥60 days19; and (3) subgroup analyses, with subgroups defined by baseline age (<65 years, ≥65 years), cardiovascular disease (congestive heart failure, ischemic heart disease, myocardial infarction, cardiac arrest, cardiac dysrhythmia, or pericarditis), and diabetes.
We calculated Kaplan–Meier survival estimates in the matched-pair cohort and subsequently fit a Cox proportional hazards regression model of survival from which we estimated the relative hazards of mortality for peritoneal dialysis versus hemodialysis during all of follow-up and during follow-up years 1, 2, 3, and 4. The Cox model included a robust covariance estimator to account for within-pair correlation.45 Estimates were derived with SAS, version 9.1.3 (Cary, NC).
Disclosures
Eric D. Weinhandl and Jon J. Snyder have no conflicts of interest. Robert N. Foley has consulted for, advised, or received honoraria from 21st Services, Affymax, Chugai, Fresenius, Merck, Novartis, Pfizer, and Pursuit Vascular. David T. Gilbertson has consulted for AMAG Pharmaceuticals and Amgen. Thomas J. Arneson has an ownership interest in Johnson & Johnson. Allan J. Collins has advised or received honoraria from AMAG Pharmaceuticals, American Renal Management, Amgen, Apothecom/Axis Healthcare, Fresenius, NxStage, Quintiles/Innovex, and United Resource Network.
Acknowledgments
This study was supported by a research contract from Baxter Healthcare, Deerfield, Illinois. The contract provides for the Minneapolis Medical Research Foundation investigators to have the final determination of the content of this manuscript. The authors wish to thank Chronic Disease Research Group colleagues Shane Nygaard, BA, for manuscript preparation and Nan Booth, MSW, MPH, for manuscript editing.
Footnotes
-
Published online ahead of print. Publication date available at www.jasn.org.
- Copyright © 2010 by the American Society of Nephrology