Related Experiment Video
Updated: Jun 4, 2025

An R-Based Landscape Validation of a Competing Risk Model
Published on: September 16, 2022
Competing Risk Modeling with Bivariate Varying Coefficients to Understand the Dynamic Impact of COVID-19
Wenbo Wu1, John D Kalbfleisch2, Jeremy M G Taylor2
1Division of Biostatistics, Department of Population Health, Division of Nephrology, Department of Medicine, Center for Data Science, New York University.
Insights
The COVID-19 pandemic significantly impacted kidney dialysis patients, with effects varying over time. A new statistical model reveals complex dynamics in readmissions and deaths for these patients.
Area of Science:
- Biostatistics
- Epidemiology
- Public Health
Background:
- The COVID-19 pandemic severely affected patients with end-stage renal disease requiring dialysis.
- Preliminary analyses showed time-varying impacts of COVID-19 on dialysis patient outcomes.
- Existing models were insufficient to capture these complex dynamics.
Purpose of the Study:
- To develop and validate a novel statistical model for analyzing competing risks in dialysis patients during the COVID-19 pandemic.
- To quantify the dynamic effects of COVID-19 on hospital readmissions and deaths in end-stage renal disease patients.
- To assess the temporal variations in COVID-19's impact relative to postdischarge time and pandemic onset.
Main Methods:
- Proposed a bivariate varying coefficient model for competing risks using tensor-product B-splines.
- Developed a proximal Newton algorithm for efficient model fitting on large Medicare dialysis patient datasets.
- Implemented difference-based anisotropic penalization and cross-validation for model stability and parameter selection.
Main Results:
- The proposed model effectively captures the complex, time-varying effects of COVID-19 on dialysis patient outcomes.
- Hypothesis testing confirmed significant variations in COVID-19 effects with postdischarge time and pandemic duration.
- Model performance was validated through applications to Medicare dialysis data and simulation studies.
Conclusions:
- The developed bivariate varying coefficient model provides a robust framework for analyzing time-dependent risks in large patient populations.
- This methodology offers crucial insights into the nuanced impact of pandemics on vulnerable patient groups like those on dialysis.
- The findings underscore the need for dynamic statistical approaches to understand and manage health crises in chronic disease populations.
Abstract:
The coronavirus disease 2019 (COVID-19) pandemic has exerted a profound impact on patients with end-stage renal disease relying on kidney dialysis to sustain their lives. A preliminary analysis of dialysis patient postdischarge hospital readmissions and deaths in 2020 revealed that the COVID-19 effect has varied significantly with postdischarge time and time since the pandemic onset. However, the complex dynamics cannot be characterized by existing varying coefficient models. To address this issue, we propose a bivariate varying coefficient model for competing risks, where tensor-product B-splines are used to estimate the surface of the COVID-19 effect. An efficient proximal Newton algorithm is developed to facilitate the fitting of the new model to the massive data for Medicare beneficiaries on dialysis. Difference-based anisotropic penalization is introduced to mitigate model overfitting and effect wiggliness; a cross-validation method is derived to determine optimal tuning parameters. Hypothesis testing procedures are designed to examine whether the COVID-19 effect varies significantly with postdischarge time and the time since the pandemic onset, either jointly or separately. Applications to Medicare dialysis patients demonstrate the real-world performance of the proposed methods. Simulation experiments are conducted to evaluate the estimation accuracy, type I error rate, statistical power, and model selection procedures. Supplementary materials for this article are available online.
Related Concept Videos
Parametric Survival Analysis: Weibull and Exponential Methods
Weibull Distribution
The Weibull distribution is a flexible model used in parametric survival analysis. It can handle both increasing and decreasing hazard rates, depending on its shape parameter...
Relative Risk
Causality in Epidemiology
Statistical Methods for Analyzing Epidemiological Data
Assumptions of Survival Analysis
Residuals and Least-Squares Property
If the observed data point lies above the line, the residual is positive, and the line underestimates the actual data value for y. If the observed data point lies below the line, the residual is negative, and the line overestimates the actual data value for y.
The process of fitting the best-fit...

