Introduction

Mortality is one of the most important outcomes to be taken into consideration when selecting dialysis modalities amongst patients entering ‘end-stage renal disease’ (ESRD). Recent observational studies suggest that the survival of peritoneal dialysis (PD) patients has improved over time, and is now comparable to the survival of hemodialysis (HD) patients1,2,3. However, given that the incidence and prevalence of ESRD in Taiwan are amongst the highest in the world, the increase has become an increasing financial burden on the Taiwan National Health Insurance (NHI) system4,5; and in deed, by 2014, the costs of dialysis were accounting for an astonishing 7.4% of the total annual NHI expenditure in Taiwan6. Since PD has a similar all-cause mortality (ACM) rate to that of HD, but with lower medical costs, the Taiwan NHI Administration has been promoting the utilization of PD as a viable alternative since 20053,7,8; and indeed, over the years, the Taiwan NHI administration has been gradually introducing a progressive program of reductions in the reimbursement rates for HD and a corresponding program of increases in the reimbursement rates for PD.

In order to promote the more widespread use of PD, from May 2008 onwards, the Taiwan NHI payment scheme not only covered continuous ambulatory PD (CAPD) but also was extended to cover the machine costs of automated PD (APD); these are the two most frequently used PD modalities. As a result, there has been a gradually increasing trend in the prevalence of PD usage in Taiwan (including APD and CAPD) from 6.5% in 2003, to 8.5% in 2007 and 10.3% in 20095,9. At the same time, the number of patients treated with PD on a global scale has also risen considerably, with a 2.5-fold increase in the prevalence of PD patients in developing countries between 1997 and 200810.

There is, however, a distinct lack of evidence on the mortality rates of both APD and CAPD modalities in eastern countries, which is clearly of importance in enabling healthcare providers to select the appropriate PD modality. Previous ESRD registries have documented racial differences in the crude mortality rates of patients on dialysis11,12, and whist the clinical outcomes of the two PD modalities have been compared in western countries, the ACM results have been somewhat inconsistent13,14,15,16,17,18. Three multicenter studies in western countries comparing ACM and technique failure (TF) rates between APD and CAPD methods reported similar outcomes13,14,15. A further study carried out in the US found that APD patients had a better technique survival rate, but no significant difference in ACM was identified between APD and CAPD patients16. Two studies undertaken in Mexico and Brazil also reported that as compared to CAPD patients, APD patients had a better technique survival rate and significantly lower ACM17,19.

Several factors may have potentially contributed to the variations in the results reported in the western studies, including geographic variations, racial and ethnic differences, patient choice bias, physician bias and differences in insurance payment scheme for APD modality; and indeed, APD has the distinct characteristic of short dwell times under automated devices, with such usage having increased over recent years. Given the increasing use of both APD and CAPD in Taiwan, an examination of the comparative assessment of the survival rates of the two PD modalities has now become a crucial issue. Our population-based study has therefore been designed to facilitate a comprehensive comparison for ACM and TF rates amongst these two PD modalities.

Results

Demographic characteristics

As shown in Fig. 1, our APD patient cohort (n = 2,346) and CAPD patient cohort (n = 7,175) were originally enrolled into the full sample between 2001 and 2010; however, as a result of the 1:1 propensity-score matching analyses, 2,287 APD and 2,287 CAPD patients were selected for subsequent analysis. The baseline patient characteristics are reported in Table 1. Following the propensity-score matching, each group, based upon age and year at cohort entry, contained the same number of APD and CAPD patients. The other baseline characteristics were found to be similar for both groups with the one exception of cirrhosis of the liver (Table 1). Further analyses were then carried out based upon the first (2001–2004), second (2005–2007), and third (2008–2010) sub-periods. Figure 2 shows the overall cohort number of PD patients. As compared to the first and third sub-periods, a much more rapid increase is discernible in the number of PD patients in the second sub-period (791 patients in 2005, reaching a peak of 1,281 patients in 2008). This sub-period coincides with the promotion of the use of PD by the Taiwan NHI Administration. Figure 2 also shows the percentage of APD patients to the overall cohort of PD patients by year at cohort entry, with a discernible gradual increase in the proportion of APD patients from 11.4% in 2001 to 18% in 2005, reaching a peak of 37.2% in 2009.

Figure 1
figure 1

Enrollment of study participants.

