Abstract
Background and objectives The relative influence of facilities and regions on the timing of dialysis initiation remains unknown. The purpose of the study is to determine the variation in eGFR at dialysis initiation across dialysis facilities and geographic regions in Canada after accounting for patient-level factors (case mix).
Design, setting, participants, & measurements In total, 33,263 dialysis patients with an eGFR measure at dialysis initiation between January of 2001 and December of 2010 representing 63 dialysis facilities and 14 geographic regions were included in the study. Multilevel models and intraclass correlation coefficients were used to evaluate the variation in timing of dialysis initiation by eGFR at the patient, facility, and geographic levels.
Results The proportion initiating dialysis with an eGFR≥10.5 ml/min per 1.73 m2 was 35.3%, varying from 20.1% to 57.2% across geographic regions and from 10% to 67% across facilities. In an unadjusted, intercept-only linear model, 90.7%, 6.6%, and 2.7% of the explained variability were attributable to patient, facility, and geography, respectively. After adjustment for patient and facility factors, 96.9% of the explained variability was attributable to patient case mix, 3.1% was attributable to the facility, and 0.0% was attributable to the geographic region. These findings were consistent when the eGFR was categorized as a binary variable (≥10.5 ml/min per 1.73 m2) or in an analysis limited to patients with >3 months of predialysis care.
Conclusions Patient characteristics accounted for the majority of the explained variation regarding the eGFR at the initiation of dialysis. There was a small amount of variation at the facility level and no variation among geographic regions that was independent of patient- and facility-level factors.
Introduction
It has been nearly four decades since the widespread availability of dialysis treatment for ESRD. However, the fundamental question of when to start dialysis to optimize survival and minimize morbidity remains unanswered. A growing body of evidence suggests that early initiation of dialysis therapy is associated with no measurable benefit, higher cost, and a potential increase in subsequent mortality (1–8).
Despite these concerns, it is clear that patients with CKD have consistently been starting dialysis at progressively higher levels of eGFR over the last two decades (9,10). These trends are concerning, and many studies have suggested that patient-related factors, such as advancing age, comorbid illness, and/or increasing fragility (case mix), may be driving earlier initiation. However, numerous regional- and facility-level factors have previously been shown to exert strong influences on dialysis-related practices, such as vascular access, dialysis adequacy, and modality selection (11–13). Local expertise, physicians opinion and knowledge, use of eGFR reporting, and reimbursement policies may all influence the initiation of dialysis at the facility or regional level. Two recent survey studies have provided considerable insight (14,15). First, a study of attitudes and opinions of physicians regarding the timing of dialysis initiation suggested that reimbursement may play a role, because physicians practicing at for-profit dialysis facilities were more likely to initiate dialysis earlier (14). Although Canada does not have for-profit dialysis facilities, physician remuneration does increase on the basis of dialysis caseload, thus possibly financially incentivizing early dialysis initiation. Second, a study presenting clinical vignettes shows that physicians were more likely to initiate dialysis earlier if the case was presented with an eGFR as opposed to a serum creatinine value (15). In Canada, eGFR reporting varies considerably by facility and geographic region. These studies suggest that nonpatient-related factors may be contributing to decision-making regarding the initiation of dialysis.
Understanding the relative contributions of patient, facility, and regional factors to the timing of dialysis initiation will aid in the development of appropriate health policy. Measures of facility variation are particularly important in this emerging era of quality improvement programs and pay-for-performance initiatives. To date, no studies have quantitated the relative variation in the timing of dialysis initiation because of regional- and facility-related factors. The purpose of this study is to determine the variation in eGFR at dialysis initiation across dialysis facilities and geographic regions in Canada after accounting for patient-level factors (case mix). We hypothesize that there will be considerable facility- and geographic-level variations across Canada.
Materials and Methods
Population and Data Sources
We studied incident dialysis patients from January of 2001 to December of 2010 captured in the Canadian Organ Replacement Registry (CORR; n=33,263). CORR is a validated registry that includes demographics, comorbidities, dialysis modality, vascular access, transplantation, and mortality on all patients with ESRD in Canada (excluding Quebec) (12,16). Data are collected by dialysis facilities and housed centrally with the Canadian Institute for Health Information. All adult patients (>18 years of age) with a serum creatinine at dialysis initiation were included in our analysis (>93%). Regional ethics board approval was obtained.
Cohort Definitions
Renal function was estimated using the Modification of Diet in Renal Disease equation on the basis of a patient’s blood-sample serum creatinine collected immediately before their first dialysis session. The eGFR was examined as a continuous variable. First visit date with a nephrologist was used to estimate the length of predialysis nephrologist care. Comorbidities (angina, chronic obstructive pulmonary disease, diabetes, malignancy, serious illness, hypertension, lung disease, coronary artery bypass grafting, pulmonary edema, peripheral vascular disease, stroke, cigarette smoking, and acute coronary syndrome) were recorded at the initiation of dialysis. Race was self-reported. Distance to facility was calculated as the direct linear distance in kilometers between a patient's primary residence (estimated from postal code at the time of dialysis initiation) to the nearest dialysis provider using the Vincenty Equation (12). Individual patients and dialysis facilities were de-identified for analytic purposes. Additional laboratory values (hemoglobin, albumin, and phosphate) were recorded in cross-section at dialysis initiation.
Facility-level variables were created to account for differences in care across centers and selected on the basis of clinical relevance, known association with outcomes, and quality of care indicators (13–15,17). The majority of centers (87.3%) provided data for the entire 10-year study period. Variables included the proportion of patients who received dialysis through a central venous catheter, whether the facility offered renal transplantation or peritoneal dialysis, the mean hemoglobin and phosphate of patients at dialysis initiation, the average distance in kilometers between a patient’s primary residence and the nearest dialysis facility, and the center size. Information on patients’ geographic regions was included and categorized as Atlantic (Newfoundland, New Brunswick, Prince Edward Island, and Nova Scotia), Ontario (North, Greater Toronto, East, and Western regions), Prairies (Manitoba, Saskatchewan, and Alberta [North and South regions]), and Pacific (Eastern British Columbia, Vancouver, and Other British Columbia regions).
Statistical Analyses
Patient, facility, and geographic characteristics were compared between patients who were started on dialysis with an eGFR>10.5 or <10.5 ml/min per 1.73 m2. This cutoff was on the basis of data from the Canada USA study, which helped to establish a recommended peritoneal dialysis target for weekly Kt/V urea of 2, translating roughly to an eGFR of 10.5 ml/min per 1.73 m2 (18). Continuous variables of interest were summarized as means or medians with SDs or interquartile ranges as appropriate. Differences in characteristics were determined by the t or the Mann–Whitney test for continuous variables and the chi-squared or the Fisher exact test for dichotomous variables.
Facility and geographic variations were examined using multilevel modeling with patients nested within facilities and facilities nested within geographic regions. A three-level linear regression model was used to assess variables associated with the eGFR at dialysis initiation. Models were adjusted for factors thought to potentially influence decisions to initiate dialysis, including facility-level factors (percentage catheter use, transplantation facility, peritoneal dialysis facility, average hemoglobin and phosphate, average distance that a patient resided from the nearest dialysis center, and number of patients), patient case mix (age, sex, body mass index, race, comorbidities, distance to facility, length of predialysis care, serum phosphate, albumin, and hemoglobin), and calendar year. Unadjusted, fully adjusted, and reduced models were created. Covariates for inclusion in our models were selected a priori, and covariates selected to be included in the reduced models were on the basis of (1) a P value <0.01 in the full model and (2) a P value <0.05 in the reduced model. Patient, facility, and geographic variations were determined by the intraclass correlation coefficient (ICC) (19). The ICCs were calculated by dividing the variance estimate at each level by the total model variance. In our study, the ICC determines the proportion of explained variation in the eGFR at dialysis initiation that is caused by being a member of a particular group, such as patient, facility, and geographic region, and is reported as a percentage. Facility-level variables were centered for the facility averages (20). R2 to determine the percentage of variation explained at each level for the full and reduced models was determined by the method by Raudenbush and Bryk (21), whereas the R2 for the total model was determined by the method by Snijders and Bosker (22). Sensitivity analyses were performed to ensure the validity of our findings. A three-level logistic regression model was created examining the eGFR≥10.5 ml/min per 1.73 m2. In a separate sensitivity analysis, an eGFR cutoff of >12 ml/min per 1.73 m2 was used. We used the SAS GLIMMIX procedure (SAS 9.2) using a logit link and the latent variable approach at the patient level (20). In this model, the ICCs were calculated for the facility and geographic regions by assuming a patient-level variance of π2/3 (23,24). Because emergent dialysis starts are not categorized separately in CORR, we also repeated our analyses by limiting our cohort to individuals who had received nephrology predialysis care for >90 days before dialysis initiation. To interrogate the patient-level variance assumption of π2/3, the ICC was further calculated using a probit link assuming a patient level variance=1. Because patients initiating peritoneal dialysis may start at higher levels of eGFR, we repeated our analyses limited to patients initiating hemodialysis only (25). Furthermore, to examine the effect of reducing the eGFR level at dialysis initiation, we excluded the top 25th percentile of facilities with the highest mean eGFR level at dialysis initiation and repeated our models.
Multiple imputation was used for missing values, with a random draw from the predictive distribution of an imputation model repeated 10 times (26). Analyses were performed using SAS, version 9.2. All hypothesis tests were two-sided, with statistical significance determined at a P value <0.05.
Results
Patient/Regional Characteristics and eGFR
The baseline patient-level (n=33,263), facility-level (n=63), and regional-level (n=14) characteristics are presented in Table 1 in relation to eGFR>10.5 or <10.5 ml/min per 1.73 m2. Patients who initiated dialysis at an eGFR≥10.5 ml/min per 1.73 m2, on average, were older, were women, were on longer predialysis care, and had more comorbid illness. They were more likely to have higher hemoglobin, have lower serum phosphate, and receive peritoneal dialysis as their initial modality. Facilities that started more patients at eGFR>10.5 ml/min per 1.73 m2 were more likely to offer peritoneal dialysis, had higher average hemoglobin, had lower serum phosphate, and had fewer total numbers of patients. Pacific regions of the country had a higher proportion of patients initiating with an eGFR≥10.5 ml/min per 1.73 m2, whereas the central region (Prairies) had the lowest.
Baseline patient, facility, and geographic characteristics of the study cohort according to individuals who initiated dialysis with an eGFR≥10.5 or <10.5 ml/min per 1.73 m2
eGFR at Dialysis Initiation across Dialysis Facilities and Geographic Regions
Overall, the mean eGFR at dialysis initiation was 9.7±4.4 ml/min per 1.73 m2, and this value increased annually from 9.1±4.1 ml/min per 1.73 m2 in 2001 to 10.1±4.5 ml/min per 1.73 m2 in 2010. The mean eGFR across geographic regions ranged from a low in Manitoba (8.2±3.8 ml/min per 1.73 m2) to a high in Other British Columbia (11.9±4.6 ml/min per 1.73 m2). Overall, the proportion of patients initiating dialysis with an eGFR≥10.5 ml/min per 1.73 m2 was 35.3%, and this value increased annually from 27.9% in 2001 to 39.5% in 2010. There was a large difference among geographic regions in the proportion of patients initiating dialysis with an eGFR≥10.5 ml/min per 1.73 m2, with the lowest proportion in Manitoba (20.1%) and the highest proportion in Other British Columbia (57.2%) (Figures 1 and 2).
Crude proportion of patients initiating dialysis with an eGFR≥10.5 ml/min per 1.73 m2 according to geographic region. N values in each region are MB (2027), North ALB (2582), NFLD/NB/NS/PEI (3338), Toronto (3818), North ON (1630), West ON (8149), South AB (1538), East ON (2832), SK (1656), Vancouver (3559), East BC (960), and Other BC (1155). ALB, Alberta; BC, British Columbia; MB, Manitoba; NFLD/NB/NS/PEI, Newfoundland/New Brunswick/Nova Scotia/Prince Edward Island; North, Northern; ON, Ontario; SK, Saskatchewan; South, Southern; West, Western.
Crude proportion of patients initiating dialysis with an eGFR≥10.5 ml/min per 1.73 m2 according to incident year (P<0.01 for trend).
Association between Patient, Dialysis Facilities, and Geographic Regions with eGFR at Initiation
The unadjusted ICCs for patient, facility, and geographic region were 90.7%, 6.6%, and 2.7%, respectively. On adjustment for case mix, facility-level factors, and calendar year, the facility and geographic ICCs reduced to 3.1% and 0.0%, respectively. The fully adjusted model had level-specific R2 values of 26.8% and 67.6% at the patient and facility levels, respectively. The total R2 for the fully adjusted model was 31%. Additional models examining eGFR≥10.5 ml/min per 1.73 m2 and limited to patients with >90 days of predialysis care yielded similar findings (Table 2). When the ICC was calculated using a probit link assuming a patient-level variance=1, the findings were similar (results not shown). In a sensitivity analysis examining eGFR>12 ml/min per 1.73 m2, the results were consistent with ICCs for patient, facility, and geographic region of 88.9%, 8.1%, and 3.0% (unadjusted) and 95.2%, 4.6%, and 0.2% (adjusted), respectively. When our cohort was limited to incident hemodialysis patients only (n=25,201), the ICCs for patient, facility, and geographic region were 89.6%, 8.1%, and 2.3% (unadjusted) and 95.9%, 4.1%, and 0.0% (adjusted), respectively. After exclusion of the top 25th percentile of facilities with the highest mean eGFR level at dialysis initiation, there remained n=25, 158 individuals, and 48 facilities. Among the reduced cohort, the proportion of patients initiating dialysis with an eGFR≥10.5 ml/min per 1.73 m2 was 29.8% (versus 35.3% with all facilities included), and the adjusted ICCs for patient, facility, and geographic region were 97.9%, 2.1%, and 0.0%, respectively.
Multilevel model analysis of the unadjusted and adjusted variations at the patient, facility, and geographic levels for the level of eGFR at dialysis initiation
Factors associated with initiation of dialysis at a higher eGFR from the reduced model are presented in Table 3. On average, women with comorbid illness (with the exception of hypertension) started dialysis earlier. Similarly, patients initiated on peritoneal dialysis and those with higher hemoglobin were more likely to initiate with an eGFR≥10.5 ml/min per 1.73 m2. In addition, lower body mass index, lower serum phosphate, lower average facility–serum phosphate, and decreased distance from a dialysis facility were associated with starting at a higher eGFR.
Variables associated with the initiation of dialysis with an eGFR≥10.5 ml/min per 1.73 m2
Discussion
In this Canadian cohort study of 33,263 patients initiating dialysis between 2001 and 2010, we observed significant variability in the eGFR at dialysis initiation across facility and geographic regions, which was largely explained by patient case mix. Patient-related factors, such as demographics, comorbid conditions, laboratory variables, and length of predialysis care, accounted for over 95% of the explained variability in eGFR at dialysis start. The treatment facility contributed 3.1% variability, whereas the influence of geographic region was minimal. These results suggest that patient factors explain much of the variability in eGFR at dialysis initiation, with a small but significant contribution on the basis of the treatment facility.
This study is the first examining facility- and geographic-level variations for the timing of dialysis initiation. Other studies have examined variability across facilities and geographic regions regarding dialysis practices (27–36). Two recent studies of 173 United States dialysis facilities reported significant facility variations for arteriovenous fistula use and dialysis adequacy (urea reduction ratio >65%) of 7.6% and 11.5%, respectively (28,29). These findings and similar findings have led to the implementation of quality improvement programs targeting these metrics and serve as a baseline measure for improvement. Our findings suggest that patient factors are integrally important regarding the decision to initiate dialysis, with nonclinical factors contributing a smaller role. The small measured degree of regional- and facility-level variation may be indicative that physicians and facilities find early dialysis initiation beneficial and practice accordingly. A recently published survey study of Canadian nephrologists found that the majority, in general, did not agree with early dialysis initiation for asymptomatic patients (37). Furthermore, it should be noted that, although we examined eGFR categorically (>10.5 or <10.5 ml/min per 1.73 m2), this finding does not imply an acceptable absolute eGFR for dialysis initiation. The most recent Kidney Disease Improving Global Outcomes guidelines suggest that an acceptable threshold for the initiation of dialysis would be the presence of more clinical criteria, such as serositis, acid-base or electrolyte abnormalities, and pruritus; inability to control volume status or BP; a progressive deterioration in nutritional status refractory to dietary intervention; and cognitive impairment (38).
We found an increase in facility-level variation when the study cohort was limited to patients on hemodialysis only. This result shows that the variability in dialysis initiation is modality-dependent and higher among patients on intermittent hemodialysis compared with patients on peritoneal dialysis. Although the variability at the facility level increases in patients on hemodialysis, a recent Canadian study reported that patients on peritoneal dialysis are more likely to initiate dialysis at higher levels of eGFR (25). Combining that finding with our observation suggests that the practice of starting patients on peritoneal dialysis at higher levels of eGFR is occurring more uniformly (or with less variability) across Canada.
Patient-level factors associated with initiating dialysis at higher eGFR levels included being a woman, residing closer to a dialysis facility, and having a lower serum phosphate. The low serum phosphate is a novel finding and may be a surrogate marker of malnutrition and decreased protein catabolic rate, which are recognized features of advancing uremia and themselves indications for elective initiation of dialysis. In this regard, other correlates of malnutrition, such as a lower body mass index, were also associated with early initiation in this study. Similar to patients with a higher burden of comorbid conditions, patients with low phosphate are also likely to have diminished protein intake and reduced muscle mass. As such, creatinine-based eGFR may be biased in this population and lead to false conclusions about the timing of dialysis initiation.
Despite the low level of facility-level variation, quality improvement interventions targeting a decrease in the eGFR level at dialysis initiation at the facility level may be warranted (39,40). Previous interventions targeted at further reducing low levels of variation at the facility or physician level have been successful. For example, a multifaceted quality improvement intervention led to improved dialysis adequacy across 196 United States dialysis centers, despite a low initial variability between facilities (33). Similar improvements have been shown in a study examining four primary care indicators (lipids, BP, patient experiences, and mammograms) all with <4% initial facility-level variation (41). On the basis of these studies, it seems that, when total patient-level variability is high but facility- or hospital-level variation is low, interventions that show only minimal improvements at the facility-level can lead to large overall reductions in variation (41,42). This finding may represent the case with our data, because explained adjusted facility variation was 3.1%. However, the unadjusted variation was high at 36.1% (Other British Columbia versus Manitoba in Figure 2). To examine the effect of a quality intervention in our population, we removed the top quartile of facilities with the highest mean eGFR at dialysis initiation; a decline in the proportion of patients initiating dialysis with an eGFR≥10.5 ml/min per 1.73 m2 by 5.5% was seen, and the adjusted variation reduced to 2.1%. The development of evidence-based guidelines, facility-level feedback with national rankings, directed education sessions, and site visits by opinion leaders may represent knowledge translation activities geared toward decreasing variations and improving care. Conversely, it still remains uncertain whether the timing of dialysis initiation is a truly modifiable quality measure, because even in the setting of a stringent randomized control trial (Initiating Dialysis Early and Late, IDEAL), it was very difficult to start patients at set, predetermined ranges of eGFR (8). Because it remains uncertain whether the timing of dialysis initiation is a modifiable metric, it would be premature to include it as a possible pay-for-performance measure.
The strengths of our study include a large nationally representative cohort capturing incident patients in Canada (excluding Quebec) during the study period. We used a multilevel model that accounted for within-facility and geographic regional clustering. We reported the ICC, a well reported method for quantitating the relative variability at differing levels in health care systems (43). Facility-level variables were created to serve as surrogates of a facility's quality of dialysis care. We examined facility-level variation as opposed to individual physicians, because many nephrology facilities in Canada use shared-care models. Our models were adjusted for a large number of covariates, and our findings were consistent across numerous sensitivity analyses.
Limitations of the study include the use of registry data, which lacks information regarding the indication for dialysis initiation or measures of frailty. Our study would not detect variability in facility and regional practices if all physicians were practicing consistently and starting dialysis at higher levels of eGFR in a uniform manner. From 2001 to 2010, the mean eGFR at dialysis initiation increased from 9.1±4.1 to 10.1±4.5 ml/min per 1.73 m2. During this period, the increasing average age and number of comorbidities make interpretation difficult. We used surrogate measures for facility quality, and they may not necessarily reflect factors at the individual physician level. We were unable to detect variation between individual clinicians, because individual physician-level data were not available. Because Canada does not have for-profit dialysis facilities, our study was unable to examine the question of whether initiation of dialysis varies on the basis of the profit status of the dialysis unit. Finally, many provinces began regional eGFR reporting during the study period, and on the basis of recent survey results, eGFR reporting may influence the timing of dialysis initiation (44).
In conclusion, we showed, in a national study of >30,000 patients with ESRD over 10 years, substantial variation in dialysis initiation in Canada that is largely attributable to patient factors, with a small but significant component at the facility level. Whether facility-level variation can be further reduced or modified requires additional investigation.
Disclosures
M.M.S. has been on advisory boards for Roche and Amgen and received speaking honoraria from Roche, Amgen, and Sanofi. S.N. has been on advisory boards for Baxter Healthcare Canada and received speaking honoraria from Baxter Healthcare Canada and Merck Frost. B.M., A.D., B.H., J.K., P.K., A.M., D.N., S.N., O.R., S.S., and M.Z. have no conflicts of interest.
Acknowledgments
We thank Bob Williams, Frank Ivis, Dolores Friesen, the Canadian Institutes of Health Research (CIHR), and the Kidney Foundation of Canada (KFoC).
M.M.S. is provided salary support by the Jindal Research Chair for the Prevention of Kidney Disease at the University of Ottawa. N.T. is supported by the Kidney Research Scientist Core Education and National Training Program New Investigator Award, a joint initiative of the KFoC, the CIHR, and the Canadian Society of Nephrology.
This work was presented as an abstract at the 50th European Renal Association-European Dialysis and Transplant Association Congress in Istanbul, Turkey on May 19, 2013.
M.M.S. 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.
Footnotes
Published online ahead of print. Publication date available at www.cjasn.org.
See related editorial, “What Drives Early Dialysis Initiation and How Do We Optimize Timing of RRT?,” on pages 1671–1673.
- Received December 5, 2013.
- Accepted May 15, 2014.
- Copyright © 2014 by the American Society of Nephrology