BACKGROUND: Identifying predictors of mortality risk in patients receiving hemodialysis in the medium term enables the development of measures aimed at reducing this risk. However, data obtained from different dialysis centers may be highly heterogeneous, which reduces the reliability of the resulting estimates. In such cases, statistical tools for aggregating results, traditionally used in meta-analyses, may be helpful. This approach can substantially increase confidence in the findings.
AIM: To identify predictors of mortality risk in patients receiving maintenance renal replacement therapy in the medium term, with adjustment for data heterogeneity in a multicenter study.
METHODS: The study included data from a retrospective, continuous November 1, 2017, to April 30, 2020, involving patients aged ≥18 years who had been receiving outpatient maintenance hemodialysis for at least 3 months and had a functioning vascular access. Censoring criteria included transition to peritoneal dialysis or kidney transplantation during the 30-month follow-up period. At the first stage, hierarchical Bayesian Cox proportional hazards models were estimated separately for four clusters of dialysis centers formed territorially. At the second stage, significant predictors of survival were selected, and generalizing models with fixed or random effects were estimated for these predictors (models were selected based on Cochran’s Q test).
RESULTS: Retrospective data over 2.5 years were collected for 2120 patients who met the inclusion and exclusion criteria across 30 dialysis centers. Over the 30-month follow-up, 468 patients died (22.08%); 7 patients (0.33%) were transferred to peritoneal dialysis (all in Ufa); and 42 patients (1.98%) underwent kidney allotransplantation. The analysis demonstrated that only patient age was a risk factor for mortality (hazards ratio [HR], 1.02; 95% confidence interval [CI], 1.01–1.03). All other factors were associated with a reduced risk of death as their values increased: duration of dialysis therapy (HR, 0.95; 95% CI, 0.92–0.99), body mass index (HR, 0.93; 95% CI, 0.87–1), duration of anti-anemic therapy in months (HR, 0.91; 95% CI, 0.86–0.97), and duration of therapy aimed at correction of calcium–phosphate metabolism in months (HR, 0.94; 95% CI, 0.92–0.96). The mean single-pool (sp) ratio of dialyzer urea clearance (K) multiplied by dialysis time (t) to the volume of urea distribution in the patient’s body (V) (spKt/V) showed no significant association with mortality risk: although the HR was 0.35, the CI was very wide (0.08–1.57), indicating that an increase in spKt/V was associated with a reduced risk of death in some patients and an increased risk in others.
CONCLUSION: Based on preliminary differences in survival among dialysis patients, four clusters belonging to different constituent entities of the Russian Federation were identified. An approach that aggregates modeling results across territorially defined clusters of dialysis centers resulted in greater confidence in the estimated survival outcomes of patients receiving outpatient maintenance hemodialysis.
Обоснование. Выявление предикторов риска смерти пациентов на гемодиализе в среднесрочной перспективе позволяет разработать мероприятия, направленные на его снижение. Однако данные, полученные из разных диализных центров, могут быть весьма неоднородными, что уменьшает надёжность получаемых оценок. В этом случае могут помочь статистические инструменты для агрегации результатов, традиционно используемые в метаанализах. Такой подход способен существенно повысить доверие к получаемым результатам.
Цель. Выявление предикторов риска смертельных исходов пациентов на постоянной заместительной почечной терапии в среднесрочной перспективе с поправкой на неоднородность данных при проведении многоцентрового исследования.
Методы. В исследование включены данные ретроспективного сплошного наблюдения с 01.11.2017 по 30.04.2020 за пациентами старше 18 лет, получавшими на постоянной основе амбулаторный гемодиализ в течение не менее 3 мес. и имевшими сформированный сосудистый доступ. Критерии цензурирования: переход пациента на перитонеальный диализ или трансплантация почки в период наблюдения 30 мес. На первом этапе оценивали иерархические байесовские модели Кокса отдельно по четырём кластерам диализных центров, сформированным по территориальному принципу; на втором этапе отбирали статистически значимые предикторы выживаемости и для них оценивали обобщающие модели с фиксированными или случайными эффектами (выбор между моделями проводили согласно Q-тесту Кохрена).
Результаты. Собраны ретроспективные данные за 2,5 года по 2120 пациентам, соответствовавшим критериям включения и исключения, в 30 диализных центрах. Всего суммарно за 30 мес. умерло 468 пациентов (22,08%), переведено на перитонеальный диализ — 7 (0,33%) (только в Уфе), получили аллотрансплантацию почки — 42 пациента (1,98%). Проведённое исследование показало, что только возраст пациента являлся фактором риска смерти (относительный риск [ОР] 1,02; 95% доверительный интервал [ДИ] 1,01–1,03), а все остальные факторы при увеличении их значения уменьшали риск смерти: длительность диализного лечения (ОР 0,95; 95% ДИ 0,92–0,99), индекс массы тела (ОР 0,93; 95% ДИ 0,87–1), длительность антианемической терапии в месяцах (ОР 0,91; 95% ДИ 0,86–0,97), длительность терапии по восстановлению фосфорно-кальциевого обмена в месяцах (ОР 0,94; 95% ДИ 0,92–0,96). Среднее значение однопулового (single-pool, sp) отношения произведения клиренса мочевины через диализатор (K) и времени диализа (t) к объёму распределения мочевины в организме пациента (V) (spKt/V) не имело статистически значимой ассоциации с риском смерти: несмотря на то, что ОР 0,35, его ДИ очень широкий (0,08–1,57), т. е. его увеличение для части пациентов приводит к снижению риска смерти, а для части, напротив, к его увеличению.
Заключение. На основании предварительно выявленных различий в выживаемости диализных пациентов было определено 4 кластера, сформированных по территориальному принципу принадлежности к разным субъектам РФ. Применение подхода, агрегирующего результаты моделирования по кластерам диализных центров, сформированным по территориальному признаку, позволило повысить доверие к получаемым оценкам выживаемости пациентов на амбулаторном гемодиализе.