Table 1 Baseline characteristics of the study patients before and after propensity-score matched cohorts.
Figure 2
figure 2

Total numbers of APD and CAPD patients (bar chart) and percentage of APD patients to total PD patients (solid line).

Univariate analyses for all-cause mortality and technique failure

The primary outcomes after matching (as reported in Table 2) show that during the 2001–2010 sample period, a total of 514 (22.5%) APD patients had died, which was slightly higher than the 470 (20.6%) deaths that had occurred amongst the CAPD patients. It is also noted from Table 2 that greater numbers of APD patients had received kidney transplants than CAPD patients during the 2001–2010 period and the first sub-period. The incident rate of ACM was higher amongst APD patients than CAPD patients in the second sub-period, with the comparable incident rate of TF also revealed a consistent result. The events of peritonitis per 1,000 patient-years were significantly lower amongst APD patients than CAPD patients during the first and third sub-periods, with the exception of the second sub-period (RR, 1.08; 95% CI, 0.93–1.24). Increase home visits of APD patients were noted from the second (22.5%) to third (24.6%) sub-period. CAPD patients also had the same trend (25.5% to 28.8%). Supplementary Table S3 shows annual peritoneal dialysis-related quantity of solution of APD and CAPD patients. The total quantities of 1.5% and 2.5% dextrose solution were significantly higher amongst APD patients than CAPD patients.

Table 2 Outcomes of the study patients after propensity-score matching.

The Kaplan-Meier curves for the patient survival rates (technique survival rates) during the overall sample period and the three sub-periods are illustrated in Fig. 3a,b), from which we can see that the survival benefit for APD patients is inferior to that of the CAPD patients during the overall sample period (p = 0.01) and the second sub-period (p < 0.01). Similar technique survival rates are discernible for both the APD and CAPD patients during the overall sample period; however, separate comparisons of each of the sub-period reveal that APD patients had an inferior technique survival rate in the second sub-period (p = 0.001).

Figure 3
figure 3

Kaplan-Meier analyses of (a) patient survival and (b) technique survival probabilities, by cohort periods.

Multivariable analyses for all-cause mortality and technique failure

We subsequently carried out Cox proportional hazard regressions and competing risks analyses to facilitate our examination of the differences on the ACM and TF rates between the APD and CAPD patients (Table 3). As compared to the CAPD patients, the APD patients were found to have a slightly higher risk of ACM and TF with HRs of 1.17 (95% CI, 1.03–1.32; p = 0.02) and 1.17 (95% CI, 1.01–1.36; p = 0.04) after Cox proportional hazard regressions. However, the results of competing risks analyses showed the APD and CAPD patients had similar risks of ACM (HR, 1.12; 95% CI: 0.98–1.27; p = 0.09) and TF (HR, 1.11; 95% CI: 0.96–1.29; p = 0.16) after multivariate adjustment. A further comparison was subsequently undertaken between the ACM and TF rates for APD vis-à-vis CAPD patients in the three sub-periods. The results on ACM, after multivariate adjustment, revealed that the APD patients had a significantly higher risk of ACM than the CAPD patients in only the second sub-period (HR, 1.51; 95% CI: 1.22–1.88; p < 0.001). As regards the results on TF, as compared to the CAPD patients, the APD patients were found to have a significantly lower TF rate in the first sub-period (HR, 0.70; 95% CI: 0.52-0.95; p = 0.02) but a significantly higher TF rate in the second sub-period (HR, 1.51; 95% CI: 1.22–1.88; p < 0.001). In addition, the results of the competing risks analyses were similar with the results of Cox proportional hazard regressions in the three sub-periods (Table 3). We also carried out sub-group Cox proportional hazard regression analyses for the ACM and TF rates of the APD and CAPD patients using the first sub-period as the reference (Supplement Table S4). The results revealed that as compared to the first sub-period, CAPD patients had a lower risk of ACM and TF in the other two sub-periods; however, no similar trend was discernible for the APD patients.

Table 3 Cox proportional hazard analysis and competing risks analysis of APD and CAPD patients after propensity-score matching.

APD-to-CAPD hazard ratios by sample period and age

Figure 4a,b illustrate the adjusted HRs of APD to CAPD after adjusted patient demographics, comorbidities, events of peritonitis and icodextrin usage for all the patients. APD patients had a significantly higher risk of both ACM and TF in the second sub-period, particularly in the first year after the commencement of PD therapy. The adjusted HRs of the TF in the first sub-period fluctuated with inconsistent risk during the follow-up period.

Figure 4
figure 4

Adjusted hazard ratios of APD to CAPD from the final multivariate model for (a) all-cause mortality and (b) technique failure, by year. The APD patients had a significantly higher risk of all-cause mortality in the 2005–2007 sub-period, particularly in the 1st and 3rd years. In contrast, a lower risk of all-cause mortality in the 2001–2004 sub-period was found, particularly in the 1st year. The adjusted hazard ratios of the technique failure fluctuated with a notably higher risk in the 1st and 2nd years of the 2005–2007 sub-period and a considerably lower risk in the 2nd and 4th years of the 2001–2004 sub-period. Patients were followed till December 31, 2011.

The HRs of APD to CAPD patients with regard to ACM and TF after multivariate adjustment by age groups are illustrated in the supplementary Fig. S1. Compared to CAPD patients, APD patients had a similar risk of ACM and TF in all age groups except aged 60 to 69 years, with APD patients having a high risk of ACM (HR, 1.68; 95% CI, 1.28–2.20; p < 0.001) and TF (HR, 1.58; 95% CI, 1.15–2.16; p = 0.01).

Sensitivity analyses for all-cause mortality and technique failure

As shown in Table 4, the results of the sub-group analyses on pure APD patients and pure CAPD patients were found to be similar to the primary analyses in the full sample period and the three sub-periods. The results after Cox proportional hazard regression revealed that APD patients still had a significantly higher risk of ACM and TF than CAPD patients in the full sample period and second sub-period. Similar results were obtained for ACM and TF risk using an alternative definition of delaying the index date to 120 and 180 days and regrouping the APD or CAPD patients. The results are shown in the supplementary tables (Tables S5 and S6). Furthermore, the results of competing risks analyses revealed similar risks of ACM and TF with the results of Cox proportional hazard regression except the result of TF in the full sample period, which showed no notable difference between APD and CAPD patients (HR, 1.15; 95% CI 0.98–1.33; p = 0.08).

Table 4 Cox proportional hazard analysis and competing risks analysis of pure APD and pure CAPD patients.

Discussion

The results of our national study of a propensity-score matched cohort of patients in receipt of APD and CAPD between 2001 and 2010 reveal that both the APD and CAPD patients had similar ACM and TF rates throughout the full sample period and the first and third sub-periods; however, APD patients were found to have a significantly higher risk of both ACM and TF in the second sub-period, particularly in the first year after the commencement of PD therapy. These results, based upon NHIRD data, remain robust under various types of sensitivity analyses. We have also found that as compared to CAPD patients, APD patients had a lower risk of developing peritonitis with the exception of the second sub-period.

There has been a considerable increase in the overall prevalence of PD patients in Taiwan ever since 2005, which was the year in which the Taiwan NHI Administration started promoting greater PD usage; furthermore, from 2008 onwards, reimbursements received from the Taiwan NHI payment scheme were also extended to cover the machine costs of APD. The impact of the policy is not only reflected in the increase in the total number of PD patients but also in the remarkable growth in the number of patients choosing APD as their preferred PD therapy, from about 18% of incident PD patients in 2005 to 30.1% in 2008, and peaking at 37.2% in 2009 (Fig. 2). This trend has also been reported on a global scale; for example, from 1997 to 2008, the use of APD in several developing countries is reported to have increased by 14.5%, whilst in the more developed countries APD usage was up by 30.3%10. In the United States and Canada, the proportion of PD patients treated under the APD modality has been increasingly rapidly, to more than 60%, essentially because the APD modality provides better quality of life and frees the patients for most of their waking hours20,21,22. In addition to coverage of APD machine costs under the NHI payment scheme, other reasons for the selection of the APD modality as the preferred PD therapy include the convenience of the APD system and better mental health quality under APD usage10,22,23.

The evidence on the mortality rates associated with APD and CAPD modalities in East Asia countries is distinctly lacking. Our findings based on the Taiwan NHIRD reveal a significantly higher risk of ACM and TF for APD patients in the second sub-period, particularly in the first year after the commencement of PD therapy (Fig. 3a). The results differ from the findings of other multi-center studies which found no significant differences in ACM and TF rates between the two modalities13,14,15. As shown in supplementary Table S4, CAPD patients had significantly better survival rates in the second and third sub-periods than in the first sub-period, although this trend was not discernible for APD patients in the second sub-period. The impact of the NHI payment scheme on PD therapy led to rapid growth in the number of patients during the second sub-period, raising concerns with regard to the lack of sufficient nursing experience to care for the growing patient numbers. Potential physician choice bias under the new payment scheme by allocating patients towards a particular dialysis modality may also occur. Under the new payment scheme, patients with insufficient self-care ability may have been encouraged to select APD therapy through such bias essentially because the peritonitis rate ratio between APD and CAPD patients was found to be at its highest during this period (Table 2).

The choice to implement CAPD or APD should not only take into account the rate of peritoneal solute transport but also patient preference and clinical condition. According to the European Best Practice Guidelines published in 2005, APD is indicated in (i) inability to achieve efficient blood purification and/or ultrafiltration by means of CAPD; (ii) need to prevent high intraperitoneal pressure; and (iii) patient preference, whereas the prescription of CAPD is well established in patients with a low rate of peritoneal solute transport24. Unfortunately, some important laboratory data regarding the choice of modality and mortality, such as peritoneal transport status, residual renal function, serum albumin level, weekly KT/V and weekly creatinine clearance were not available in this study because the NHIRD does not contain laboratory results. We suggested that the increased number of APD patients during the second sub-period could well be due to high transporters essentially because APD modality is more suitable for these patients25. A meta-analysis study has shown that there was a higher risk of death in high transporters treated with APD as compared to those treated with CAPD (relative risk, 1.15; 95% CI: 1.07–1.23; p < 0.001)26. As such, in the case of limitations of the NHIRD, some potential impact factors could not be adjusted and these patients may well have been over-represented in the APD patient cohort.

The major difference between APD and CAPD procedures include the methods and frequency of setting up manual connections and disconnections from the PD catheter and the dialysate bags, with APD having fewer dialysate exchanges than CAPD; however, the results of our study show that there was a higher risk of TF in the APD modality during the second sub-period, a result which is not consistent with the other two sub-periods or the results of prior multi-center studies in different parts of the world13,14,15. Ever since the Y-set (twin-bag) connecting system was introduced into Taiwan in 1997, it has become standard practice for CAPD therapy, which means that its usage is encompassed in the full period examined in the present study. It has also been the main reason for the prevention of peritonitis in CAPD27. A prior single-center study in Taiwan comparing APD with the twin-bag system revealed lower incidences of peritonitis with the use of the APD (146.4 events/1,000 patient-years vs. 273.6 events/1,000 patient-years, p < 0.001)28. Our national cohort study also reveals that as compared to CAPD patients, APD patients had a lower peritonitis incidence rate ratio in all of the sample periods, with the exception of the second sub-period (Table 2), which could be the major cause of the higher TF rate found in APD patients during the same period. Our sub-group analysis also revealed that APD patients had a higher risk of TF rate in the first two years (Fig. 4b), which indicates that the rapid increase in the numbers of PD patients choosing the APD modality in the second sub-period also resulted in them returning to the HD modality in the first two years. The higher risk of TF rate in APD patients returned to a similar risk level in the third sub-period once the rapid growth in PD patients had subsided. Given the increasing use of APD in Taiwan, it is particularly important to clarify the outcome differences between APD and CAPD patients. We therefore suggest that the significant differences in the risk of ACM and TF between APD and CAPD patients in the second sub-period could be attributable to the impact of the Taiwan NHI payment scheme and the rapid growth in the total numbers of PD patients. The trend of rapid increase in the number of incident PD patient flattened since 2008 probably owing to the detected higher TF rate during the second sub-period by physicians, less promotion of PD therapy by the Taiwan NHI Administration, and the intervention of increased home visits by better experienced nurses.

The main strengths of this study include the use of a national population-based sample and propensity-score matched cohort study involving an East Asian population. However, despite these strengths, these data sources also have some limitations, which largely coincide with those commonly found in administrative database research studies. First, the underlying diseases causing ESRD were unknown and the assignment of patients to the two PD therapies was not random; thus, caution is urged in attempting to infer causality. Second, detailed laboratory results were not included in the NHIRD. So we could not use the assigned ICD-9-CM codes to identify the albumin level, residual renal function, nutritional status, ultrafiltration rate, peritoneal function test and dialysis clearance (weekly KT/V and weekly creatinine clearance), which might be associated with mortality. Third, the CAPD patients selected for this study were those who matched the APD patients; however, some characteristics of the selected CAPD patients may differ from the general population of CAPD patients. Other residual confounding may have been present, even after adjustment for the most relevant covariables; for example, we could not identify from the Taiwan NHIRD how patients or physicians chose the APD or CAPD modality at the time of treatment-modality decision, the self-care ability of PD patients after starting PD, the support system within their family, the educational process of PD or the therapeutic compliance of the PD patients. These biases were inherent to PD modality prescription and could not be overcome without a randomized controlled trial.

In conclusion, this study clarifies the overall survival benefits of APD and CAPD in an East Asian population based upon the use of the Taiwan NHIRD. The APD patients had similar ACM and TF rates as CAPD patients, except in the second sub-period, particularly in the first year after the commencement of PD therapy, which is clearly associated with the rapid growth in the overall number of PD patients. The findings of the study do not suggest the presence of a clear advantage of CAPD over APD. Differences observed between modalities might be attributed to specials circumstances of sample periods. Prospective randomized studies are needed to verify these findings.

Methods

Study design and data sources

We conducted a nationwide retrospective cohort study by retrieving data from the Taiwan NHI Research Database (NHIRD) on all patients in receipt of PD between January 1, 1999 and December 31, 2011. The NHIRD contains healthcare data on over 99% of the entire population of Taiwan (23 million in population); indeed, this database is one of the largest and most comprehensive databases in the world, encompassing information on virtually all patients in receipt of renal replacement therapy in Taiwan between 1995 and 2011. The dataset adopted for this study has been widely utilized for epidemiologic research with the results having been validated for both chronic kidney disease and ESRD29. Comprehensive details on the files used from the NHIRD have already been provided in our previous work30. Each of the diseases examined in this study was defined using the International Classification of Diseases, 9th Revision, Clinical Modification (ICD-9-CM) codes and the corresponding types of treatment received. The Joint Institutional Review Board of Taipei Medical University approved this study, and the informed consent was waived as a result of the personal information obtained from the NHIRD having been de-identified.

Study cohort and participants

All ESRD patients with long-term PD were identified from the registry for catastrophic illnesses; those ESRD patients in need of long-term dialysis in Taiwan, who are required to be confirmed by two nephrologists as suffering from a catastrophic illness according to the NHI Administration, are subsequently exempted from co-payments under the NHI system. Long-term PD patients were defined as those who had received PD therapy as their dialysis modality 90 days after the first dialysis commencement and continued for at least 90 days between January 1, 2001 and December 31, 2010. Baseline characteristics were collected on the sample patients, with these patients subsequently being followed up by referring to the NHIRD database covering the period from January 1, 1999 to December 31, 2011. Patients were excluded from the sample if: (1) they were less than 18 years of age; (2) they had received HD more than 90 days after their initial dialysis treatment; or (3) they were in receipt of a kidney transplant prior to commencing PD treatment. The patients were then categorized as APD or CAPD patients according to whether or not they were in receipt of PD therapy via a cycler for at least 90 days after PD commencement. The types of APD patients enrolled into the present study included: (1) nocturnal intermittent PD with an empty cavity during the day; (2) continuous optimized PD with one or two exchanges before the night session; (3) continuous cycling PD with a full peritoneal cavity during the day; and (4) tidal PD. The index date was defined as the first day on which the corresponding patient first received either APD or CAPD therapy for a period lasting at least 90 days. A total of 2,346 APD and 7,175 CAPD therapy patients were included in this study. Figure 1 provides a schematic diagram of the study sample along with the exclusion criteria.

Propensity score computation and matching

The computation of the propensity scores was based upon the variables listed in Table 1, including patient age and gender, number of hospitalizations before enrollment, Charlson comorbidity index score, comorbidities, medications received at study enrollment, and premium wage classes. The comorbidities were defined as at least three outpatient visit claims or one claim for incident hospitalization. All baseline covariates were extracted from the assigned period during which the subjects were identified (within 3 months to 2 years prior to the index date) (Table 1).

The logistic regression model used in this study for both the calculation and distribution of the propensity scores amongst APD and CAPD patients are described in the supplementary tables (Tables S1 and S2). Although the full sample size ratio of APD to CAPD differed before matching (about 1:3), the ratio increased gradually every year (from near 1:8 in 2001 to less than 1:2 in 2010) (Fig. 2). In order to match by both age group and year at cohort entry, the APD and CAPD patients were 1:1 matched according to the ‘nearest available matching without replacement’ on the estimated propensity score and based upon a ≤0.1 difference in their propensity scores (Fig. 1 and Table 1)31.

Outcome measures

The main outcomes of interest to this study are the ACM rate and the TF rate during the five-year period after the index date. For all of the sample patients under APD or CAPD therapy, the censoring criteria established for this study were that the patients were followed-up for up to 5 years until the last day of follow-up (December 31, 2011), transferred to maintenance HD for at least 60 days, in receipt of a kidney transplant or death (whichever happened first). Given that switching between APD and CAPD modalities during the follow-up period is attributed to the original PD therapy, it was not subjected to censoring.

Since the NHI in Taiwan is a compulsory program, those patients on whom follow-up had failed had invariably died, or in very rare cases, moved abroad. Cases of the death in this study were therefore identified as: (1) the date of death was available and obtained from the NHIRD; or (2) the patients had been withdrawn from the NHI program and had not been enrolled in the NHI beneficiary registry files. If the death of a patient occurred within 60 days after switching from PD to maintenance HD, then the death was attributed to the original PD modality; conversely, if the death of a patient occurred within 60 days of a kidney transplant, then the death was attributed to the transplant.

TF was defined in the present study as a transfer from PD to maintenance HD for at least 60 days, with the examination being undertaken without taking into account any occurrences of death. We also calculated the cumulated episodes of PD-related peritonitis for each patient during the follow-up periods, with peritonitis being defined using ICD-9-CM codes 567.9 and 996.59 and concurrent antibiotics treatment.

Statistical analysis

Patients undergoing APD and CAPD were compared at the baseline, with the patient cohorts being divided into three sub-periods, 2001–2004, 2005–2007, and 2008–2010 based upon the index date, in order to observe any changes over time. We used the independent sample Student t-test to examine normally-distributed continuous variables and the Wilcoxon rank sum test to examine non-normally-distributed continuous variables, whereas categorical variables were analyzed using the Pearson chi-squared (χ2) test. Secondly, the incidence rates for ACM and TF were subsequently calculated based upon the time between the index date and the date of censor, death, or end of follow-up as the patient-months of follow-up. The log-rank test was also used to examine the therapeutic effects, with the survival curves then being charted based upon the Kaplan-Meier method. Both Cox proportional hazard regression and competing risks analysis were subsequently applied to examine the outcomes of ACM and TF between APD and CAPD patients, with the proportional hazard assumption being checked using the Schoenfeld residuals test. In the competing risks analysis, the hazard of death is analyzed regarding kidney transplant and TF as censoring and the analysis of TF hazard is performed with kidney transplant and death regarded as censoring. Loss of follow-up or recovery of residual renal function was extremely rare in this study because that the NHI in Taiwan is a compulsory program and the catastrophic illness of ESRD is required to be confirmed by two nephrologists. The hazard ratio (HR) estimates and 95% confidence intervals (CIs) were calculated according to time. The Cox proportional hazard model also measured the age-specific HRs (decades) for ACM and TF. Differences between groups were considered significant if the two-sided p-value was <0.05. All analyses were performed using SAS 9.3 software (SAS Institute Inc., Cary, NC, USA).

Sensitivity analyses

In order to assess the robustness of our findings, we restricted the analysis to the sub-groups of patients so as to compare the therapeutic outcomes between those patients treated entirely with APD modality (pure APD patients, switching to CAPD modality for less than 60 days) and those treated entirely with CAPD modality (pure CAPD patients, switching to APD modality for less than 60 days). We also conducted a series of analyses defining APD/CAPD usage at intervals of 120 and 180 days after the index date to minimize any misclassification bias, with the APD or CAPD patients being regrouped according to whether or not they were in receipt of PD therapy via a cycler for at least 30 days after the intervals.

Additional Information

How to cite this article: Tang, C.-H. et al. Do Automated Peritoneal Dialysis and Continuous Ambulatory Peritoneal Dialysis Have the Same Clinical Outcomes? A Ten-year Cohort Study in Taiwan. Sci. Rep. 6, 29276; doi: 10.1038/srep29276 (2016).